clc;
clear;
close;
n=0;
for x = -0.5:0.01:0.5,
n=n+1;
m=0;
for y = -0.5:0.01:0.5,
m=m+1;
z(n,m) = abs(sqrt(1 - x.^2 - y.^2));
end
end
surf(z), title('3-D Plot'), grid
hold on
contour(z,30)
axis('square');
grid on