Comment by btilly

The general case always has a solution. At least until the point where two of the three bodies meet (which is a singularity). We can approximate that solution numerically.

The problem is that the solutions very strongly tend to be chaotic. Meaning that small differences in initial conditions, tend to grow exponentially with time. Which means that if you measure everything to 3 digits of precision, in finite time it will stop looking like the actual solution. Every additional digit of precision adds a similar finite time to how long the approximation is good for.

So when finally found, say, the 1953 BC conjunction described in https://en.wikipedia.org/wiki/Conjunction_%28astronomy%29?#N... - that was a very good stress test for our estimated planetary data. Because surprisingly small errors in modern data would have kept that conjunction from happening.