Numerical methods for ordinary differential equations by C. Vuik