$G$ $\epsilon$ $s$
Cząstka 1 $M_1$ $r_{01}$ $\phi_{01}$ $v_{01}$ $\alpha_{01}$
Cząstka 2 $M_2$ $r_{02}$ $\phi_{02}$ $v_{02}$ $\alpha_{02}$
Cząstka 3 $M_3$ $r_{03}$ $\phi_{03}$ $v_{03}$ $\alpha_{03}$
10x wolniej 5x wolniej 1x normalny 5x szybciej 10x szybciej

Równania ruchu dla $i$-tej cząstki

$\textbf{a}_i \approx G \sum\limits_{j = 1}^{2} \frac{M_j \textbf{r}_{ij}}{(||\textbf{r}||^2 + \epsilon^2)^{3/2}}$

Stała $\epsilon^2 > 0$ jest pewnym czynnikiem zmiękczającym. W przypadku jego braku siła między dwoma zbliżającymi się obiektami rosła by bez ograniczeń co w konsekwencji nikorzystnie wpłynęłoby na obliczenia numeryczne.