intg - intégration numérique adaptative

Calling Sequence

[v,err]=intg(a,b,f [,ea [,er])

Parameters

• a,b : nombres réels
• f : fonction externe (fonction Scilab ou chaîne de caractères ou liste).
• ea, er : nombres réels
• ea : erreur absolue désirée (0 par défaut)
• er : erreur relative (1.d-8 par défaut)
• err : estimation de l'erreur absolue sur le résultat

Description

intg(a,b,f) approche l'intégrale de a à b de f(t)dt . L'évaluation satisfait si possible abs(I-v)<= max(ea,er*abs(I)) où I représente la valeur exacte de l'intégrale.

f est une fonction externe :

Si f est une fonction Scilab elle doit avoir la liste d'appel : y = f(t)

Si f est une liste, cette liste doit avoir la structure suivante : list(f,x1,x2,...) où f est une fonction Scilab avec la liste d'appel : f(t,x1,x2,...) .

Si f f est une chaîne de caractères, elle désigne le nom d'une subroutine Fortran ou C (voir le code source de fintg.f )

Examples

deff('[y]=f(x)','y=x*sin(30*x)/sqrt(1-((x/(2*%pi))^2))')
exact=-2.5432596188;
abs(exact-intg(0,2*%pi,f))
// voir routines/default/Ex-intg.f 
abs(exact-intg(0,2*%pi,'intgex'))
 

See Also

intc ,   intl ,   inttrap ,   intsplin ,   ode ,