x=[1.0,1.5,2.0,2.5,3.0,3.5,4.0]; y=[1.1,1.3,1.6,2.0,2.7,3.4,4.1]; p=polyfit(x,y,2); y_fit=polyval(p,x); plot(x,y,'-bo',x,y_fit,'-md') disp(p)
