& gt& gty=x.^5+x.^4+4*random('normal',0,1,1,length(x));
& gt& gtp=polyfit(x,y,5);
& gt& gtf=polyval(p,x);
& gt& gtplot(x,y,':r ',x,f,' b ');
polyfit和polyval就是這麽做的。如果要在計算方法導論的課程中使用,恐怕不太合適
附圖的壹部分似乎不清楚。