============================spa
clear allblog
close allget
[x1,y1,z1]=meshgrid(-26:1:26,-100:5:100,-100:5:100);it
v1=(25*y1-25*sqrt(26^2-x1.^2)+8.5)./(sqrt(26^2-x1.^2)-7.25)-z1;class
[x2,y2,z2]=meshgrid(-26:1:26,-50:5:70,0:5:60);grid
v2=sqrt(3600-120*sqrt(26^2-x2.^2)+2*sqrt(26^2-x2.^2).*y2-y2.^2)-z2;程序
v2(abs(imag(v2))>eps)=0;im
p1=patch(isosurface(x1,y1,z1,v1,0));hold onnw
p2=patch(isosurface(x2,y2,z2,v2,0));hold offimg
set(p1,'FaceColor','r','FaceAlpha',0.3,'EdgeColor','r');
set(p2,'FaceColor','b','FaceAlpha',0.3,'EdgeColor','b');
%daspect([1 1 1])
view(3);
axis([-30 30 -100 100 -100 100]);
grid on;
xlabel('x');ylabel('y');zlabel('z');