MATLAB 半球面

这个题我先用的两个参数画的两个图形，但是我真不会用相同的参数把交线表示出来，问了问老师最后是暴力跑一遍点画的

 

r=linspace(0,1,100);
t=0:pi/50:2*pi;
[R T]=meshgrid(r,t);
x=R.*cos(T);
y=R.*sin(T);
z=sqrt((1-x.*x - y.*y)+eps);
mesh(x,y,z);
axis equal
hold on;
z1=R;
x1=cos(T).^2;
y1=sin(T).*cos(T);
mesh(x1,y1,z1);
shading flat
[x y]=meshgrid(-1:.002:1);
z=sqrt((1-x.*x - y.*y)+eps);
for i=1:1:1000
    for j=1:1:1000
        if((1-x(i,j).*x(i,j) - y(i,j).*y(i,j))>=0)
            if(abs(x(i,j)*x(i,j)+y(i,j)*y(i,j)-x(i,j))<=1e-3)
                plot3(x(i,j),y(i,j),sqrt((1-x(i,j).*x(i,j) - y(i,j).*y(i,j))),'b*');
            end
        end
        end
end