30 . PXI, PYI ,PZI ,PXJ, PYJ, PZJ, HIJ ,
31 . FXX, FYY ,FZZ ,FXY, FYZ, FXZ)
35#include "implicit_f.inc"
42 . pxi(*), pxj(*),pyi(*), pyj(*), pzi(*), pzj(*),
43 . fxx(*),fyy(*),fzz(*),fxy(*),fxz(*),fyz(*),hij(*)
51 hij(i)=fxx(i)*pxi(i)*pxj(i)+fyy(i)*pyi(i)*pyj(i)+
53 . fxy(i)*(pxi(i)*pyj(i)+pyi(i)*pxj(i))+
54 . fyz(i)*(pyi(i)*pzj(i)+pzi(i)*pyj(i))+
55 . fxz(i)*(pxi(i)*pzj(i)+pzi(i)*pxj(i))
subroutine s8zkeh3(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, hij, fxx, fyy, fzz, fxy, fyz, fxz)