Poincaré sections with solve_ivp: is a star's orbit regular or chaotic?

Afterwards you can compute a Poincaré section with solve_ivp events, launch orbits on an energy surface, and measure how much of the section is chaotic.