! ! PES for ^7Li_3 ! subroutine vli3(r,v) implicit none real*8 r(3),v,vli2,vli3_3b ! ! Compute the two body part ! All this is in a.u. ! v = 0d0 v = vli2(r(1))+vli2(r(2))+vli2(r(3)) v = v + vli3_3b(r) return end