Guiding Center 4d Deeply Trapped with Gauss-Legendre Runge-Kutta Methods

The time step and number of time steps are those of this orbit, set in src/guiding-center-4d.jl; the toroidal momentum is shown alongside the energy, as the second invariant of the guiding centre dynamics.