[x,y] = meshgrid(-3:0.1:3,-2:0.1:2);
z=(x.^2-2*x).*exp(-x.^2-y.^2-x.*y); mesh(x,y,z)

axis([-3,3,-2,2,-0.7,1.5])

ezmesh('(x^2-2*x)*exp(-x^2-y^2-x*y)')
figure; ezsurf('(x^2-2*x)*exp(-x^2-y^2-x*y)')