AR

Lennard-Jones potential with $\sigma$ appropriate for Ar, i.e. $3.4\,$Å. Hence, input is expected in Angstrom and energies are calculated for $\epsilon=1$. Subroutines used: ljdiff.