cp = C1 + C2*T + C3*T^-2
function integral_cp_dT input Real T; output Real Int_cp_T; end integral_cp_dT;