I have to get my math skills back up to scratch, so I'm trying to (re)derive Newton's solution to the two body problem.
I managed to get:
x'' = -GMx / ||x||³
with '' being the second time derivative, x being the vector distance between the two masses as a function of time, and M being the sum of the masses.
Which is what I expected to get, but now I have to find x. As far as I know, this is done by guessing forms for x and finding one that fits, and the important clue is that x'' is dependent on -x.
So I tried x(t) = K.exp(ivt+b), with K, v and b are unknown constant vectors and i² = -1, because I think this will have enough freedom for circular, elliptical, parabolic and hyperbolic orbits.
x'' = -v².K.exp(ivt+b) = -v².x (with v² being v with each element squared; not v dot v. Similarly, the full stops are per element multiplication, not vector multiplication)
-v².x = -GMx / ||x||³
So either x = 0, or
-v² = -GM / ||x||³
||x||³ = GM/ v² (So the distance between the masses remains constant, and v must have the same value for all elements)
I think that bringing complex numbers in might have been a mistake, I'll try ordinary trig functions next.
P.S. I mostly wrote this thing so that I could analyze my thought processes, but am posting it in case someone else enjoys reading it.