O movimento dos corpos é baseado na segunda lei de Newton e na sua lei da gravitação universal.
A força resultante no corpo de massa $m_i$ é dada por
$$
\vec{F}_i = \sum_{j \neq i}^N \dfrac{Gm_im_j}{r_{ij}^2}\hat{r}_{ij},
$$
onde
- $N$: número de corpos
- $G$: constante da gravitação universal
O sistema de EDO's é resolvido pelo método de Runge-Kutta de quarta ordem.
Outras considerações:
- O centro de massa do sistema está mantido fixo.
- As velocidades dos corpos bem como suas posições são relativas ao centro de massa do sistema.
Leia mais: