[X,Y] = meshgrid(-2:0.2:2);
Z = X.*Y.*exp(-X.^2-Y.^2);
figure(1)
colormap(hot);
pcolor(X,Y,Z)
shading interp
colorbar

figure(2)
colormap(jet);
surf(X,Y,Z)
colorbar










