# On a Derivative-Free Variant of King’s Family with Memory

Authors: M. Sharifi; S. Karimi Vanani; F. Khaksar Haghani; M. Arab; S. Shateyi
Journal: The Scientific World Journal
 (2015)

---

## Abstract

The aim of this paper is to construct a method with memory according to King’s family of methods without memory for nonlinear equations. It is proved that the proposed method possesses higherR-order of convergence using the same number of functional evaluations as King’s family. Numerical experiments are given to illustrate the performance of the constructed scheme.


---

## Body

## 1. Introduction

Many problems arising in diverse disciplines of mathematical sciences can be described by a nonlinear equation of the following form (see, e.g., [1]): (1)
f
(
x
)
=
0
,where f
:
D
⊆
R
→
R is a sufficiently differentiable function in a neighborhood D of a simple zero α of (1). If we are interested in approximating the root α, we can do it by means of an iterative fixed-point method in the following form: (2)
x
k
+
1
=
ψ
(
x
k
)
,
k
≥
0
,provided that the starting point x
0 is given.In this work, we are concerned with the fixed-point methods that generate sequences presumably convergent to the true solution of a given single smooth equation. These schemes can be divided into one-point and multipoint schemes. We remark that the one-point methods can possess high order by using higher derivatives of the function, which is expensive from a computational point of view. On the other hand, the multipoint methods are allowing the user not to waste information that had already been used. This approach provides the construction of efficient iterative root-finding methods [2].In such circumstance, special attention is devoted to multipoint methods with memory that use already computed information to considerably increase convergence rate without additional computational costs. This would be the focus of this paper.Traub in [2] proposed the following method with memory (TM):(3)
w
k
=
x
k
+
β
k
f
x
k
,
β
k
=
-
1
f
x
k
,
x
k
-
1
,
k
=
0,1
,
2
,
…
,
x
k
+
1
=
x
k
-
f
x
k
f
x
k
,
w
k
,with the order of convergence 1
+
2.The iterative methods with memory can improve the order of convergence of the without memory method without any additional functional calculations, and this results in a higher computational efficiency index. We remark that it is assumed that an initial approximationx
0 close enough to the sought simple zero and β
0 are given for iterative methods of type (3).Recently, authors in [3] designed an approach to make derivative-free families with low complexity out of optimal methods. In fact, they conjectured that every time that one applies the approximation of the derivative f
′
(
x
n
)
≈
f
[
x
n
,
w
n
], with w
n
=
x
n
+
β
f
(
x
n
)
l, on an optimal order 2
q, we will need l
≥
q for preserving the order of convergence.For instance, choosing the well-known optimal two-step family of King (KM) [4],(4)
y
k
=
x
k
-
f
x
k
f
′
x
k
,
k
=
0,1
,
2
,
…
,
x
k
+
1
=
y
k
-
f
(
y
k
)
f
′
(
x
k
)
f
(
x
k
)
+
γ
f
(
y
k
)
f
(
x
k
)
+
(
γ
-
2
)
f
(
y
k
)
,
γ
∈
R
,and the conjecture of Cordero-Torregrosa, one may propose the following method (DKM): (5)
y
k
=
x
k
-
f
x
k
F
D
,
k
=
0,1
,
2
,
…
,
x
k
+
1
=
y
k
-
f
(
y
k
)
F
D
f
(
x
k
)
+
γ
f
(
y
k
)
f
(
x
k
)
+
(
γ
-
2
)
f
(
y
k
)
,
γ
∈
R
,wherein (6)
F
D
=
f
(
x
k
)
-
f
(
w
k
)
x
k
-
w
k
,
w
k
=
x
k
+
β
f
x
k
2
,
β
∈
R
∖
0
.In this work, we propose a two-step method with memory possessing a high efficiency index according to the well-known family of King’s methods (5).Our inspiration and motivation for constructing a higher-order method are linked in a direct manner with the fundamental concept of numerical analysis that any numerical method should give as accurate as possible output results with minimal computational cost. To state the matter differently, it is necessary to pursue methods of higher computational efficiency.For more background concerning this topic, one may refer to [5, 6].The paper is organized as follows. In Section2, the aim of this paper is presented by contributing an iterative method with memory based on (5) for solving nonlinear equations. The proposed scheme is an extension over (4) and has a simple structure with an increased computational efficiency. In Section 3, we compare the theoretical results by applying the definition of efficiency index and further supports are furnished whereas numerical reports are stated. Some concluding remarks will be drawn in Section 4 to end the paper.


