The three-body problem only states the problem to solve, it doesn’t itself state anything about the existence or non-existence of solutions. It has been proven that there is no general closed-form solution. And there are obvious solutions for trivial special cases, such as three equal masses in an equilateral triangle rotating around each other.

Further reading: https://en.wikipedia.org/wiki/Three-body_problem#Solutions

There is no closed-form solution for finding the roots of >4th degree polynomials in general, but that doesn’t preclude many families of >4th degree polynomials from having closed-form solutions. As a trivial example, x^5 - 1. The exact same thing with the three-body problem.

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.

In the case of a planet going around the sun we know that the planet travels in an ellipse, more generally a conic section including interstellar comets. Orbital periods and everything else can be computed straightforwardly with formulas.

In the three body problem you can always do a numerical integration (e.g. simulation) and this is valid for a certain amount of time but will not be valid forever because of: (1) chaotic motion which amplifies small errors exponentially over time and (2) celestial mechanics is symplectic which means it conserves certain geometric properties and most integrators are not symplectic and must have different long term dynamics. There are symplectic integrators but they don’t work as well overall as ordinary integrators.

We do not know, for instance, if the solar system is stable. In the short term the planets seem to be basically doing their own thing in their orbits with just minor perturbations. We know the orbits vary a bit over millions of years. We aren’t sure which side of the sun Pluto will be on in 30 million years. It’s very believable that the planets are doing the same thing in 4 billion years but we can’t rule out that the orbits could change in a big way or one could get ejected.

I'm pretty sure there is always a unique solution to the equations of motions (safe for some pathological edge cases perhaps). Classical mechanics is deterministic, after all. But for more than two bodies, there is in general no solution in closed form, and it's often chaotic, so not even computeable for arbitrary time frames.

The "about" info states that all of these are computed numerically.

Yes. It's kind of like the halting problem: You cannot write a general computer program that will analyze the source code of any random other computer program and tell you if it will halt.

You can write a program that will analyze the code of a few specific other programs and tell you if they will halt. You just can't do it in general.

The 3-body problem is like that. Except it's much harder to find stable 3-body problems than computer programs that are predictable.

That's exactly the case

[dead]