This paper describes a package of computer programs for the unified treatment of initial-value problems for systems of ordinary differential equations. The programs implement a numerical method which is efficient for a general class of differential equations. The user may determine the solutions over finite or infinite intervals. The solutions may have singularities at the end-points of the interval for which the solution is sought. Besides giving the initial values and the analytical expression for the differential equations to be solved the user needs to specify the nature of the singularities and give some other analytical information as described in the paper in order to take advantage of the speed and accuracy of the package described here.
['Frank Stenger', 'Sven-Åke Gustafson', 'B. M. Keyes', "M. O'Reilly", 'K. Parker']
ODE-IVP-PACK via Sinc indefinite integration and Newton's method