public class VirialH2PYCorrection
Computes corrections to Percus-Yevick approximations of B4 and B5 for a hydrogen pair potential.
Select the compressibility or virial route via boolean variable compressibility.
Use semiclassical boolean to change pair potential to the quadratic Feynman-Hibbs potential.
Use calcApprox boolean to calculate quantities with the approximate potential.
Use subtractApprox boolean to compute the correction to that approximate potential.
If determining which option is most efficient via short calculations to estimate standard error,
maintain a 50-50 split of steps between reference and target during data collection.