[c,r2,eqn] = lsqcurvefit_approx(x,y,'linear') does the exact same thing as the syntax above. Only the linear and polynomial fits are true linear least squares fits. The nonlinear fits (power, ...