solve Redlich-Kwong-Soave for density
function density_pT_Soave input Temperature T; input AbsolutePressure p; input AbsolutePressure psat = Ancillary.saturationPressure_T(T = T); output Density d; end density_pT_Soave;