returns specifc internal energy u from EoS
function u input HelmholtzDerivs f; output SpecificEnergy u; end u;