m5_1.asv

来自「matlab 相关的数值积分 例如高斯积分 龙贝格积分等等 可参考。」· 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));