## 2. A New Method with Memory

In this section, we propose the following iterative method with memory based on (5):(7)
y
k
=
x
k
-
f
x
k
F
D
,
w
k
=
x
k
+
β
k
f
x
k
2
,
k
=
0,1
,
2
,
…
,
x
k
+
1
=
y
k
-
f
(
y
k
)
F
D
f
(
x
k
)
+
γ
f
(
y
k
)
f
(
x
k
)
+
(
γ
-
2
)
f
(
y
k
)
,
γ
=
-
1
2
,wherein the self-accelerating parameter is β
k. The error equation of (5) is (γ
=
-
1
/
2) (8)
e
k
+
1
=
-
c
2
f
′
α
2
β
c
2
+
c
3
e
k
4
+
O
e
k
5
,where c
j
=
(
1
/
j
!
)
(
f
(
j
)
(
α
)
/
f
′
(
α
)
). We now must find a way so as to vanish the asymptotic error constant η
=
-
c
2
f
′
(
α
)
2
β
c
2
+
c
3.Toward this goal, one can increase theR-order by considering the following substitution: (9)
β
=
-
c
3
f
′
α
2
c
2
.Since the zero is not known, relation (9) cannot be used in its exact form and we must approximate it recursively. This builds a variant with memory for King’s family by using (10)
β
k
≈
-
c
-
3
f
-
′
α
2
c
-
2
,where c
-
j
≈
c
j. Now if we consider N
3
(
t
) to be Newton’s interpolation polynomial of third degree set through four available approximations x
k
,
x
k
-
1
,
y
k
-
1
,
w
k
-
1 at the end of each cycle, we can propose the following new method with memory:(11)
β
k
=
-
N
3
′′′
x
k
3
N
3
′
x
k
2
N
3
′′
(
x
k
)
,
y
k
=
x
k
-
f
x
k
F
D
,
w
k
=
x
k
+
β
k
f
x
k
2
,
k
=
0,1
,
2
,
…
,
x
k
+
1
=
y
k
-
f
(
y
k
)
F
D
f
(
x
k
)
-
1
/
2
f
(
y
k
)
f
(
x
k
)
-
5
/
2
f
(
y
k
)
.Note that, for example, we have the following formulation for the interpolating polynomial:(12)
N
3
′
x
k
=
d
d
t
N
3
(
t
)
t
=
x
k
=
f
[
x
k
,
x
k
-
1
]
+
f
[
x
k
,
x
k
-
1
,
y
k
-
1
]
(
x
k
-
x
k
-
1
)
+
f
[
x
k
,
x
k
-
1
,
y
k
-
1
,
w
k
-
1
]
(
x
k
-
x
k
-
1
)
(
x
k
-
y
k
-
1
)
.Acceleration in convergence for (11) is based on the use of a variation of one free nonzero parameter in each iterative step. This parameter is calculated using information from the current and previous iteration(s) so that the developed method may be regarded as method with memory according to Traub’s classification [2].We are at the time to write about the theoretical aspects of our proposed solver (11).Theorem 1.
Let the functionf
(
x
) be sufficiently differentiable in a neighborhood of its simple zero α. If an initial approximation x
0 is sufficiently close to α, then, the R-order of convergence of the two-step method (11) with memory is at least 4.23607.Proof.
Let{
x
k
} be a sequence of approximations generated by an iterative method. The error relations with the self-accelerating parameter β
=
β
k for (11) are in what follows: (13)
e
^
k
=
w
k
-
α
~
c
k
,
1
e
k
,
(14)
e
~
k
=
y
k
-
α
~
c
k
,
2
e
k
2
,
(15)
e
k
+
1
=
x
k
+
1
-
α
~
c
k
,
4
e
k
4
.Using a symbolic computations, we attain that (16)
-
c
2
f
′
α
2
β
c
2
+
c
3
~
e
k
-
1
.Substituting the value of -
c
2
f
′
(
α
)
2
β
c
2
+
c
3 from (16) in (15), one may obtain (17)
e
k
+
1
~
c
k
,
4
e
k
-
1
e
k
4
.Note that in general we know that the error equation should read e
k
+
1
~
A
e
k
p, where A and p are to be determined. Hence, one has e
k
~
A
e
k
-
1
p, and subsequently (18)
e
k
-
1
~
A
-
1
/
p
e
k
1
/
p
.Thus, it is easy to obtain (19)
e
k
p
~
A
-
1
/
p
C
e
k
4
+
1
/
p
,wherein C is a constant. This results in(20)
p
=
4
+
1
p
,with two solutions {
-
0.236068,4.23607
}. Clearly the value for p
=
4.23607 is acceptable and would be the convergence R-order of method (11) with memory. The proof is complete.The increase ofR-order is attained without any (new) additional function calculations so that the novel method with memory possesses a high computational efficiency index. This technique is an extension over scheme (5) to increase the R-order from 4 to 4.23607.The accelerating method (11) is new, simple, and useful, providing considerable improvement of convergence rate without any additional function evaluations in contrast to the optimal two-step methods without memory.We also remark that an alternative form of our proposed method with memory could be deduced using backward finite difference formula at the beginning of the first substep and a minor modification in the accelerators; that is to say, we have the following alternative method with memory possessing 4.23607 as itsR-order (APM) as well: (21)
β
k
=
N
3
′′′
(
x
k
)
3
N
3
′
x
k
2
N
3
′′
(
x
k
)
,
y
k
=
x
k
-
f
x
k
F
D
,
w
k
=
x
k
-
β
k
f
x
k
2
,
x
k
+
1
=
y
k
-
f
(
y
k
)
F
D
f
(
x
k
)
-
1
/
2
f
(
y
k
)
f
(
x
k
)
-
5
/
2
f
(
y
k
)
.Theorem 2.
Let the functionf
(
x
) be sufficiently differentiable in a neighborhood of its simple zero α. If an initial approximation x
0 is sufficiently close to α, then, the R-order of convergence of the two-step method (21) with memory is at least 4.23607.Proof.
The proof of this theorem is similar to Theorem1. It is hence omitted.


