龙贝格算法的matlab实现.doc

《龙贝格算法的matlab实现.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、作业三——龙贝格算法的matlab实现程序流程图：输入a,b,e,k输出“所求次数不够或不可积”输出计算结果abs(T(i+1,i+1)-T(i,i))<=e&i>=4i=1,2,3,···k程序源代码：文件f.mfunctionfx=f(x)ifx==0fx=1;elsefx=sin(x)/x;endend文件longbeige.mclcclearall;formatlonga=input('请输入你要求得积分的下限:');b=input('请输入你要求得积分的上限:');e=input('请输入你要求得积分的结束精度:');k=input('请输入你要求得积

2、分的最大次数:');fx=@(x)sin(x)/x;lbg(@f,a,b,k,e)文件lbg.mfunctionlbg(fx,a,b,k,e)T=zeros(k,k);T(1,1)=(b-a)*(1+fx(b))/2;fori=1:km=0;forj=1:2^(i-1)m=m+fx(a+(2*j-1)*(b-a)/(2^i));endT(i+1,1)=0.5*T(i,1)+(b-a)*m/2^i;forn=1:iT(i+1,n+1)=(4^n*T(i+1,n)-T(i,n))/(4^n-1);endifabs(T(i+1,i+1)-T(i,i))<=e&i>=4

3、break;else;endendfori=1:kifT(i,1)==0j=i;break;else;endendifj==kerror('所求次数不够或不可积')else;endT=T(1:j-1,1:j-1)disp('所求的积分值为：')disp(T(j-1,1))