数值计算方法实验报告--龙贝格积分法的上机实验.docx

数值计算方法实验报告实验报告题目：龙贝格积分法的上机实验实验要求龙贝格积分法的上机实现二、 实验分析（包括数学原理，小组分析讨论后确定实验方案和实现思路）例子：用程序计算。WTkdx的值。精确到小数点后第10位。三、 实脸步骤(过程)(包括程序及上机的实现的结果)首先要定义龙贝格积分函数，再调用该函数进行定积分求解。龙贝格积分函数程序如下：functionR,quad,err,h=romber(f,a,b,n,de1ta)为f是被积函数；a,b分别是积分的上下限；n+1是T数表的列数；de1ta是允许误差；R是T数表；quad是所求积分值；M=I;h=b-a;err=1J=O;R=zeros(4,4);R(1,1)=h*(feva1(,f,a)+feva1(,f,b)/2whi1e(err>de1ta)&(J<n)(J<4)J=J+1;h=h2;s=0;forp=1:Mx=a+h*(2*p-1);s=s+feva1('f,x)endR(J+1,1)=R(Jr1)2+h*s;M=2*M;forK=I:JR(J+1,K+1)=R(J+1,K)+(R(J+1,K)-R(JzK)/(4K-1)；enderr=abs(R(J,J)-R(J+1,K+1);endquad=R(J+1,J+1)在命令窗口输入：>>romber(,sqrt(4*x-x2,z0,2,7,0.5*(10(-10)回车得到：quad=3.1413ans=2.000000000002.73212.97610000002.99573.08363.0908000003.08983.12123.12373.124200003.12333.13443.13533.13553.13550003.13513.13913.13943.13943.13943.1394003.13933.14073.14083.14083.14083.14083.140803.14083.14133.14133.14133.14133.14133.14133.1413故4x-x2dx3.1413四、 总结（包括实险过程遇到的情况等，组长总结组员在整个过程的参与情况）采用龙贝格求积公式，即逐次对分积分区间的方法，可以把前面计算的结果作为一个整体代入对分后的计算公式中，只需增加新的分点的函数值。