Guiding Center 4d Deeply Trapped with Lobatto 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.
LobattoIIIA(2)

LobattoIIIA(3)

LobattoIIIA(4)

LobattoIIIB(2)

LobattoIIIB(3)

LobattoIIIB(4)

LobattoIIIC(2)

LobattoIIIC(3)

LobattoIIIC(4)

LobattoIIID(2)

LobattoIIID(3)

LobattoIIID(4)

LobattoIIIE(2)

LobattoIIIE(3)

LobattoIIIE(4)

LobattoIIIF(2)

LobattoIIIF(3)

LobattoIIIF(4)

LobattoIIIG(2)

LobattoIIIG(3)

LobattoIIIG(4)
