圆周率的计算C程序（数值法）梯形积分公式

圆周率的计算C程序(数值法)梯形积分公式，能运行！！！

#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