double precision d_ssxm(1:3),d_ljxm(1:3),d_ssm(1:3),d_ljm(1:3)
double precision f1,f2,h1,h2,hd1,hd2
double precision omega,delta_inv,deltasq_inv,fac1,fac2
+ double precision xi,yi,zi
c-------FIRST METHOD
double precision xm,d_xm(1:3)
c-------END FIRST METHOD