407.sagews

x, y = var('x'),var('y');
W = plot3d(sqrt(1-x^2),(x,-1,1),(y,-1,1), frame=False, color='red');
X = plot3d(-sqrt(1-x^2),(x,-1,1),(y,-1,1), frame=False, color='red'); 
Y = plot3d(sqrt(1-y^2),(x,-1,1),(y,-1,1), frame=False, color='blue'); 
Z = plot3d(-sqrt(1-y^2),(x,-1,1),(y,-1,1), frame=False, color='blue');
show(W+X+Y+Z)

u,v = var('u'),var('v')
E1= parametric_plot3d([u+v-u*v,v,sqrt(1-(u+v-u*v)^2)],(u,0,1),(v,0,1), frame=False, color='red')
E2= parametric_plot3d([u+v-u*v,v,-sqrt(1-(u+v-u*v)^2)],(u,0,1),(v,0,1), frame=False, color='red')
E3= parametric_plot3d([u+v-u*v,-v,sqrt(1-(u+v-u*v)^2)],(u,0,1),(v,0,1), frame=False, color='red')
E4= parametric_plot3d([u+v-u*v,-v,-sqrt(1-(u+v-u*v)^2)],(u,0,1),(v,0,1), frame=False, color='red')
N1 = parametric_plot3d([u*v,v,sqrt(1-(v)^2)],(u,-1,1),(v,0,1), frame=False, color='blue')
N2 = parametric_plot3d([u*v,v,-sqrt(1-(v)^2)],(u,-1,1),(v,0,1), frame=False, color='blue')
W1= parametric_plot3d([-(u+v-u*v),v,sqrt(1-(u+v-u*v)^2)],(u,0,1),(v,0,1), frame=False, color='red')
W2= parametric_plot3d([-(u+v-u*v),v,-sqrt(1-(u+v-u*v)^2)],(u,0,1),(v,0,1), frame=False, color='red')
W3= parametric_plot3d([-(u+v-u*v),-v,sqrt(1-(u+v-u*v)^2)],(u,0,1),(v,0,1), frame=False, color='red')
W4= parametric_plot3d([-(u+v-u*v),-v,-sqrt(1-(u+v-u*v)^2)],(u,0,1),(v,0,1), frame=False, color='red')
S1 = parametric_plot3d([u*v,-v,sqrt(1-(v)^2)],(u,-1,1),(v,0,1), frame=False, color='blue')
S2 = parametric_plot3d([u*v,-v,-sqrt(1-(v)^2)],(u,-1,1),(v,0,1), frame=False, color='blue')
show(E1+E2+E3+E4+N1+N2+W1+W2+W3+W4+S1+S2)