A self-starting and inherently stable method to solve ordinary differential equation numerically. The numerical solution is a table of the dependent variable and its derivatives over only a limited portion of the range of the independent variable