График контура продолжает появляться как точка на 0, но когда я вычислял контуры вручную, они были различны и существенны.Не уверен, что с моим кодом или моей математикой что-то не так?:
Re = 6371000
rho = np.linspace(-20*Re, 20*Re, 2000)
z1 = np.linspace(-10*Re, 10*Re, 2000)
z = z1[:, np.newaxis]
Be = 3120*10**(-9)
F = (rho**2)*Be*(Re/np.sqrt(rho**2 + z**2))**3
plt.contour(rho, z1, F)