%例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)]'