This function calculates the derivative of liquid density of R134a in the two-phase region with regard to the state variable p (absolute pressure). The non-derivative function is bubbleDensity.
It is only valid in the two-phase region (i.e., ptriple ≤ p ≤ pcrit ).