例87用不同的方法求函数(x)的数值导数,并在同一个坐 标系中做出(x)的图像。 程序如下: f= inline('qr(x^3+2*x^2-x+12)+(x+5).^(1/6)+5x+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,(x),5);%用5次多项式p拟合x) dp=polder(p); %对拟合多项式p求导数dp dpx=polyval(dp, x); %求d在假设点的函数值 dx=df(o30)0.01;%直接对f(x)求数值导数 gX=g(x); %求函数导函数g在假设点的导数 plot(x,dpx,x,dx,,x,gx,'");%作图例8-7 用不同的方法求函数f(x)的数值导数，并在同一个坐 标系中做出f'(x)的图像。 程序如下： 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,'-'); %作图