Solving equation of motion of two massive particle exerting the gravitational force each other I'm trying to analyze the motion of the particles which exert the gravitational force each other. Let $M_1$, $M_2$ be the masses of the particles, and the equation of motion of particle $M_1$
$$
F=G\frac{M_1M_2}{r^2}=M_1\ddot{r} \\
\ddot{r}r^2-GM_2=0
$$
I find it doubtful that the simple motion has difficult equation of motion.
Should I just use some form of successive approximation?
 A: I suggest to see the scenario physically first, and proceed from there.
Here, the gravitational force does not change the distance from the center, but it provides the centripetal force for the circular motion (or any conic section for that matter, based on initial values). 
Assuming it is a circular motion for now, we can just get right into the mathematical description of the motion.
For such systems, masses move in a circle with their COM at the center.

Because there is no other force in system, the COM does not accelerate. It is free to move at a constant velocity though, but this is (usually) irrelevant when analyzing the system.
Because they move in a circle about the COM, just get the distances from their COM, and equate the gravitational force to the centripetal force. This way you have fully described the motion.
Example for $M_2:$
$$
\dfrac{GM_1}{\underbrace{R}_{\text{distance between bodies}}} = \dfrac{{v_2^2}}{\underbrace{R_{2_{cm}}}_{\text{distance of 2 from COM}}}
$$
Similar for $M_1$.
