integrate - intégration numérique
calcule :
/x1
[
x = I f(v)dv
]
/x0
L'approximation satisfait (on l'espère) l'affirmation suivante sur l'erreur : abs(I-x)<= max(ea,er*abs(I)) où I représente la valeur exacte de l'intégrale.
integrate('sin(x)','x',0,%pi)
integrate(['if x==0 then 1,';
'else sin(x)/x,end'],'x',0,%pi)