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

《复化辛普森求积公式的Matlab实现.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

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

2、mel(c)-1;t=t+feval('f',d);ends=h/6*(feval('f',a)+feval('f',b)+2*s+4*t);