function y = Cp_water(T)

% Script for calculation of Cp(T) for water at 1 MPa and T in the range
% of 0.04 - 179.91 C.

% Parameters
a = 4.211744250047330;
b = -0.001332574855572;
c = 0.000013326290865;

y = a + (b*T)+ (c*T*T);