Next Article in Journal
Improved Direct Linear Transformation for Parameter Decoupling in Camera Calibration
Previous Article in Journal
An Improved Dynamic Joint Resource Allocation Algorithm Based on SFR
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Comment

Comment on: On the Kung-Traub Conjecture for Iterative Methods for Solving Quadratic Equations. Algorithms 2016, 9, 1

1
Dipartimento di Scienza e Alta Tecnologia, Universita dell’Insubria, Via Valleggio 11, Como 22100, Italy
2
Departament de Física i Enginyeria Nuclear, Universitat Politècnica de Catalunya, Comte d’Urgell 187, 08036 Barcelona, Spain
Submission received: 21 January 2016 / Revised: 12 April 2016 / Accepted: 13 April 2016 / Published: 26 April 2016

Abstract

:
Kung-Traub conjecture states that an iterative method without memory for finding the simple zero of a scalar equation could achieve convergence order 2 d 1 , and d is the total number of function evaluations. In an article “Babajee, D.K.R. On the Kung-Traub Conjecture for Iterative Methods for Solving Quadratic Equations, Algorithms 2016, 9, 1, doi:10.3390/a9010001”, the author has shown that Kung-Traub conjecture is not valid for the quadratic equation and proposed an iterative method for the scalar and vector quadratic equations. In this comment, we have shown that we first reported the aforementioned iterative method.

1. Iterative Methods for Solving Quadratic Equations Presented in [1]

According to Kung-Traub conjecture (KTC) [2], an iterative method without memory for solving nonlinear equations in the case of simple zeros, could achieve a maximum convergence order of 2 d 1 , where d is the number of function evaluations. Recently, an article was published for solving quadratic equations [1] with arbitrary order of convergence by using one function and two derivatives evaluations. The details of the proposed formulation in [1] is given as follows. Let f ( x ) = κ 2 x 2 + κ 1 x + κ 0 be a quadratic function where κ 2 0 and κ 1 , κ 0 are constants. The proposed iteration function ψ ( r + 2 ) th BQIM ( x ) in [1] is
u ( x ) = f ( x ) f ( x ) τ = f x 2 / 3 u ( x ) f ( x ) H ( τ , r ) = 1 + i = 1 r a i ( τ 1 ) i ψ ( r + 2 ) th BQIM ( x ) = x u ( x ) H ( τ , r )
The error equation of Equation (1) is given as ψ ( r + 2 ) th BQIM ( x ) x * = C r c 2 r + 1 e ( k ) r + 2 + O e ( k ) r + 3 where C r is asymptotic error constant, x * is the simple root of quadratic equation and c 2 = 1 / 2 f ( x * ) 1 f ( x * ) . The error equation clearly shows that KTC is not valid in the case of quadratic equations. By using binomial expansion, the weight function H ( τ , r ) can be written as
H ( τ , r ) = 1 + i = 1 r j = 0 i a i i i j ( 1 ) ( i j ) τ j H ( τ , r ) = 1 + i = 0 r b i τ i
where b i is constant and can be computed by comparing two expressions of H ( τ , r ) in Equation (2). The powers of τ can be computed recursively and hence the iteration function Equation (1) is written as
u ( x ) = f ( x ) f ( x ) y = x 2 / 3 u ( x ) τ = f ( y ) f ( x ) ϕ 0 = u ( x ) for i = 1 , r ϕ i = τ ϕ i 1 end ψ ( r + 2 ) th BQIM ( x ) = x ϕ 0 i = 0 r b i ϕ i .
The computationally efficient vector version of iteration function Equation (3) is
F ( x ) ϕ 0 = F ( x ) y = x 2 / 3 ϕ 0 for i = 1 , r F ( x ) ϕ i = F ( y ) ϕ i 1 end ψ ( r + 2 ) th BQIM ( x ) = x ϕ 0 i = 0 r b i ϕ i .

2. Iterative Methods for Solving Matrix-Vector Quadratic Equations Presented in [3]

In this direction, a manuscript [3] was posted on 4 May 2015 on Researchgate in which the author provided models of three iterative methods, with their respective convergence orders, for computing the solution of matrix-vector quadratic equations. As the proposed iterative methods are valid for solving systems of nonlinear equations with quadratic nonlinearity, they are also valid for scalar quadratic equations. In the article [3] the model of iterative Method II can be written as:
x 0 = initial guess F ( x 0 ) ϕ 1 = F ( x 0 ) x 1 = x 0 + α 1 , 1 ϕ 1 for i = 1 , m F ( x 0 ) ϕ i + 1 = F ( x 1 ) ϕ i end x 2 = x 0 + j = 1 m α 2 , j ϕ j x 2 = x 0
where F ( x ) = 0 is a quadratic equation i.e., F ( x ) = B is a bilinear form and F ( s ) ( x ) = O (zero tensor) for s 3 and the convergence order is m + 2 . The scalar version of iterative method Equation (5) can be written as
x 0 = initial guess f ( x 0 ) ϕ 1 = f ( x 0 ) x 1 = x 0 + α 1 , 1 ϕ 1 for i = 1 , m f ( x 0 ) ϕ i + 1 = f ( x 1 ) ϕ i end x 2 = x 0 + j = 1 m α 2 , j ϕ j x 0 = x 2
where f ( x ) = a x 2 + b x + c with a 0 and convergence order is m + 2 . The convergence proofs of different iterative methods are established in Figure 1, Figure 2, Figure 3 and Figure 4. We can see that the error equations in all cases are the same. The Figure 4 shows that the iterative method (1) is a particular case of iterative method Equation (6) for α 1 , 1 = 2 / 3 . Finally we provide the convergence proof of the iterative method Equation (5).
Theorem 2.1. 
Let F : D R n R n be a function with all continuous Fréchet derivatives, F ( x ) = B and F ( j ) ( x ) = O for j 3 , where B is a bilinear form and D is convex open subset of R n . If we take x 0 in the vicinity of a simple root x * of F ( x ) = 0 then the sequence of successive approximation generated by iterative method Equation (5) for m = 6 and α 1 , 1 = 2 / 3 converges to x * with convergence order eight.
Proof. 
We denote e k = x k x * , C 1 = F ( x * ) and C 2 = F ( x * ) 1 F ( x * ) / 2 . By expanding F ( x 0 ) around x * we get
F ( x 0 ) = C 1 e 0 + C 2 e 0 2
The Fréchet derivative of Equation (7) with respective to e 0 is
F ( x 0 ) = C 1 I + 2 C 2 e 0
The inverse of Equation (8) is
F ( x 0 ) 1 = ( I 2 C 2 e 0 + 4 C 2 2 e 0 2 8 C 2 3 e 0 3 + 16 C 2 4 e 0 4 32 C 2 5 e 0 5 + 64 C 2 6 e 0 6 128 C 2 7 e 0 7 + 256 C 2 8 e 0 8 ) C 1 1 + O e 0 9
We compute ϕ 1 using Equations (9) and (7)
ϕ 1 = e 0 C 2 e 0 2 + 2 C 2 2 e 0 3 4 C 2 3 e 0 4 + 8 C 2 4 e 0 5 16 C 2 5 e 0 6 + 32 C 2 6 e 0 7 64 C 2 7 e 0 8 + O e 0 9 .
The expression for e 1 is
e 1 = 1 / 3 e 0 + 2 C 2 e 0 2 4 C 2 2 e 0 3 + 8 C 2 3 e 0 4 16 C 2 4 e 0 5 + 32 C 2 5 e 0 6 + 64 C 2 6 e 0 7 + 128 C 2 7 e 0 8 + O e 0 9
First order Fréchet derivative of F ( x ) at x 1 is
F ( x 1 ) = C 1 ( I + 1 / 3 ( 2 C 2 e 0 + 4 C 2 2 e 0 2 8 C 2 3 e 0 3 + 16 C 2 4 e 0 4 32 C 2 5 e 0 5 + 64 C 2 6 e 0 6 128 C 2 7 e 0 7 + 256 C 2 8 e 0 8 ) ) + O e 0 9
Next we compute M = F ( x 0 ) 1 F ( x 1 )
M = I 4 / 3 C 2 e 0 + 4 C 2 2 e 0 2 32 / 3 C 2 3 e 0 3 + 80 / 3 C 2 4 e 0 4 64 C 2 5 e 0 5 + 448 / 3 C 2 6 e 0 6 1024 / 3 C 2 7 e 0 7 + 768 C 2 8 e 0 8 + O e 0 9
with help of Equation (13), we obtain the expressions for ϕ i for i { 2 , 3 , 4 , 5 , 6 , 7 }
ϕ 2 = e 0 7 / 3 C 2 e 0 2 + 22 / 3 C 2 2 e 0 3 64 / 3 C 2 3 e 0 4 + 176 / 3 C 2 4 e 0 5 464 / 3 C 2 5 e 0 6 + 1184 / 3 C 2 6 e 0 7 2944 / 3 C 2 7 e 0 8 + O e 0 9 ϕ 3 = e 0 11 / 3 C 2 e 0 2 + 130 / 9 C 2 2 e 0 3 460 / 9 C 2 3 e 0 4 + 168 C 2 4 e 0 5 1568 / 3 C 2 5 e 0 6 + 4672 / 3 C 2 6 e 0 7 4480 C 2 7 e 0 8 + O e 0 9 ϕ 4 = e 0 5 C 2 e 0 2 + 70 / 3 C 2 2 e 0 3 2584 / 27 C 2 3 e 0 4 + 9712 / 27 C 2 4 e 0 5 34208 / 27 C 2 5 e 0 6 + 114496 / 27 C 2 6 e 0 7 367616 / 27 C 2 7 e 0 8 + O e 0 9 ϕ 5 = e 0 19 / 3 C 2 e 0 2 + 34 C 2 2 e 0 3 4252 / 27 C 2 3 e 0 4 + 53512 / 81 C 2 4 e 0 5 208624 / 81 C 2 5 e 0 6 + 767968 / 81 C 2 6 e 0 7 2697920 / 81 C 2 7 e 0 8 + O e 0 9 ϕ 6 = e 0 23 / 3 C 2 e 0 2 + 418 / 9 C 2 2 e 0 3 6448 / 27 C 2 3 e 0 4 + 89168 / 81 C 2 4 e 0 5 1137712 / 243 C 2 5 e 0 6 + 4543840 / 243 C 2 6 e 0 7 5747840 / 81 C 2 7 e 0 8 + O e 0 9 ϕ 7 = e 0 9 C 2 e 0 2 + 182 / 3 C 2 2 e 0 3 9236 / 27 C 2 3 e 0 4 + 46264 / 27 C 2 4 e 0 5 637376 / 81 C 2 5 e 0 6 + 24618880 / 729 C 2 6 e 0 7 100011520 / 729 C 2 7 e 0 8 + O e 0 9 .
We obtain expression for e 2 by using Equation (14)
e 2 = 1 + α 2 , 1 + α 2 , 2 + α 2 , 3 + α 2 , 4 + α 2 , 5 + α 2 , 6 + α 2 , 7 e 0 + α 2 , 1 7 / 3 α 2 , 2 11 / 3 α 2 , 3 5 α 2 , 4 19 / 3 α 2 , 5 23 / 3 α 2 , 6 9 α 2 , 7 C 2 e 0 2 + ( 2 α 2 , 1 + 22 / 3 α 2 , 2 + 130 / 9 α 2 , 3 + 70 / 3 α 2 , 4 + 34 α 2 , 5 C 2 2 + 418 / 9 α 2 , 6 + 182 / 3 α 2 , 7 ) C 2 2 e 0 3 + ( 4 α 2 , 1 64 / 3 α 2 , 2 460 / 9 α 2 , 3 2584 / 27 α 2 , 4 4252 / 27 α 2 , 5 6448 / 27 α 2 , 6 9236 / 27 α 2 , 7 ) C 2 3 e 0 4 + ( 8 α 2 , 1 + 176 / 3 α 2 , 2 + 168 α 2 , 3 + 9712 / 27 α 2 , 4 + 53512 / 81 α 2 , 5 + 89168 / 81 α 2 , 6 + 46264 / 27 α 2 , 7 ) C 2 4 e 0 5 + ( 16 α 2 , 1 464 / 3 α 2 , 2 1568 / 3 α 2 , 3 34208 / 27 α 2 , 4 208624 / 81 α 2 , 5 1137712 / 243 α 2 , 6 637376 / 81 α 2 , 7 ) C 2 5 e 0 6 + ( 32 α 2 , 1 + 1184 / 3 α 2 , 2 + 4672 / 3 α 2 , 3 + 114496 / 27 α 2 , 4 + 767968 / 81 α 2 , 5 + 4543840 / 243 α 2 , 6 + 24618880 / 729 α 2 , 7 ) C 2 6 e 0 7 + ( 64 α 2 , 1 2944 / 3 α 2 , 2 4480 α 2 , 3 367616 / 27 α 2 , 4 2697920 / 81 α 2 , 5 5747840 / 81 α 2 , 6 100011520 / 729 α 2 , 7 ) C 2 7 e 0 8 + O e 0 9 .
By equating the coefficients of powers of e 0 in Equation (15), we get system of seven equations
e q 1 : = 1 + α 2 , 1 + α 2 , 2 + α 2 , 3 + α 2 , 4 + α 2 , 5 + α 2 , 6 + α 2 , 7 = 0 e q 2 : = α 2 , 1 7 / 3 α 2 , 2 11 / 3 α 2 , 3 5 α 2 , 4 19 / 3 α 2 , 5 23 / 3 α 2 , 6 9 α 2 , 7 = 0 e q 3 : = 2 α 2 , 1 + 22 / 3 α 2 , 2 + 130 / 9 α 2 , 3 + 70 / 3 α 2 , 4 + 34 α 2 , 5 + 418 / 9 α 2 , 6 + 182 / 3 α 2 , 7 = 0 e q 4 : = 4 α 2 , 1 64 / 3 α 2 , 2 460 / 9 α 2 , 3 2584 / 27 α 2 , 4 4252 / 27 α 2 , 5 6448 / 27 α 2 , 6 9236 / 27 α 2 , 7 = 0 e q 5 : = 8 α 2 , 1 + 176 / 3 α 2 , 2 + 168 α 2 , 3 + 9712 / 27 α 2 , 4 + 53512 / 81 α 2 , 5 + 89168 / 81 α 2 , 6 + 46264 / 27 α 2 , 7 = 0 e q 6 : = 16 α 2 , 1 464 / 3 α 2 , 2 1568 / 3 α 2 , 3 34208 / 27 α 2 , 4 208624 / 81 α 2 , 5 1137712 / 243 α 2 , 6 637376 / 81 α 2 , 7 = 0 e q 7 : = 32 α 2 , 1 + 1184 / 3 α 2 , 2 + 4672 / 3 α 2 , 3 + 114496 / 27 α 2 , 4 + 767968 / 81 α 2 , 5 + 4543840 / 243 α 2 , 6 + 24618880 / 729 α 2 , 7 = 0 .
The solution set of seven equations is
s o l = { α 2 , 1 = 43903 1024 , α 2 , 2 = 55767 256 , α 2 , 3 = 497763 1024 , α 2 , 4 = 37719 64 , α 2 , 5 = 416421 1024 , α 2 , 6 = 38637 256 , α 2 , 7 = 24057 1024 } .
After simplification the error equation we get
e 2 = 429 C 2 7 e 0 8 + O e 0 9 .
As we have seen, the convergence order of the iterative method Equation (5) is eight for m = 6 and this confirms our claimed order of convergence which is s + 2 for m = s . Now we provide the proof of convergence order via mathematical induction.
Theorem 2.2. 
Let F : D R n R n be a function with all continuous Fréchet derivatives, F ( x ) = B and F ( j ) ( x ) = O for j 3 , where B is a bilinear form and D is a convex open subset of R n . If we take x 0 in the vicinity of a simple root x * of F ( x ) = 0 then the sequence of successive approximation generated by the iterative method Equation (5) for m = s and α 1 , 1 = 2 / 3 converges to x * with convergence order s + 2 .
Proof. 
We suppose that our claim about the convergence order of the iterative method Equation (5) is true for m = s , which means we have
e 2 ( s ) = H s C 2 s 1 e 0 s + 2 + O e 0 s + 3
where H s is asymptotic error constant and superscript “ ( s ) ” means the value of e 2 when m = s . We can write e 2 ( s )
e 2 ( s ) = e 0 2 / 3 ϕ 1 + α 2 , 1 + α 2 , 2 M + + α 2 , s 1 M s ϕ 1 + H s C 2 s 1 e 0 s + 2 + O e 0 s + 3
It is convenient to express the combination
α 2 , 1 + α 2 , 2 M + α 2 , 3 M 2 + + α 2 , s + 1 M s
in the powers of I M [1], we establish the following identity
α 2 , 1 + α 2 , 2 M + + α 2 , s 1 M s β 2 , 1 + β 2 , 2 I M + + β 2 , s 1 I M s
By comparing the same powers of M on both sides we can easily compute the value of β i ’s. By using Equation (21), error Equation (20) can be written as
e 2 ( s ) = e 0 2 / 3 ϕ 1 + β 2 , 1 + β 2 , 2 I M + + β 2 , s 1 I M s ϕ 1 + H s C 2 ( s + 1 ) e 0 ( s + 2 ) + O e 0 ( s + 3 )
However, according to our assumption we can find the value of unknowns to make the following expression equal to zero
e 0 2 / 3 ϕ 1 + β 2 , 1 + β 2 , 2 I M + + β 2 , s 1 I M s ϕ 1 = 0
We can notice from Equation (13) that I M s = 3 4 s C 2 s e 0 s + O e 0 ( s + 1 ) and ϕ 1 = O ( e 0 ) .
Now we consider
e 2 ( s + 1 ) = e 2 ( s ) + β 2 , s I M s + 1 ϕ 1 = H s C 2 s + 1 e 0 s + 2 + β 2 , s 3 4 s + 1 C 2 s + 1 e 0 ( s + 2 ) + O e 0 ( s + 3 ) = H s + 3 4 s + 1 β 2 , s C 2 s + 1 e 0 ( s + 2 ) + O e 0 ( s + 3 )
As we know the value of H s we can find β 2 , s to make the coefficient of e 0 ( s + 2 ) equals to zero. Hence we get
e 2 ( s + 1 ) = O e 0 ( s + 3 )
which completes the proof. ☐

3. Numerical Testing

We adopt the following definition of computational convergence order (COC)
COC = l o g | | F ( x k + 1 ) | | / | | F ( x k ) | | l o g | | F ( x k ) | | / | | F ( x k 1 ) | |
To verify the claimed convergence order of our proposed iterative method Equation (5), we study the following system of quadratic equations
F ( x ) = x 2 x 3 + x 4 ( x 2 + x 3 ) = 0 x 1 x 3 + x 4 ( x 4 + x 3 ) = 0 x 1 x 2 + x 4 ( x 1 + x 4 ) = 0 x 1 x 2 + x 3 ( x 1 + x 2 ) = 1
In Table 1, we listed the norm of the residue of F ( x ) and COC against the sequence of iterations for different values of parameters α 11 . The Table 1 confirms the claimed convergence order. For two different values of the parameter α 11 , we obtained the record of the norm of the residue of F ( x ) equal to the system of quadratic Equations (26). The possible reason for his could be the same error equation of iterative method for different values of parameter α 11 in the iterative method Equation (5).

4. Conclusions

We conclude that the iterative structure of iterative method Equation (1) was first reported in article [3] as a particular case and our proposed iterative method [3] is general because α 1 , 1 is a free parameter.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Babajee, D.K.R. On the Kung-Traub Conjecture for Iterative Methods for Solving Quadratic Equations. Algorithms 2016, 9, 1. [Google Scholar] [CrossRef]
  2. Traub, J.F. Iterative Methods for the Solution of Equations; Prentice-Hall: Englewood Cliffs, NJ, USA, 1964. [Google Scholar]
  3. Ahmad, F. Higher Order Iterative Methods for Solving Matrix Vector Equations. Researchgate 2015. [Google Scholar] [CrossRef]
