m5_1.asv

这是我从一本基于MATLAB实现的数值分析算法的书中自带的光盘程序

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