x = r0 * np.cos(u)
y = r0 * np.sin(u)
z = h0 * np.ones(len(x))
XYZ = np.vstack([x, y, z])
xp, yp, zp = euler_rot(XYZ, phi, theta, psi)
ax.plot(xp, yp, zs=zp)
создание 3D-конуса или диска и постоянное обновление его оси симметрии с помощью matplotlib