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