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