fig = orbit.plot_3d(alpha=0.5, color='k')