A classic. This
code
contains (a) the physics engine to solve the two body problem
using the universal Kepler equation and Stumpff functions - also
a nice example of converting from a state vector to orbital
elements and (b) the lick mongo subroutine graphics
calls to make the plot and movie below. See the README in the
code.
Two equal mass bodies about the center of mass.
Same as the pdf to the left, but in the time domain.
See the circular and elliptical restricted three body problem
here.
Please cite the relevant references if you publish a piece of work
that use these codes, pieces of these codes, or modified versions of them.
Offer co-authorship as appropriate.