import numpy as np
import matplotlib.pyplot as plt
from matplotlib import cm

ax = plt.figure().add_subplot(projection='3d')

R = 0.2
pi = np.pi
x1 = np.arange(0,1,0.01)
x2 = np.arange(-0.5,0.5,0.05)
L = len(x1)
X1, X2 = np.meshgrid( x1, x2)
# print(X2)

# definition du cylindre
v1_cyl = R*np.outer( np.cos(2*pi*X1), np.ones(np.size(X2)) )
v2_cyl = R*np.outer( np.sin(2*pi*X1), np.ones(np.size(X2)) )
v3_cyl =   np.outer( np.ones(np.size(X1)), X2 )


surf = ax.plot_surface(v1_cyl, v2_cyl, v3_cyl, cmap=cm.coolwarm,
                       linewidth=0, antialiased=False)

#surf_cyl = ax.plot_wireframe( v1_cyl, v2_cyl, v3_cyl )
plt.show()