clear; clc; close all;

[X, Y] = meshgrid(linspace(-10, 10, 51));
R = sqrt(X .^ 2 + Y .^ 2);
R(R == 0) = eps; % avoid the singularity
Z = sin(R) ./ R;
surf(X, Y, Z);