Semilocal convergence analysis of an eighth order iterative method for solving nonlinear systems

: In this paper, the semilocal convergence of the eighth order iterative method is proved in Banach space by using the recursive relation, and the proof process does not need high order derivative. By selecting the appropriate initial point and applying the Lipschitz condition to the ﬁrst order Fr´echet derivative in the whole region, the existence and uniqueness domain are obtained. In addition, the theoretical results of semilocal convergence are applied to two nonlinear systems, and satisfactory results are obtained


Introduction
In mathematics and various sciences, the changes of output and input of nonlinear systems are out of proportion. Most of the systems involved in life are essentially nonlinear, so solving nonlinear problems has attracted various scientists. Scholars have proposed more efficient iterative methods for solving nonlinear systems. One of the most famous iterative methods for solving nonlinear systems is Newton's method [1], for k = 0, 1, 2, ..., x 0 is the starting point. The Newton's method is second order convergent and effective in solving some nonlinear systems.
The theoretical results of local convergence and semilocal convergence of the iterative method are also important in the study. Local convergence requires the existence of the assumed solution and the initial value is close enough to the solution. Semilocal convergence does not require the existence of an assumed solution, but the selection of initial values also needs to meet certain conditions (see [14][15][16][17][18]). Therefore, for some systems that cannot be analyzed and solved, the results of semilocal convergence cannot only prove the convergence of iterative sequences, but also prove the existence of solutions of these systems, so as to obtain the existence domain and uniqueness domain of system solutions; for further study (see [19][20][21][22]). Based on this, we perform a semilocal convergence analysis on the method (1.2).
This paper consists of five sections. In Section 2 of the paper, the recurrence relation is explained. The semilocal convergence of the iterative method (1.2) is proved in Section 3. In Section 4, the numerical experiments of two nonlinear systems are completed. Finally, the conclusion of this paper is made.

Recurrence relations
In this section, let X and Y be Banach spaces and let F : Ω ⊆ X → Y be a twice differetiable nonlinear Fréchet operator in an open Ω [23]. Let us assume that the inverse of the Jacobian matrix of the system in the iteration (1.2) is Γ 0 ∈ L(Y, X), which is the set of linear operation from Y to X.

(2.4)
These functions play a key role in the analysis that will be performed next.
Preliminary results. In order to get the difference of the first two elements in the iterative method (1.2), we have The Taylor series expansion of F around x 0 evaluated in y 0 is where the term F(x 0 ) + F (x 0 )(y 0 − x 0 ) is equal to zero, since it comes from a Newton's step. With the change Then, (2.8) Taking norms and applying Lipschitz condition, we get so that, Using a method similar to (2.5), we get w 0 − y 0 Next, the next step analysis Using Taylor's expansion of F(w 0 ) around x 0 and applying Lipschitz condition, we obtain (2.14) By applying Banach's lemma, one has Then, as far as a 0 (h(a 0 )) < 1 (by taking a 0 < 0.45807), Banach's lemma guarantees that .
Based on the above analysis, we can obtain the following theorem.

Semilocal convergence analysis
According to the convergence property of x k sequence in Banach space, we need to prove that this sequence is a Cauchy sequence. Based on the auxiliary function, we can obtain the following results. Lemma 1. According h(x), f (x) and g(x), we have: i. f (x) is inceasing and f (x) > 1 for x ∈ (0, 0.45807), ii. h(x) and g(x) are increasing for x ∈ (0, 0.45807). The above lemma can be calculated from the Section 2, and the process is omitted.
Proof. It is straightforword that i, ii are satisfied. As f (a 0 ) 2 g(a 0 ) < 1, then by construction of a k , it is a dereasing sequence. So a k < a 0 ≤ 0.21715, for all k ≥ 1.
1− f (a 0 )g(a 0 ) . In the case, the iterates {x k } and {y k } are contained in B e (x 0 , Rη) and x * ∈ B e (x 0 , Rη). In addition, the x * is the only solution of equation F(x) = 0 in B n (x 0 , 2 Kβ − Rη) ∩ Ω. Proof. By recursively applying (O4 k ), we can write Then, f (a i )g(a i ) As h(x) is increasing and a k dreasing, it can be stated that Moreover, according Lemmas 1 and 2, by using the expression for the partial sum of a geometrical series, So, the Cauchy sequence if and only if f (a 0 )g(a 0 ) < 1 (Lemma 2).
when m → ∞, we get the radius od convergence Rη = h(a 0 ) 1− f (a 0 )g(a 0 ) η. Let's prove that x * is the solution of F(x) = 0 starting from the boundary of F (x n ) , as h(x), f (x) and g(x) are increasing and a k is the decreasing sequence, Taking into account that F (x k ) is bounded and ( f (a 0 )g(a 0 )) k tends to zero when k → ∞, we conclude that F(x k ) → 0. As F is continuous in Ω, then F(x * ) = 0. Finally, the uniqueness of (3.9) In order to guarantee that y * − x * = 0 it is necesssary to prove that operator 1 0 F (x * + t(y * − x * ))dt is invertible. Applying hypothesis (M 3 ), (3.10) By the Banach lemma, the intergal operator is invertible and hence y * = x * .

