qumian:=plot3d([a*sin(phi)*cos(theta),a*sin(phi)*sin(theta),a*cos(phi)],phi=0..Pi,theta=0..2*Pi,color=blue,style=wireframe): 整个球面
qumian2:=plot3d([r*cos(theta),r*sin(theta),sqrt(a^2-r^2)],theta=-Pi/2..Pi/2,r=0..a*cos(theta),color=green,style=patchnogrid): 含在圆柱面内的球面（上半部分）
qumian3:=plot3d([r*cos(theta),r*sin(theta),-sqrt(a^2-r^2)],theta=-Pi/2..Pi/2,r=0..a*cos(theta),color=green,style=patchnogrid): 含在圆柱面内的球面（下半部分）
zhumian:=plot3d([a*(cos(theta))^2,a*sin(theta)*cos(theta),z],theta=-Pi/2..Pi/2,z=-a*abs(sin(theta))..a*abs(sin(theta)),color=yellow): 含在球面内的圆柱面
quxian:=spacecurve([a*(cos(theta))^2,a*sin(theta)*cos(theta),a*sin(theta)],theta=0..2*Pi,color=red,thickness=5):
维维安尼曲线
x_axis:=plot3d([u,0,0],u=-1.2*a..1.5*a,v=0..0.01,thickness=3): （x轴）
y_axis:=plot3d([0,u,0],u=-1.2*a..1.5*a,v=0..0.01,thickness=3): （y轴）
z_axis:=plot3d([0,0,u],u=-1.2*a..1.5*a,v=0..0.01,thickness=3): （z轴）
display(qumian,qumian2,qumian3,zhumian,quxian,x_axis,y_axis,z_axis,scaling=constrained);
评论