matlab画直交圆柱面,三个直交圆柱面所围成的立体

with(plots): R:=1:

zhumian1:=plot3d([R*cos(t),R*sin(t),z],z=-2*R..2*R, t=0..2*Pi,color=yellow):

zhumian2:=plot3d([R*cos(t),y,R*sin(t)],y=-2*R..2*R, t=0..2*Pi,color=green):

zhumian3:=plot3d([x,R*cos(t),R*sin(t)],x=-2*R..2*R, t=0..2*Pi,color=red):

display(zhumian1,zhumian2,zhumian3,scaling=constrained,orientation=[40,65]);

三个圆柱面所围成的立体：

with(plots): R:=1:

zhumian_a1:=plot3d([R*cos(t),R*sin(t),z],z=-R*sin(t)..R*sin(t),t=-Pi/4..Pi/4,color=yellow):

zhumian_a2:=plot3d([R*cos(t),R*sin(t),z],z=-R*cos(t)..R*cos(