Guiding Center 4d Barely Passing with Radau IIA VPRK 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.