Guiding Center 4d Barely Passing 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)

Plots

Guiding Center 4d Barely Passing

LobattoIIIA(3)

Plots

Guiding Center 4d Barely Passing

LobattoIIIA(4)

Plots

Guiding Center 4d Barely Passing

LobattoIIIB(2)

Plots

Guiding Center 4d Barely Passing

LobattoIIIB(3)

Plots

Guiding Center 4d Barely Passing

LobattoIIIB(4)

Plots

Guiding Center 4d Barely Passing

LobattoIIIC(2)

Plots

Guiding Center 4d Barely Passing

LobattoIIIC(3)

Plots

Guiding Center 4d Barely Passing

LobattoIIIC(4)

Plots

Guiding Center 4d Barely Passing

LobattoIIID(2)

Plots

Guiding Center 4d Barely Passing

LobattoIIID(3)

Plots

Guiding Center 4d Barely Passing

LobattoIIID(4)

Plots

Guiding Center 4d Barely Passing

LobattoIIIE(2)

Plots

Guiding Center 4d Barely Passing

LobattoIIIE(3)

Plots

Guiding Center 4d Barely Passing

LobattoIIIE(4)

Plots

Guiding Center 4d Barely Passing

LobattoIIIF(2)

Plots

Guiding Center 4d Barely Passing

LobattoIIIF(3)

Plots

Guiding Center 4d Barely Passing

LobattoIIIF(4)

Plots

Guiding Center 4d Barely Passing

LobattoIIIG(2)

Plots

Guiding Center 4d Barely Passing

LobattoIIIG(3)

Plots

Guiding Center 4d Barely Passing

LobattoIIIG(4)

Plots

Guiding Center 4d Barely Passing