[szj1123szj]

clear,clc
r=5;
[x,y]=meshgrid(-r-r/100:r/115:r+r/100);
z=(sqrt(r^2-x.^2-y.^2));
n=size(z);
for i=1:n
for j=1:n
if ~isreal(z(i,j)) & imag(z(i,j))>.75
z(i,j)=nan;
z1(i,j)=nan;
else
z(i,j)=real(z(i,j));
z1(i,j)=-real(z(i,j));
end
end
end
mesh(x,y,z)
hold on
%surf(x,y,z)
mesh(x,y,z1)
%surf(x,y,z1)
axis tight
shading interp
alpha(.9)


转载自 作者：cxf208