typedef float real;
typedef real rvec[3];
void calc_dx2dx(real *, real *);
void phi_sr(int nj,rvec x[],int k)
{
  int i,j;
  for(i=0; (i<nj-1); i++)
    for(j=i+1; (j<nj); j++)
      if (k)
	calc_dx2dx(x[i],x[j]);
}