%var theta, z fx = cos(theta) fy = sin(theta) fz = z gx = cos(theta)*(0.75 + 0.25*cos((z-4.5)*pi)) gy = sin(theta)*(0.75 + 0.25*cos((z-4.5)*pi)) gz = z hx = 0.5*cos(theta) hy = 0.5*sin(theta) hz = z cylinder = parametric_plot3d([fx, fy, fz], (theta, 0, 2*pi), (z, 0, 4.5), color="brown", opacity=.9, mesh=1, spin=10, plot_points = 25) cylindersticker = parametric_plot3d([fx, fy, fz], (theta, 0, pi), (z, 0.75, 3.75), color="blue", opacity=.9, mesh=1, spin=5, plot_points = 25) arc = parametric_plot3d([gx, gy, gz], (theta, 0, 2*pi), (z, 4.5, 5.5), color="brown", opacity=.9, mesh=11, spin=20, plot_points = 12.5) cylinder2 = parametric_plot3d([hx, hy, hz], (theta, 0, 2*pi), (z, 5.5, 9), color="brown", opacity=0.9, mesh=1, spin=10, plot_points = 25) cylinder2sticker = parametric_plot3d([hx, hy, hz], (theta, 0, pi), (z, 6.75, 7.75), color="blue", opacity=0.9, mesh=1, spin=10, plot_points = 25) show(cylinder + cylindersticker + arc + cylinder2 + cylinder2sticker)
3D rendering not yet implemented