龙贝格算法的matlab实现

作业三——龙贝格算法的matlab实现程序流程图：输入a,b,e,k)]()([2)1(1bfafabT1221)(0)1(11]2)12([221iiiliiabiafabTTi=1,2,3,···k144)(1)1(1)1(1ininiiniTTTabs(T(i+1,i+1)-T(i,i))=e&i=4输出计算结果输出“所求次数不够或不可积”-欢迎下载2程序源代码：文件f.mfunctionfx=f(x)ifx==0fx=1;elsefx=sin(x)/x;endend文件longbeige.mclcclearall;formatlonga=input('请输入你要求得积分的下限:');b=input('请输入你要求得积分的上限:');e=input('请输入你要求得积分的结束精度:');k=input('请输入你要求得积分的最大次数:');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=4break;else;endend-欢迎下载3fori=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))