d = dsat / [1 - c×ln{(b + p) / (b + psat)}] ; where
b = åBi+1×t i,
c = åCi+1×t i, and
t = (T-Tcenter)/100, where Tcenter is a constant parameter.
Summations are from i = 0 to nTerms - 1, where nTerms is the number of B or C terms.
dsat = the density of the saturated liquid and psat = the saturation vapor pressure.