simp.asv

来自「基于matlab的数值计算试验程序。equ_direct」· ASV 代码 · 共 14 行

function I=simp(f,a,b,n)
%n为偶数,复化simpson公式
h=(b-a)/n;
s=feval(f,a);
for i=1:2:n-1
x(i)=a+h*i;
s=s+4*feval(f,x(i));
end
for i=2:2:n-2
    x(i)=a+h*i;
    s=s+2*feval(f,x(i));
end
s=s+feval(f,b);
I=h*s/3;