We discuss numerical solution of a system of nonlinear ordinary differential equations with a time-varying delay. This system describes the mutual influence between mass bodies in the case of a finite speed of the gravitational interaction. We show that this fact has to be taken into account in...