import math


def pressure_lrc ( density, r_cut ):
    """Calculates long-range correction for Lennard-Jones pressure."""

    
    # density, r_cut, and the results, are in LJ units where sigma = 1, epsilon = 1
    sr3 = 1.0 / r_cut**3 
    return math.pi * ( (32.0/9.0) * sr3**3  - (16.0/3.0) * sr3 ) * density**2