The Newton algorithm based on the “continuation” method may be written as being governed by the equation xj(t) + Bij-1 Fi(xj) = 0, where Fi(xj) = 0, i, j = 1, …n are nonlinear algebraic equations (NAEs) to be solved, and Bij = ∂Fi /∂xj is the corresponding Jacobian matrix. It is known that the Newton's algorithm is quadratically convergent; however, it has some drawbacks, such as being sensitive to the initial guess of solution, and being expensive in the computation of the inverse of Bij at each iterative step. How to preserve the convergence speed, and to remove the drawbacks is a very important issue in the solutions of NAEs. In this paper we discretize the above equation being written as Bij xj(t) + Fi(xj) = 0, by a backward difference scheme in a new time scale of s = 1–e -t , and an ODEs system is derived by introducing a fictitious time-like variable. The new algorithm is obtained by applying a numerical integration scheme to the resultant ODEs. The new algorithm does not need the inverse of Bij, and is thus resulting in a significant reduction in computational time than the Newton's algorithm. A similar technique is also used to modify the homotopy method. Numerical examples given confirm that the modified Newton method is highly efficient, insensitive to the initial condition, to find the solutions with a very small the residual error.
Atluri, Satya N.; Liu, Chein-Shan; and Kuo, Chung-Lun
"A MODIFIED NEWTON METHOD FOR SOLVING NON-LINEAR ALGEBRAIC EQUATIONS,"
Journal of Marine Science and Technology: Vol. 17
, Article 9.
Available at: https://jmstt.ntou.edu.tw/journal/vol17/iss3/9