This function uses the package poliastro, which is incompatible with Python 3.13. Poliastro is an undocumented dependency (so there's no issues regarding this other than when calling sdf.plot_orbits()). This function entirely relies on the poliastro package, so I'm not sure what the solution would be.