Poincaré section (R, v_R, z = 0, v_z > 0) of the two orbits displayed in Figure 10 (top panel). The left orbit from Figure 10 is shown in blue and the right orbit is displayed in yellow. The code used to generate this surface of section is given in the bottom panel; this code requires the code displayed at the bottom of Figure 10 to be run first.