x=-8:.5:8;
y=x';
X=ones(size(y))*x;
Y=y*ones(size(x));
R=sqrt(X.^2+Y.^2)+eps;
z=sin(R)./R;
figure(gcf);
mesh(z);