We derive a family of eighth-order multipoint methods for the solution of nonlinear equations. In terms of computational cost, the family requires evaluations of only three functions and one first derivative per iteration. This implies that the efficiency index of the present methods is 1.682. Kung and Traub (1974) conjectured that multipoint iteration methods without memory based on n evaluations have optimal order 2n-1. Thus, the family agrees with Kung-Traub conjecture for the case n=4. Computational results demonstrate that the developed methods are efficient and robust as compared with many well-known methods.
1. Introduction
Solving nonlinear equations is one of the most important problems in science and engineering [1, 2]. The boundary value problems arising in kinetic theory of gases, vibration analysis, design of electric circuits, and many applied fields are reduced to solving such equations. In the present era of advance computers, this problem has gained much importance than ever before.
In this paper, we consider iterative methods to find a simple root r of the nonlinear equation f(x)=0, where f:R→R be the continuously differentiable real function. Newton’s method [1] is probably the most widely used algorithm for solving such equations, which starts with an initial approximation x0 closer to the root r and generates a sequence of successive iterates {xi}0∞ converging quadratically to the root. It is given by the following:
(1.1)xi+1=xi-f(xi)f′(xi),i=0,1,2,3,….
In order to improve the local order of convergence, a number of ways are considered by many researchers, see [3–26] and references therein. In particular, King [3] developed a one-parameter family of fourth-order methods defined by
(1.2)wi=xi-f(xi)f′(xi),xi+1=wi-f(xi)+βf(wi)f(xi)+(β-2)f(wi)f(wi)f′(xi),
where wi is the Newton point and β is a constant.
This family requires two evaluations of the function f and one evaluation of first derivative f′ per iteration. The famous Ostrowski’s method [4, 5] is a member of this family for the case β=0. From practical point of view, the methods (1.2) are important because of higher efficiency than Newton’s method (1.1).
Traub [5] has divided iterative methods into two classes, namely, one-point methods and multipoint methods. Each class is further divided into two subclasses, namely, one-point methods with and without memory, and multipoint methods with and without memory. The important aspects related to these classes of methods are order of convergence and computational efficiency. Order of convergence shows the speed with which a given sequence of iterates converges to the root while the computational efficiency concerns with the economy of the entire process. Investigation of one-point methods with and without memory, has demonstrated theoretical restrictions on the order and efficiency of these two categories (see [5]). However, Kung and Traub [6] have conjectured that multipoint iteration methods without memory based on n evaluations have optimal order 2n-1. In particular, with three evaluations a method of fourth-order can be constructed. The King’s method (1.2) is a well-known example of fourth-order multipoint methods without memory.
Recently, based on Ostrowski’s or King’s methods some higher-order multipoint methods have been proposed and analyzed for solving nonlinear equations. For example, Grau and Díaz-Barrero [10], Sharma and Guha [11], and Chun and Ham [12] have developed sixth-order modified Ostrowski’s methods each requires three f and one f′ evaluations per iteration. Kou et al. [15] presented a family of variants of Ostrowski’s method with seventh-order convergence requiring three fand one f′ evaluations. With same number of evaluations, Bi et al. [18] developed a seventh-order family of modified King’s methods. Bi et al. [19] also presented an eighth-order family of modified King’s methods requiring four evaluations which agrees with the Kung-Traub conjecture.
In this paper, we present a new family of eighth-order methods without using second and higher derivatives. In terms of computational cost, it requires the evaluations of three functions and one first derivative per iteration. Thus the present methods provide a new example of multipoint methods without memory that with four evaluations a method of optimum order eight can be achieved as conjectured by Kung and Traub. The performance and effectiveness of the developed family of methods is tested and compared through some test functions.
Contents of the paper are summarized as follows. Some basic definitions relevant to the present work are presented in Section 2. In Section 3, we obtain new methods. Convergence analysis, for establishing eighth-order convergence, is carried out in Section 4. In Section 5, we provide some particular cases of the family. In Section 6, the method is tested and compared with other well-known methods on a number of problems. Concluding remarks are given in Section 7.
2. Basic DefinitionsDefinition 2.1.
Let f(x) be a real function with a simple root r and let {xi}i∈N be a sequence of real numbers that converges towards r. Then, we say that the order of convergence of the sequence is p, if there exits a number p∈ℝ+ such that
(2.1)limi→∞xi+1-r(xi-r)p=C,
for some C≠0, C is known as the asymptotic error constant.
If p=1, 2 or 3, the sequence is said to have linear convergence, quadratic convergence or cubic convergence, respectively.
Definition 2.2.
Let ei=xi-r be the error in the ith iteration, we call the relation
(2.2)ei+1=Ceip+O(eip+1),
the error equation.
Definition 2.3.
Let n be the number of new pieces of information required by a method. A “piece of information” typically is any evaluation of a function or one of its derivatives. The efficiency of the method is measured by the concept of efficiency index [27] and is defined by the following:
(2.3)E=p1/n,
where p is the order of the method.
Definition 2.4.
Suppose that xi+1,xi and xi-1 are three successive iterations closer to the root r. Then, the computational order of convergence ρ (see [24, 25, 28]) is approximated by using (2.2) as follows:
(2.4)ρ≅ln|(xi+1-r)/xi-r|ln|xi-r/xi-1-r|.
3. The Method
We consider the iteration scheme of the form
(3.1)wi=xi-f(xi)f′(xi),zi=wi-ω(λi)f(wi)f′(xi),xi+1=zi-W(μi)f(zi)f′(zi),
where λi=f(wi)/f(xi), μi=f(zi)/f(xi), and ω(t) and W(t) represent the real-valued functions (here onwards called weight functions). This scheme consists of three steps in which the first step represents Newton’s method and last two are weighted-Newton steps. It is quite obvious that formula (3.1) requires five evaluations per iteration. However, we can reduce the number of evaluations to four by using some suitable approximation of the derivative f′(zi). We obtain this approximation by considering the approximation of f(x) by a rational linear function of the form
(3.2)y(x)-y(xi)=(x-xi)+ab(x-xi)+c,
where the parameters a,b, and c are determined by the condition that f and y coincide at xi, wi and zi. That means y(x) satisfies the conditions
(3.3)y(xi)=f(xi),y(wi)=f(wi),y(zi)=f(zi).
From (3.2) and first condition of (3.3), it is easy to show that
(3.4)a=0.
Substituting the value of a into (3.2) then using the last two conditions of (3.3), after some simple calculations we obtain
(3.5)b(wi-xi)+c=1f[xi,wi],b(zi-xi)+c=1f[xi,zi],
where f[xi,wi]=(f(wi)-f(xi))/(wi-xi) and f[xi,zi]=(f(zi)-f(xi))/(zi-xi) are first-divided differences.
Solving these equations, we can obtain b and c as follows:
(3.6)b=1wi-zi(1f[xi,wi]-1f[xi,zi]),c=1wi-zi(xi-zif[xi,wi]-xi-wif[xi,zi]).
Differentiation of (3.2) gives
(3.7)y′(x)=c[b(x-xi)+c]2.
We can now approximate the derivative f'(x) with the derivative y'(x) of rational function (3.2) and obtain
(3.8)f′(zi)≈y′(zi).
Substituting the values of b and c obtained in (3.6) into (3.7) then using (3.8), we get after simplifications
(3.9)f′(zi)=f[xi,zi]f[wi,zi]f[xi,wi].
Then the iteration scheme (3.1) in its final form is given by the following:
(3.10)wi=xi-f(xi)f′(xi),zi=wi-ω(λi)f(wi)f′(xi),xi+1=zi-W(μi)f[xi,wi]f(zi)f[xi,zi]f[wi,zi],
where λi=f(wi)/f(xi), μi=f(zi)/f(xi), and ω(t) and W(t) are the weight functions.
Thus the scheme (3.10) defines a new family of multipoint methods with two weight functions ω(t) and W(t). In the next section, we will see that both of these functions play an important role in establishing eighth-order convergence of the methods.
4. Convergence of the Method
In order to examine the convergence property of the family (3.10), we prove the following theorem.
Theorem 4.1.
Let the function f:R→R be sufficiently smooth in R. If f(x) has a simple root r in R and x0 is sufficiently close to r, then the sequence {xi} generated by any method of the family (3.10) converges to r with convergence order eight, provided the weight functions ω(t) and W(t) satisfy the conditions ω(0)=1,ω′(0)=2,ω′′(0)=8,W(0)=1, W'(0)=1, and |ω′′′(0)|<∞.
Proof.
Let ei=xi-r be the error in the iterate xi. Using Taylor’s series expansion, we get
(4.1)f(xi)=f′(r)[ei+A2ei2+A3ei3+A4ei4+A5ei5+A6ei6+A7ei7+A8ei8+O(ei9)],f′(xi)=f′(r)[1+2A2ei+3A3ei2+4A4ei3+5A5ei4+6A6ei5+7A7ei6+8A8ei7+O(ei8)],
where Ak=f(k)(r)/k!f′(r) for k∈N, N is the set of natural numbers.
Now,
(4.2)f(xi)f′(xi)=ei-A2ei2-2(-A22+A3)ei3-(4A23-7A2A3+3A4)ei4-K1ei5-K2ei6-K3ei7-K4ei8+O(ei9),
Following are the expressions of Kn(n=1,2,3,4)(4.3)K1=-8A24+20A22A3-6A32-10A2A4+4A5,K2=16A25-52A23A3+33A2A32+28A22A4-17A3A4-13A2A5+5A6,K3=-32A26+128A24A3-126A22A32+18A33-72A23A4+92A2A3A4-12A42+36A22A5-22A3A5-16A2A6+6A7,K4=64A27-304A25A3+408A23A32-135A2A33+176A24A4-348A22A3A4+75A32A4+64A2A42-92A23A5+118A2A3A5-31A4A5+44A22A6-27A3A6-19A2A7+7A8.
For the sake of brevity, we omit their specific forms. We will use the same means in the following.
For
(4.4)wi=xi-f(xi)f′(xi)=r+A2ei2+2(-A22+A3)ei3+(4A23-7A2A3+3A4)ei4+K1ei5+K2ei6+K3ei7+K4ei8+O(ei9).e~i=A2ei2+2(-A22+A3)ei3+(4A23-7A2A3+3A4)ei4+K1ei5+K2ei6+K3ei7+K4ei8+O(ei9).
Using Taylor’s series expansion, we get
(4.5)f(wi)=f′(r)[e~i+A2e~i2+A3e~i3+A4e~i4+O(ei9)],
therefore,
(4.6)f(wi)f′(xi)=e~i-2A2e~iei+[(4A22-3A3)e~iei2+A2e~i2]+[-2A22e~i2ei+(-8A23+12A2A3-4A4)e~iei3]+M1ei6+M2ei7+M3ei8+O(ei9),M1=(4A23-3A2A3)e~i2ei4+(16A24-36A22A3+9A32+16A2A4-5A5)e~iei2+A3e~i3ei6,M2=-2A2A3e~i3ei6+(-8A24+12A22A3-4A2A4)e~i2ei4+(-32A25+96A23A3-54A2A32-48A22A4+24A3A4+20A2A5-6A6)e~iei2,M3=(4A22A3-3A32)e~i3ei6+(16A25-36A23A3+9A2A32+16A22A4-5A2A5)e~i2ei4+(64A26-240A24A3+216A22A32-27A33+128A23A4-144A2A3A4+16A42-60A22A5+30A3A5+24A2A6-7A7)e~iei2+A3e~i4ei8.
Also
(4.7)λi=f(wi)f(xi)=e~iei-A2e~i+[A2e~i2ei+(A22-A3)e~iei]+[(-A23+2A2A3-A4)e~iei2-A22e~i2]+[(A24-3A22A3+A32+2A2A4-A5)e~iei3+(A23-A2A3)e~i2ei+A3e~i3ei]+[(-A25+4A23A3-3A2A32-3A22A4+2A3A4+2A2A5-A6)e~iei4+(-A24+2A22A3-A2A4)e~i2ei2-A2A3e~i3]+O(ei7).
Thus, using the Taylor expansion, we get
(4.8)ω(λi)=ω(0)+ω′(0)λi+12!ω′′(0)λi2+13!ω′′′(0)λi3+O(λi4)=ω(0)+ω′(0)e~iei+[12ω′′(0)e~i2ei2-A2ω′(0)e~i]+[16ω′′′(0)e~i3ei3+A2(ω′(0)-ω′′(0))e~i2ei+(A22-A3)ω′(0)e~iei]+L1ei4+L2ei5+L3ei6+O(ei7),
where
(4.9)L1=124ωiν(0)e~i4ei8+12A2(2ω′′(0)-ω′′′(0))e~i3ei6+(-A3ω′′(0)+A22(-ω′(0)+32ω′′(0))e~i2ei4-(A23+2A2A3-A4)ω′(0)e~iei2,L2=1120ων(0)e~i5ei10+16A2(3ω′′′(0)-ω(iν)(0))×(A3(ω′(0)-12ω′′′(0))+A22(-2ω′′(0)+ω′′′(0)))e~i3ei6+(-A2A3(ω′(0)-3ω′′(0))+A23(ω′(0)-2ω′′(0))-A4ω′′(0))e~i2ei4+(A24-3A22A3+A32+2A2A4-A5)ω′(0)e~iei2,L3=1720ωνi(0)e~i6ei12+124A2(4ωiν(0)-ων(0))e~i5ei10+112(2A3(6ω′′(0)-ωiν(0))+A22(6ω′′(0)-18ω′′′(0)+5ωiν(0)))e~i4ei8+(-A2A3(ω′(0)+2ω′′(0)-2ω′′′(0))+A23(3ω′′(0)-53ω′′′(0))-12A4ω′′′(0))e~i3+(2A22A3(ω′(0)-3ω′′(0))-A2A4(ω′(0)-3ω′′(0))+12(3A32-2A5)ω′′(0)+A24(-ω′(0)+52ω′′(0)))e~i2ei4+(-A25+4A23A3-3A2A32-3A22A4+2A3A4+2A2A5-A6)ω′(0)e~iei2.
Using (4.6) and (4.8), we have
(4.10)zi=wi-ω(λi)f(wi)f'(xi)=r+[1-ω(0)]e~i-[ω′(0)e~i2ei-2ω(0)A2e~iei]-[12ω′′(0)e~i3ei2+A2(ω(0)-3ω′(0))e~i2+ω(0)(4A22-3A3)e~iei2]-[12ω′′(0)e~i3ei2+A2(ω(0)-3ω′(0))e~i2+ω(0)(4A22-3A3)e~iei2]ei5-M6ei6-M7ei7-M8ei8+O(ei9),
where Mn(n=6,7,8) are the expression about An(n=2,3,…,8).
If (0)=1, ω′(0)=2 and substituting the value of e~i from (4.4), we get
(4.11)e^i=zi-r=[(5-12ω′′(0))A23-A2A3]ei4+(-36+5ω′′(0)-16ω′′′(0))A24+(32-3ω′′(0))A22A3-2A32-2A2A4]ei5+M6ei6+M7ei7+M8ei8+O(ei9).
Using Taylor’s series expansion, we get
(4.12)f(zi)=f′(r)[e^i+A2e^i2+O(ei9)],
furthermore,
(4.13)f[xi,wi]=f(wi)-f(xi)wi-xi=f′(r)[1+A2ei+(A22+A3)ei2+(-2A23+3A2A3+A4)ei3+(4A24-8A22A3+2A32+4A2A4+A5)ei4+O(ei5)].f[xi,zi]=f(zi)-f(xi)zi-xi=f′(r)[((5-12ω′′(0))A24-A22A3+A5)1+A2ei+A3ei2+A4ei3+((5-12ω′′(0))A24-A22A3+A5)ei4+O(ei5)].f[wi,zi]=f(zi)-f(wi)zi-wi=f′(r)[((9-12ω′′(0))A24-7A22A3+3A2A4)1+A22ei2-2(A23-A2A3)ei3+((9-12ω′′(0))A24-7A22A3+3A2A4)ei4+O(ei5)].
Using the above results, we obtain
(4.14)f[xi,wi]f[wi,zi]f[xi,zi]=1f′(r)[((-7+ω′′(0))A24-4A22A3+2A32+A2A4)1+(-A23+A2A3)ei3+((-7+ω′′(0))A24-4A22A3+2A32+A2A4)ei4+O(ei5)].
Also
(4.15)μi=f(zi)f(xi)=e^iei-A2e^i+O(ei5).
Thus, using the Taylor expansion and |W′′(0)|<∞,we get
(4.16)W(μi)=W(0)+W′(0)μi+O(μi2)=W(0)+W′(0)(e^iei-A2e^i)+O(ei5).
Using these results in
(4.17)xi+1=zi-W(μi)f[xi,wi]f(zi)f[wi,zi]f[xi,zi],
we obtain
(4.18)ei+1=e^i-[W(0)+W′(0)e^iei-W′(0)A2e^i+O(ei5)][e^i+A2e^i2+O(e^i3)]×[1+(-A23+A2A3)ei3+((-7+ω′′(0))A24-4A22A3+2A32+A2A4)ei4+O(ei5)]=e^i-[W(0)+W′(0)e^iei-W′(0)A2e^i](e^i+A2e^i2)×[1+(-A23+A2A3)ei3+((-7+ω′′(0))A24-4A22A3+2A32+A2A4)ei4]+O(ei9)=[1-W(0)]e^i-[W(0)(-A23+A2A3)ei3+W′(0)e^iei]e^i-[W(0)-W′(0)]A2e^i2-[(-7+ω′′(0))A24-4A22A3+2A32+A2A4]W(0)ei4e^i+O(ei9)=[1-W(0)]e^i-[W(0)(-A23+A2A3)+W′(0)(5-12ω′′(0)A23-A2A3)]ei3e^i+[W′(0)-W(0)]A2e^i2+[W(0)((7-ω′′(0))A24+4A22A3-2A32-A2A4)-W′(0)((-36+5ω′′(0)-16ω′′′(0))A24+(32-3ω′′(0))A22A3-2A32-2A2A4)16ω′′′(0)]×ei4e^i+O(ei9).
This means that convergence order of the family (3.10) is seventh-order with W(0)=1 and the error equation is
(4.19)ei+1=[(A23-A2A3)((5-12ω′′(0)A23-A2A3)-W′(0)((5-12ω′′(0)A23-A2A3)2]ei7+O(ei8),
and if W is any function with W(0)=1,W′(0)=1, and ω′′(0)=8, then the convergence order of any method of the family (3.10) arrives to eight, and the error equation is
(4.20)ei+1=A22(A22-A3)[(-5+16ω′′′(0))A23-4A2A3+A4]ei8+O(ei9).
Thus if ω and Ware any functions with (0)=1,ω′(0)=2,ω′′(0)=8,W(0)=1, and W′(0)=1, then the eighth-order convergence is established. This completes the proof of the theorem.
Note that per iteration every method of the family (3.10) uses four pieces of information, namely, f(xi),f′(xi),f(wi),f(zi) and has eighth-order convergence with the conditions ω(0)=1,ω′(0)=2,ω′′(0)=8,W(0)=1, andW′(0)=1, which is in accordance with Kung-Traub conjecture for 4 evaluations.
5. Some Particular Forms
Here, we consider some forms of the functions ω(t) and W(t) satisfying the conditions of the Theorem 4.1. Based on these forms some methods of the family (3.10) are also presented.
5.1. Forms of ω(t)Form 1.
For the function ω given by the following:
(5.1)ω1(t)=1+2t+4t2+αt3,
where α∈R is a constant, it is clear that the conditions of Theorem 4.1 are satisfied.
Form 2.
For the function ω defined by the following:
(5.2)ω2(t)=1+2t1-2t+αt2,
where α∈R, it can be easily seen that this function satisfies the conditions of Theorem 4.1.
Form 3.
For the function ω defined by the following:
(5.3)ω3(t)=(1+2αt+βt2)1/α,β=2α(α+1),
where α∈R-{0}, it can be seen that this function also satisfies the conditions of Theorem 4.1.
5.2. Forms of W(t)Form 1.
For the function W given by the following:
(5.4)W1(t)=1+t+γt2,
where γ∈R, it can be seen the function W1(t) satisfies the conditions of Theorem 4.1.
Form 2.
For the function W defined by the following:
(5.5)W2(t)=1+t1+γt,
where γ∈R, it is simple to see that W2(t) satisfies the conditions of Theorem 4.1
Form 3.
For the function Wdefined by the following:
(5.6)W3(t)=(1+γt)1/γ,
where γ∈R-{0}, again it can be seen that W3(t) satisfies the conditions of Theorem 4.1.
5.3. Forms of Methods
To form a concrete method we can take any combination of the above defined ω(t) and W(t). For simplicity, we consider only three such combinations. For example, by taking ω2(t) with Wi(t),i=1,2,3 the following methods can be formed
Method 1.
Taking ω2(t) and W1(t), we get a new two-parameter family of eighth-order methods
(5.7)wi=xi-f(xi)f′(xi),zi=wi-f2(xi)+αf2(wi)f2(xi)-2f(xi)f(wi)+αf2(wi)f(wi)f′(xi),xi+1=zi-[1+f(zi)f(xi)+γf2(zi)f2(xi)]f[xi,wi]f(zi)f[wi,zi]f[xi,zi].
Method 2.
Considering ω2(t) and W2(t), we get another new two-parameter family of eighth-order methods
(5.8)wi=xi-f(xi)f′(xi),zi=wi-f2(xi)+αf2(wi)f2(xi)-2f(xi)f(wi)+αf2(wi)f(wi)f′(xi),xi+1=zi-[f(xi)+(1+γ)f(zi)f(xi)+γf(zi)]f[xi,wi]f(zi)f[wi,zi]f[xi,zi].
Method 3.
Considering now ω2(t) and W3(t), we get another new two-parameter family of eighth-order methods
(5.9)wi=xi-f(xi)f′(xi),zi=wi-f2(xi)+αf2(wi)f2(xi)-2f(xi)f(wi)+αf2(wi)f(wi)f′(xi),xi+1=zi-[1+γf(zi)f(xi)]1/γf[xi,wi]f(zi)f[wi,zi]f[xi,zi].
The proposed families require three evaluations of the function f and one evaluation of first derivative f′per iteration, and achieve eighth-order convergence. Thus the efficiency index (E) defined by (2.3) of the present methods (3.10) is E=84≈1.682 which is better than E=2≈1.414 of Newton’s method, E=43≈1.587 of King’s [3] and Ostrowski’s [4] methods, E=64≈1.565 of sixth-order methods [10–12] and E=74≈1.627 of seventh-order methods [15, 18].
6. Numerical Examples
We employ the present methods (4.1), and (4.4) denoted by M81, M82 and M83, respectively to solve some nonlinear equations and compare with Newton’s method (NM) defined by (1.1), the eighth-order method developed by Cordero et al. [23] denoted by C8 and defined as follows:
(6.1)wi=xi-f(xi)f′(xi),zi=xi-f(xi)-f(wi)f(xi)-2f(wi)f(xi)f'(xi),ui=zi-(f(xi)-f(wi)f(xi)-2f(wi)+12f(zi)f(wi)-2f(zi))f(zi)f′(xi),xi+1=ui-3(β2+β3)(ui-zi)β1(ui-zi)+β2(wi-xi)+β3(zi-xi)f(zi)f′(xi),
where β1,β2,β3∈R and β2+β3≠0, eighth-order method developed by Liu and Wang [22] denoted by L8 and defined as follows:
(6.2)wi=xi-f(xi)f′(xi),zi=wi-f(xi)f(xi)-2f(wi)f(wi)f′(xi),xi+1=zi-[(f(xi)-f(wi)f(xi)-2f(wi))2+f(zi)f(wi)-α1f(zi)+4f(zi)f(xi)+α2f(zi)]f(zi)f′(xi),
where α1,α2∈R, eighth-order method developed by Petković et al. [20] denoted by P8 and defined as follows:
(6.3)wi=xi-f(xi)f′(xi),zi=wi-f(xi)f(xi)-2f(wi)f(wi)f′(xi),xi+1=zi-(1+a4(zi-xi))2f(zi)a2-a1a4+a3(zi-xi)(2+a4(zi-xi)),
where a1=f(xi), a3=(f′(xi)f[wi,zi]-f[xi,wi]f[xi,zi])/((xif[wi,zi]+(wif(zi)-zif(wi))/(wi-zi))-f(xi)),
(6.4)a4=a3f[xi,wi]+f′(xi)-f[xi,wi](wi-xi)f[xi,wi],a2=f′(xi)+a4f(xi),
eighth-order method developed by Thukral and Petković [21] denoted by T8 and defined as follows:
(6.5)wi=xi-f(xi)f′(xi),zi=wi-f(xi)f(xi)-2f(wi)f(wi)f′(xi),xi+1=zi-[f2(xi)f2(xi)-2f(xi)f(wi)-f2(wi)+f(zi)f(wi)-αf(zi)+4f(zi)f(xi)]f(zi)f′(xi),
where α∈R, eighth-order methods presented in Section 3 of [19] by Bi et al. denoted by B81, B82 and B83.
The test functions and root r correct up to 16 decimal places are displayed in Table 1. The first eight functions we have selected are same as in [19]. The last function is selected from [18]. In Table 2, we exhibit the absolute values of the difference of root r and its approximation xi, where r is computed with 350 significant digits and xi is calculated by costing the same total number of function evaluations (TFE) for each method. The TFE is counted as sum of the number of evaluations of the function itself plus the number of evaluations of the derivatives. In the calculations, 12 TFE are used by each method. That means 6 iterations are used for NM and 3 iterations for the remaining methods. The absolute values of the function |f(xi)| and the computational order of convergence (ρ) are also displayed in Table 2. It can be observed that the computed results, displayed in Table 2, overwhelmingly support the theory of convergence and efficiency analyses discussed in the previous sections. From the results, it can be concluded that the proposed methods are competitive with existing methods and possess quick convergence for good initial approximations. Among the eighth-order methods, we are not able to select one as the best. For some initial guess one is better while for other initial guess the another one would be appropriate. Thus the present methods can be of practical interest.
Test functions.
f(x)
r
f1(x)=x5+x4+4x2-15
1.3474280989683050
f2(x)=sin(x)-x/3
2.2788626600758283
f3(x)=10xe-x2-1
1.6796306104284499
f4(x)=cos(x)-x
0.7390851332151606
f5(x)=e-x2+x+2-1
−1.0000000000000000
f6(x)=e-x+cos(x)
1.7461395304080124
f7(x)=ln(x2+x+2)-x+1
4.1525907367571583
f8(x)=arcsin(x2-1)-x/2+1
0.5948109683983692
f9(x)=xex2-sin2x+3cosx+5
−1.2076478271309189
Comparison of methods using same total number of function evaluations for all methods (TFE = 12).
NM
C8
L8
P8
T8
B81
B82
B83
M81
M82
M83
β1=β3=0,β2=1
α1=α2=1
α=0
α=1
α=γ=1
α=1
α=γ=1
α=γ=1
α=γ=1
f1, x0 = 1.6
|xi-r|
9.75e-41
4.35e-338
4.51e-300
0.00e+00
3.67e-305
6.55e-304
8.24e-304
9.60e-304
0.00e+00
1.40e-349
1.70e-350
|f(xi)|
3.61e-39
1.61e-336
1.67e-298
0.00e+00
1.36e-303
2.43e-302
3.05e-302
3.55e-302
6.96e-350
5.19e-348
6.28e-349
ρ
2.0000000
7.9999999
7.9999999
8.0000000
7.9999999
7.9999989
7.9999989
7.9999989
7.9999997
7.9999997
7.9999997
f2, x0 = 2.0
|xi-r|
4.27e-57
0.00e+00
4.76e-334
0.00e+00
2.88e-333
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
|f(xi)|
4.20e-57
0.00e+00
4.68e-334
0.00e+00
2.83e-333
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
ρ
2.0000000
8.0000000
7.9999993
8.0000000
7.9999994
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
f3, x0 = 1.8
|xi-r|
4.41e-58
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
|f(xi)|
1.22e-57
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
ρ
2.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
7.9999997
7.9999997
7.9999997
f4, x0 = 1.0
|xi-r|
1.80e-83
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
|f(xi)|
3.00e-83
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
ρ
2.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
f5, x0 = −0.5
|xi-r|
3.46e-27
7.71e-202
8.70e-163
5.99e-246
3.75e-162
3.12e-222
1.20e-221
2.98e-221
4.57e-195
7.97e-194
1.98e-194
|f(xi)|
1.04e-26
2.31e-201
2.61e-162
1.80e-245
1.13e-161
9.36e-222
3.61e-221
8.93e-221
1.37e-194
2.39e-193
5.95e-194
ρ
2.0000000
7.9999339
7.9996776
7.9999847
7.9996400
8.0000580
8.0000594
8.0000603
7.9995556
7.9995321
7.9995437
f6, x0 = 2.0
|xi-r|
7.97e-85
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
|f(xi)|
9.24e-85
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
ρ
2.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
f7, x0 = 3.2
|xi-r|
4.66e-74
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
|f(xi)|
2.81e-74
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
0.00e+00
ρ
2.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
f8, x0 = 1.0
|xi-r|
1.68e-54
4.42e-342
4.21e-309
0.00e+00
1.20e-309
0.00e+00
5.83e-346
1.68e-341
0.00e+00
0.00e+00
0.00e+00
|f(xi)|
1.78e-54
4.68e-342
4.46e-309
0.00e+00
1.27e-309
0.00e+00
6.17e-346
1.78e-341
0.00e+00
0.00e+00
0.00e+00
ρ
2.0000000
8.0000004
8.0000018
8.0000001
8.0000021
8.0000012
8.0000009
8.0000014
8.0000001
8.0000001
8.0000001
f9, x0 = −1
|xi-r|
8.63e-33
0.00e+00
1.76e-339
0.00e+00
0.00e+00
4.78e-218
3.12e-214
8.95e-212
1.83e-216
3.85e-218
2.63e-217
|f(xi)|
1.75e-31
0.00e+00
3.58e-338
0.00e+00
0.00e+00
9.71e-217
6.34e-213
1.82e-210
3.72e-215
7.81e-217
5.34e-216
ρ
2.0000000
8.0000001
8.0000002
8.0000000
8.0000000
7.9999987
7.9999985
7.9999984
8.0001475
8.0001379
8.0001426
7. Conclusions
In this work, we have obtained a new simple and elegant family of eighth-order multipoint methods for solving nonlinear equations. Thus, one requires three evaluations of the function f and one of its first-derivative f′ per full step and therefore, the efficiency index of the present methods is 1.682 which is better than the efficiency index of Newton method, fourth-order methods, sixth-order methods, and seventh-order methods.
Many numerical applications use higher precision in their computations. In these types of applications, numerical methods of higher-order are important. The numerical results show that the methods associated with a multiprecision arithmetic floating point are very useful, because these methods yield a clear reduction in number of iterations. Finally, we conclude that the methods presented in this paper are preferable to other recognized efficient methods, namely, Newton’s method, King’s methods, sixth-order methods [10–12], seventh-order methods [15, 18], etc.
OrtegaJ. M.RheinboldtW. C.1970New York, NY, USAAcademic Press0273810ZBL0241.65046ChapraS. C.CanaleR. P.1988New York, NY, USAMcGraw-Hill Book CompanyKingR. F.A family of fourth order methods for nonlinear equations197310876879034358510.1137/0710072ZBL0266.65040OstrowskiA. M.1960New York, NY, USAAcademic PressTraubJ. F.1964Englewood Cliffs, NJ, USAPrentice-Hall0169356KungH. T.TraubJ. F.Optimal order of one-point and multipoint iteration197421643651035365710.1145/321850.321860ZBL0289.65023JarrattP.Some efficient fourth order multipoint methods for solving equations196991191240246505ZBL0188.22101EzquerroJ. A.HernándezM. A.SalanovaM. A.Construction of iterative processes with high order of convergence1998691-219120110.1080/002071698088047171649487ZBL0908.65031GutiérrezJ. M.HernándezM. A.An acceleration of Newton's method: super-Halley method20011172-322323910.1016/S0096-3003(99)00175-71806048ZBL1023.65051GrauM.Díaz-BarreroJ. L.An improvement to Ostrowski root-finding method2006173145045610.1016/j.amc.2005.04.0432203399ZBL1090.65053SharmaJ. R.GuhaR. K.A family of modified Ostrowski methods with accelerated sixth order convergence2007190111111510.1016/j.amc.2007.01.0092335433ZBL1126.65046ChunC.HamY.Some sixth-order variants of Ostrowski root-finding methods2007193238939410.1016/j.amc.2007.03.0742385796ZBL1193.65055KouJ.The improvements of modified Newton's method2007189160260910.1016/j.amc.2006.11.1152330238ZBL1122.65332KouJ.LiY.An improvement of the Jarratt method200718921816182110.1016/j.amc.2006.12.0622332135ZBL1122.65338KouJ.LiY.WangX.Some variants of Ostrowski's method with seventh-order convergence2007209215315910.1016/j.cam.2006.10.0732387122ZBL1130.41006ChunC.Some improvements of Jarratt's method with sixth-order convergence200719021432143710.1016/j.amc.2007.02.0232339734ZBL1122.65329ParhiS. K.GuptaD. K.A sixth order method for nonlinear equations20082031505510.1016/j.amc.2008.03.0372451537ZBL1154.65327BiW.RenH.WuQ.New family of seventh-order methods for nonlinear equations2008203140841210.1016/j.amc.2008.04.0482452299ZBL1154.65323BiW.RenH.WuQ.Three-step iterative methods with eighth-order convergence for solving nonlinear equations2009225110511210.1016/j.cam.2008.07.0042490175ZBL1161.65039PetkovićL. D.PetkovićM. S.DžunićJ.A class of three-point root-solvers of optimal order of convergence2010216267167610.1016/j.amc.2010.01.1232601535ZBL1188.65068ThukralR.PetkovićM. S.A family of three-point methods of optimal order for solving nonlinear equations201023392278228410.1016/j.cam.2009.10.0122577765ZBL1180.65058LiuL.WangX.Eighth-order methods with high efficiency index for solving nonlinear equations201021593449345410.1016/j.amc.2009.10.0402576835ZBL1183.65051CorderoA.TorregrosaJ. R.VassilevaM. P.Three-step iterative methods with optimal eighth-order convergence2011235103189319410.1016/j.cam.2011.01.0042773303ZBL1215.65091KhattriS. K.LogT.Derivative free algorithm for solving nonlinear equations201192216917910.1007/s00607-010-0135-72794922ZBL1232.65073KhattriS. K.ArgyrosI. K.Sixth order derivative free family of iterative methods2011217125500550710.1016/j.amc.2010.12.0212770167ZBL1229.65080GeumY. H.KimY. I.A biparametric family of optimally convergent sixteenth-order multipoint methods with their fourth-step weighting function as a sum of a rational and a generic two-variable function2011235103178318810.1016/j.cam.2011.01.0032773302ZBL1215.65093GautschiW.1997Boston, Mass, USABirkhäuser1454125WeerakoonS.FernandoT. G. I.A variant of Newton's method with accelerated third-order convergence2000138879310.1016/S0893-9659(00)00100-21791767ZBL0973.65037