Numerical results
In this section, we provide some numerical examples to illustrate the theoretical results introduced earlier. [24], which is given as follows:

Example 1. Hammerstein equation is a kind of important nonlinear integral equation
where x ∈ C[0, 1], s, t ∈ [0, 1], with the kernel N is To solve (4.1) we transform it into a syste of nonlinear equations through a discretization process. We approximate the integral appearing in Eq (4.1) by using Gauss-Legendre quadrature, being t j and w j the nodes and the weights of the Gauss-Legendre polynomial. Denoting the approximation of x(t j ) as x i , i = 1, ..., 7, then we estimate (4.1) with the nonlinear system of equations where So,the system can be rewritten as where F if a nonlinear operator in the Banach space R L ,and F is its Fréchet derivative in L(R L , R L ). According the method (1.2), we will use it to solve the nonlinear systems. Taking x 0 = (1.8, 1.8, ..., 1.8) T , L = 7 and the infinity norm, we get The above results satisfy the semilocal convergence condition, so this method can be applied to the system. Thus, we guarantee the existence of the solution in B e (x 0 , 0.5646), and the uniqueness in B n (x 0 , 22.4874). Table 1 shows the the radius of the existence domain and the radius of the unique domain under different initial values. For x 0i > 1.87, i = 1, 2, ..., 7, convergence conditions are not satisfied and, therefore, the convergence is not guaranteed. Example 2. Let X = Y = R 2 be equipped with the max-norm. Choose: x 0 = (0.9, 0.9) T , s ∈ [ 0, 1 2 ). Let s = 0.49, define function F by The fréchet-derivative of operator F is given by Taking x 0 = (0.9, 0.9) T and the infinity norm, we get The convergence conditions are met and consequently the method can be applied to the system. The existence domain of the solution is B e (x 0 , 0.9101), and the uniqueness domain is B n (x 0 , 1.2192). Taking x 0 = (0.73, 0.73) T and the infinity norm, then The existence domain of the solution is B e (x 0 , 0.5095), and the uniqueness domain is B n (x 0 , 0.943534). When the initial value satisfies the Kantorovich condition and the range of a 0 obtained, the initial value within that range is taken to solve the system. Iterative method (1.2) for solving nonlinear (4.4) with roots of x * = (0.7884, 0.7884) T .
Similar results can be obtained in Tables 2 and 3, that is, under the Kantorovich condition, by selecting different initial values, we can converge to a unique solution. When the initial value is closer to the root, the error estimate is lower. This semilocal convergence that can prove the existence and uniqueness of solutions under certain assumptions is very valuable.

Conclusions
In this paper, the semilocal convergence of the eighth order iterative method (1.2) is studied. By analyzing the behavior of the iterative method under the Kantorovich condition, the Lipschitz condition is applied to the first derivative, and the theory of semilocal convergence of the iterative method is obtained by using the recurrence relation. The existence and uniqueness domain of the solution of the nonlinear system is obtained. In the experimental part, a classical Hammerstein nonlinear integral equation and a matrix function are solved. The experimental results are consistent with expectations, and the high-precision approximation of the system solution also proves the effectiveness of the method numerically.

Use of AI tools declaration
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.