Figure 1. Iterative method Equation (1).
Figure 1. Iterative method Equation (1).
Algorithms 09 00030 g001
Figure 2. Iterative method Equation (6) for arbitrary finite α 1 , 1 .
Figure 2. Iterative method Equation (6) for arbitrary finite α 1 , 1 .
Algorithms 09 00030 g002
Figure 3. Iterative method (6) for α 1 , 1 = 1 .
Figure 3. Iterative method (6) for α 1 , 1 = 1 .
Algorithms 09 00030 g003
Figure 4. Iterative method (6) for α 1 , 1 = 2 / 3 .
Figure 4. Iterative method (6) for α 1 , 1 = 2 / 3 .
Algorithms 09 00030 g004
Table 1. Computational convergence order of iterative method Equation (5), initial guess = [ 0 . 5 , 0 . 5 , 0 . 5 , 0 . 25 ] .
Table 1. Computational convergence order of iterative method Equation (5), initial guess = [ 0 . 5 , 0 . 5 , 0 . 5 , 0 . 25 ] .
iter | | F ( x ) | | COC | | F ( x ) | | COC
α 11 = 1 α 11 = 2 / 3
02.50e1 2.50e1
12.71e–6 2.71e–6
27.56e–478.177.56e–478.17
32.80e–3718.002.80e–3718.00
49.79e–29678.009.79e–29678.00

Share and Cite

MDPI and ACS Style

Ahmad, F. Comment on: On the Kung-Traub Conjecture for Iterative Methods for Solving Quadratic Equations. Algorithms 2016, 9, 1. Algorithms 2016, 9, 30. https://0-doi-org.brum.beds.ac.uk/10.3390/a9020030

AMA Style

Ahmad F. Comment on: On the Kung-Traub Conjecture for Iterative Methods for Solving Quadratic Equations. Algorithms 2016, 9, 1. Algorithms. 2016; 9(2):30. https://0-doi-org.brum.beds.ac.uk/10.3390/a9020030

Chicago/Turabian Style

Ahmad, Fayyaz. 2016. "Comment on: On the Kung-Traub Conjecture for Iterative Methods for Solving Quadratic Equations. Algorithms 2016, 9, 1" Algorithms 9, no. 2: 30. https://0-doi-org.brum.beds.ac.uk/10.3390/a9020030

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