Added calculation and plotting of heat capacity polynomials