二、 实验步骤（过程）：（一） 龙贝格积分以书本例4.6来说明clc,clear;syms x;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）;if a=0 p=1;endT（1,1）=（b-a）/2*（p+q）;for i=2:n f=0; for j=1:2（i-2） t=a+（2*j-1）/2（i-1）*（b-a）; z=subs（y,x,t）; f=f+z; end T（i,1）=0.5*T（i-1,1）+f*（b-a）/2（i-1）;for j=2: n=n-1; for i=1: T（i,j）=（4（j-1）*T（i+1,j-1）-T（i,j-1）/（4（j-1）-1）; T=vpa（T,7）T=eval（T）;w=ones（1,7）; for j=2: if abs（T（1,j）-T（1,j-1）=R break; enda=T（1,j）；a=vpa（a,7） %a即为所求答案三、 结论：（一）龙贝格积分
升级会员 