Next Article in Journal
Numerical Non-Linear Modelling Algorithm Using Radial Kernels on Local Mesh Support
Previous Article in Journal
Interpolative Reich–Rus–Ćirić and Hardy–Rogers Contraction on Quasi-Partial b-Metric Space and Related Fixed Point Results
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

General Local Convergence Theorems about the Picard Iteration in Arbitrary Normed Fields with Applications to Super–Halley Method for Multiple Polynomial Zeros

Faculty of Physics and Technology, University of Plovdiv Paisii Hilendarski, 24 Tzar Asen, 4000 Plovdiv, Bulgaria
Mathematics 2020, 8(9), 1599; https://doi.org/10.3390/math8091599
Submission received: 23 June 2020 / Revised: 5 August 2020 / Accepted: 14 September 2020 / Published: 17 September 2020

Abstract

:
In this paper, we prove two general convergence theorems with error estimates that give sufficient conditions to guarantee the local convergence of the Picard iteration in arbitrary normed fields. Thus, we provide a unified approach for investigating the local convergence of Picard-type iterative methods for simple and multiple roots of nonlinear equations. As an application, we prove two new convergence theorems with a priori and a posteriori error estimates about the Super-Halley method for multiple polynomial zeros.

1. Introduction and Preliminaries

Let ( K , | · | ) be a normed field and f : K K be an arbitrary function. It is well known that one of the most commonly used tools for finding the zeros of f is the Picard iteration
x k + 1 = T x k , k = 0 , 1 , 2 , ,
where T : D K K is some iteration function (see, e.g., Traub [1]). Although the Picard iteration is an old method, it still draws strong interest as a powerful tool for solving a wide range of mathematical problems (see e.g., [2,3,4] and references therein). Nowadays, there are a plethora of convergence results about the iterative methods of the type (1) but the proofs are often based on Taylor’s expansions or other methods that require the existence of higher-order derivatives of f and ensure only asymptotic error constants where available. Moreover, the most of the existing convergence theorems do not provide exact information about the sets of initial approximations that guarantee the convergence of the iteration (1) (see e.g., [5,6,7] and the references therein).
In the last decade, a general convergence theory of the Picard iteration (1) in cone metric spaces and in n-dimensional vector spaces has been developed by Proinov [8,9,10,11]. The main role in this theory is played by a real-valued function called function of initial conditions of the iteration function T. Following the ideas of Proinov [11], in this paper we make two contributions. First, in Section 2, using two different functions of initial conditions we prove two general local convergence theorems (Theorem 1 and Theorem 2) that provide domains of initial approximations (convergence domains) to guarantee the convergence of the Picard iteration (1) with a priori and a posteriori error estimates right from the first step. In this manner, we reduce the convergence analysis of (1) up to studying of some simple properties of the iteration function T. Second, applying this general approach, in Section 3 we obtain two new local convergence theorems with error estimates (Theorem 3 and Theorem 4) about the Super-Halley method for approximation of polynomial zeros with known multiplicity m. The obtained theorems are new even in the case of simple zeros ( m = 1 ).
In the remainder of this section, we recall some definitions and theorems of Proinov [9,12] which are crucial for the proof of the upcoming results.
Furthermore, J always denotes an interval on R + containing 0 and we assume that 0 0 1 . Also, we denote by S k ( r ) the sum of the first k terms of geometric sequence 1 , r , r 2 , , i.e., for all k N , we have
S k ( r ) = 1 + r + + r k 1 .
In the case k = 0 , we set S 0 ( r ) = 0 .
Definition 1
([9] Definition 2.1). A function φ : J R R + is called quasi-homogeneous of degree p 0 if it is such that φ ( λ t ) λ p φ ( t ) for all λ [ 0 , 1 ) and t J .
Recall some simple but useful properties of the quasi-homogeneous functions.
(P1)
A function f is quasi-homogeneous of degree p = 0 on some interval J if and only if f is nondecreasing on J;
(P2)
If a function f is quasi-homogeneous of degree p on J, then there exists a nondecreasing function F : J R + such that f ( t ) = t p F ( t ) ;
(P3)
If g and h are quasi-homogeneous of degree p 0 and q 0 on J, then g h is quasi-homogeneous of degree p + q on J.
Definition 2
([9] Definition 2.4). A function φ : J R + is said to be a gauge function of order p 1 on J if it satisfies the following conditions:
(i) 
φ is a quasi-homogeneous function of degree p on J;
(ii) 
φ ( t ) t for all t J .
A gauge function φ of order p on J is said to be a strict gauge function if the inequality (ii) holds strictly whenever t J \ { 0 } .
Definition 3
([9] Definition 3.1). Let T : D X X be a map on an arbitrary set X. A function E : D R + is said to be a function of initial conditions of T (with gauge function φ on an interval J) if there exists a function φ : J J such that
E ( T x ) φ ( E ( x ) ) f o r a l l x D w i t h T x D a n d E ( x ) J .
Definition 4
([9] Definition 3.2). Let T : D X X be a map on an arbitrary set X. Suppose E : D R + is a function of initial conditions of T (with gauge function on an interval J). Then a point x D is said to be an initial point of T if E ( x ) J and all of the iterates T n ( x ) ( n = 0 , 1 , ) are well defined and belong to D.
The next theorem provides a sufficient condition for initial points.
Proposition 1
([9] Proposition 4.1). Let T : D X X be a map on a set X, and let E : D R + be a function of initial conditions of T with a gauge function φ on J. Suppose that
x D w i t h E ( x ) J i m p l i e s T x D .
Then every point x 0 D such that E ( x 0 ) J is an initial point of T.
Proposition 2
([12] Lemma 3.6). For every initial point x 0 D of T and every k 0 , we have
E ( x k + 1 ) φ ( E ( x k ) ) a n d E ( x k ) φ k ( E ( x 0 ) ) .
In addition, if φ is a gauge function of order r 1 , then
E ( x k ) E ( x 0 ) λ S k ( r ) a n d ϕ ( E ( x k ) ) λ r k ,
where λ = ϕ ( E ( x 0 ) ) and ϕ is a non-negative and nondecreasing function on J such that φ ( t ) = t ϕ ( t ) for all t J .

2. General Local Convergence Theorems

Let T : D K K be an arbitrary iteration function. In this section, we present two general theorems about the local convergence of the Picard iteration (1). The first one (Theorem 1) is a case of Theorem 3.1 of Proinov [11] but provides an extra error estimate that ensures the Q-order of convergence of the Picard iteration (1) (see, e.g., [13]). The second theorem (Theorem 2) is of bit more practical importance than Theorem 1 since there is less unknown information involved into the initial conditions and error estimates.
Let f : D K K be an arbitrary function which has at least one zero in D . For two points x , ξ K , we define the functions E : K R + and ε : D K R + by
E ( x ) = | x ξ | d and ε ( x ) = | x ξ | ρ ( x ) ,
where d denotes the distance between ξ and the nearest zero of f and ρ ( x ) is the radius of a disk with center x that not contains zeros of f. If ξ is a zero of f and it is unique, then we set E ( x ) 0 . In this case, we assume that d = . Also, we observe that the domain of ε is the set
D = { x K : ρ ( x ) > 0 } .
We have to note that usually ξ is a zero of f but we do not assume this in Theorem 1 and Theorem 2.
The following theorem is our first main result.
Theorem 1.
Let T : D K K be an iteration function, ξ K and E : K R + be defined by (2). Suppose ϕ : J R + is a quasi-homogeneous function of degree p 0 such that for each x K with E ( x ) J , the following two conditions are satisfied:
(a) 
x D
(b) 
| T x ξ | ϕ ( E ( x ) ) | x ξ | .
Also, let x 0 K be an initial approximation such that
E ( x 0 ) J a n d ϕ ( E ( x 0 ) ) < 1 ,
then the following statements hold:
(i) 
The Picard iteration (1) is well defined and converges to ξ with order r = p + 1 .
(ii) 
For all k 0 , we have the following error estimates:
| x k + 1 ξ | λ r k | x k ξ | a n d | x k ξ | λ S k ( r ) | x 0 ξ | ,
where λ = ϕ ( E ( x 0 ) ) .
(iii) 
The Picard iteration (1) converges to ξ with Q-order r = p + 1 and with the following error estimate:
| x k + 1 ξ | ( R d ) 1 r | x k ξ | r f o r a l l k 0 ,
where R is the minimal solution of the equation ϕ ( t ) = 1 in the interval J \ { 0 } .
Proof. 
The claims (i) and (ii) follow from Theorem 3.1 of Proinov [11] in the one dimensional case. So, it remains to prove the estimate (iii).
By (i) and (ii), for all k 0 we have x k U , where U = { x K : E ( x ) J } . Since, ϕ is a quasi-homogeneous function of degree p 0 , then by (P2) it follows that there exists a nondecreasing function Φ : J R + such that ϕ ( t ) = t p Φ ( t ) . This together with ϕ ( R ) = 1 implies that Φ ( R ) = 1 / R p . Hence, from conditions (a) and (b), for all x U we get
| T x ξ | ϕ ( E ( x ) ) | x ξ | = Φ ( E ( x ) ) d p | x ξ | p + 1 1 ( R d ) p | x ξ | p + 1 .
Now, setting x = x k and p = r 1 in the last inequality, we get the estimate of (iii) which implies that the Picard iteration (1) converges to ξ with Q-order r = p + 1 and completes the proof. □
In what follows, for a nondecreasing function β : J R + , we define the function ψ : J R by
ψ ( t ) = 1 t ( 1 + β ( t ) ) .
If the function ψ is positive, then we can define the function ϕ : J R + by
ϕ ( t ) = β ( t ) / ψ ( t ) .
The second main result of this paper is represented by the following theorem:
Theorem 2.
Let T : D K K be an iteration function, ξ K and ε : D K R + be defined by (2). Suppose β : J R + is a nonzero quasi-homogeneous function of degree p 0 and for each x K with ε ( x ) J , the following two conditions are satisfied:
(a) 
x D ;
(b) 
| T x ξ | β ( ε ( x ) ) | x ξ | .
Let also, x 0 K be an initial guess such that
ε ( x 0 ) J a n d β ( ε ( x 0 ) ) ψ ( ε ( x 0 ) ) ,
where the function ψ is defined by (4). Then the Picard iteration (1) is well defined and converges to ξ with the following error estimates:
| x k + 1 ξ | θ λ r k | x k ξ | a n d | x k ξ | θ k λ S k ( r ) | x 0 ξ | f o r a l l k 0 ,
where λ = ϕ ( ε ( x 0 ) ) and θ = ψ ( ε ( x 0 ) ) with ϕ and ψ defined by (5) and (4), respectively. In addition, if the second inequality in (6) is strict, then the order of convergence of Picard iteration (1) is at least r = p + 1 .
Proof. 
The initial condition (6) can be written in the form ε ( x 0 ) Δ , where Δ = { t J : β ( t ) ψ ( t ) } . It is not hard to verify that for all t Δ , the following inequalities hold (see, e.g., ([11] Lemma 7.3)):
0 < ψ ( t ) 1 , 0 ϕ ( t ) 1 and 0 β ( t ) < 1 .
According to the first of these inequalities and the fact that ρ ( x ) > 0 , we get
ρ ( T x ) ψ ( ε ( x ) ) ρ ( x ) > 0 .
Indeed, if η D is a zero of f : D K K , then by the triangle inequality, the definition of ρ ( x ) and the condition (b), we obtain
| T x η | | x η | | x ξ | | T x ξ | ρ ( x ) ( 1 + β ( ε ( x ) ) | x ξ | = [ 1 ε ( x ) ( 1 + β ( ε ( x ) ) ) ] ρ ( x ) = ψ ( ε ( x ) ) ρ ( x ) .
This implies (9) which in turn means that T x D . So, dividing both sides of (b) by ρ ( T x ) , we reach the inequality
ε ( T x ) φ ( ε ( x ) ) ,
where φ ( t ) = t ϕ ( t ) with ϕ defined by (5). Since the inequality (10) implies that ε ( T x ) J , then we have both
T x D and ε ( T x ) J
which according to (a) means that T x D . This, together with (10) purports that ε is a function of initial conditions of T with gauge function φ on the interval J (see Definition 3).
Now let U be the set U = { x D : ε ( x ) Δ } . Then the condition (6) implies that x 0 U . On the other hand (10) and (11) imply that T ( U ) U which means that starting from x 0 the Picard iteration (1) is well defined and remains in the set U. The convergence to ξ follows from the second estimate of (7). So, it remains to prove the estimates (7). Since x 0 is an initial point of T with respect to ε , and φ is a gauge function of order r = p + 1 , then from Proposition 2 it follows that
ε ( x k ) λ S k ( r ) ε ( x 0 ) ,
where λ = ϕ ( ε ( x 0 ) ) . From this, taking into account that β is quasi-homogeneous function of degree p 0 , and that β = ϕ ψ , we get
β ( ε ( x k ) ) β ( λ S k ( r ) ε ( x 0 ) ) λ p S k ( r ) β ( ε ( x 0 ) ) = λ p S k ( r ) + 1 ψ ( ε ( x 0 ) ) = θ λ r k .
From this and condition (b) applied to x k , we obtain the first estimate in (7). From the first estimate in (7) one can easily obtain the second one which in turn implies the convergence of the Picard iteration (1) due to S k ( r ) k and 0 θ λ < 1 . □

3. Local Convergence of the Super-Halley Method for Multiple Polynomial Zeros

In 2008, Osada [5] derived a modification of the Chebyshev–Halley family of iterative methods for computation of multiple zeros of known multiplicity m. The most famous iterative methods, namely Newton’s, Halley’s and Chebyshev’s methods for multiple zeros are members of this family. A detailed local convergence analysis of these illustrious methods applied to multiple polynomial zeros can be found in the papers [8,14,15,16]. Another member of the Chebyshev–Halley family that has been rarely studied is the so called Super-Halley method for multiple zeros, which can be defined by the following iteration (see, e.g., [17,18] and references therein):
x k + 1 = S x k ,
where the iteration function S : D K K is defined by
S x = x N ( x ) 2 m + 1 1 L ( x ) if f ( x ) 0 , x if f ( x ) = 0
with N ( x ) and L ( x ) defined as follows
N ( x ) = f ( x ) f ( x ) and L ( x ) = N ( x ) f ( x ) f ( x ) .
It is seen that the domain of the function (13) is the set
D = x K : f ( x ) 0 1 L ( x ) 0 .
In this section, we implement Theorem 1 and Theorem 2 to the Super-Halley method (12) applied for computation of polynomial zeros of known multiplicity m. Thus, we obtain two new local convergence theorems (Theorem 3 and Theorem 4) that provide sufficient conditions to guarantee the cubic convergence of the method (12) right from the first step. A priory and a posteriori error estimates at any iteration are also provided. The obtained results are new even in the case of simple zeros. Note that in this case, i.e., when m = 1 the Super-Halley method (12) was presented in 1992 by Hernández-Verón [19] and later investigated in [20,21], etc.
Before proceeding further, we give two useful technical lemmas.
Lemma 1.
Let K be an arbitrary field, and let f K [ x ] be polynomial of degree n 2 which splits over K. Suppose that ξ 1 , , ξ s are all distinct zeros of f with multiplicities m 1 , , m s j = 1 s m j = n , respectively.
(i) 
If x K is not a zero of f, then for any i = 1 , , s we have
f ( x ) f ( x ) = m i + a i x ξ i , w h e r e a i = ( x ξ i ) j i m j x ξ j .
(ii) 
If x K is not a zero of both f and f , then for any i = 1 , , s we have
f ( x ) f ( x ) = ( m i + a i ) 2 m i b i ( x ξ i ) ( m i + a i ) , w h e r e a i i s d e f i n e d i n ( i ) a n d b i = ( x ξ i ) 2 j i m j ( x ξ j ) 2 .
Proof. 
(i) Using a well known identity (see, e.g., ([15] Lemma 4.2(i))), we obtain
f ( x ) f ( x ) = j = 1 s m j x ξ j = m i x ξ i + j i m j x ξ j = m i + a i x ξ i , where a i = ( x ξ i ) j i m j x ξ j
which proves the first claim.
The second claim follows from (i) and the following identities (see, e.g., ([15] Lemma 4.2(ii)) and ([14] Equation (15))):
f ( x ) f ( x ) = f ( x ) f ( x ) f ( x ) f ( x ) j = 1 s m j ( x ξ j ) 2 and j = 1 s m j ( x ξ j ) 2 = m i + b i ( x ξ i ) 2 .
 □
Lemma 2.
Let x , ξ K and ξ 1 , , ξ s K be all zeros of f which are not equal to ξ, then for any j = 1 , , s the following inequality holds:
| x ξ j | ( 1 E ( x ) ) d ,
where E : K R + is defined by (2).
Proof. 
Since d | ξ ξ j | for all j = 1 , , s , then by the triangle inequality in K , we get
| x ξ j | = | ξ ξ j + x ξ | | ξ ξ j | | x ξ | ( 1 E ( x ) ) d .
 □

3.1. Local Convergence of the First Type

From now on K [ x ] denotes the ring of polynomials over the normed field K . Let f K [ x ] be a polynomial of degree n 2 . We shall study the convergence of the Super-Halley method (12) with respect to the function of initial conditions E : K R + defined by (2).
Onwards, for n m 1 , we define the real functions g and h by
g ( t ) = ( n m ) 2 ( n m ) t 1 t + n t 2 and h ( t ) = 2 ( m n t ) ( ( 2 m n ) t 2 2 m t + m ) 1 t .
Obviously, the function g is positive and quasi-homogeneous of second degree on the interval [ 0 , 1 ) and the function h is decreasing and positive on the interval [ 0 , τ ) , where
τ = m / n if n 2 m , m / ( m + m ( n m ) ) if n < 2 m .
Hence, we can define the function ϕ : [ 0 , τ ) R + by
ϕ ( t ) = g ( t ) h ( t ) = ( n m ) ( n + ( n 2 m ) t ) t 2 2 ( m n t ) ( ( 2 m n ) t 2 2 m t + m ) .
Note that the function ϕ is quasi-homogeneous of second degree on the interval [ 0 , τ ) owing to the properties (P1) and (P3).
The next lemma shows that the function ϕ defined by (19) satisfies the conditions (a) and (b) of Theorem 1 with the Super-Halley iteration function and the function E defined by (2).
Lemma 3.
Let f K [ x ] be a polynomial of degree n 2 which splits over K , and let ξ K be a zero of f with multiplicity m. Let x K be such that
E ( x ) < τ ,
where the function E is defined by (2) and τ is the number (18). Then:
(i) 
x belongs to the set D defined by (15).
(ii) 
| S x ξ | ϕ ( E ( x ) ) | x ξ | , where the function ϕ is defined by (19).
Proof. 
Let x K satisfy (20). If m = n or x = ξ , then S x = ξ and the statements of the lemma hold. Let m < n and x ξ . Let also ξ 1 , , ξ s be all distinct zeros of f with respective multiplicities m 1 , , m s and the quantities a i and b i be defined as in Lemma 1. Put ξ = ξ i , m = m i , a = a i and b = b i for some 1 i s .
(i) According to (15), we have to prove that f ( x ) 0 implies 1 L ( x ) 0 , where L ( x ) is defined by (14). From E ( x ) < τ < 1 and Lemma 2, we get
| x ξ j | ( 1 E ( x ) ) d > 0
for each j i which means that f ( x ) 0 . Then, from (14) and Lemma 1 (i), we get
1 N ( x ) = f ( x ) f ( x ) = m + a x ξ , where a = ( x ξ ) j i m j x ξ j .
Now, from the triangle inequality and (21), we obtain the following estimate:
| a | | x ξ | j i m j | x ξ j | | x ξ | ( 1 E ( x ) ) d j i m j = ( n m ) E ( x ) 1 E ( x ) .
From this, using the triangle inequality and E ( x ) < τ m / n , we get
| m + a | m | a | m ( n m ) E ( x ) 1 E ( x ) = m n E ( x ) 1 E ( x ) > 0 .
Hence, m + a 0 which implies f ( x ) 0 . Then, from Lemma 1 (ii), we have
f ( x ) f ( x ) = ( m + a ) 2 m b ( x ξ ) ( m + a ) , where b = ( x ξ ) 2 j i m j ( x ξ j ) 2 .
Now, from (14), (22) and (25), we obtain
1 L ( x ) = 1 N ( x ) f ( x ) f ( x ) = 1 ( m + a ) 2 m b ( m + a ) 2 = m + b ( m + a ) 2 .
Further, by means of the triangle inequality, (21) and E ( x ) < τ , we reach the estimates (see e.g., ([14] Equation (11)))
| b | ( n m ) E ( x ) 2 ( 1 E ( x ) ) 2 and | m + b | m | b | h ( E ( x ) ) 2 ( m n E ( x ) ) ( 1 E ( x ) ) > 0 .
The last estimate implies that | 1 L ( x ) | > 0 and so x D .
(ii) From the definition of the Super-Halley iteration function (13), (22) and (25), we easily get
S x ξ = x ξ x ξ 2 ( m + a ) m + ( m + a ) 2 m + b = σ ( x ξ ) ,
where
σ = 1 2 m 2 + m b + 2 m a + a 2 2 ( m + a ) ( m + b ) = 2 a b + m b a 2 2 ( m + a ) ( m + b ) .
Hence, to complete the proof it is sufficient to estimate | σ | from above. With this aim, using the estimates (23), (24) and (26), we get
| σ | 2 | a | | b | + m | b | + | a | 2 2 | m + a | | m + b | 2 ( n m ) E ( x ) 1 E ( x ) ( n m ) E ( x ) 2 ( 1 E ( x ) ) 2 + m ( n m ) E ( x ) 2 ( 1 E ( x ) ) 2 + ( n m ) 2 E ( x ) 2 ( 1 E ( x ) ) 2 2 m n E ( x ) 1 E ( x ) h ( E ( x ) ) 2 ( m n E ( x ) ) ( 1 E ( x ) ) = g ( E ( x ) ) h ( E ( x ) ) = ϕ ( E ( x ) ) .
This completes the proof of the lemma. □
The following theorem is the first convergence result about the Super-Halley method (12).
Theorem 3.
Let f K [ z ] be a polynomial of degree n 2 which splits over K , and let ξ K be a zero of f with multiplicity m. Suppose x 0 K n is an initial approximation satisfying
E ( x 0 ) < R = 2 m n + m + 3 ( n m ) ( n + m ) ,
where the function E is defined by (2). Then the Super-Halley iteration (12) is well defined and converges cubically to ξ with the following error estimates for all k 0 :
| x k + 1 ξ | λ 3 k | x k ξ | a n d | x k ξ | λ ( 3 k 1 ) / 2 | x 0 ξ | ,
where λ = ϕ ( E ( x 0 ) ) and the function ϕ is defined by (19). Besides, the iteration (12) converges Q-cubically to ξ with the following error estimate
| x k + 1 ξ | 1 ( R d ) 2 | x k ξ | 3 f o r a l l k 0 .
Proof. 
It is easy to verify that R τ , where τ is defined by (18) and ϕ ( t ) < 1 for all t [ 0 , R ) . Hence, the initial condition (29) is equivalent to (3) with J = [ 0 , τ ) and ϕ defined by (19). Therefore, the proof follows from Lemma 3 and Theorem 1. □
Remark 1.
It is interesting to note that in the case n > 2 m , Theorem 3 gives larger convergence domain and better error estimates for the Super-Halley method (12) than Theorem 4.5 of [15] for the Halley’s method.

3.2. Local Convergence of the Second Type

Here, we study the convergence of the Super-Halley method (12) with respect to the function of initial conditions ε : D K R + defined by (2).
For n > m 1 , we define the real function β : [ 0 , m / ( n m ) ) R + by
β ( t ) = ( n m ) ( n + 2 ( n m ) t ) t 2 2 ( m ( n m ) t ) ( m ( n m ) t 2 ) .
Apparently, β is a quasi-homogeneous function of second degree on [ 0 , m / ( n m ) ) .
The role played by the following lemma is the same as the one of Lemma 3.
Lemma 4.
Let f K [ x ] be a polynomial of degree n 2 which splits over K , and let ξ K be a zero of f with multiplicity m. Let x K be such that
ε ( x ) < m / ( n m ) ,
where the function ε is defined by (2). Then:
(i) 
x D , where D is the set (15).
(ii) 
| S x ξ | β ( ε ( x ) ) | x ξ | , where the function β is defined by (32).
Proof. 
The proof is the same as those of Lemma 3. One just should use the estimate | x ξ i | ρ ( x ) > 0 instead of (21). □
Now we can state the second main result of this section.
Theorem 4.
Let f K [ x ] be a polynomial of degree n 2 which splits over K , and let ξ K be a zero of f with multiplicity m. Suppose an element x 0 K satisfies the following initial conditions
ε ( x 0 ) < m / ( n m ) a n d β ( ε ( x 0 ) ) ψ ( ε ( x 0 ) ) ,
where the function ε is defined by (2) and the real functions β and ψ are defined by (32) and (4), respectively. Then the Super-Halley iteration (12) is well defined and converges to ξ with error estimates
| x k + 1 ξ | θ λ 3 k | x k ξ | a n d | x k ξ | θ k λ ( 3 k 1 ) / 2 | x 0 ξ | f o r a l l k 0 ,
where θ = ψ ( ε ( x 0 ) ) and λ = β ( ε ( x 0 ) ) / ψ ( ε ( x 0 ) ) . In addition, if the second inequality in (34) is strict, then the iteration (12) converges at least cubically to ξ.
Proof. 
The proof follows immediately from Lemma 4 and Theorem 2. □

4. Conclusions

Two general theorems (Theorem 1 and Theorem 2) that give an easy algorithm for investigating the local convergence of Picard-type iterative processes in arbitrary normed fields are proven in this paper. Any of the presented theorems is supplied with a priori and a posteriori error estimates. Moreover, the initial conditions of Theorem 1 guarantee the Q-order convergence of the Picard iteration right from the first step. Furthermore, the presented general theorems are applied to study the local convergence of the Super-Halley method for multiple polynomial zeros. Thus, two types of local convergence theorems (Theorem 3 and Theorem 4) with error estimates for this method are proven. The obtained theorems are new even in the case of simple zeros.
Finally, the ideas of this paper can be further developed to obtain general convergence theorems with computationally verifiable initial conditions and error estimates that are of significant practical importance.

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Traub, J. Iterative Methods for the Solution of Equations, 2nd ed.; Chelsea Publishing Company: New York, NY, USA, 1982. [Google Scholar]
  2. Tisdell, C.C. On Picard’s iteration method to solve differential equations and a pedagogical space for otherness. Int. J. Math. Educ. Sci. Technol. 2019, 50, 788–799. [Google Scholar] [CrossRef]
  3. Proinov, P.D. Fixed point theorems for generalized contractive mappings in metric spaces. J. Fixed Point Theory Appl. 2020, 22, 21. [Google Scholar] [CrossRef]
  4. Secelean, N. A New Kind of Nonlinear Quasicontractions in Metric Spaces. Mathematics 2020, 8, 661. [Google Scholar] [CrossRef]
  5. Osada, N. Chebyshev–Halley methods for analytic functions. J. Comput. Appl. Math. 2008, 216, 585–599. [Google Scholar] [CrossRef] [Green Version]
  6. Behl, R.; Kanwar, V.; Kim, Y.I. Higher-order families of multiple root finding methods suitable for non-convergent cases and their dynamics. Math. Model. Anal. 2019, 24, 422–444. [Google Scholar] [CrossRef] [Green Version]
  7. Dehghan, R. A new iterative method for finding the multiple roots of nonlinear equations. Afr. Mat. 2019, 30, 747–753. [Google Scholar] [CrossRef]
  8. Proinov, P.D. General local convergence theory for a class of iterative processes and its applications to Newton’s method. J. Complex. 2009, 25, 38–62. [Google Scholar] [CrossRef] [Green Version]
  9. Proinov, P.D. New general convergence theory for iterative processes and its applications to Newton–Kantorovich type theorems. J. Complex. 2010, 26, 3–42. [Google Scholar] [CrossRef] [Green Version]
  10. Proinov, P.D. General convergence theorems for iterative processes and applications to the Weierstrass root-finding method. J. Complex. 2016, 33, 118–144. [Google Scholar] [CrossRef]
  11. Proinov, P.D. Unified convergence analysis for Picard iteration in n-dimensional vector spaces. Calcolo 2018, 55, 6. [Google Scholar] [CrossRef]
  12. Proinov, P.D. A generalization of the Banach contraction principle with high order of convergence of successive approximations. Nonlinear Anal. 2007, 67, 2361–2369. [Google Scholar] [CrossRef]
  13. Jay, L.O. A note on Q-order of convergence. BIT 2001, 41, 422–429. [Google Scholar] [CrossRef]
  14. Ivanov, S.I. On the convergence of Chebyshev’s method for multiple polynomial zeros. Results Math. 2016, 69, 93–103. [Google Scholar] [CrossRef]
  15. Proinov, P.D.; Ivanov, S.I. On the convergence of Halley’s method for multiple polynomial zeros. Mediterr. J. Math. 2015, 12, 555–572. [Google Scholar] [CrossRef]
  16. Kyncheva, V.K.; Yotov, V.V.; Ivanov, S.I. Convergence of Newton, Halley and Chebyshev iterative methods as methods for simultaneous determination of multiple polynomial zeros. Appl. Numer. Math. 2017, 112, 146–154. [Google Scholar] [CrossRef]
  17. Kanwar, V.; Sharma, K.K.; Behl, R. A new family of Schröder’s method and its variants based on power means for multiple roots of nonlinear equations. Int. J. Math. Educ. Sci. Technol. 2010, 41, 209–226. [Google Scholar] [CrossRef]
  18. Kanwar, V.; Tomar, S.K.; Singh, S.; Kumar, S. Note on Super-Halley Method and its Variants. Tamsui Oxf. J. Inf. Math. Sci 2012, 28, 191–216. [Google Scholar]
  19. Hernández-Verón, M.A. Newton-Raphson’s method and convexity. Univ. Novom Sadu Zb. Rad. Period. Mat. Fak. Ser. Mat. 1992, 22, 159–166. [Google Scholar]
  20. Gutiérrez, J.M.; Hernàndez, M.A. An acceleration of Newton’s method: Super-Halley method. Appl. Math. Comput. 2001, 117, 223–239. [Google Scholar] [CrossRef]
  21. Amat, S.; Busquier, S.; Gutiérrez, J.M. Geometric constructions of iterative functions to solve nonlinear equations. J. Comput. Appl. Math. 2003, 157, 197–205. [Google Scholar] [CrossRef] [Green Version]

Share and Cite

MDPI and ACS Style

Ivanov, S.I. General Local Convergence Theorems about the Picard Iteration in Arbitrary Normed Fields with Applications to Super–Halley Method for Multiple Polynomial Zeros. Mathematics 2020, 8, 1599. https://doi.org/10.3390/math8091599

AMA Style

Ivanov SI. General Local Convergence Theorems about the Picard Iteration in Arbitrary Normed Fields with Applications to Super–Halley Method for Multiple Polynomial Zeros. Mathematics. 2020; 8(9):1599. https://doi.org/10.3390/math8091599

Chicago/Turabian Style

Ivanov, Stoil I. 2020. "General Local Convergence Theorems about the Picard Iteration in Arbitrary Normed Fields with Applications to Super–Halley Method for Multiple Polynomial Zeros" Mathematics 8, no. 9: 1599. https://doi.org/10.3390/math8091599

APA Style

Ivanov, S. I. (2020). General Local Convergence Theorems about the Picard Iteration in Arbitrary Normed Fields with Applications to Super–Halley Method for Multiple Polynomial Zeros. Mathematics, 8(9), 1599. https://doi.org/10.3390/math8091599

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop