A Preconditioned Iterative Method for Solving Systems of Nonlinear Equations Having Unknown Multiplicity

A modification to an existing iterative method for computing zeros with unknown multiplicities of nonlinear equations or a system of nonlinear equations is presented. We introduce preconditioners to nonlinear equations or a system of nonlinear equations and their corresponding Jacobians. The inclusion of preconditioners provides numerical stability and accuracy. The different selection of preconditioner offers a family of iterative methods. We modified an existing method in a way that we do not alter its inherited quadratic convergence. Numerical simulations confirm the quadratic convergence of the preconditioned iterative method. The influence of preconditioners is clearly reflected in the numerically achieved accuracy of computed solutions.


Introduction
The design of an iterative method for solving nonlinear equations and systems of nonlinear equations is an active area of research. Many researchers have proposed iterative methods for solving nonlinear and systems of nonlinear equations for finding simple zeros or zeros with multiplicity greater than one [1][2][3][4][5][6][7][8][9][10][11][12][13][14][15]. The classical iterative method for solving nonlinear and systems of nonlinear equations to find simple zeros is the Newton method, which offers quadratic convergence [16,17] under certain conditions. When we are talking about the iterative method for solving nonlinear equations or systems of nonlinear equations to find zeros with multiplicities greater than one, the classical Newton all z. The first order Fréchet derivative of (7) can be computed as If we apply the Newton method to (7), we obtain The convergence order of (10) is two. The iterative method (6) with a general preconditioner can be written as The convergence order of (11) is also two. The modified Newton method [17,21,22] for solving nonlinear equations with unknown multiplicity can be developed in this way. We define a new function The application of the Newton method to (12) gives The order of convergence of (13) is two. Noor and his co-researchers [23] have constructed a family of iterative methods for solving nonlinear equations with unknown multiplicity by introducing a preconditioner. They defined a new function and application of the Newton method to (14) gives where λ(z) is a non-zero function. The order of convergence of (15) is two.

Proposed Method
When we observe (14), we can notice that the preconditioner is only introduced for φ(z), and not for φ (z). We will also introduce a preconditioner for φ (z), and will show that the convergence order of (14) is still quadratic. We define a new function and after applying the Newton method, we obtain where ω(z) is a non-zero function. For the purpose of generalization of the iterative method (17) to a system of nonlinear equations, we define a new function Q(z) The first order Fréchet derivative of (18) can be written as Further simplification of Q (z) −1 Q(z) gives If compare the underlined expressions in (17) and (20), they are different. Generally, it is not possible to commute ( However, we , and get the following iterative method. It can be seen that the iterative method (21) is not the application of the Newton method to (18). The iterative method (17) for solving scalar nonlinear equations with unknown multiplicity and vector version (21) are exactly the same. We will only provide the proof of quadratic convergence for (21), and it is automatically applicable to scalar version (17). An iterative method was proposed in [24] to compute the zeros with multiplicity of system of nonlinear equations that used preconditioners for a system of nonlinear equations, but not for the Jacobian of the system of nonlinear equations. Notice that in this article we are introducing preconditioners for the system of nonlinear equations as well as the Jacobian of the system of nonlinear equations.

Convergence
In the following theorem, we established the proof of quadratic convergence of (21).

Numerical Testing
The two preconditioners ω(z) and λ(z) produce families of iterative methods. If we define ω(z) = exp( z) and λ(z) = exp(ϑ z) , we get the following two-parameter family of iterative methods for solving nonlinear equations that have zeros with unknown multiplicity.
We only conducted numerical testing for the system of nonlinear equations, and the cases for the nonlinear equations are similar. It is important to test the computational convergence order (CCO) of the proposed iterative methods. In all our simulations, we adopted the following definition of CCO: For numerical simulations, three problems were selected with different multiplicities. The performance of iterative method (11) is not better, comparatively. The various choices for the preconditioners are made in Tables 1-3 for all three problems. In Table 1, we have shown that the selection of preconditioners has an influence on the numerical accuracy of computed zeros with multiplicities. Moreover, the computational cost of performing the different operation is reasonable, because in all cases, we selected preconditioners in a way that their first-and second-order Fréchet derivatives are diagonal matrices. When we selected Λ Λ Λ(z) = 6 + cos(z)/10 and Ω Ω Ω(z) = 6 + cos(z)/10 for Problem 1, we achieved the best accuracy in computed zeros with different multiplicities. For the second problem, Table 2 shows that the selection of Λ Λ Λ(z) produces good accuracy. In Table 3, again the selection of both preconditioners provides the best accuracy, comparatively.

Conclusions
The inclusion of preconditioners in the existing iterative methods for finding zeros with multiplicities for solving a system of nonlinear equations gives benefits in numerical stability and numerical accuracy. The proposed methodology is equally effective for nonlinear and systems of nonlinear equations. It is assumed in all cases that the preconditioners should be non-zero, because in this way, it does not affect the zeros of nonlinear or systems of nonlinear equations. The different selections of preconditioners provide different families of iterative methods. The claimed order of convergence is also verified by computing the computational order of convergence in all numerical simulations. Study of the dynamics of nonlinear preconditioners for finding zeros with multiplicities of nonlinear equations and systems of nonlinear equations could be an interesting topic for research.
Author Contributions: First author established the idea and all other authors contributed equally in the article.

Conflicts of Interest:
The authors declare no conflict of interest.