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