Calcul intégral: matlab

Bonjour,

J'ai un exercice à réaliser pour mon cours.

Je doit calculer le moment d'ordre 1 d'une gaussienne.

Soit y(x) ma gaussienne et x= [-3*std,+3*std]
si je calcule le moment d'ordre 1, normalement je dois trouver le centre de ma gaussienne ????
Le moment d'odre1:
m1=intégrale (y(x)*x dx)/intégrale(y(x)dx) ??
donc
en matlab:

Code:

---

1 2	A=y*x%dénominateur m1=trapz(A)/trapz(y);

---

Mais en faisant ça, je ne trouve pas le centre de ma gaussienne.

Je ne comprends pas mon erreur.

Merci d'avance.
Marc