iterative calculation of saturation properties from EoS with Newton-Raphson algorithm
function setSat_d input Density d; output SaturationProperties sat; end setSat_d;