A new multistage technique for approximate analytical solution of nonlinear differential equations

The article introduces a new multistage technique for solving a polynomial system of nonlinear initial and boundary value problems of differential equations. The radius of convergence R of the series solution to the problem is derived a-priorly in terms of the parameters of the polynomial system. Then guided by the convergence-control parameter [Formula: see text] , the domain of the problem is split into subintervals. By stepping out in a multistage manner, corresponding subproblems are defined which are then subsequently solved with conventional Parker-Sochacki method to get a piecewise continuous solution with very high accuracy. The method is applied to SIR epidemic model, stiff differential equation modelling combustion, Lorenz chaotic problem, and the Troesch's boundary value problem. The results obtained showed a remarkable accuracy when compared with Runge-Kutta Method of order 4. The article showcased the proposed method as a simple, yet accurate approximate analytical technique for nonlinear differential equations.