x1=-4:0.3:4;
y1=-4:0.3:4;
[x,y]=meshgrid(x1,y1);
surf(x,y,sin(x).*cos(y).*exp(-(x.^2+y.^2).^0.5))