复化辛普森求积公式的Matlab实现1页.doc

%复化辛普森求积公式的Matlab实现function s=simpson(f,a,b,n,x,x1)n=input(please input n=);a=input(please input a=);b=input(please input b=);% f是被积函数，调用时表达式加单引号% a,b是积分上下线% n是子区间的个数% s是梯形总面积h=(b-a)/n);s=0;t=0;for k=1:n-1x=a+k*h;s=s+feval(f,x);endfor k=1:n-1 c=a+(1/2+k)*h; d=numel(c)-1; t=t+feval(f,d);ends=h/6*(feval(f,a)+feval(f,b)+2*s+4*t);