function f ( r ) implicit double precision (a-h,o-z) e = 1.d0 s = 1.d0 div = (s/r)**6 f = 4.d0*e*(div**2-div) end