var('t,s')
C = parametric_plot3d([sin(t), cos(t), (sin(t))^2], (t,0,2*pi), color= "blue", thickness = 4)
Cyl = parametric_plot3d([cos(t), sin(t), s], (t,0,2*pi), (s, -1, 2), color="green", opacity = 0.25)
Cyl2 = parametric_plot3d([t,s,t^2],(t,-1,1),(s,-1.5,1.5), color="red", opacity = 0.25)
show(C+Cyl + Cyl2)