Error in the angle coordinates as a function of integration time for the same orbit and the same action–angle methods as in Figure 15. The error is computed as Δθ = θ(t) − θ(t = 0) − Ω t, where the frequency Ω is calculated as the mean frequency over the orbit. The errors on the frequencies for both methods are ≈2 × 10⁻⁹.