Calculating the classical virial expansion using automated algebra

Using schematic model potentials, we calculate exactly the virial coefficients of a classical gas up to sixth order and use them to calculate the virial expansion of basic thermodynamic quantities such as pressure, density, and compressibility. At sufficiently strong couplings, as expected, the virial expansion fails to converge. However, at least for the interactions and parameter ranges we explored, we find that Padé-Borel resummation methods are very effective in improving the convergence of the expansion.