var('x,y,u,v')
L = plot3d(lambda x,y: 2, (-5,5), (-5,5), color="lightblue", opacity=0.8)
P = plot3d(4 - x^3 - y^2, (-2,2), (-2,2), color='green')
Q = plot3d(x^3 + y^2 - 4, (-2,2), (-2,2), color='orange')
I = icosahedron(color='red').translate(0,0,6).scale(2)
M = parametric_plot3d([cos(u)*(1+v*cos(u/2)), sin(u)*(1+v*cos(u/2)), 0.2*v*sin(u/2)],
(u,0, 4*pi+0.5), (v,0, 0.3),plot_points=[50,50],opacity=0.7).scale(6)
L + P + Q + I + M