Solving forward Lorentz-Dirac-like equations

It is shown that successive approximations can be used to implement a numerical method to integrate forward the Lorentz-Dirac equation, as well as other equations with the same singular structure. The code automatically selects the physical solution and avoids the so-called 'runaway solutions'. The method's convergence is analytically discussed in a particular but illustrative case. The convergence is also numerically studied in a capture motion and a chaotic scattering process.