1  typedef float real;
       2  typedef real rvec[3];
       3  void calc_dx2dx(real *, real *);
       4  void phi_sr(int nj,rvec x[],int k)
       5  {
       6    int i,j;
       7    for(i=0; (i<nj-1); i++)
       8      for(j=i+1; (j<nj); j++)
       9        if (k)
      10  	calc_dx2dx(x[i],x[j]);
      11  }