求助，请进来看一个简单的MPI程序

是《MPI与OpenMP并行程序设计－C语言版》的作业；91页里面习题4.11：

通过计算在曲线下图形的面积来计算pi，一个简单的方法是使用矩形规则，将区间【0，1】分成k个等长子区间。用这些高度可以构造k个矩形。这些矩形的面积逼近曲线下图形的面积。

如下是课文中的C程序 

#define  intervals 1000000

int  main( int  argc,  char   * argv[])
... {
   double area;
   double ysum;
   double xi;
   int i;   
   ysum=0.0;
   for(i=0;i<intervals;i++)
   ...{
      xi=(1.0/intervals)*(i+0.5);
      ysum+=4.0/(1.0+xi*xi);
   }
   area=ysum*(1.0/intervals);
   printf("Area is %f ",area); 
   return 0;
}

 现要求编写一个使用矩形规则在100000个区间的并行程序计算pi；

我的程序如下，

#include  < mpi.h >
#include  < stdio.h >
#include  < math.h >
                   
#define  intervals 1000000
double  area,aa;
int  main( int  argc,  char   * argv[])
{   int  id;                                /* Process rank */
    int  p;                                 /* Number of Processes */
    int  i;
    void  pi( int , int );

   MPI_Init( & argc,  & argv);
   MPI_Comm_rank(MPI_COMM_WORLD, & id);
   MPI_Comm_size(MPI_COMM_WORLD, & p);

    for  (i = id;i < intervals;i += p)
      {pi(id,i);
      aa += area;}
      printf( " Area is %f " ,aa);
      MPI_Finalize();
       return   0 ;
}

void  pi( int  id, int  z){
    double  ysum;
    double  xi;
   ysum = 0.0 ;
   xi = ( 1.0 / intervals) * (z + 0.5 );
   ysum = 4.0 / ( 1.0 + xi * xi);
   area = ysum * ( 1.0 / intervals);
   }

编译之后执行不了，菜鸟的sb程序不知道是哪里出错了呢、