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.