m5_1.asv

数值微分和数值积分:Newton-Cotes公式、复化公式和区间逐次半分法、外推法和Romberg积分、自适应Simpson积分法 [ReMap.rar] - EASYARM2200教学实验平台上

%%%%%%%%%%%%%%%%  Trapezoidal or Simpson  %%%%%%%%%%%%%%%%%%%%%%%%%%%
function Simpson=Simpson(n,a,b)
h=(b-a)/n;
x=[a:h:b];
for i=1:n
     xx(i)=x(i)+h/2;
end
yy=f(xx);
y=f(x);
Result=h/6*(4*sum(yy)+2*sum(y)-f(a)-f(b));