a=1x=-1:0.1:1;y=-1:0.1:1;[Y,X]=meshgrid(x,y);Z=arrayfun(@(x,y) fzero(@(z) x.^2+y.^2+(z-a).^2-a.^2,1),X,Y);
surf(X,Y,Z)
