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łaby bez ograniczeń, co w konsekwencji niekorzystnie wpłynęłoby na obliczenia numeryczne.