Explicit algorithmic regularization in the few-body problem for velocity-dependent perturbations

A new algorithm is presented for the numerical integration of second-order ordinary differential equations with perturbations that depend on the first derivative of the dependent variables with respect to the independent variable; it is especially designed for few-body problems with velocity-dependent perturbations. The algorithm can be used within extrapolation methods for enhanced accuracy, and it is fully explicit, which makes it a competitive alternative to standard discretization methods.