Llegue tres años y medio tarde, pero si alguna otra persona necesita hacer lo mismo aquí les dejo el código que desarrolle, primeramente uso coordenadas polares y después lo pase a cartesianas; en mi caso, necesito generar el cilindro completo, si solo requieren la tapa inferior basta con eliminar la lista `[X_t, Y_t, Z1_t]` en el `return`. Notas: 1. He desarrollado una función, la cual retorna una lista de listas; las listas interiores contienen los parámetros que necesita `plot_surface()` para graficar, primero se debe crear el objeto donde se graficara y después mediante un bucle se van anexando los gráficos de la lista retornada. 2. El cilindro es paralelo al eje Z 3. `x0, y0` son el centro del cilindro 4. `r` es el radio del cilindro 5. `z0, z1` definen la altura del cilindro ``` def cyl(x0, y0, r, z0, z1): edge = 11 # Número de vértices->caras-1 para graficar la superficie la superficie #-------Radio, ángulo phi y altura de la ecuación en cilíndricas; e: ecuación r_e = np.linspace(0, r, edge) p_e = np.linspace(0, 2*np.pi, edge) z_e = np.linspace(z0, z1, edge) #-------Mallados para las tapas y el cilindro; t: tapa, c: cilindro R_t, P_t = np.meshgrid(r_e, p_e) R_t, Z_c = np.meshgrid(r_e, z_e) #-------Coordenadas cartesianas del cilindro para graficar X_c = r*np.cos(p_e)+x0 Y_c = r*np.sin(p_e)+y0 #-------Asignamos la altura de los puntos a graficar en Z z0_lambda = lambda x: z0*(x+1)/(x+1) # x+1 para que no dé problemas con la evaluación en 0/0 z1_lambda = lambda x: z1*(x+1)/(x+1) #-------Puntos en Z Z0_t = z0_lambda(R_t) Z1_t = z1_lambda(R_t) #-------Pasamos el grid a cartesianas para las tapas X_t, Y_t = R_t*np.cos(P_t)+x0, R_t*np.sin(P_t)+y0 return [[X_c, Y_c, Z_c], [X_t, Y_t, Z0_t], [X_t, Y_t, Z1_t]] ``` La forma en que lo gráfico es la siguiente ``` fig = plt.figure() ax = fig.add_subplot(111, projection='3d') cilindro = cyl(3,3,4,-6,2) for i in range(len(cilindro)): ax.plot_surface(cilindro[i][0], cilindro[i][1], cilindro[i][2], color='blue') plt.show() ``` Esto genera un cilindro centrado en (3,3) de radio 3, base en -6 y altura en 2, como en el siguiente gráfico: [![Cilindro completo][1]][1] Si le quitamos la lista `[X_t, Y_t, Z1_t]` en el `return` (que señale en el primer párrafo), se obtiene: [![Cilindro sin tapa superior][2]][2] [1]: https://i.sstatic.net/LZQaI.png [2]: https://i.sstatic.net/pbeVT.png Espero les sea de utilidad.