這個比較容易。
clc;clear
x=linspace(-pi,pi,50);
[X,Y]=meshgrid(x,x);
z=(X.^2-2*X).*exp(-X.^2-Y.^2-X.*Y);
contourf(X,Y,z)
axis equal
figure(2)
hx=2*pi/(50-1);hy=hx;
[px,py]=gradient(z,hx,hy);
quiver(X,Y,px,py)