gusucode.com > 《精通matlab 7》包括配套光盘 > 配套光盘里的程序附件/Ex-35/Ex3512.m
%例35-12:求方向导数。 syms x y z f=x*y^2-x*y*z+z^3; s=jacobian(f) sx=subs(s,'x','1'); sy=subs(sx,'y','1'); sz=subs(sy,'z','2'); g=vpa(sz) a=pi/3; b=pi/3; c=pi/3; L=g*[cos(a),cos(b),cos(c)]'