clear;
[x,y]=meshgrid([-2:0.1:2]);
z=x.*exp(-x.^2-y.^2);
contour3(x,y,z,30);
xlabel('x');ylabel('y');
zlabel('x*exp(-x^2-y^2');
axis('equal');
view(3);
