This function calculates the vapor phase density of R1234yf from the state variable p (absolute pressure). It is modelled by cubic splines which are fitted with non-equidistant grid points derived from the fundamental equation of state of Richter et al (2011) and the Maxwell criteria.
It is only valid in the two-phase region (i.e., ptriple ≤ p ≤ pcrit ).