x=-10:10; y=-10:10; [X,Y]=meshgrid(x,y); Z=(sin(sqrt(X.^2+Y.^2))+eps)./(sqrt(X.^2+Y.^2)+eps); mesh(X,Y,Z)