z = [-1:0.1:1];
theta = [0:pi/18:2*pi];
[R, THETA] = meshgrid(z, theta);
X = R.*cos(THETA);
Y = R.*sin(THETA);
Z = sqrt(1-R.^2);

mesh(X, Y, Z);
hold on
mesh(X, Y, -Z);