2012-2013-2实验9 数值积分实验.docVIP

实验9 数值积分实验 成绩 专业班级 数学112 学号 201112010 姓名 陈贵红 报告日期 5月27日 . 实验类型：●验证性实验 ○综合性实验 ○设计性实验 实验目的：进一步熟练掌握变步长数值积分算法，提高编程能力和解决定积分问题的实践技能。 实验内容：用龙贝格积分算法计算 实验原理 龙贝格积分算法：生成j=k 的近似积分结果逼近表，并以R(j+1,j+1)为最终解来逼近积分。 R(j,0)=T(j), j=0,T(j)为区间逐次减半递推梯形求积分公式算出的结果； R(j,1)=S(j), j=1,S(j)为区间逐次减半递推辛普森求积分公式算出的结果； R(j,2)=B(j), j=2,B(j)为递推布尔求积分公式算出的结果； 实验步骤 1 要求上机实验前先编写出程序代码 2 编辑录入程序 3 调试程序并记录调试过程中出现的问题及修改程序的过程 4 经反复调试后，运行程序并验证程序运行是否正确。 5 记录运行时的输入和输出。 实验总结 实验报告：根据实验情况和结果撰写并递交实验报告。 参考程序 在编辑窗口中输入： function [R,quad,err,h]=romber(f,a,b,n,tol) M=1; h=b-a; err=1; J=0; R=zeros(4,4) R(1,1)=h*(feval(f,a)+feval(f,b))/2; while