Using Runge-Kutta-Fehlberg 7(8) method for solving numerical problems in astro-dynamics, for example, two-body and three-body problems.

Modify the parameters at the beginning of twobody.cpp, and type following command in your *nix terminal:

make
make run_twobody
make plot_twobody

Solve the differential equation required by central configuration

make
make run_centconf
make plot_centconf

Plot Hill surface using python