计算方法上机程序之龙贝格求积公式

用c++实现龙贝格求积公式的数学算法，是数值计算的上机程序之一

#include
#include
#define eps  0.000001
void romberg（double adouble b）;
double  f（double y）;
void main（）
{
	double a=1.0;
	double b=3.0;
	romberg（1.03.0）;

}
double  f（double x）
{
	double y;
	y=1/x;
	return y;
}
void romberg（double adouble b）
{
	
	double hT2S2C2R2T1C1S1R1;
	h=（b-a）/2;
	int n=1;
	int k=0;
	T2=h*（f（a）+f（b））;
	S2=0;
	C2=0;
	R2=0;
while （fabs（R2-R1）>eps）
{
	R1=R2;
	T1=T2;
	S1=S2;
	C1=C2;

	double sum=0.0;
	for（int i=1;i<=n;i++）
	{
		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.0;
}
//cout<//cout<cout<<“隆贝格求解的结果：“<
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       8440  2008-12-22 22:53  龙贝格求积\Debug\longbeige.obj

     文件      61440  2008-12-22 22:53  龙贝格求积\Debug\vc60.pdb

     文件     225347  2008-12-22 22:53  龙贝格求积\Debug\龙贝格求积.exe

     文件     517120  2008-12-22 22:53  龙贝格求积\Debug\龙贝格求积.pdb

     文件        730  2008-12-22 22:53  龙贝格求积\longbeige.cpp

     文件       4331  2008-12-16 18:57  龙贝格求积\龙贝格求积.dsp

     文件        528  2008-12-16 18:50  龙贝格求积\龙贝格求积.dsw

     文件      41984  2008-12-22 22:53  龙贝格求积\龙贝格求积.ncb

     文件      48640  2008-12-22 22:53  龙贝格求积\龙贝格求积.opt

     文件        906  2008-12-22 22:53  龙贝格求积\龙贝格求积.plg

     目录          0  2008-12-23 08:38  龙贝格求积\Debug

     目录          0  2009-01-02 19:41  龙贝格求积

----------- ---------  ---------- -----  ----

               909466                    12