interval=[-3, 3, -3, 3];
ezsurfc('x.*y.*exp(-x.^2-y.^2)',interval);







