# An Efficient Family of Root-Finding Methods with Optimal Eighth-Order Convergence

Authors: Rajni Sharma; Janak Raj Sharma
Journal: Advances in Numerical Analysis
 (2012)

---

## Abstract

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 order2n-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.


---

## Body

## 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 rootr 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 functionf 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 Section2. 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 Definitions

Definition 2.1.
Letf(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.
Ifp=1, 2 or 3, the sequence is said to have linear convergence, quadratic convergence or cubic convergence, respectively.Definition 2.2.
Letei=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.
Letn 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 thatxi+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 ofa 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 obtainb 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 derivativef'(x) with the derivative y'(x) of rational function (3.2) and obtain
(3.8)f′(zi)≈y′(zi).Substituting the values ofb 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 functionf: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.
Letei=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 ofW(t)

Form 1.
For the functionW 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 functionW 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.1Form 3.
For the functionWdefined 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 formedMethod 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].


## 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 ofW(t)

Form 1.
For the functionW 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 functionW 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.1Form 3.
For the functionWdefined 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 formedMethod 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 rootr 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.Table 1
Test functions.
f
(
x
)
r
f
1
(
x
)
=
x
5
+
x
4
+
4
x
2
-
15
1
.
3474280989683050
f
2
(
x
)
=
sin
(
x
)
-
x
/
3
2.2788626600758283
f
3
(
x
)
=
10
x
e
-
x
2
-
1
1
.
6796306104284499
f
4
(
x
)
=
cos
(
x
)
-
x
0.7390851332151606
f
5
(
x
)
=
e
-
x
2
+
x
+
2
-
1
−1.0000000000000000
f
6
(
x
)
=
e
-
x
+
cos
(
x
)
1.7461395304080124
f
7
(
x
)
=
ln
(
x
2
+
x
+
2
)
-
x
+
1
4.1525907367571583
f
8
(
x
)
=
arcsin
(
x
2
-
1
)
-
x
/
2
+
1
0.5948109683983692
f
9
(
x
)
=
x
e
x
2
-
si
n
2
x
+
3
cos
x
+
5
−1.2076478271309189Table 2
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
f
1, x0 = 1.6
|
x
i
-
r
|
9.75
e
-
41
4.35
e
-
338
4.51
e
-
300
0.00
e
+
00
3.67
e
-
305
6.55
e
-
304
8.24
e
-
304
9.60
e
-
304
0.00
e
+
00
1.40
e
-
349
1.70
e
-
350
|
f
(
x
i
)
|
3.61
e
-
39
1.61
e
-
336
1.67
e
-
298
0.00
e
+
00
1.36
e
-
303
2.43
e
-
302
3.05
e
-
302
3.55
e
-
302
6.96
e
-
350
5.19
e
-
348
6.28
e
-
349
ρ
2.0000000
7.9999999
7.9999999
8.0000000
7.9999999
7.9999989
7.9999989
7.9999989
7.9999997
7.9999997
7.9999997
f
2, x0 = 2.0
|
x
i
-
r
|
4.27
e
-
57
0.00
e
+
00
4.76
e
-
334
0.00
e
+
00
2.88
e
-
333
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
|
f
(
x
i
)
|
4.20
e
-
57
0.00
e
+
00
4.68
e
-
334
0.00
e
+
00
2.83
e
-
333
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
ρ
2.0000000
8.0000000
7.9999993
8.0000000
7.9999994
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
f
3, x0 = 1.8
|
x
i
-
r
|
4.41
e
-
58
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
|
f
(
x
i
)
|
1.22
e
-
57
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
ρ
2.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
7.9999997
7.9999997
7.9999997
f
4, x0 = 1.0
|
x
i
-
r
|
1.80
e
-
83
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
|
f
(
x
i
)
|
3.00
e
-
83
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
ρ
2.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
f
5, x0 = −0.5
|
x
i
-
r
|
3.46
e
-
27
7.71
e
-
202
8.70
e
-
163
5.99
e
-
246
3.75
e
-
162
3.12
e
-
222
1.20
e
-
221
2.98
e
-
221
4.57
e
-
195
7.97
e
-
194
1.98
e
-
194
|
f
(
x
i
)
|
1.04
e
-
26
2.31
e
-
201
2.61
e
-
162
1.80
e
-
245
1.13
e
-
161
9.36
e
-
222
3.61
e
-
221
8.93
e
-
221
1.37
e
-
194
2.39
e
-
193
5.95
e
-
194
ρ
2.0000000
7.9999339
7.9996776
7.9999847
7.9996400
8.0000580
8.0000594
8.0000603
7.9995556
7.9995321
7.9995437
f
6, x0 = 2.0
|
x
i
-
r
|
7.97
e
-
85
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
|
f
(
x
i
)
|
9.24
e
-
85
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
ρ
2.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
f
7, x0 = 3.2
|
x
i
-
r
|
4.66
e
-
74
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
|
f
(
x
i
)
|
2.81
e
-
74
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
ρ
2.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
8.0000000
f
8, x0 = 1.0
|
x
i
-
r
|
1.68
e
-
54
4.42
e
-
342
4.21
e
-
309
0.00
e
+
00
1.20
e
-
309
0.00
e
+
00
5.83
e
-
346
1.68
e
-
341
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
|
f
(
x
i
)
|
1.78
e
-
54
4.68
e
-
342
4.46
e
-
309
0.00
e
+
00
1.27
e
-
309
0.00
e
+
00
6.17
e
-
346
1.78
e
-
341
0.00
e
+
00
0.00
e
+
00
0.00
e
+
00
ρ
2.0000000
8.0000004
8.0000018
8.0000001
8.0000021
8.0000012
8.0000009
8.0000014
8.0000001
8.0000001
8.0000001
f
9, x0 = −1
|
x
i
-
r
|
8.63
e
-
33
0.00
e
+
00
1.76
e
-
339
0.00
e
+
00
0.00
e
+
00
4.78
e
-
218
3.12
e
-
214
8.95
e
-
212
1.83
e
-
216
3.85
e
-
218
2.63
e
-
217
|
f
(
x
i
)
|
1.75
e
-
31
0.00
e
+
00
3.58
e
-
338
0.00
e
+
00
0.00
e
+
00
9.71
e
-
217
6.34
e
-
213
1.82
e
-
210
3.72
e
-
215
7.81
e
-
217
5.34
e
-
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 functionf 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.



---

*Source: 346420-2012-10-18.xml*