#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