Integrators#
Why physicskit.classical defaults every conservative system to a symplectic
integrator (RK4’s energy drifts monotonically, Verlet/Yoshida4 stay
bounded), and automated timestep selection with
estimate_dt().
Why symplectic integration: RK4 leaks energy, Yoshida4 does not
Why symplectic integration: RK4 leaks energy, Yoshida4 does not
