Using Runge-Kutta-Fehlberg 7(8) method for solving numerical problems in celestial dynamics, for example, two-body and three-body problems, mainly on C++ with some python scripts for plotting.
Modify the parameters at the beginning of twobody.cpp
, and
type following command in your *nix terminal:
make
make run_twobody
make plot_twobody
Planar circular restricted 3-body problem
make
make run_gaussperturb
make plot_gaussperturb