Mean orbital eccentricities and inclinations of Miranda obtained from the numerical simulations (only those with initial e₄ = 0.0001 and Q₂ ≤ 250) after the disruption of the 2:1 mean motion resonance, both in the case of 4:2:1 resonant crossing (top) and 3:2:1 resonant crossing (bottom). The values are taken when the system reaches n₂/n₃ = 1.94 (see text in Sect. 4.1). The dashed lines indicate the current values of the mean orbital elements. The red points show the orbital elements of the simulations presented in Figs. 7 and 9.