Starting step size for an ODE solver

Abstract One of the more critical issues in solving ordinary differential equations by a step-by-step process occurs in the starting phase. Somehow the procedure must be supplied with an initial step size that is on scale for the problem at hand. It must be small enough to yield a reliable solution by the process, but not so small as to significantly affect the efficiency of solution. In this paper, we discuss an algorithm for obtaining a good starting step size and present a subroutine which can be readily used in most ODE solvers.