MATLAB: Numerical Integration of Equally Spaced Data

equally spaced dataMATLABnumerical integrationsimpsons ruletrapezoid rule

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.

Best Answer

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
Related Question