Bonjour, je ne comprends pas , après maintes vérifications, pourquoi les 2 programmes suivants ne donnent pas, par la méthodes trapèzes le même résultat (environ 450) pour le calcul de l'intégrale de -10 à 8 de x-> x^2+ 2*x -1.
La première façon donne 453 si je me rappelle et la deuxième en faisant : I=trapez_n('f',-10,8,18) en mode interactif, on obtient 3224,5. Je ne vois pas pourquoi après recherche des fautes. Merci de votre aide. Amicalement.
(J'ai choisi 18 comme dernier argument de trapez_n afin de retrouver h=1).
1ère façon :
suivi en mode interactif dans MATLAB(mode commande) des lignes :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2 function I=trapez_v(g,h) I=(sum(f)-(f(1)+f(length(f)))/2)*h;
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3 >>x=-10:1:8; >>f=x.^2+2*x-1; >>h=1;
2ème façon :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5 function I=trapez_n(fonction_f,a,b,n) h=(b-a)/n; x=a+(0:n)*h; f=feval(fonction_f,x); I=trapez_v(f,h)
Partager