function Cp = Cp_O2(T)
Cp = (3.630 + -1.794e-3*T + 0.658e-5*T^2 + -0.601e-8*T^3 + 0.179e-11*T^4)*8.3145;
end