%gradient_example.m
%計算二維高斯函數的梯度場
v=-2:0.25:2;                                 
[x,y]=meshgrid(v,v);                    %產生自變量x,y
z= exp(-(x.^2+y.^2+0.5*x.*y));          %二維高斯函數
[px py]=gradient(z,0.25);               %梯度場
contour(v,v,z,4);                       %畫4條等高線
hold on;                               
quiver(v,v,px,py);                      %畫出梯度場