Same as Fig. 20, except that the integration is performed using the unaveraged equations of the restricted three-body problem including Saturn’s oblateness. The output is then digitally filtered to remove its short-period component. Due to the chaotic divergence of trajectories, we ran several simulations with different migration velocity and chose two of them that resemble the trajectories in Fig. 20.