Next Article in Journal
Ostrowski-Type Fractional Integral Inequalities: A Survey
Previous Article in Journal
The Magical “Born Rule” and Quantum “Measurement”: Implications for Physics
Previous Article in Special Issue
Improved Higher Order Compositions for Nonlinear Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison between Two Competing Newton-Type High Convergence Order Schemes for Equations on Banach Spaces

1
Department of Computing and Mathematical Sciences, Cameron University, Lawton, OK 73505, USA
2
Mangalmay Institute of Engineering and Technology, Plot No. 8, Knowledge Park II, Greater Noida 201310, India
3
Department of Mathematics, University of Houston, Houston, TX 77204, USA
*
Author to whom correspondence should be addressed.
Foundations 2023, 3(4), 643-659; https://doi.org/10.3390/foundations3040039
Submission received: 11 September 2023 / Revised: 15 October 2023 / Accepted: 20 October 2023 / Published: 30 October 2023
(This article belongs to the Special Issue Iterative Methods with Applications in Mathematical Sciences II)

Abstract

:
We carried out a local comparison between two ninth convergence order schemes for solving nonlinear equations, relying on first-order Fréchet derivatives. Earlier investigations require the existence as well as the boundedness of derivatives of a high order to prove the convergence of these schemes. However, these derivatives are not in the schemes. These assumptions restrict the applicability of the schemes, which may converge. Numerical results along with a boundary value problem are given to examine the theoretical results. Both schemes are symmetrical not only in the theoretical results (formation and convergence order), but the numerical and dynamical results are also similar. We calculated the convergence radii of the nonlinear schemes. Moreover, we obtained the extraneous fixed points for the proposed schemes, which are repulsive and are not part of the solution space. Lastly, the theoretical and numerical results are supported by the dynamic results, where we plotted basins of attraction for a selected test function.

1. Introduction

Let F : Ω E 1 E 2 be an operator, where E 1 and E 2 are Banach spaces and Ω ϕ is an open convex subset of E 1 . F is a Fréchet differentiable operator at each point of Ω . We use this scheme to find a solution x * of the nonlinear operator equation in the form
F ( x ) = 0 .
Newton’s scheme is most commonly used for solving such equations. However, it is only of order two under some conditions ([1,2,3]). It is defined as follows
x n + 1 = x n F ( x n ) 1 F ( x n ) , n = 0 , 1 , 2 , .
There are several papers on the variation or modification of Newton’s scheme in real ([1,4]) and in Banach space ([5,6]) to achieve a higher order.
A plethora of iterative scheme of convergence orders three or higher [7,8,9,10,11,12,13,14] have been developed in the literature provided that E 1 = E 2 = R J . Third convergence order schemes, each involving two linear mapping inversions, are given by Cordero and Torregrosa [15], Homeier [9], Grau- Sanchez et al. [16] and Noor and Waseem [17].
Cordero et al. [18] developed a fourth-order converging scheme involving the computation of two operators, two derivatives (first) and one linear mapping inversion. Sharma and Gupta provided in [19] a fifth convergence order scheme requiring the evaluation of two operators, two derivative (first) and two linear operator inversions. Other such schemes can be found in [20] and the references therein. The convergence orders can be extended further if multistep schemes are introduced involving more than two or three steps. Let us revisit two such schemes as follows:
Let k be a natural number; H 1 and H 2 denote the iteration function of order q. Define the k step schemes for each n = 0 , 1 , 2 , ,
x 0 Ω , y n = x n F ( x n ) 1 F ( x n ) , z n = H 1 ( x n , y n ) , z n ( 1 ) = z n ψ 1 ( x n , y n ) F ( z n ) , . . . . . . . . . . z n ( k 1 ) = z n ( k 2 ) ψ 1 ( x n , y n ) F ( z n ( k 2 ) ) , x n + 1 = z n ( k ) = z n ( k 1 ) ψ 1 ( x n , y n ) F ( z n ( k 1 ) ) ,
and
x 0 Ω , y n = x n F ( x n ) 1 F ( x n ) , z n = H 2 ( x n , y n ) , z n ( 1 ) = z n ψ 2 ( x n , y n ) F ( z n ) , . . . . . . . . . . z n ( k 1 ) = z n ( k 2 ) ψ 2 ( x n , y n ) F ( z n ) ( k 2 ) , x n + 1 = z n ( k ) = z n ( k 1 ) ψ 2 ( x n , y n ) F ( z n ) ( k 1 ) ,
where we define A n = 3 F ( y n ) F ( x n ) , B n = F ( y n ) 1 F ( x n ) F ( y n ) 1 + F ( x n ) 1 , ψ 1 ( x n , y n ) = 1 3 ( 4 A n 1 + F ( x n ) 1 ) and ψ 2 ( x n , y n ) = 1 2 B n .
Notice that both the schemes require k function evaluations and two inverses per complete step. The convergence order q + 3 k is established by Xiao and Yin in [21,22], respectively. However, there exist limitations restricting the utilizations of the aforementioned schemes. Below is a list.
[ l 1 ] —The boundedness and existence of at least F , F ( 2 ) , F ( 3 ) , F ( 4 ) and even higher orders is assumed although only F is on these schemes. There even exist equations on the real line, such as the results in the previous references, which cannot assure the convergence of the schemes to a solution x * Ω of the equation F ( x ) = 0 . Let Ω = [ 0.45 , 1.45 ] .
Define the function F : Ω as
F ( t ) = t 4 + t 5 + 2 t 3 l n t , f o r t 0 , 0 , f o r t = 0
The definition of the function F gives that F ( t ) = 22 24 t + 60 t 2 + 12 l n t . The function F ( t ) is not continuous at t = 0 and t * = 1 Ω solves the equation F ( t ) = 0 . Thus, all the results requiring the existence of at least the third derivative of F cannot assure the convergence of the afforementioned scheme to t * , although they may converge.
[ l 2 ] —A priori estimates of ( x n x * ) are not provided. Thus, it is not known in advance the number of iterations to be performed to satisfy a certain predecided error tolerance.
[ l 3 ] —There is no information about if there are solutions other than x * Ω .
[ l 4 ] —The results are restricted in j .
Notice, however, that the schemes (1) and (2) are extended in the setting of Banach space. Limitations ( l 1 ) ( l 4 ) constitute the motivation for writing this article. Moreover, the novelty of this article is that the items ( l 1 ) ( l 4 ) are addressed positively as follows:
[ l 1 ] —The convergence is established using only F , which is in these schemes, and the idea of generalized continuity [5,18,20].
[ l 2 ] —The number of iterations to reach a predecided error tolerance is calculated.
[ l 3 ] —The isolation of x * is addressed.
[ l 4 ] —The results are valid in Banach space.
The techniques to achieve the aforementioned objectives are demonstrated for some specializations of schemes (1) and (2). However, the technique can be analogously used on the rest of the schemes.
Let k = 2 ;
H 1 ( x n , y n ) = x n 2 3 ( A n 1 + F ( x n ) 1 ) F ( x n ) ; and
H 2 ( x n , y n ) = x n 1 2 ( F ( y n ) 1 + F ( x n ) 1 ) F ( x n ) .
These iteration functions are of convergence order three. Therefore, (1) and (2) specialize to schemes of convergence order 3 + 3 × 2 = 9 , since q = 3 and k = 2 . In particular, these specializations are of schemes (1) and (2), respectively.
The convergence order can be found using the following formula
μ ( C O C ) = ln ( x n + 1 x * ) / ( x n x * ) ln ( x n x * ) / ( x n 1 x * ) ,
or
μ 1 ( A C O C ) = ln ( x n + 1 x n ) / ( x n x n 1 ) ln ( x n x n 1 ) / ( x n 1 x n 2 ) .
These computations do not require the F or x * (in the case of Formula (6)).
We study the ball convergence comparison of the two ninth-order iterative schemes
x 0 Ω y n = x n F ( x n ) 1 F ( x n ) , z n = x n 2 / 3 ( A n 1 + F ( x n ) 1 ) F ( x n ) , w n = z n 1 / 3 ( 4 A n 1 + F ( x n ) 1 ) F ( z n ) , x n + 1 = w n 1 / 3 ( 4 A n 1 + F ( x n ) 1 ) F ( w n ) ,
and
x 0 Ω y n = x n F ( x n ) 1 F ( x n ) , z n = x n 1 2 ( F ( y n ) 1 + F ( x n ) 1 ) F ( x n ) , w n = z n 1 2 B n F ( z n ) , x n + 1 = w n 1 2 B n F ( w n ) ,
where F ( x ) is the Fréchet derivative of operator F at the point x Ω .
Numerical results consist of the comparative study of the proposed schemes along with Newton’s scheme by using the some test functions for nonlinear equations, systems of equations and a boundary value problem. One important characteristic of this work is the comparability of the dynamics of the proposed schemes along with Newton’s scheme for the solution of nonlinear equations.

