& 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就是這麽做的。如果要在計算方法導論的課程中使用，恐怕不太合適
附圖的壹部分似乎不清楚。