Based on the numerical solution of initial problem of ordinary differential equation is applied to nonlinear dynamics analysis of structure, the solution of the state equations for nonlinear dynamics system governed by the equation is discussed. The numerical solutions of which include single-step method such as Runge-Kutta method and multi-step method such as Adams-Bashforth-Moulton's predictor-corrector method, Milne-Simpson's predictor-corrector method and Milne-Hamming's predictor- corrector method. The stability and accuracy of those numerical solutions are dis...