f=inline('sqrt(x.^3+2*x.^2-x+12)+(x+5).^(1/6)+5*x+2');
g=inline('(3*x.^2+4*x-1)./sqrt(x.^3+2*x.^2-x+12)/2+1/6./(x+5).^(5/6)+5');
x=-3:0.01:3;
p=polyfit(x,f(x),5);          %用5次多項式p擬合f(x)
dp=polyder(p);                 %對擬合多項式p求導數dp
dpx=polyval(dp,x);            %求dp在假設點的函數值
dx=diff(f([x,3.01]))/0.01;   %直接對f(x)求數值導數
gx=g(x);                         %求函數f的導函數g在假設點的導數
plot(x,dpx,x,dx,'.',x,gx,'-');   %作圖