Matlab - Calcul d'intégrale

quad : algorithme de Simpson

 

quad8 : algorithme de Newton-Cote (version 4.2.)

quadl (version 6.5.)

 

Il faut commencer par créer le fichier .m de la fonction à intégrer :

>> fplot( 'f4', [ 0 1 ])

>> grid on

>> quad( 'f4' , 0 , 1 )

ans =

3.1416

>> format long e

>> quad( 'f4' , 0 , 1 )

ans =

3.141592502458707e+000

 

La valeur exacte de cette intégrale est pi.

>> pi

ans =

3.141592653589793e+000

 

Il y a un petit écart ...

Améliorons la précision avec une tolérance de calcul de 1e-10 (la tolérance par défaut est 1e-3) :

>> quad( 'f4' , 0 , 1 , 1e-10 )

ans =

3.141592653589792e+000

 

... ce qui est beaucoup mieux !

>> quad8( 'f4' , 0 , 1 , 1e-10)

ans =

3.141592653589793e+000

 

© Fabrice Sincère

http://perso.orange.fr/fabrice.sincere/