摘要

In this paper, we propose a new modification of the Adomian decomposition method for solution of higher-order nonlinear initial value problems with variable system coefficients and solutions of systems of coupled nonlinear initial value problems. We consider various algorithms for the Adomian decomposition series and the series of Adomian polynomials to calculate the solutions of canonical first- and second-order nonlinear initial value problems in order to derive a systematic algorithm for the general case of higher-order nonlinear initial value problems and systems of coupled higher-order nonlinear initial value problems. Our new modified recursion scheme is designed to decelerate the Adomian decomposition series so as to always calculate the solution's Taylor expansion series using easy-to-integrate terms. The corresponding nonlinear recurrence relations for the solution coefficients are deduced. Next we consider convergence acceleration and error analysis for the sequence of solution approximations. Multistage decomposition and numeric algorithms are designed and we debut efficient MATHEMATICA routines PSSOL and NSOL that implement our new algorithms. Finally we investigate several expository examples in order to demonstrate the rapid convergence of our new approach.