## 3. Numerical Computations

Computational efficiency of different iterative methods with and without memory can be measured in a prosperous manner by applying the definition of efficiency index. For an iterative method with convergence (R-)order r that requires θ functional evaluations, the efficiency index (also named computational efficiency) is calculated by Ostrowski-Traub’s formula [2]: (22)
E
=
r
1
/
θ
.According to this, we find (23)
E
SM
≈
1.4142
<
E
(
3
)
≈
1.5737
=
E
(
4
)
≈
1.5874
=
E
(
5
)
≈
1.5874
<
E
(
11
)
≈
1.6180
,where SM is the quadratically convergent method of Steffensen without memory [7].It should be remarked that Džunić in [8] designed an efficient one-step Steffensen-type method with memory possessing (
1
/
2
)
(
3
+
17
)  R-order of convergence as follows: (24)
w
k
=
x
k
+
β
k
f
(
x
k
)
,
β
k
=
-
1
N
2
′
(
x
k
)
,
p
k
=
-
N
3
′′
(
w
k
)
2
N
3
′
(
w
k
)
,
x
k
+
1
=
x
k
-
f
x
k
f
x
k
,
w
k
+
p
k
f
w
k
,and Cordero et al. in [9] presented a two-step biparametric Steffensen-type iterative method with memory possessing seventh R-order of convergence:(25)
w
k
=
x
k
+
β
k
f
(
x
k
)
,
β
k
=
-
1
N
3
′
(
x
k
)
,
p
k
=
-
N
4
′′
(
w
k
)
2
N
4
′
(
w
k
)
,
y
k
=
x
k
-
f
(
x
k
)
f
[
x
k
,
w
k
]
+
p
k
f
(
w
k
)
,
x
k
+
1
=
y
k
-
f
y
k
f
x
k
,
y
k
+
f
w
k
,
x
k
,
y
k
y
k
-
x
k
.Note that our main aim was to develop King’s family in terms of efficiencies index and was not to achieve the highest possible efficiency index.Although these methods possess higher computational efficiency indices than our proposed method (11), we exclude them from numerical comparisons since our method is not a Steffensen-type method and it is a Newton-type method with memory. For more refer to [10].Now, we apply and compare the behavior of different methods for finding the simple zeros of some different nonlinear test functions in the programming package Mathematica [11] using multiple precision arithmetic to clearly reveal the high R-order of PM and APM. We compare methods with the same number of functional evaluations per cycle.We notice that, by applying any root solver with local convergence, a special attention must be paid to the choice of initial approximations. If initial values are sufficiently close to the sought roots, then the expected (theoretical) convergence speed is obtainable in practice; otherwise, the iterative methods show slower convergence, especially at the beginning of the iterative process.In this section, the computational order of convergence (coc) has been computed by(26)
coc
=
ln
⁡
f
(
x
k
)
/
f
(
x
k
-
1
)
ln
⁡
f
(
x
k
-
1
)
/
f
(
x
k
-
2
)
.The calculated value coc estimates the theoretical order of convergence well when pathological behavior of the iterative method (i.e., slow convergence at the beginning of the implemented iterative method, oscillating behavior of approximations, etc.) does not exist.Here the results of comparisons for the test functions are given by applying 1000 fixed floating point arithmetic using the stop termination|
f
(
x
k
)
|
≤
1
0
-
100.Example 3.
We consider the following nonlinear test function in the intervalD
=
[
1.5,2.5
]: (27)
f
(
x
)
=
x
-
2
tan
x
x
3
-
8
,using the initial approximation x
0
=
1.7. The results are provided in Table 1.Table 1
Results of comparisons for Example3 and to find α
=
2.
Methods
|
f
(
x
1
)
|
|
f
(
x
2
)
|
|
f
(
x
3
)
|
|
f
(
x
4
)
|
coc
KM
18.577
64890.
7.2226
×
1
0
10
3.2493
×
1
0
9
—
OM
4.7484
0.0023129
1.3928
×
1
0
-
16
1.8313
×
1
0
-
69
4.00000
DKM
0.53362
5.3207
×
1
0
-
7
5.2711
×
1
0
-
31
5.0774
×
1
0
-
127
4.00000
PM
0.53362
1.9202
×
1
0
-
6
3.6106
×
1
0
-
30
1.6392
×
1
0
-
130
4.22928In this section, we have usedβ
0
=
0.0001 whenever required. Furthermore, for DKM we considered γ
=
-
1
/
2.Example 4.
We compare the behavior of different methods for finding the complex solution of the following nonlinear equation:(28)
g
(
x
)
=
(
-
1
+
2
I
)
+
1
x
+
x
+
sin
(
x
)
,using the initial approximation x
0
=
1
-
3
I where α
=
0.28860
⋯
-
1.24220
⋯
I. The results for this test are given in Table 2.Table 2
Results of comparisons for Example4.
Methods
|
f
(
x
1
)
|
|
f
(
x
2
)
|
|
f
(
x
3
)
|
|
f
(
x
4
)
|
coc
KM
2.0873
0.0095650
7.7971
×
1
0
-
12
3.4597
×
1
0
-
48
4.00000
OM
0.81344
0.0010884
1.5476
×
1
0
-
15
6.3280
×
1
0
-
63
4.00000
DKM
2.1909
0.013379
2.9909
×
1
0
-
11
7.5008
×
1
0
-
46
4.00000
PM
2.1909
0.0011772
7.0556
×
1
0
-
16
8.4197
×
1
0
-
68
4.23539
APM
1.9861
0.00089226
2.3251
×
1
0
-
16
7.5243
×
1
0
-
70
4.23526It is evident from Tables1 and 2 that approximations to the roots possess great accuracy when the proposed method with memory is applied. Results of the fourth iterate in Tables 1 and 2 are given only for demonstration of convergence speed of the tested methods and in most cases they are not required for practical problems at present.We also incorporated and applied the developed methods with memory (11) and (21) for different test examples and obtained results with the same behavior as above. Hence, we could mention that the theoretical results are upheld by numerical experiments and thus the new method is good with a high computational efficiency index.


## 4. Summary

In this paper, we have proposed a new two-step Steffensen-type iterative method with memory for solving nonlinear scalar equations. Using one self-correcting parameter calculated by Newton interpolatory polynomial, theR-order of convergence of the constructed method was increased from 4 to 4.23607 without any additional calculations.The new method was compared in performance and computational efficiency with some existing methods by numerical examples. We have observed that the computational efficiency index of the presented method with memory is better than those of other existing two-step King-type methods.



---

*Source: 514075-2015-03-25.xml*