第壹種：首先用壹個4次多項式p（x）擬合y=f（x），並對p（x）求壹般意義下的導數dp（x），然後求出dp（x）在假設點的值；
第二種：直接求y=f（x）在假設點的數值導數；
第三種：先手工算出y=f（x）導函數，然後求在假設點的值。
程序如下：
f=inline('sqrt(x.^3+2*x.^2-x+12)+7*x.^2+36');
x=-3:0.01:3;
p=polyfit(x,f(x),4);
dp=polyder(p);
dpx=polyval(dp,x);
dx=diff(f([x,3.01]))/0.01;
g=inline('(3*x.^2+4*x-1)./sqrt(x.^3+2*x.^2-x+12)/2+14*x');
gx=g(x);
plot(x,dpx,x,dx,'.',x,gx,'-');