1 typedef float real;
2 typedef real rvec[3];
3 void calc_dx2dx(real *, real *);
phi_sr(int nj,rvec x[],int k)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 }
12