>> [x,y]=meshgrid(-8:0.5:8,-10:0.5:10);
>> R=sqrt(x.^2+y.^2);
>> R=sqrt(x.^2+y.^2)+eps;
>> z=sin(R)./R;
>> surf(x,y,z)
