%% Data
T = 50:50:700;
CpCO = [0.02918 0.02942 0.02968 0.02996 0.03024 0.03053 0.03083 0.03114 ...
    0.03145 0.03176 0.03208 0.03239 0.03271 0.03301];
CpH2 = [0.02889 0.02895 0.02903 0.02912 0.02922 0.02934 0.02947 0.02962 ...
    0.02977 0.02994 0.03011 0.03029 0.03049 0.03069];
CpCH3OH = [0.04703  0.05104 0.05493 0.05872 0.06239 0.06593 0.06935 ...
    0.07263 0.07263 0.07876 0.08159 0.08427 0.08678 0.08912];
%% Fit and Estimations
CpCOOrder3 = polyfit(T,CpCO,3);
CpH2Order3 = polyfit(T,CpH2,3);
CpCH3OHOrder3 = polyfit(T,CpCH3OH,3);
EstCpCOOrder3 = polyval(CpCOOrder3,T);
EstCpH2Order3 = polyval(CpH2Order3,T);
EstCpCH3OHOrder3 = polyval(CpCH3OHOrder3,T);

𝐶𝑂 = 𝑏𝐶𝑂(𝑇 − 25) + 𝑚1𝐶𝑂(𝑇2 − 252)/2 + 𝑚2𝐶𝑂(𝑇3 − 253)/3 + 𝑚3𝐶𝑂(𝑇4 − 254)/4