simp.asv

基于matlab的数值计算试验程序。equ_direct

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;