[X,Y]=meshgrid(-2:.2:2, -2:.2:2);
Z=X .* exp(-X.^2 - Y.^2);surf(X,Y,Z);
xlabel('X')
ylabel('Y')
zlabel('z')
title('Z(X,Y)=X*exp(-X^2-Y^2)')
