%Energy loss heating
function Cp = Cp_N2(T)
Cp = (3.539 + -0.261e-3*T + 0.007e-5*T^2 + 0.157e-8*T^3 + -0.099e-11*T^4)*8.3145;
end