“…Examples of multi-stage methods that have been developed recently to solve IVPs for chaotic and nonchaotic systems include the, multi-stage homotopy analysis method [2,4,5], piecewise homotopy perturbation methods [9,25], multi-stage Adomian decomposition method [1,20], multi-stage differential transformation method, [3,15,24], multi-stage variational iteration method [14,19]. Because these methods attempt to obtain analytical solutions at each interval they involve timeconsuming and tedious computational operations and if too many small intervals are considered, as may be the case when dealing with highly oscillatory systems, the analytical integration process will be too much to handle even with the use of symbolic scientific software.…”