• Purpose
• Contents
  • Lorenz attractor
  • Two-body problem
  • Three-body problems
    • Restricted three-body problem
      • PCR3BP
      • Poincaré Section
      • Sitnikov problem
  • Gaussian perturbation
  • Laplace’s parameter
  • Spin-orbit resonance

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