x=[0 1 2 3 4 5]; 

 y=[0 20 60 68 77 110]; 

 newx=0:0.05:5; % 以新陣列形成更小增量以利回歸計算及繪圖 

for n=2:9 

f(:,n)=polyval(polyfit(x,y,n),newx)'; 

 plot(newx,f(:,n),x,y,'o') 

title(['Poly. regression, deg=',int2str(n)]) 

 xlabel('Time'), ylabel('Temp'), grid 

 pause % 每次要暫停看清楚圖再執行下一步驟 

 end