Error correction: The reference in the function enthalpyOfWater and its derivative was 0°C rather than 20°C. That was causing the error in the T_phX function