returns internal energy derivative (du/dT)@d=const
function duTd input EoS.HelmholtzDerivs f; output DerEnergyByTemperature duTd; end duTd;