Numerical solution of ordinary differential equations - Atkinson