m5_1.asv

来自「数值微分和数值积分:Newton-Cotes公式、复化公式和区间逐次半分法、外推」· ASV 代码 · 共 13 行

%%%%%%%%%%%%%%%%  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));