Discussion :

[saad.pizzicato]

Bonjour tout le monde,

j'ai une fonction F définie comme vecteur de valeurs V=(F(i)).
Je voudrai calculer une l'intégrale de cette fonction.
J'ai déjà essayé la commande TRAPZ avec un pas d'intégration de 1/120 entre 0 et 1...Mais je ne suis pas très satisfait de la précision des résultats...

N'y a t-il pas des méthodes plus précises déjà présente sur MATLAB: Méthode de Simpson, Romberg etc...


Merci...

[membreComplexe12]

bonjour,

1) méthode de simpson adapatative :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

q = quad(fun,a,b,tol)

2) méthode de newton cotes adapatative :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

q = quadl(fun,a,b,tol)

tu peux fixer la précision en remplissant donnant la tolérance "tol"

la méthode de Romberg n'est pas présente à ma connaissance, il faut l'implémenter soit même (si tu le fais ça m'intéresse de voir comment du là programmée)