#include
#include
#include
#define?N????20
#define?Max??10
#define??a????0
#define??b????1
#define?eps???1e-6
#define?pi???3.14159

double?f（double?x）???????????????????????????????????????/*計算f（x）在各點值的函數*/
{double?y;
	y=4/（1+x*x）;
	return（y）;
}
double?trapz（double?a1double?b1long?int?n）??????????????????/*梯形公式求小區間上的積分函數*/
{int?i;
??double?sumh=（b1-a1）/n;
??for（i=1;i??sum+=f（a1+i*h）;
??sum+=（f（a1）+f（b1））/2;
??return（h*sum）;
}
void?main（）????????????????????????????????????????????/*Romberg積分程序主函數*/
{int?ik;
?long?int?n=Nm=0;
?double?T[Max+1][2];
?T[0][1]=trapz（abn）;
?n*=2;
?for（m=1;m??{
???for（i=0;i??????{T[i][0]=T[i][1];}
???T[0][1]=trapz（abn）;
???n*=2;
???for（i=1;i<=m;i++）
????T[i][1]=T[i-1][1]+（T[i-1][1]-T[i-1][0]）/（pow（22*m）-1）;
??if（fabs（T[m-1][1]-T[m][1]）????{
?????printf（“pi的估計值為：%lf\n“T[m][1]/2）;
?????break;
	?if（?k?!=?1?）?exit（0）;?}
????}
}

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1056??2009-05-24?14:44??Text1.c

-----------?---------??----------?-----??----

?????????????????1056????????????????????1