复化辛普森求积公式.doc

《复化辛普森求积公式.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、f=input('请输入函数:f(x)=','s')%输入函数例如：1/(1+x^4)a=input('请输入:a=')b=input('请输入:b=')%区间[ab]n=input('请输入节点数目（包括两个端点）n=')%取不同的n值比较误差h=(b-a)/(n-1)fori=1:n-1x(i)=a+h*(i-1)%计算节点X(i)=x(i)+0.5*h%计算节点中点yn=eval(subs(f,'x','b'))%y(1)=ay(i)=eval(subs(f,'x','x(i)'))%计算节点函数值Y(i)=eval(subs(f,'x','X(i)'))%计算节点中点函数值S1

2、(1)=0S1(i+1)=S1(i)+y(i)%求f(xK)的和S2(1)=0S2(i+1)=S2(i)+Y(i)%求f（xk+0.5）Sn=h/6*(-y(1)+4*S2(i+1)+2*S1(i+1)+yn)%复化辛普森求积公式end