龙贝格积分MATLAB编程实验报告(共2页).doc

精选优质文档-倾情为你奉上龙贝格积分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