[x,y] = meshgrid(-8:0.5:8);
z = sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);
subplot(2,2,1);
mesh(x,y,z);
title('mesh(x,y,z)');
subplot(2,2,2);
meshc(x,y,z);
title('meshc(x,y,z)');
subplot(2,2,3);
meshz(x,y,z)
title('meshz(x,y,z)');
subplot(2,2,4);
surf(x,y,z)
title('surf(x,y,z)');