function y = Cp_model(par,T_list)

alfa = par(1);
beta = par(2);
gamma = par(3);

n = length(T_list);
y = zeros(n,1);
for i = 1:n
    y(i) = alfa + (beta*T_list(i)) + (gamma*(T_list(i))^2);
end