The present plots, tables and programs give the location and estimated strength
in relative units of the three-body resonances between a massless particle and two planets.
The estimated strength Delta rho was calculated following Gallardo (2014, Icarus 231, 273-286).
It depends on the eccentricity,
inclination and the argument of perihelion (w) of the particle's orbit.
For more details see preprint here.
From 0 to 1000 au. Code colors: resonances involving Mercury as the most inner planet = red, Venus=green, Earth=blue, Mars=pink, Jupiter=black, Saturn=red, Uranus = green.
From 2 au to 3.5 au. Blue: two-body resonances, red: three-body resonances. Two-body and three-body resonances are not to the same scale.
Or run the program ATLAS3BR and compute the strengths of all resonances near
the semimajor axis of the particle according
to the values of (e,i,w) of the particle's orbit. This is the best choice. Then compute and look at the time evolution of
the critical angles corresponding to the strongest resonances in the interval.
The semimajor axis of the actual resonances can be something different from this theory because we have not taken into account
the shift in "a" generated by the time evolution of the perihelion and node of the particle's orbit. For larger semimajor
axes we have larger errors in "a".