figure(1)
[X,Y]=meshgrid(-10:0.1:10);
R=sqrt(X.^2+Y.^2)+eps;
Z = X.*sin(R)./R;
figure(1)
imagesc(Z)
colormap(gray)
colorbar

figure(2)
surf(X,Y,Z)
colorbar
shading interp
axis square














