In this study, the step size strategies are obtained such that the local error is smaller than the desired error level in the numerical integration of a type of nonlinear equation system in interval [t 0,T]. The algorithms are given for calculating step sizes and numerical solutions according to these strategies. The algorithms are supported by the numerical examples. © 2011 Elsevier B.V. All rights reserved.