def lennard_jones(r, sigma, epsilon):
    """LJ Potential. https://en.wikipedia.org/wiki/Lennard-Jones_potential"""

    V = 4 * epsilon * (
        (sigma / r)**12 - (sigma / r)**6
    )

    return V