以下是我使用plot_surface、meshgrid和三角函数的一些尝试：

import numpy as np
import matplotlib.pyplot as plt

a2, b2 = 11.39, 6.26

X = np.linspace(-np.sqrt(a2), np.sqrt(a2), 100)
Theta = np.linspace(0, 2.1*np.pi, 1000)
X, Theta = np.meshgrid(X, Theta)

Y0 = np.sqrt(b2 * (1 - X**2 / a2))
Y = Y0 * np.cos(Theta)
Z = Y0 * np.sin(Theta)

fig,ax = plt.subplots(subplot_kw={'projection':'3d'})
ax.set_box_aspect((np.ptp(X), np.ptp(Y), np.ptp(Z)))
ax.plot_surface(X, Y, Z)
ax.set_xlabel('x')
plt.show()