Convergent regions of the newton homotopy method for nonlinear systems: Theory and computational applications

This paper introduces the concept of the convergent region of a solution of a general nonlinear equation using the Newton homotopy method. The questio n of whether an initial guess converges to the solution of our interest usi ng the Newton homotopy method is investigated. It is shown that convergent regions of the Newton homotopy method are equal to stability regions of a c orresponding Newton dynamic system. A necessary and sufficient condition fo r the adjacency of two solutions using the Newton homotopy method is derive d. An algebraic characterization of a convergent region and its boundary fo r a large class of nonlinear systems is derived. This characterization is e xplicit and computationally feasible. A numerical method to determine the c onvergent region and to establish simple criteria to avoid revisits of the same solutions from different initial guesses is developed. It is shown tha t for general nonlinear systems or gradient systems, it is computationally infeasible to construct a set of initial guesses which converge to the set of all type-one equilibrium points on the stability boundary of a stable eq uilibrium point x(s) from a finite number of function values and derivative s near x(s) using the Newton homotopy method. Several examples are applied to illustrate the theoretical developments.