数值分析实验报告-C语言编程-龙贝格.doc

《数值分析实验报告-C语言编程-龙贝格.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、班级：*****学号：******姓名：成绩：________数值分析实验报告实验3数值积分3.1实验目的通过本实验理解数值积分与微分的基本原理。掌握数值积分中常见的复合求积公式的编程实现。掌握龙贝格算法的基本思路和迭代步骤；培养编程与上机调试能力。3.2算法描述3.2.1龙贝格算法基本思路先算出他T0（0），从而计算T0（1）以此类推，在判断

2、Tn（0）-Tn-1（0）

3、

4、数据积分近似值。利用Romberg方法计算函数3.3实验内容用龙贝格算法计算：3.4实验步骤3.4.1代码#include"iostream.h"#include"math.h"#definee0.001doublef(doublex){doubley;if(x==0){returny=1.0;}elsey=sin(x)/x;returny;}voidromberg(doublea,doubleb){intn=1,k=0;doubleh,T2,S2=0,C2=0,R2=0,T1,C1,S1,R1;h=(b-a)/2

5、;T2=h*(f(a)+f(b));while(fabs((R2-R1))>e){R1=R2;T1=T2;S1=S2;C1=C2;doublesum=0;inti;for(i=1;i<=n;i++){sum=sum+f(a+(2*i-1)*h);}T2=T1/2+sum*h;S2=(4*T2-T1)/3;C2=(16*S2-S1)/15;R2=(64*C2-C1)/63;n=n*2;k++;h=h/2;}cout<<"*****最后结果为:"<<"I="<

6、b;cout<<"***输入上下限a,b的值用空格隔开***"<>a>>b;cout<<"***下限a="<

7、************************************************************************************************