This is a PDF version of a preprint submitted to Springer. The definitive version was published in the Numerical algorithms. The original publication is available at www.springerlink.comThis preprint discusses a method for a numerical solution of a nonlinear fractional differential equation, which can be seen as a generalisation of the Adams–Bashforth–Moulton scheme