Specific heat capacity at constant pressure | turns inifity in two-phase region! | use setState_phX function for input