Define the Lennard Jones potential function

LJ[r_, eps_, sig_]:= 4eps ((sig/r)^12 - (sig/r)^6)

LJ[0.4, 0.6502, 0.3166]