.GNU_ScientificLibrary.Examples.specfunc.CopperSpecHeat

But the main feature here -- namely, the temperature-dependent heat capacity -- is still thereby revealed. Note the c_v~T³-behavior for low temperatures, accurately predicted in the Debye model (for high-T, this correctly turns over to an approach to a constant). Note also the relatively small contribution from the conduction electrons (well approximated by a gas of degenerate fermions): 

See also the Physics.DebyeEoS and Physics.FermiDiracEoS blocks.