In the two phase region this function returns the interpolated heat capacity between the liquid and vapour state heat capacities.