Discussion :

[phil971]

Bonsoir a tous je fais actuellement un TP sur l'intégration numérique,
et je cherche a calculer grâce a la formule du point milieu .
voici la formule :
Qm=((b-a)/n)somme (de j=0 à n-1) f((aj+aj+1)/n)

voici mon code:

Code :

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

1
2
3
4
5
6
7
function Qm = mil(f,a,b,n)
Qm=0;
x=[0:n-1];
feval(f,x);
for k=0:n-1
    Qm=Qm+f(a+k*(b-a)/n);
end

Et ca ne fonctionne pas j'ai une erreur en ligne 9
Qui pourrait m'aider svp ?
Merci d'avance pour votre aide !!

[christophe_halgand]

Salut,

Peux tu faire un copier coller de toute l'erreur donné par matlab ?
De plus donne nous cette ligne 9...

A bientôt

Christophe

[phil971]

Voici le code final de la méthode composite du point milieu

Code :

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

1
2
3
4
5
6
7
8
9
10
%%%%%%%%%% CODE %%%%%%%%%
 
function Qm = mil(f,a,b,n)
H=(b-a)/n;
x = linspace (a+H/2,b-H/2,n);
fmp = feval (f,x);
Qm =H*sum (fmp );
 
 
%%%%%%%%%% CODE %%%%%%%%%%