>>?x=[6,7,8,9,10,11];
>>?y=[0,3.3,6.3,10.9,26.7,108.8];
>>?xx=6:.2:11;
>>?yy=interp1(x,y,xx,'cpline');%插值
>>?p=polyfit(xx,yy,4);%多項式擬合
>>?py=polyval(p,xx);
>>?dp=polyder(p);
>>?dpy=polyval(dp,xx);%多項式求導
>>?plot(x,y,xx,py,xx,dpy)
>>?legend('y-x(原始曲線）','py-xx（擬合曲線）','dpy-xx（求導曲線）')