syms x y
z=(x^2-2*x)*exp(-x^2-y^2-x*y);
zx=simple(diff(z,x))

zy=diff(z,y)

[x,y]=meshgrid(-3:.2:3,-2:.2:2);
z=(x.^2-2*x).*exp(-x.^2-y.^2-x.*y);
surf(x,y,z), axis([-3 3 -2 2 -0.7 1.5]) % 直接繪制三維曲面

contour(x,y,z,30), hold on   % 繪制等值線
zx=-exp(-x.^2-y.^2-x.*y).*(-2*x+2+2*x.^3+x.^2.*y-4*x.^2-2*x.*y);
zy=-x.*(x-2).*(2*y+x).*exp(-x.^2-y.^2-x.*y);    % 偏導的數值解
quiver(x,y,zx,zy)  % 繪制引力線