In order to predict astronomical events, it is necessary to know the positions of planets and moons at particular points in time. Johannes Kepler determined that orbits are elliptical through observation and geometry. Sir Isaac Newton refined the definition of orbits by describing gravity as an attractive force between masses. A description of Kepler's laws cat be found here.
Actually, the orbits of planets and moons are not true ellipses. This is because the gravitational effects of the planets are constantly changing the shape of the orbits of the other planets. To get more accurate predictions better models are required.
The VSOP87 data and ELP2000 data provide means of calculating the positions of the planets and the Moon respectively.
The NASA JPL Ephemerides use numerical integration to generate accurate positions for the Moon and planets from 1550 to 2050. There is a New Horizons web site. There is also ASCII data which take the form of Chebyshev polynomials for expert users. The data sets are often refered to as DE4XX. I started with DE430. These data sets get refined using space probe data. The latest data set is DE438. The computations required to use DE4XX data can be found here.
The orbits of Neptune and Pluto cross each other. When Pluto is close to perihelion, it is inside the orbit of Neptune. If Pluto got too close to Neptune it would probably be ejected from the solar system. This currently can't happen. The reason for this being that Pluto and Neptune are in a 2:3 resonance. This means that for every 2 orbits of Pluto, Neptune completes 3 orbits. This means that when Pluto is inside Neptune's orbit, Neptune is always distant on average a quarter of an orbit away
Jupiter's moons Io, Europa and Ganymede are in a 1:2:4 resonance.
When one or more bodies are in orbit around a planet or a star, they don't orbit around the centre of the parent body unless it is significantly more massive than the orbiting bodies. In fact all of the bodies in the system orbit around the centre of mass of the system. This is called the barycentre. The barycentre is usually in constant motion due to the fact that the centre of mass of a system constantly changes due to orbital motion.
The Earth doesn't orbit around the Sun. Rather the Sun and all of the planets orbit about the centre of mass of the Solar System which is called the Solar System Barycentre (SSB). The graph shows the position of the barycentre relative to the centre of the Sun using NASA DE430 data from 1938 to 2052. It shows just how much the SSB moves as the relative positions of the planets change.
The units on the axes are one solar radius which is 696,000km.
The graph shows the path of the barycentre in the plane of the ecliptic and perpendicular to the plane of the ecliptic. The planets' orbits are a few degrees inclined to the ecliptic which means that the barycentre actually moves in three dimensions.
In 1951 the EMB was close to the centre of the Sun. In 2022 it will be orbiting around a point some 600,000km above the surface of the Sun, around 1,000,000km from the centre of the Sun. All of this is mainly due to Jupiter and Saturn.
In 1982 all of the planets were aligned on the same side of the Sun. It pulled the barycentre way off centre as seen in the top right quadrant of the graph.