2. Convergence: Scheme 1

We describe the ball convergence of proposed schemes (3) and (4), which are based on some real functions and positive parameters. Let T = [ 0 , ) .
Suppose:
(a) There exists a nondecreasing and continuous function (NDCF) h 0 from T T such that equation
h 0 ( t ) 1 = 0
has a minimal positive (MP) zero R 0 . Set T 0 = [ 0 , R 0 ) . Consider (NDCF) function h from T 0 T .
(b) Equation
v 1 ( t ) 1 = 0
has an MP zero r 1 ( 0 , R 0 ) , where
v 1 ( t ) = 0 1 h ( ( 1 τ ) t ) d τ 1 h 0 ( t ) .
(c) Equation
p ( t ) 1 = 0 ,
has an MP zero denoted by r p , where
p ( t ) = 1 2 ( 3 h 0 ( v 1 ( t ) t ) + h 0 ( t ) ) .
Set R 1 = m i n { R 0 , R p } and T 1 = [ 0 , R 1 ) . Consider (NDCF) h 1 from T 1 T nondecreasing and continuous.
(d) Equation
v 2 ( t ) 1 = 0 ,
has an MP zero denoted by r 2 ( 0 , R 1 ) , where
v 2 ( t ) = v 1 ( t ) + h 0 ( t ) + h 0 ( v 1 ( t ) t ) 0 1 h 1 ( τ t ) d τ 2 ( 1 h 0 ( t ) ) ( 1 p ( t ) ) .
(e) Equation
h 0 ( v 2 ( t ) t ) 1 = 0
has an MP zero denoted by R 2 . Set R 3 = m i n { R 1 , R 2 } and T 2 = [ 0 , R 3 ) .
(f) Equation
v 3 ( t ) 1 = 0
has an MP zero denoted by r 3 ( 0 , R 3 ) , where
v 3 ( t ) = [ v 1 ( v 2 ( t ) t ) + h 0 ( t ) + h 0 ( v 2 ( t ) t ) 0 1 h 1 ( τ v 2 ( t ) t ) d τ ( 1 h 0 ( t ) ) ( 1 h 0 ( v 2 ( t ) t ) ) + h 0 ( t ) + h 0 ( v 1 ( t ) t ) 0 1 h 1 ( τ v 2 ( t ) t ) d τ ( 1 h 0 ( t ) ) ( 1 p ( t ) ) ] v 2 ( t ) .
(g) Equation
h 0 ( v 3 ( t ) t ) 1 = 0
has an MP zero denoted by R 4 . Set R 5 = m i n { R 3 , R 4 } and T 3 = [ 0 , R 5 ) .
(h) Equation
v 4 ( t ) 1 = 0
has a minimal positive zero denoted by r 4 ( 0 , R 5 ) , where
v 4 ( t ) = [ v 1 ( v 3 ( t ) t ) + h 0 ( t ) + h 0 ( v 3 ( t ) t ) 0 1 h 1 ( τ v 3 ( t ) t ) d τ ( 1 h 0 ( t ) ) ( 1 h 0 ( v 3 ( t ) t ) ) + h 0 ( t ) + h 0 ( v 1 ( t ) t ) 0 1 h 1 ( τ v 3 ( t ) t ) d τ ( 1 h 0 ( t ) ) ( 1 p ( t ) ) ] v 3 ( t ) .
We shall show that parameter ρ , given by
r = m i n { r k } , k = 1 , 2 , 3 , 4 ,
is a convergence radius for scheme (3). By this definition, it follows for all t [ 0 , r ) that
0 h 0 ( t ) < 1 ,
0 h 0 ( v 2 ( t ) t ) < 1 ,
0 h 0 ( v 3 ( t ) t ) < 1 ,
0 p ( t ) < 1
and
0 v k ( t ) < 1 , i = 1 , 2 , 3 , 4 .
Let S ( x * , r ) and S ¯ ( x * , r ) denote the open and closed balls in B 1 , respectively, of center x * and radius r > 0 .
Next, we list the hypotheses in hypothesis ( A ) used in our convergence analysis. They relate to the functions as described previously.
Suppose:
( a 1 ) F : Ω B 2 is continuously differentiable and there exists a simple solution x * of equation F ( x ) = 0 .
( a 2 ) For all x Ω ,
F ( x * ) 1 ( F ( x ) F ( x * ) h 0 ( x x * ) .
Set Ω 0 = Ω S ( x * , ρ 0 ) .
( a 3 ) For all x , y Ω 0 ,
F ( x * ) 1 ( F ( y ) F ( x * ) h ( y x * ) .
F ( x * ) 1 F ( x ) h 1 ( x x * ) .
( a 4 ) S ¯ ( x * , r ˜ ) Ω , where r ˜ > 0 is to be determined.
( a 5 ) There exists r * r ˜ such that
0 1 h 0 ( τ r * ) d τ < 1 .
Set
Ω 1 = Ω S ¯ ( x * , r * ) .
Based on hypothesis ( A ) and the developed notation, we show the local convergence result for scheme (3).
Theorem 1.
Under the hypotheses in ( A ) for r ˜ = r , further suppose x 0 S ( x * , r ) { x * } . Then, the sequence generated by scheme (3) starting at x 0 is well defined, remains in S ( x * , r ) { x * } and converges at x * . Moreover, the following estimates hold true
y n x * v 1 ( x n x * ) x n x * x n x * < r ,
z n x * v 2 ( x n x * ) x n x * x n x * ,
w n x * v 3 ( x n x * ) x n x * x n x * ,
x n + 1 x * v 4 ( x n x * ) x n x * x n x * ,
where functions v k are defined earlier and r is given by (16). Furthermore, x * is unique as a solution of equation F ( x ) = 0 in the domain Ω 1 , given in ( a 5 ) .
Proof. 
We shall prove items (22)–(25) using Induction. Consider x S ( x * , r ) x * . Then, using ( a 1 ) and ( a 2 ) , we have
F ( x * ) 1 ( F ( x ) F ( x * ) ) h 0 x x * h 0 ( r ) < 1 .
Using the Banach lemma on invertible operators [12] and the preceding inequality, F ( x ) 1 exists so that
F ( x ) 1 F ( x * ) 1 1 h 0 ( x x * ) .
In particular, if x = x 0 , F ( x 0 ) 1 exists, since x 0 S ( x * , r ) . Then, the iterate y 0 exists by the first substep of scheme (3). We can write
y 0 x * = x 0 x * F ( x 0 ) 1 F ( x 0 ) F ( x 0 ) 1 F ( x * ) × 0 1 F ( x * ) 1 F ( x * + τ ( x 0 x * ) ) F ( x 0 ) ( x 0 x * ) d τ
So, by conditions ( a 1 ) , ( a 3 ) and (16) and (26),
y 0 x * 0 1 h ( ( 1 τ ) ) x 0 x * d τ x 0 x * 1 h 0 ( x 0 x * ) v 1 ( x 0 x * ) x 0 x * x 0 x * < r .
Hence, the iterate y 0 S ( x * , r ) and (22) is valid for n = 0. We should show A 0 is invertible, so z 0 , v 0 and x 1 exist by scheme (3) for n = 0. Indeed, we have by ( a 2 ) , (17), (26) and (27)
( 2 F ( x * ) ) 1 ( A 0 2 F ( x * ) ) ( 2 F ( x * ) ) 1 ( 3 ( F ( y 0 ) F ( x * ) ) + ( F ( x * ) F ( x 0 ) ) 1 2 [ 3 F ( x * ) ) 1 ( F ( y 0 ) F ( x * ) ) + F ( x * ) ) 1 ( F ( x 0 ) F ( x * ) ) ] 1 2 ( h 0 ( x 0 x * ) + 3 h 0 ( y 0 x * ) ) 1 2 ( h 0 ( x 0 x * ) + 3 h 0 ( v 1 ( x 0 x * ) x 0 x * ) ) p ( x 0 x * ) p ( r ) < 1 .
Thus, the linear operator A 0 is invertible and
A 0 1 F ( x * ) 1 2 ( 1 p ( x 0 x * ) ) .
Then, using the second substep of schemes (3), (11), (21) (for i = 2), (23) (for x = x 0 ) , (26) and (28), we first have
z 0 x * = x 0 x * F ( x 0 ) 1 F ( x 0 ) + ( F ( x 0 ) 1 2 3 A 0 1 2 3 F ( x 0 ) 1 ) F ( x 0 ) = x 0 x * F ( x 0 ) 1 F ( x 0 ) + 1 3 F ( x 0 ) 1 ( 3 F ( y 0 ) F ( x 0 ) 2 F ( x 0 ) ) A 0 1 F ( x 0 ) = ( x 0 x * F ( x 0 ) 1 F ( x 0 ) ) + ( F ( x 0 ) 1 ( F ( y 0 ) F ( x 0 ) ) A 0 1 F ( x 0 ) .
So, we obtain, by using also the triangle inequality,
z 0 x * [ v 1 ( x 0 x * ) + ( h 0 ( y 0 x * ) + h 0 ( x 0 x * ) ) 0 1 h 1 ( τ x 0 x * ) d τ 2 ( 1 h 0 ( x 0 x * ) ) ( 1 p ( x 0 x * ) ) ] x 0 x * v 2 ( x 0 x * ) x 0 x * x 0 x * .
Hence, the iterate z 0 S ( x * , r ) and (23) is valid for n = 0 . By the third substep of scheme (3) for n = 0 , we write
w 0 x * = z 0 x * F ( z 0 ) 1 F ( z 0 ) + ( F ( z 0 ) 1 4 3 A 0 1 1 3 F ( x 0 ) 1 ) F ( z 0 ) = ( z 0 x * F ( z 0 ) 1 F ( z 0 ) ) + F ( z 0 ) 1 ( F ( x 0 ) F ( z 0 ) ) F ( x 0 ) 1 F ( z 0 ) + 2 F ( x 0 ) 1 × ( F ( y 0 ) F ( x 0 ) ) A 0 1 F ( z 0 ) .
Then, using (12), (13), (16) (for m = 3 ), (21) (for x = z 0 ) and (24)–(30), we have
w 0 x * [ v 1 ( z 0 x * ) + ( h 0 ( z 0 x * ) + h 0 ( x 0 x * ) ) 0 1 h 1 ( τ z 0 x * ) d τ ( 1 h 0 ( z 0 x * ) ) ( 1 h 0 ( x 0 x * ) ) + ( h 0 ( y 0 x * ) + h 0 ( x 0 x * ) ) 0 1 h 1 ( τ z 0 x * ) d τ ( 1 h 0 ( x 0 x * ) ) ( 1 p ( x 0 x * ) ) ] ( z 0 x * ) v 3 ( x 0 x * ) x 0 x * x 0 x * < r .
So, the iterate w 0 S ( x * , r ) and (24) holds true for n = 0. Similarly, if we exchange the role of z 0 with w 0 we first obtain
x 1 x * = w 0 x * F ( w 0 ) 1 F ( w 0 ) + F ( w 0 ) 1 ( F ( x 0 ) F ( w 0 ) ) F ( x 0 ) 1 + 2 F ( x 0 ) 1 ( F y 0 ) F ( x 0 ) A 0 1 F ( w 0 ) .
Hence, we see that
x 1 x * [ v 1 ( w 0 x * ) + ( h 0 ( w 0 x * ) + h 0 ( ( x 0 x * ) ) 0 1 h 1 τ w 0 x * ) d τ ( 1 h 0 ( x 0 x * ) ) ( 1 h 0 ( w 0 x * ) ) + ( h 0 ( y 0 x * ) + h 0 ( ( x 0 x * ) 0 1 h 1 ( τ w 0 x * ) d τ ( 1 h 0 ( x 0 x * ) ) ( 1 p 0 ( ( x 0 x * ) ) ) ] w 0 x * v 4 ( x 0 x * ) x 0 x * x 0 x * < r .
Replace x 0 , y 0 , z 0 , v 0 and x 1 with x m , y m , z m , w m and x m + 1 in the previous calculations to complete the induction for items (22)–(25). It then follows by the estimation
x m + 1 x * q ( x m x * ) < r ˜
with q = v 4 ( x 0 x * ) [ 0 , 1 ) , that lim m x m = x * and x m + 1 S ( x * , r ˜ ) . Set M = 0 1 F ( x * + τ ( x * * x * ) ) d τ for some x * * Ω 1 with F ( x * * ) = 0 . Then, by ( a 1 ) , ( a 2 ) and ( a 5 ) , we see in turn that, in view of ( a 2 ) and ( a 6 ) ,
F ( x * ) 1 ( M F ( x * ) ) 0 1 h 0 ( ( 1 τ ) x * x * * ) d τ 0 1 h 0 ( τ r * ) d τ < 1 .
Therefore, from the invertability of M and the estimate 0 = F ( x * * ) F ( x * ) = M ( x * * x * ) , we conclude that x * * = x * . □

3. Convergence: Scheme 2

In a similar way, we provide the local convergence analysis for scheme (4). In this case, the functions v 1 ¯ , v 2 ¯ , v 3 ¯ and v 4 ¯ are defined as follows
v 1 ¯ = v 1 ,
v 2 ¯ ( t ) = v 1 ¯ ( t ) + h 0 ( t ) + h 0 ( v 1 ¯ ( t ) t ) 0 1 h 1 ( τ t ) d τ 2 ( 1 h 0 ( t ) ) ( 1 h 0 ( v 1 ¯ ( t ) t ) ) ,
v 3 ¯ ( t ) = [ v 1 ¯ ( v 2 ¯ ( t ) t ) + h 0 ( v 2 ¯ ( t ) t ) + h 0 ( v 1 ¯ ( t ) t ) 0 1 h 1 ( τ v 2 ¯ ( t ) t ) d τ ( 1 h 0 ( v 2 ¯ ( t ) t ) ) ( 1 h 0 ( v 1 ¯ ( t ) t ) ) + 1 2 h 0 ( v 1 ¯ ( t ) t ) + h 0 ( t ) 0 1 h 1 τ v 2 ¯ ( t ) t d τ ( 1 h 0 ( v 1 ¯ ( t ) t ) ) ( 1 h 0 ( t ) ) + 1 2 h 0 ( v 1 ¯ ( t ) t ) + h 0 ( t ) 0 1 h 1 τ v 2 ¯ ( t ) t d τ ( 1 h 0 ( v 1 ¯ ( t ) t ) ) 2 ] v 2 ¯ ( t ) ,
v 4 ¯ ( t ) = [ v 1 ¯ ( v 3 ¯ ( t ) t ) + ( h 0 ( v 1 ¯ ( t ) t ) + h 0 ( v 3 ¯ ( t ) t ) ) 0 1 h 1 ( τ v 3 ¯ ( t ) t ) d τ ( 1 h 0 ( v 3 ¯ ( t ) t ) ) ( 1 h 0 ( v 1 ¯ ( t ) t ) ) + 1 2 ( h 0 ( t ) + h 0 ( v 1 ¯ ( t ) t ) ) 0 1 h 1 ( τ v 3 ¯ ( t ) t ) d τ ( 1 h 0 ( t ) ) ( 1 h 0 ( v 1 ¯ ( t ) t ) ) + 1 2 ( h 0 ( t ) + h 0 ( v 1 ¯ ( t ) t ) ) 0 1 h 1 ( τ v 3 ¯ ( t ) t ) d τ ( 1 h 0 ( v 1 ¯ ( t ) t ) ) 2 ] v 3 ¯ ( t ) .
Define r ¯ by
r ¯ = m i n { r k ¯ } ,
where we suppose that the MP zero r k ¯ exits for equations
v k ¯ ( t ) 1 = 0
Then, we have, as in scheme (4),
z n x * = x n x * F ( x n ) 1 F ( x n ) 1 2 ( F ( y n ) 1 F ( x n ) 1 ) F ( x n ) = x n x * F ( x n ) 1 F ( x n ) 1 2 F ( y n ) 1 ( F ( x n ) F ( y n ) ) F ( x n ) 1 F ( x n ) .
Therefore, we have
z n x * [ v 1 ¯ ( x n x * ) + ( h 0 ( y n x * ) + h 0 ( x n x * ) ) 0 1 h 1 ( τ ( x n x * ) ) d τ 2 ( 1 h 0 ( x n x * ) ) ( 1 h 0 ( y n x * ) ) ] x n x * v 2 ¯ ( x n x * ) x n x * x n x * < r ¯ .
Moreover, we can write
w n x * = z n x * F ( z n ) 1 F ( z n ) + F ( z n ) 1 1 2 F ( y n ) 1 F ( x n ) F ( y n ) 1 1 2 F ( x n ) 1 F ( z n ) = z n x * F ( z n ) 1 F ( z n ) + [ ( F ( z n ) 1 + F ( y n ) 1 ) + 1 2 ( F ( y n ) 1 F ( x n ) 1 ) + 1 2 F ( y n ) 1 ( F ( y n ) F ( x n ) ) F ( y n ) 1 ] F ( z n ) ,
leading to
w n x * [ v 1 ¯ ( z n x * ) + ( h 0 ( z n x * ) + h 0 ( y n x * ) ) 0 1 w 1 ( τ ( z n x * ) ) d τ ( 1 h 0 ( z n x * ) ) ( 1 h 0 ( y n x * ) ) + 1 2 ( h 0 ( y n x * ) + h 0 ( x n x * ) ) 0 1 w 1 ( τ ( z n x * ) ) d τ ( 1 h 0 ( y n x * ) ) ( 1 h 0 ( x n x * ) ) + 1 2 ( h 0 ( y n x * ) + h 0 ( x n x * ) ) 0 1 h 1 ( τ ( z n x * ) ) d τ ( 1 h 0 ( y n x * ) ) 2 ] z n x * v 3 ¯ ( x n x * ) x n x * x n x * < r ¯ .
Next, by the third substep of scheme (4),
x n + 1 x * = w n x * F ( w n ) 1 F ( w n ) + F ( w n ) 1 1 2 F ( y n ) 1 F ( x n ) F ( y n ) 1 1 2 F ( x n ) 1 F ( w n ) = w n x * F ( w n ) 1 F ( w n ) + [ ( F ( w n ) 1 F ( y n ) 1 ) + 1 2 ( F ( y n ) 1 F ( x n ) 1 ) + 1 2 F ( y n ) 1 ( F ( y n ) F ( x n ) ) F ( y n ) 1 ] F ( w n ) .
Hence, we have
x n + 1 x * [ v 1 ¯ ( w n x * ) + ( h 0 ( w n x * ) + h 0 ( y n x * ) ) 0 1 h 1 ( τ ( w n x * ) ) d τ ( 1 h 0 ( w n x * ) ) ( 1 h 0 ( y n x * ) ) + 1 2 ( h 0 ( y n x * ) + h 0 ( x n x * ) ) 0 1 h 1 ( τ ( w n x * ) ) d τ ( 1 h 0 ( y n x * ) ) ( 1 h 0 ( x n x * ) ) + 1 2 ( h 0 ( y n x * ) + h 0 ( x n x * ) ) 0 1 h 1 ( τ ( w n x * ) ) d τ ( 1 h 0 ( y n x * ) ) 2 ] × ( w n x * ) v 4 ¯ ( x n x * ) x n x * x n x * < r ¯ .
Hence, we arrive at the local convergence result for scheme (4) corresponding to scheme (3).
Theorem 2.
Under hypotheses A with r ˜ = r ¯ , choose x 0 S ( x * , r ¯ ) x * . Then, the conclusions of Theorem 1 hold for scheme (4) but with v k ¯ and r ¯ replacing v k and r, respectively.

4. Numerical Results

In this section, we study the efficiency of iterative proposed schemes (3) and (4). We performed a comparative study of schemes (3) and (4) along with the classical Newton’s scheme.
Example 1.
Let E 1 , E 2 = R , Ω = ( 1 , 1 ) and F : Ω be a function defined by
F ( x ) = e x 1 , x Ω .
Then, F is Fréchet differentiable and its Fréchet derivative F ( x ) at any point x Ω is given by
F ( x ) = e x .
We computed the numerical results with the help of MATLAB 2007 and the stopping criterion used for the computation is | x n + 1 x * | + | f ( x n + 1 ) | < 10 14 . The initial approximation is 1.0 and approximate solution is 0. The numerical solutions for example 1 using the second-order Newton scheme and the proposed ninth-order scheme are given in Table 1. The numerical results in Table 1 reveal that the proposed schemes (3) and (4) perform with the same number of iterations with a little advantage to proposed scheme (3).
Example 2.
Let E 1 = E 2 = , Ω = ( 2 , 2 ) and F : Ω be an operator defined by
F ( x ) = x 3 1 , x Ω .
Then, F is Fréchet differentiable and its Fréchet derivative F ( x ) at any point x Ω is given by
F ( x ) = 3 x 2 .
We computed the numerical results with the help of MATLAB 2007 and the stopping criterion is | x n + 1 x * | + | f ( x n + 1 ) | < 10 14 . The initial approximation is 3.5 and the approximate solution is 1.0. The numerical solutions for example 2 using different schemes are given in Table 2. The numerical results in Table 2 reveal that the proposed schemes (3) and (4) perform with the same number of iterations with a little advantage to proposed scheme (4).
Example 3.
Let Ω = E 1 = E 2 = 2 . Consider an operator F : 2 2 defined by
F ( x ) = x 2 y 0.2 , x + y 2 0.3 , x = ( x , y ) 2 .
The starting vector is [0.1, 0.1] and the approximate solution is [−0.2860321636288604, −0.11818560136979284]. The numerical solutions for example 3 using different methods are given in Table 3. The numerical results show that the proposed scheme (3) converges at the solution in fewer iterations in comparison with scheme (4).
Example 4.
Consider the following boundary problem
x + 3 x x = 0 , x ( 0 ) = 0 , x ( 2 ) = 1 .
We take t 0 = 0 < t 1 < t 2 < t 3 < , < t n 1 < t n = 2 , t i + 1 = t i + h , and h = 2 n . Here, x 0 = x ( t 0 ) = 0 , x 1 = x ( t 1 ) , x 2 = x ( t 2 ) , x 3 = x ( t 3 ) , , x n 1 = x ( t n 1 ) and x n = x ( t n ) = 1 .
We discretize the above problem by using the central difference schemes for the first and second-order derivatives, i.e.,
x i = x i 1 2 x i + x i + 1 h 2 , i = 1 , 2 , 3 , , n 1 ,
x i = x i + 1 x i 1 2 h , i = 1 , 2 , 3 , , n 1 ,
x i = x i + 1 + x i 1 2 , i = 1 , 2 , 3 , , n 1 .
Thus, we have an ( n 1 ) × ( n 1 ) nonlinear system:
4 ( x i 1 2 x i + x i + 1 ) + 3 h ( x i + 1 2 x i 1 2 ) = 0 , i = 1 , 2 , 3 , , n 1 .
Next, we solve the above problem for n = 3 with the proposed scheme along with Newton’s scheme using the initial approximations x 0 = [ 0.1 , 0.1 ] . The solution [0.7321436796857499, 0.9820632479169275] of the problem is shown in Table 4 with x = [ x 1 , x 2 ] and F = [ f , g ] . From Table 4, it is confirmed that scheme (4) converges at the root in fewer iterations and, hence, scheme (4) is better than (3). Notice that accelerated methods are vital for multivariable problems (see, e.g., [23]).

5. Extraneous Fixed Points

The Newton-like iterative schemes described in earlier sections may be viewed as fixed-point iteration:
x n + 1 = x n E f ( x n ) f ( x n ) f ( x n ) , n = 0 , 1 , 2 , .
Clearly, the root x * of f ( x ) = 0 is a fixed point in the scheme. If the right side of (37) also vanishes at some points ξ x * for E f ( ξ ) = 0 , then ξ is also a fixed point in the scheme. These fixed points are known as extraneous fixed points (see [2]). Now, we describe the extraneous fixed points of some Newton-like scheme for z 3 1 .
Remark 1.
Newton’s scheme does not have any extraneous fixed points as for Newton’s scheme, E f ( x n ) = 1 .
Theorem 3.
The proposed scheme given by Equation (3) has 106 extraneous fixed points.
Proof. 
For proposed scheme (3), E f ( x n ) is given by the following equation 1 / ( 1 + z 3 ) z 2 ( 6 z + ( 2 + 6 z 3 8 z 9 ) / ( 3 ( z 2 + 4 z 5 + z 8 ) ) ( ( 1 + z 3 ) 3 ( 1 + 4 z 3 + 13 z 6 ) ( 8 204 z 3 1641 z 6 5438 z 9 6672 z 1 2 618 z 1 5 + z 1 8 ) ) / ( 2187 z 8 ( 1 + 4 z 3 + z 6 ) 4 ) ( ( 1 + 2 z 3 ) 2 ( 1 + ( 8 212 z 3 419 z 6 + 20,760 z 9 + 218,580 z 1 2 + 1,011,516 z 1 5 + 2,454,006 z 1 8 + 2,937,180 z 2 1 + 1,389,756 z 2 4 + 428,120 z 2 7 + 43,061 z 3 0 + 716 z 3 3 ) 3 / ( 282,429,536,481 z 2 4 ( 1 + 4 z 3 + z 6 ) 1 2 ) ) ) / ( 3 z 2 ( 1 + 4 z 3 + z 6 ) ) . In this equation, the numerator is of degree 106 and, hence, proposed scheme (3) has 106 extraneous fixed points.
z = 2.0563466004961662355811893 0.3095651829791553623908596 I , z = 2.0563466004961662355811893 + 0.3095651829791553623908596 I , z = 1.983769642614411418048416 0.116719374008501542977146 I , z = 1.983769642614411418048416 + 0.116719374008501542977146 I , z = 1.908911702325388970327643 , z = 1.2324703408496331888011102 + 1.1027252320147195003269012 I , z = 1.2962646128352097743260829 1.6260658035258724274991578 I , z = 1.2962646128352097743260829 + 1.6260658035258724274991578 I .
These fixed points are repulsive since the magnitude of the derivative at these points is > 1 . □
Theorem 4.
The proposed scheme given by Equation (4) has 129 extraneous fixed points.
Proof. 
For proposed scheme (4), E f ( x n ) given by the following equation. ( 1 / ( 2 ( 1 + z 3 ) ) ) 3 z 2 ( 4 z + ( 1 + 3 z 3 + 9 z 6 13 z 9 ) / ( 3 ( z + 2 z 4 ) 2 ) ( 1 / ( 3 z 2 ) + ( 27 z 1 0 ) / ( 1 + 2 z 3 ) 4 ) ( 1 + ( 1 + 9 z 3 + 33 z 6 + 11 z 9 ) 3 / ( 216 ( z + 2 z 4 ) 6 ) ) + ( 1 / ( 3 z 2 ) + ( 27 z 1 0 ) / ( 1 + 2 z 3 ) 4 ) ( 1 + ( z + ( 1 + 3 z 3 + 9 z 6 13 z 9 ) / ( 6 ( z + 2 z 4 ) 2 ) 1 / 2 ( 1 / ( 3 z 2 ) + ( 27 z 1 0 ) / ( 1 + 2 z 3 ) 4 ) ( 1 + ( 1 + 9 z 3 + 33 z 6 + 11 z 9 ) 3 / ( 216 ( z + 2 z 4 ) 6 ) ) ) 3 ) ) . In this equation, the numerator is of degree 129 and, hence, proposed scheme (4) has 129 extraneous fixed points.
z = 1.3906494250828819032610349 , z = 1.3089391963407268999831171 0.3564555150226694666922243 I , z = 1.3089391963407268999831171 + 0.3564555150226694666922243 I , z = 1.11019566319196274717 , z = 1.09582782332284367467 0.00666195819164429163 I , z = 1.09582782332284367467 + 0.00666195819164429163 I z = 0.7824279535043961805671 0.6482932690512660338217 I , z = 0.7824279535043961805671 + 0.6482932690512660338217 I , z = 0.9631691294990608108553360 0.9553468385289218973982227 I , z = 0.9631691294990608108553360 + 0.9553468385289218973982227 I .
Since the magnitude of the derivative at these points is > 1 , these fixed points are repulsive. □
Remark 2.
As the magnitude of the derivative at these points is > 1 , these fixed points are repulsive. These fixed points can be seen in the basins of attraction plot for example 3 ( z 3 1 ), Figure 2 (see Section 6.2).

6. Dynamics of Scheme

We studied the dynamics and fractal patterns of the functions in example 1 ( F ( z ) = e x p ( z ) 1 ) and example 2 ( F ( z ) = z 3 1 ) by using proposed iterative schemes (3) and (4) along with Newton’s scheme. The dynamical analysis help us to study the convergence and stability of the schemes (see [6]).

6.1. For Example 1

We considered a square R × R = [ 4.0 , 4.0 ] × [ 4.0 , 4.0 ] of 500 × 500 points with a tolerance- | f ( z n ) | < 5 × 10 2 and a maximum of 11 iterations to study the dynamics of the function F ( z ) = e x p ( z ) 1 . We described the basins of attraction with a fixed color for the second-order Newton scheme, the ninth-order proposed scheme (3) and the ninth-order proposed scheme (4) for finding complex roots of the above-mentioned functions (Figure 1).
1.
The basins for all the iterative schemes contain a fractal Julia set and the basins of all the schemes look almost similar.
2.
The basins of attraction of the second-order Newton scheme contain a higher number of orbits and are less dark in comparison with the ninth-order schemes.
3.
Again, the Fatou set with blue color shows the basins of the schemes. The blue-colored area shows that the proposed scheme (4) contains the Fatou set with bigger and darker orbits.

6.2. For Example 2

We plotted the fractal pattern graph of example 2 ( F ( z ) = z 3 1 ) for the different iterative schemes under the same previous conditions with a different color fixed to each root of the basins of attraction.
The basins of attraction for schemes to find the complex roots of example 2 ( F ( z ) = z 3 1 ) are shown in Figure 2. We can see that there are no extraneous fixed points for the second-order Newton scheme, which is in agreement with the findings of the section Extraneous Fixed Points. Again, there are 106 extraneous fixed points for the proposed ninth-order scheme (3) and 129 extraneous fixed points for the proposed ninth-order scheme (4). Since, the magnitude of the derivative at these points is > 1 , these fixed points are repulsive. Thus, we see that scheme (3) is better in terms of extraneous fixed points.

7. Convergence Radii

In the next two examples, we compute the convergence radii for proposed schemes (3) and (4).
Example 5.
By the example in the Section 1, conditions (A) are satisfied if we choose h 0 ( t ) = h ( t ) = 97 t and h 1 ( t ) = 2 . Then, by solving the convergence radii, we have
r 1 = 0.010256678622334536
r 2 = 0.36382236284223446 , 0.0103087361577125
r 3 = 0.10153752845355392 , 0.22527857007361768 , 1.3993484775051468
r 4 = 0.11929997098554257 , 0.21910981171003618
and
r 1 ¯ = 0.010256678622334536
r 2 ¯ = 0.08299300112444563
r 3 ¯ = 0.5291548166116622
r 4 ¯ = 0.5420088529012402
Example 6.
Let B 1 and B 2 be the space of continuous functions on the interval [ 0 , 1 ] with the max-norm and Ω = S ¯ ( 0 , 1 ) . Define F : Ω B 2 by
F ( ψ ) ( x ) = ψ ( x ) 0 1 x τ ψ ( τ ) 3 d τ .
Then, the Fréchet derivative is given as
F ( ψ ( μ ) ) ( x ) = μ ( x ) 3 0 1 x τ ψ ( τ ) 2 μ d τ
for all μ Ω . Then, conditions (A) are satisfied if we choose h 0 ( t ) = 1.5 t , h ( t ) = 3 t and h 1 ( t ) = 2 .
Then, we have
r 1 = 0.4574271077563381
r 2 = 0.7496405298576232 , 0.46028836028846015
r 3 = 0.7444569229567777 , 0.5894382155567577 , 0.501131729398928
r 4 = 0.6180132678614112
and
r 1 ¯ = 0.4574271077563381
r 2 ¯ = 0.8074036383761246 , 0.7071105546988605 , 0.5879931572730634 , 0.486927560161037
r 3 ¯ = 0.5291548166116622
r 4 ¯ = 0.5420088529012402

8. Conclusions

We developed two ninth-order Newton-like schemes for solving nonlinear equations in Banach space and discussed the ball convergence analysis for both of them. We performed a local convergence comparison with the use of only the first-order derivative. The study is used to prove the convergence for scheme (3) and scheme (4) under weak conditions, extending the usage of these schemes. Earlier work relies on the existence and boundedness of F 4 , which is not in these schemes. Thus, these results are not applicable in cases where these hypothesis are violated. However, these schemes may converge. We checked the theoretical results by using the numerical examples along with the boundary value problem. We also examined the numerical results with the basins of attraction for some selected examples. All the results (theoretical, numerical, dynamical) are generative for the advanced study of higher-order Newton-like schemes. The new approach is applicable in other schemes. This is revealing for our future research.

Author Contributions

Conceptualization, I.K.A., M.K.S. and S.R.; methodology, I.K.A., M.K.S. and S.R.; software, I.K.A., M.K.S. and S.R.; validation, I.K.A., M.K.S. and S.R.; formal analysis, I.K.A., M.K.S. and S.R.; investigation, I.K.A., M.K.S. and S.R.; resources, I.K.A., M.K.S. and S.R.; data curation, I.K.A., M.K.S. and S.R.; writing—original draft preparation, I.K.A., M.K.S. and S.R.; writing—review and editing, I.K.A., M.K.S. and S.R.; visualization, I.K.A., M.K.S. and S.R.; supervision, I.K.A., M.K.S. and S.R.; project administration, I.K.A., M.K.S. and S.R.; funding acquisition, I.K.A., M.K.S. and S.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Singh, M.K. A Six-order variant of Newton’s method for solving non linear equations. Comput. Meth. Sci. Technol. 2009, 15, 185–193. [Google Scholar] [CrossRef]
  2. Vrscay, E.R.; Gilbert, W.J. Extraneous fixed points, basin boundaries and chaotic dynamics for Schroder and Konig rational iteration functions. Numer. Math. 1987, 52, 1–16. [Google Scholar] [CrossRef]
  3. Wang, K.; Kou, J.; Gu, C. Semilocal convergence of a sixth-order jarrat method in Banach spaces. Numer. Algorithms 2011, 57, 441456. [Google Scholar] [CrossRef]
  4. Singh, M.K.; Singh, A.K. Variant of Newton’s Method Using Simpson’s 3/8th Rule. Int. J. Appl. Comput. 2020, 6, 20. [Google Scholar] [CrossRef]
  5. Argyros, I.K. Convergence and Applications of Newton-Type Iterations; Springer: Berlin, Germany, 2008. [Google Scholar]
  6. Argyros, I.K.; Magrenan, A.A. Iterative Methods and Their Dynamics with Applications: A Contemporary Study; Taylor and Francis: Boca Raton, FL, USA, 2017. [Google Scholar]
  7. Cordero, A.; Torregrosa, J.R. Variants of Newton’s method using fifth-order quadrature formulas. Appl. Math. Comput. 2007, 190, 686–698. [Google Scholar] [CrossRef]
  8. Grau-Sánchez, M.; Grau, A.; Noguera, M. On the computational efficiency index and some iterative methods for solving systems of nonlinear equations. J. Comput. Appl. Math. 2011, 236, 1259–1266. [Google Scholar] [CrossRef]
  9. Homeier, H.H.H. A modified Newton method with cubic convergence: The multivariable case. J. Comput. Appl. Math. 2004, 169, 161–169. [Google Scholar] [CrossRef]
  10. Homeier, H.H.H. On Newton-type methods with cubic convergence. J. Comput. Appl. Math. 2005, 176, 425–432. [Google Scholar] [CrossRef]
  11. Jarratt, P. Some fourth order multipoint iterative methods for solving equations. Math. Comput. 1996, 20, 434–437. [Google Scholar] [CrossRef]
  12. Kantorovich, L.V.; Akilov, G.P. Funtional Analysis; Pergamon Press: Oxford, UK, 1982. [Google Scholar]
  13. Xiao, X.Y.; Yin, H.W. Accelerating the convergence speed of iterative methods for solving nonlinear systems. Appl. Math. Compt. 2018, 333, 8–19. [Google Scholar] [CrossRef]
  14. Xiao, X.Y.; Yin, H. Achieving higher order of convergence for solving systems of nonlinear equations. Appl. Math. Compt. 2017, 311, 251–261. [Google Scholar] [CrossRef]
  15. Cordero, A.; Torregrosa, J.R. Variants of Newton’s method for functions of several variables. Appl. Math. Comput. 2006, 183, 199–208. [Google Scholar] [CrossRef]
  16. Sharma, J.R.; Guha, R.K.; Sharma, R. An efficient fourth order weighted-Newton method for systems of nonlinear equations. Numer. Algorithms 2013, 62, 307–323. [Google Scholar] [CrossRef]
  17. Noor, M.A.; Waseem, M. Some iterative methods for solving a system of nonlinear equations. Comput. Math. Appl. 2009, 57, 101–106. [Google Scholar] [CrossRef]
  18. Cordero, A.; Martínez, E.; Torregrosa, J.R. Iterative methods of order four and five for systems of nonlinear equations. J. Comput. Appl. Math. 2009, 231, 541–551. [Google Scholar] [CrossRef]
  19. Sharma, J.R.; Gupta, P. An efficient fifth order method for solving systems of nonlinear equations. Comput. Math. Appl. 2014, 67, 591–601. [Google Scholar] [CrossRef]
  20. Argyros, I.K. Computational Theory of Iterative Methods, Series: Studies in Computational Mathematics; Elsevier Publishing Company: New York, NY, USA, 2007. [Google Scholar]
  21. Xiao, X.Y.; Yin, H.W. A new class of methods with higher order of convergence for solving systems of nonlinear equations. Appl. Math. Comput. 2015, 264, 300–309. [Google Scholar] [CrossRef]
  22. Xiao, X.Y.; Yin, H.W. Increasing the order of convergence for iterative methods to solve nonlinear systems. Calcolo 2016, 53, 285–300. [Google Scholar] [CrossRef]
  23. Sandor, K. Optimizing single Slater determinant for electronic Hamiltonian with Lagrange multipliers and Newton-Raphson methods as an alternative to ground state calculations via Hartree-Fock self consistent field. AIP Conf. Proc. 2019, 2116, 450030. [Google Scholar] [CrossRef]
Figure 1. Basin of attraction for e x p ( z ) 1 with different schemes.
Figure 1. Basin of attraction for e x p ( z ) 1 with different schemes.
Foundations 03 00039 g001
Figure 2. Basin of attraction for f 2 = z 3 1 with different schemes.
Figure 2. Basin of attraction for f 2 = z 3 1 with different schemes.
Foundations 03 00039 g002
Table 1. Comparison of different schemes for example 1.
Table 1. Comparison of different schemes for example 1.
SchemeNx f ( x )
11.000000001.71828182845905
20.367879441171440.44466786100977
30.060080068726790.06192156984951
Newton scheme40.001769199442640.00177076539934
5 1.564110789898428 × 10 6 1.564112013019425  × 10 6
61.223321565989411  × 10 12 1.223243728531998  × 10 12
77.783745890945912  × 10 17 0.00000000
11.000000001.71828182845905
Proposed scheme (3)2−8.566001524658931  × 10 4 −8.562333752899498  × 10 4
38.017605522905244  × 10 18 0.00000000
11.000000001.71828182845905
Proposed scheme (4)20.001806631404570.00180826434631
33.652768952695270  × 10 17 0.00000000
Table 2. Comparison of different schemes for example 2.
Table 2. Comparison of different schemes for example 2.
SchemeNx f ( x )
13.5000000041.87500000000000
22.3605442176870712.15335132155504
31.633517254842433.35884252127395
41.213931306812980.78888464195259
Newton scheme51.035486455037460.11028191827017
61.001202239852960.00361105743855
71.000001443067224.329207893061238  × 10 6
81.000000000002086.247669048775606  × 10 12
91.00000000000.00000000
13.5000000041.87500000000000
30.80025523102213−0.48750980007799
Proposed scheme (3)31.000003822894831.146872833501789  × 10 5
41.00000000000.00000000
13.5000000041.87500000000000
31.252062089857640.96280700078081
Proposed scheme (4)31.000002649459617.948399889379232  × 10 6
41.00000000000.00000000
Table 3. Comparison of different schemes for example 3.
Table 3. Comparison of different schemes for example 3.
SchemeNxy f ( x , y ) g ( x , y )
Newton scheme10.150000.15000−0.3275−0.4275
2−0.427747−0.3508240.3337920.250825
3−0.38617−0.05260160.001728690.0889367
4−0.289566−0.1254840.009332370.00531187
5−0.286091−0.1181640.00001207410.0000535827
6−0.286032−0.118186 3.44119 × 10 9 4.61865 × 10 10
7−0.286032−0.1181860.0000.000
Proposed scheme (3)10.150000.15000−0.3275−0.4275
2−0.182026−0.2408760.0740096−0.0599525
3−0.285858−0.1190560.0007704760.0000322954
4−0.286032−0.1181860.0000.000
Proposed scheme (4)10.150000.15000−0.3275−0.4275
2−3.351613.604747.4285716.0457
3−0.692548−0.0175140.2971370.39285
4−0.283154−0.117859−0.00196551−0.00295
5−0.286032−0.1181860.0000.000
Table 4. Solution for example 4 (B V P) using the proposed scheme.
Table 4. Solution for example 4 (B V P) using the proposed scheme.
SchemeNxy f ( x , y ) g ( x , y )
N S10.1000000.1000000.38000005.580000
20.5466777408637870.9985049833887041.614622410348670−0.3990420083663530
30.7431654008582190.9931608473034420.0000571195809957458−0.077214801060175
40.7321430794201050.9820935404168350.0002449705634444132−0.0002429831397698922
50.7321436794213340.9820632478816041.835275931227897  × 10 9 −7.20090653771876  × 10 13
60.7321436796857490.982063247916927−2.22044604925031  × 10 16 0.000000
P S (3)10.1000000.1000000.3800005.580000
20.7227652639465430.9787818530599810.04903313242663990.01602697777088990
30.7321436934712140.9820632699956406.47619073923522  × 10 8 −1.618596074948186  × 10 7
40.7321436796857490.982063247916927−2.22044604925031  × 10 16 0.000000
P S (4)10.1000000.1000000.3800005.580000
20.7321718116237600.9820897901242680.0000146207603477499−0.0001821979719456301
30.7321436796857490.982063247916927−2.22044604925031  × 10 16 0.000000
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Argyros, I.K.; Singh, M.K.; Regmi, S. Comparison between Two Competing Newton-Type High Convergence Order Schemes for Equations on Banach Spaces. Foundations 2023, 3, 643-659. https://doi.org/10.3390/foundations3040039

AMA Style

Argyros IK, Singh MK, Regmi S. Comparison between Two Competing Newton-Type High Convergence Order Schemes for Equations on Banach Spaces. Foundations. 2023; 3(4):643-659. https://doi.org/10.3390/foundations3040039

Chicago/Turabian Style

Argyros, Ioannis K., Manoj K. Singh, and Samundra Regmi. 2023. "Comparison between Two Competing Newton-Type High Convergence Order Schemes for Equations on Banach Spaces" Foundations 3, no. 4: 643-659. https://doi.org/10.3390/foundations3040039

APA Style

Argyros, I. K., Singh, M. K., & Regmi, S. (2023). Comparison between Two Competing Newton-Type High Convergence Order Schemes for Equations on Banach Spaces. Foundations, 3(4), 643-659. https://doi.org/10.3390/foundations3040039

Article Metrics

Back to TopTop