Code to fit a linear and cubic polynomial for the Cp data

Objective: clear allclose allclc % Preparing the data cp_data = load(\'data\');linear = cp_data(:,1);cp = cp_data(:,2); % Curve fit % cp = a*l^3+b*l^2+c*l+dco_effs = polyfit(linear,cp,3);predicted_cp = polyval(co_effs,linear); % Compare my curve fit with the original data plot(linear,cp,\'linewidth\',3)hold onplot(linear,predicted_cp,\'linewidth\',3,\'color\',\'r\')xlabel(\'linear…