溫度梯度矢量場的分佈

clear;

[x,y] = meshgrid(-2:.2:2,-2:.2:2);

z = x.*exp(-x.^2-y.^2);

[px,py] = gradient(z,.1,.1);

contour(z);

hold on;

quiver(px,py,1);

hold off;