Numerical Integration of Equally Spaced Data

I have been using trapz to integrate a vector of equally spaced data. Is there a more accurate method available in script form that carries out the Simpson's method or perhaps the Newton-Cotes method? Thanks.

Answer 1

[Youssef Khmou]

hi David,

I think using trapz integrator with the x/t axis gives good result, in the mean wile try to interpolate your data to get better accuracy , as this example :

 t=0:1/100:pi/2;
 y=sin(t);
 f=trapz(t,y);

 t2=interp(t,4);
 y2=sin(t2);
 f2=trapz(t2,y2)      %  Df=1.0064-0.9992=0.0072

[David Koenig]

Thanks, Youssef. I will give your idea a try.

[Youssef Khmou]

ok, do not accept the answer until it works to get other suggestions from other contributors