This is just a main method for PercusYevick.java. It creates a discretization of the Lennard-Jones Mayer function, which PercusYevick then employs to compute Percus-Yevick
approximations of the virial coefficients up to mth order. The second and third coefficients are fully accurate (the
PY approximation is exact).
Use the setRoute() method in PercusYevick to select whether the compressibility-route or virial-route approximation is used.