龙贝格积分MATLAB编程实验报告 一、 问题描述：编写函数实现龙贝格积分。二、 实验步骤（过程）：（1） 龙贝格积分以书本例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);endfor j=2:n n=n-1; for i=1:n T(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); for j=2:n if abs(T(1,j)-T(1,j-1)=R break; end enda=T(1,j)；a=vpa(a,7) %a即为所求答案三、 结论：（一）龙贝格积分