龙贝格积分MATLAB编程实验报告.doc

《龙贝格积分MATLAB编程实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、龙贝格积分MATLAB编程实验报告一、问题描述：编写函数实现龙贝格积分。二、实验步骤（过程）：（一）龙贝格积分——以书本例4.6来说明clc,clear;symsx;n=8;a=0;b=1;R=0.5*10^(-6);%取足够大的0矩阵来构建Romberg积分表T=zeros(n,n);y=sin(x)/x;p=subs(y,x,a);q=subs(y,x,b);ifa==0p=1;endT(1,1)=(b-a)/2*(p+q);fori=2:nf=0;forj=1:2^(i-2)t=a+((2*j

2、-1)/2^(i-1))*(b-a);z=subs(y,x,t);f=f+z;endT(i,1)=0.5*T(i-1,1)+f*(b-a)/2^(i-1);endforj=2:nn=n-1;fori=1:nT(i,j)=(4^(j-1)*T(i+1,j-1)-T(i,j-1))/(4^(j-1)-1);endendT=vpa(T,7)T=eval(T);n=8;w=ones(1,7);forj=2:nifabs(T(1,j)-T(1,j-1))<=Rbreak;endenda=T(1,j)；a=vp

3、a(a,7)%a即为所求答案一、结论：（一）龙贝格积分