[MATLAB] probleme de fonction

j'ai fait une fonction de simpson 1/3 pour un devoir et je me demande ce qu'il manque pour quelle fonctionne


à l'entree je marque disons

fonct = sin(x)
a=0
b=pi/2
n=2
Simpson_Un_Tiers (fonct,a, b, n)


Voici la fonction:

function I = Simpson_Un_Tiers (fonct,a, b, n)

h = (b - a)/2*n;
I = 0;
for i = 1:n
fgauche = fonct(a);
fmilieu = fonct(a+n);
fdroite = fonct(a+2*n);
I = I + h*(1/3*fgauche + 4/3*fmilieu + 1/3*fdroite);
end

Merci Keven