returns ideal gas heat capacity from EoS
function cp0 input HelmholtzDerivs f; output SpecificHeatCapacity cp0; end cp0;