[x,y] = meshgrid([-2:0.2:2]);
R = (x.^2+y.^2)+eps;
z=sin(R)./R;
fh = figure('Color','w');
ah = axes('Color', [1 1 1], ...
'XTick',[-2 -1 0 1 2],...
'YTick',[-2 -1 0 1 2], ...
'XGrid','on','YGrid','on', ... 
'ZGrid','on');
sh = surface('XData',x, ...
'YData',y,'ZData',z,...
'FaceColor', [1 1 1], ... 
'EdgeColor','k','LineWidth',2);
view(3)












