function Cp = Cp_CO2(T)
Cp = (3.259 + 1.356e-3*T + 1.502e-5*T^2 + -2.374e-8*T^3 + 1.056e-11*T^4)*8.3145;
end