# A New Family of Iterative Methods Based on an Exponential Model for Solving Nonlinear Equations

Authors: Tianbao Liu; Hengyan Li; Zaixiang Pang
Journal: Journal of Applied Mathematics
 (2013)

---

## Abstract

We present two new families of iterative methods for obtaining simple
roots of nonlinear equations. The first family is developed by fitting the modelm(x)=epx(Ax2+Bx+C) to the function f(x) and its derivative f′(x), f″(x) at a point xn. In order to remove the second derivative of the first methods, we construct the second
family of iterative methods by approximating the equation f(x)=0 around the point (xn,f(xn)) by the quadratic equation. Analysis of convergence shows that the new
methods have third-order or higher convergence. Numerical experiments show that
new iterative methods are effective and comparable to those of the well-known existing
methods.


---

## Body

## 1. Introduction

In this paper, we consider iterative methods to find a simple rootα, that is,f(α)=0andf′(α)≠0, of a nonlinear equation:
(1)f(x)=0,
wheref:I⊂R→Rfor an open intervalIis a scalar function.Many of the complex problems in science and engineering contains the function of nonlinear and transcendental nature in (1), so finding the simple roots of the nonlinear equation is one of the most important problems in numerical analysis. Numerical iterative methods are often used to obtain the approximate solution of such problems because it is not always possible to obtain its exact solution by usual algebraic process. We all know that Newton’s method is an important and basic approach for solving nonlinear equations [1, 2], its formulation is given by
(2)xn+1=xn-f(xn)f′(xn),
and this method converges quadratically. Earlier, many investigations [3–21] have been made to explain the root of nonlinear algebraic and transcendental equations.The outline of the paper is as follows. In Section2, we firstly describe a one-parameter family of third-order methods by fitting the modelm(x)=epx(Ax2+Bx+C)to the functionf(x)and its derivativef′(x),f′′(x)at a pointxn, and then we use a quadratic equation for approximating the equationf(x)=0to obtain a four-parameter family of second-derivative-free iterative methods. In Section 3, we obtain some different iterative methods by taking several parameters. In Section 4, different numerical tests confirm the theoretical results, and the new methods are comparable with other known methods and give better results in many cases. Finally, we infer some conclusions.


## 2. Development of Methods and Convergence Analysis

Consider the exponential model:(3)m(x)=epx(Ax2+Bx+C),
whereA,B,C, andpare parameters. We construct a new iteration scheme by fitting model (3) to the functionf(x) and its derivativef′(x) andf′′(x)at a pointxn. Imposing the conditions as follows
(4)m(xn)=f(xn),m′(xn)=f′(xn),m′′(xn)=f′′(xn)
at the point (xn,m(xn)) and then solving (3) and (4) forA,B, and C, we have
(5)A=12[f′′(xn)+p2f(xn)-2pf′(xn)]e-pxn,(6)B=[(f′′(xn)+p2f(xn)-2pf′(xn))f′(xn)-pf(xn)-(f′′(xn)+p2f(xn)-2pf′(xn))xn]e-pxn,(7)C=[f(xn)-12[f′′(xn)+p2f(xn)-2pf′(xn)]xn2-[f′(xn)-pf(xn)-(f′′(xn)+p2f(xn)-2pf′(xn))xn]xn[f(xn)-12]e-pxn.
From (3), (4), and (5), we take the values ofA,B, andCinto (3) and give
(8)m(x)-f(xn)epx-pxn=12[f′′(xn)+p2f(xn)-2pf′(xn)](x-xn)2epx-pxn+[f′(xn)-pf(xn)](x-xn)epx-pxn.
At the root estimatexn+1, it follows thatm(xn+1)=0. We consider reducing (8) as follows
(9)12[f′′(xn)+p2f(xn)-2pf′(xn)](x-xn)2+[f′(xn)-pf(xn)](x-xn)+f(xn)=0.
Thus from (9) we develop the iteration formula:
(10)xn+1=xn-[f′(xn)-pf(xn)]-f′(xn)1-Lf,p(xn)f′′(xn)+p2f(xn)-2pf′(xn),
where
(11)Lf,p(xn)=2f′′(xn)f(xn)+p2f2(xn)-2pf′(xn)f(xn)f′2(xn).
The square root is required in (10); however, this may cost expensively and even fail in the case1-Lf,p(xn)<0. In order to avoid the calculation of the square roots, we will derive some forms free from square roots by Taylor approximation [5].It is easy to know that Taylor approximation of1-Lf,p(xn)is
(12)1-Lf,p(xn)=∑k≥0(12k)(-Lf,p(xn))k.Using (12) in (10), we can obtain the following form:
(13)xn+1=xn-((∑k≥0(12k)[f′(xn)-pf(xn)]-f′(xn)∑k≥0m(12k)(-Lf,p(xn))k)×(f′′(xn)+p2f(xn)-2pf′(xn))-1∑k≥0(12k)).We have the convergence analysis of the methods by (13).Theorem 1.
Letα∈Ibe a simple zero of sufficiently differentiable functionf:I⊂R→Rfor an open intervalI. Ifx0is sufficiently close toα, form≥1, the methods defined by (13) are cubically convergent.Proof.
Leten=xn-α, and we use the following Taylor expansions:
(14)f(xn)=f′(α)[en+c2en2+c3en3+c4en4+O(en5)],
whereck=(1/k!)(f(k)(α)/f′(α)); furthermore, we have
(15)f′(xn)=f′(α)[1+2c2en+3c3en2+4c4en3+5c5en4+O(en5)],(16)f′′(xn)=f′(α)[2c2+6c3en+12c4en2+20c5en3+30c6en4+O(en5)].
Since (15), we obtain
(17)f′2(xn)=f′2(α)[1+4c2en+(6c3+4c22)en2+(8c4+12c2c3)en3+(10c5+16c2c4+9c32)en4+O(en5)].
Since (14), we get
(18)f2(xn)=f′2(α)[en2+2c2en3+(2c3+c22)en4+O(en5)].
From (14) and (15), we also easily have
(19)f(xn)f′(xn)=f′2(α)[en+3c2en2+(2c22+4c3)en3+(5c4+5c2c3)en4+O(en5)].
We obtain the following expression by taking into account (14) and (16):
(20)f(xn)f′′(xn)=f′2(α)[2c2en+(2c22+6c3)en2+(8c2c3+12c4)en3+(14c2c4+6c32+20c5)en4+O(en5)].
From (18), (19), and (20), we obtain
(21)2f′′(xn)f(xn)+p2f2(xn)-2pf′(xn)f(xn)=f′2(α)[(4c2-2p)en+(4c22+12c3+p2-6pc2)en2+(16c2c3+24c4+2p2c2-4pc22-8pc3)en3+(28c2c4+12c32+40c5+2p2c3+p2c22-10pc4-10pc2c3(28c2c4+12c32+40c5+2p2c3+p2c22)en4+O(en5)].
From (17) and (21), we have
(22)Lf,p(xn)=f′(α)[O(en4)(4c2-2p)en+(12c3-12c22+p2+2pc2)en2+(24c4-56c2c3-2p2c2-4pc22+4pc3+32c23(24c4-56c2c3-2p2c2-4pc22)en3+O(en4)].
From (14) and (15), we have
(23)f′(xn)-pf(xn)=f′(α)[1+(2c2-p)en+(3c3-pc2)en2+(4c4-pc3)en3+(5c5-pc4)en4+O(en5)].
Using (14), (15), and (16), we have
(24)f′′(xn)+p2f(xn)-2pf′(xn)=f′(α)[(2c2-2p)+(6c3+p2-4pc2)en+(12c4+p2c2-6pc3)en2+(20c5+p2c3-8pc4)en3+(30c6+p2c4-10pc5)en4+O(en5)].
We know that
(25)∑k≥0m(12k)(-Lf,p(xn))k=1-12Lf,p(xn)-18Lf,p(xn)2-116Lf,p(xn)3+⋯=1-(2c2-p)en+(pc2+4c22-6c3-p2)en2-(2p2c2-16c2c3-4pc3+12c4+8c23-p3)en3+O(en4).
From (15) and (25), we obtain
(26)f′(xn)∑k≥0m(12k)(-Lf,p(xn))k=f′(α)[1+pen+(3pc2-p2-3c3)en2+(7pc3-4p2c2-2c2c3+p3-8c4+2pc22)en3+O(en4)].
From (23) and (26), we have
(27)[f′(xn)-pf(xn)]-f′(xn)∑k≥0m(12k)(-Lf,p(xn))k=2(c2-p)en-(4pc2-6c3-p2)en2-(8pc3-12c4-4p2c2-2c2c3+p3+2pc22)en3+O(en4).
Using (13), (24), and (27), we have
(28)xn+1=xn-((∑k≥0m(12k)[f′(xn)-pf(xn)]-f′(xn)∑k≥0m(12k)(-Lf,p(xn))k)×(f′′(xn)+p2f(xn)-2pf′(xn))-1∑k≥0m(12k))=xn-en-(12p2+c3-pc2)en3+O(en4).
Fromen+1=xn+1-α, we have
(29)en+1=(pc2-12p2-c3)en3+O(en4),
which completes the proof.The family methods given by (13) are novel third-order methods, but the methods depend on the second derivatives in computing process, and therefore their practical applications are restricted in some cases. In recent years, several methods with free second derivatives have been developed; see [4–15] and references therein.In order to avoid the calculation of the second derivatives, we consider approximating the equationf(x)=0 around the point(xn,f(xn))by the quadratic equation inxandyin the following form [8]:
(30)a1x2+a2y2+a3xy+a4x+a5y+a6=0,
whereai∈R, i=1,2,…,6, are parameters. We impose the tangency conditions
(31)y(xn)=f(xn),y′(xn)=f′(xn),y(wn)=f(wn),
wherexnisnth iterate and
(32)wn=xn-f(xn)f′(xn).
From (30) and (31), we have
(33)f′′(xn)≈y′′(xn)=Lai(xn,wn)=([2a1+2a2f′2(xn)+2a3f′(xn)]f′2(xn)f(wn))×(a1f2(xn)+a2f′2(xn)[f(xn)-f(wn)]2+a3f(xn)f′(xn)[f(xn)-f(wn)][2a1+2a2f′2(xn)+2a3f′(xn)])-1,
wherei=1,2,3.From (33) we can approximate
(34)Lai,f,p(xn,wn)=2Lai(xn,wn)f(xn)+p2f2(xn)-2pf′(xn)f(xn)f′2(xn).
UsingLai,f,p(xn,wn)instead ofLf,p(xn)(11), we obtain a new four-parameter family of methods free from second derivative:
(35)xn+1=xn-((∑k≥0m(12k)[f′(xn)-pf(xn)]-f′(xn)∑k≥0m(12k)(-Lai,f,p(xn,wn))k)×(Lai(xn,wn)+p2f(xn)-2pf′(xn))-1∑k≥0m(12k)),
where ai (i=1,2,3), p∈R,m≥0.We also have the convergence analysis of the methods by (35).Theorem 2.
Letα∈Ibe a simple zero of sufficiently differentiable functionf:I⊂R→Rfor an open intervalI. Ifx0is sufficiently close toα, form≥1, ai (i=1,2,3), p∈R, the methods defined by (35) are at least cubically convergent; as particular cases, ifm≥2, a2=a3=p=0, a1∈R and the methods have convergence order four.Proof.
Leten=xn-α, and we use the following Taylor expansions:
(36)f(xn)=f′(α)[en+c2en2+c3en3+c4en4+O(en5)],
whereck=(1/k!)(f(k)(α)/f′(α)); furthermore, we have
(37)f2(xn)=f′2(α)[en2+2c2en3+(2c3+c22)en4+O(en5)],(38)f′(xn)=f′(α)[1+2c2en+3c3en2+4c4en3+5c5en4+O(en5)].
Dividing (36) by (38)
(39)f(xn)f′(xn)=en-c2en2+2(c22-c3)en3+(7c2c3-4c23-3c4)en4+O(en5).
From (39), we get
(40)wn=xn-f(xn)f′(xn)=α+c2en2-2(c22-c3)en3-(7c2c3-4c23-3c4)en4+O(en5).
Expandingf(wn)in Taylor’s Series aboutαand using (40), we get
(41)f(wn)=f′(α)[wn-α+c2(wn-α)2+c3(wn-α)3+c4(wn-α)4+⋯]=f′(α)[c2en2+2(c3-c22)en3+(5c23+3c4-7c2c3)en4+O(en5)].
From (36) and (41), we have
(42)f(xn)f(wn)=f′2(α)[c2en3+(2c3-c22)en4+O(en5)].
From (38), we obtain
(43)f′2(xn)=f′2(α)[1+4c2en+(6c3+4c22)en2+(8c4+12c2c3)en3+(10c5+16c2c4+9c32)en4+O(en5)].
From (38), (41), and (43) we also easily obtain
(44)[2a1+2a2f′2(xn)+2a3f′(xn)]f′2(xn)f(wn)=f′3(α)[2(a2f′2(α)+a3f′(α)+a1)c2en2+4(a2f′2(α)c3+3a2f′2(α)c22+a3f′(α)c3+2f′(α)a3c22+a1c3+a1c22)en3+2(3a2f′2(α)c4+21f′2(α)a2c2c3+5f′(α)a3c23+3a3f′(α)c4+14f′(α)a3c2c3+7a1c3c2)en4+O(en5)].
Substituting (36), (37), (38), (41), and (43) in the denominator ofLai(xn,wn), we obtain
(45)a1f2(xn)+a2f′2(xn)[f(xn)-f(wn)]2+a3f(xn)f′(xn)[f(xn)-f(wn)]=f′2(α)[+4a3c22f′(α)+3a3c3f′(α))(a1+a2f′2(α)+a3f′(α))en2+(2a1c2+4a2c2f′2(α)+3a3c2f′(α))en3+(+4a3c22f′(α)+3a3c3f′(α))2a1c3+a1c22+8a2c22f′2(α)+4a2c3f′2(α)+4a3c22f′(α)+3a3c3f′(α))en4+O(en5)].
Using (44) and (45), we have
(46)Lai(xn,wn)=f′(α)[(a2f′2(α)+a3f′(α)+a1)22c2+2(2a2f′2(α)c3+2a2f′2(α)c22+2a3f′(α)c3+f′(α)a3c22+2a1c3f′2)en/(a2f′2(α)+a3f′(α)+a1)-(+16a22f′4(α)c23-14f′3(α)a2c2c3a3-6a2f′3(α)c4a3-9f′4(α)a22c2c3+17f′3(α)a3c23a2-5f′2(α)a32c2c3-3a2f′2(α)c4a1+2f′(α)a3c23a1-3a3f′(α)c4a1+a12c23-3a22f′(α)4c4+2f′(α)2a32c23-3a32f′(α)2c4-a12c3c2-10f′2(α)a2c2c3a1-6f′(α)a3c2c3a1+13a1c23a2f′2(α))en2/(a2f′2(α)+a3f′(α)+a1)2+O(en3)].
From (36), (37), and (46), we have
(47)2Lai(xn,wn)f(xn)+p2f2(xn)-2pf′(xn)f(xn)=f′2(α)[(a2f′2(α)+a3f′(α)+a1)22(2c2-p)en+(12a2f′2(α)c22+8f′(α)a3c22+4a1c22+p2a2f′(α)2+p2a2f′(α)+p2a1+8a2f′2(α)c3+8a3f′(α)c3+8a1c3-6pc2a2f′2(α)-6pc2a3f′(α)-6pc2a1a2f′2(α)c22)en2/(a2f′2(α)+a3f′(α)+a1)-2(2pc22c2f′2(α)+2pc22a22f′4(α)-6f′2(α)a32c4-6f′4(α)a22c4+28f′4(α)a22c23+2a12c23-p2c2a22f′4(α)-p2c2a32f′2(α)-2p2c2a2f′3(α)a3-2p2c2a2f′2(α)a1-2p2c2a3f′(α)a1+2pc22a12+4pc3a12+2f′(α)2a32c23-32f′2(α)a2c2c3a1-16f′2(α)c2c2c3-40f′3(α)a2c2c3a3-6f′2(α)a2c4a1+2f′(α)a3c23a1-4f′2(α)a1c23a2-12f′3(α)a2c4a3-6f′(α)a3c4a1-6f′(α)3a3c23a2-24f′(α)a3c2c3a1-8a12c3c2-24f′4(α)a22c2c3+4pc22a2f′(α)3a3+4pc22a2f′2(α)a1+4pc22a3f′(α)a1+8pc3a2f′3(α)a3+8pc3a2f′2(α)a1+8pc3a3f′(α)a1-p2c2a12+4pc3a32f′2(α)+4pc3a22f′(α)4)en2/(a2f′2(α)+a3f′(α)+a1)2+O(en3)].
Since (43) and (47), we get
(48)Lai,f,p(xn,wn)=[(a2f′2(α)+a3f′(α)+a1)22(2c2-p)en-(4a2f′2(α)c22+8a3f′(α)c22+12a1c22-p2a2f′2(α)-p2a3f′(α)-p2a1-8a2f′2(α)c3-8a3f′(α)c3-8a1c3-2pc2a2f′2(α)-2pc2a3f′(α)-2pc2a1f′2(α)c22)en2/(a2f′2(α)+a3f′(α)+a1)+2(-2pc22a32f′2(α)-2pc22a22f′4(α)+6f′2(α)a32c4+6f′4(α)a22c4-28f′4(α)a22c23+14a12c23-p2c2a22f′4(α)-p2c2a32f′2(α)-2p2c2a2f′3(α)a3-2p2c2a2f′2(α)a1-2p2c2a3f′(α)a1-2pc22a12+2pc3a12+6f′2(α)a32c23-24f′2(α)a2c2c3a1-12f′2(α)a32c2c3-16f′3(α)a2c2c3a3+6f′2(α)a2c4a1+22f′(α)a3c23a1+20f′2(α)a1c23a2+12f′3(α)a2c4a3+6f′(α)a3c4a1+14f′3(α)a3c23a2-32f′(α)a3c2c3a1-20a12c3c2-4f′4(α)a22c2c3-4pc22a2f′3(α)a3-4pc22a2f′2(α)a1-4pc22a3f′(α)a1+4pc3a2f′3(α)a3+4pc3a2f′2(α)a1+4pc3a3f′(α)a1-p2c2a12+2pc3a32f′2(α)+2pc3a22f′4(α))en3/(a2f′2(α)+a3f′(α)+a1)2+O(en4)].
Using (48), we write
(49)∑k≥0m(12k)(-Lai,f,p(xn,wn))k=1-12Lai,f,p(xn,wn)-18Lai,f,p(xn,wn)2-116Lai,f,p(xn,wn)3+⋯=1-(2c2-p)en+12(4a2f′2(α)c22+8cf′(α)c22+12a1c22-p2a2f′2(α)-p2a3f′(α)-p2a1-8a2f′2(α)c3-8a3f′(α)c3-8a1c3-2pc2a2f′2(α)-2pc2a3f′(α)-2pc2a14a2f′2(α)c22)en2/(a2f′2(α)+a3f′(α)+a1)+O(en3).
Using (36), (38), and (46), we obtain
(50)Lai(xn,wn)+p2f(xn)-pf′(xn)=f′(α)[(a2f′2(α)+a3f′(α)+a1)22(c2-p)+(4a2f′2(α)c3+4a2f′2(α)c22+4a3f′(α)c3+2a3f′(α)c22+4a1c3+p2a2f′2(α)+p2a3f′(α)+p2a1-4pc2a2f′2(α)-4pc2a3f′(α)-4pc2a1)en/(a2f′2(α)+a3f′(α)+a1)-(6a3f′(α)c4a16pc3a12-p2c2a22f′4(α)-p2c2a12+6pf′4(α)c3a22+6pf′2(α)c3a32-2p2c2a3f′(α)a1-2p2c2a2f′2(α)a1-2p2c2a2f′3(α)a3-p2c2a32f′2(α)+12f′3(α)pc3a2a3+12pf′2(α)c3a2a1+12pf′(α)c3a3a1-18a22f′4(α)c2c3+26a1c23a2f′2(α)-20c2a2f′2(α)c3a1-12c2a3f′(α)c3a1-12f′3(α)a2c4a3-6a22f′4(α)c4+4a32f′2(α)c23-6a32f′2(α)c4-2a12c3c2+2a12c23+32a22f′4(α)c23-28c2a2f′3(α)c3a3+34f′3(α)a3c23a2-10a32f′2(α)c2c3-6a2f′2(α)c4a1+4a3f′(α)c23a1-6a3f′(α)c4a1)en2/(a2f′2(α)+a3f′(α)+a1)2+O(en3)].
Taking into account (36), (38), (49), and (50), we finally obtain
(51)xn+1=xn-((∑k≥0m(12k)(-Lai,f,p(xn,wn))k[f′(xn)-pf(xn)]-f′(xn)×∑k≥0m(12k)(-Lai,f,p(xn,wn))k)×(Lai(xn,wn)+p2f(xn)-2pf′(xn))-1∑k≥0m(12k)(-Lai,f,p(xn,wn))k)=xn-en+12(-40f′3(α)a3c23a2-30a1c23a2f′2(α)-2a3f′(α)c23a1-2a32f′2(α)c23-4a22f′4(α)c23+2c22pa12+6pf′(α)c22a3a1+8pf′2(α)c22a2a1+4pf′2(α)c22a32+10f′3(α)pc22a2a3+6pf′4(α)c22a22-3p2c2a32f′2(α)-6p2c2a2f′2(α)a1-6p2c2a3f′(α)a1-3p2c2a12-3p2c2a22f′4(α)-6p2c2a2f′3(α)a3+2f′2(α)p3a1a2+2f′3(α)p3a2a3+p3a12+2f′(α)p3a1a3+f′4(α)p3a22+f′2(α)p3a32)en3/(c2-p)/(a2f′2(α)+a3f′(α)+a1)2+O(en4).
Taking into account the last expression (51) anden+1=xn+1-α, we have
(52)en+1=12(-40f′3(α)a3c23a2-30a1c23a2f′2(α)-2a3f′(α)c23a1-2a32f′2(α)c23-4a22f′4(α)c23+2c22pa12+6pf′(α)c22a3a1+8pf′2(α)c22a2a1+4pf′2(α)c22a32+10f′3(α)pc22a2a3+6pf′4(α)c22a22-3p2c2a32f′2(α)-6p2c2a2f′2(α)a1-6p2c2a3f′(α)a1-3p2c2a12-3p2c2a22f′4(α)-6p2c2a2f′3(α)a3+2f′2(α)p3a1a2+2f′3(α)p3a2a3+p3a12+2f′(α)p3a1a3+f′4(α)p3a22+f′2(α)p3a32)en3/(c2-p)/(a2f′2(α)+a3f′(α)+a1)2+O(en4).
This means that the methods defined by (35) are at least of order three for any ai (i=1,2,3),  p∈R. Furthermore, we consider that if m≥2, a2=a3=p=0, then the methods defined by (35) are shown to converge the order four.


## 3. Some Special Cases

From (33)–(35), we have
(53)La1,a2,a3(xn,wn)=([2a1+2a2f′2(xn)+2a3f′(xn)]f′2(xn)f(wn))×(a1f2(xn)+a2f′2(xn)[f(xn)-f(wn)]2+a3f(xn)f′(xn)[f(xn)-f(wn)f′2(xn)]f′2(xn)f(wn))-1,
wherea1,a2,a3∈R. From (33) we can approximate
(54)La1,a2,a3,f,p(xn,wn)=(2La1,a2,a3(xn,wn)f(xn)+p2f2(xn)-2pf′(xn)f(xn))×(f′2(xn))-1.
Let
(55)Ka1,a2,a3,f,p(xn,wn)=(∑k≥0m‍(12k)[f′(xn)-pf(xn)]-f′(xn)∑k≥0m‍(12k)(-La1,a2,a3,f,p(xn,wn))k)×(La1,a2,a3(xn,wn)+p2f(xn)-2pf′(xn))-1,
wherea1,a2,a3,p∈R,m≥0.10:
Ifa1=a2=a3=p=1, and m=3, from (55), we obtain a third-order method (LM1):
(56)xn+1=xn-K1,1,1,f,1(xn,wn).20:
Ifa1=a2=0, a3=p=1, and m=3, from (55) we also obtain a third-order method (LM2):
(57)xn+1=xn-K0,0,1,f,1(xn,wn).30:
Ifa1=1, a2=0, a3=p=1, and m=3, from (55) we obtain a new third-order method (LM3):
(58)xn+1=xn-K1,0,1,f,1(xn,wn).40:
Ifa1=-1, a2=a3=0, p=1, and m=3, from (55) we obtain a third-order method (LM4):
(59)xn+1=xn-K-1,0,0,f,1(xn,wn).50:
Ifa1=a2=0, a3=-1, p=1, and m=3, from (55) we obtain a new third-order method (LM5):
(60)xn+1=xn-K0,0,-1,f,1(xn,wn).


## 4. Numerical Examples

In this section, some numerical examples commonly used in the literature are presented in Table1 to check the effectiveness of the new methods. The following methods were compared: Newton method (NM), the method of Weerakoon and Fernando [10] (WF), the method of Potra and Pták (PP) [11], Chebyshev’s method (CHM) [12, 13], Halley’s method (HM) [12], and our new methods (56) (LM1), (57) (LM2), (58) (LM3), (59) (LM4), and (60) (LM5). Displayed in Table 1 are the number of iterations (IT), the number of function evaluations (NFE) counted as the sum of the number of evaluations of the function itself plus the number of evaluations of the derivative, the valuef(xn+1), the computing time (TIME, the unit of time is one second), and the distance of two consecutive approximationsδ=|xn+1-xn|. All computations were done using Matlab 7.1 environment with a ADM athlon (tm) II X2 250-3.01 GHz based PC. We accept an approximate solution rather than the exact root, depending on the precisionϵof the computer. We use the following stopping criteria for computer programs:|f(xn+1)|<ϵ, we used the fixed stopping criterionϵ=10-15.Table 1
Comparison of various third-order methods and Newton’s method.
IT
NFE
f
(
x
n
+
1
)
Time
δ
f
1
:
x
0
=
1
NM
5
10
0
0.062577
2.126987475037367
e
-
011
WF
3
9
0
0.018731
2.284722713019605
e
-
006
PP
4
12
0
0.044520
1.558753126573720
e
-
013
CHM
4
12
0
0.040102
1.643130076445232
e
-
014
HM
3
9
0
0.019876
3.698649917449615
e
-
007
LM1
4
12
0
0.029137
1.043609643147647
e
-
014
LM2
3
9
0
0.022242
1.788683664960544
e
-
006
LM3
3
9
0
0.023550
1.295199573814188
e
-
006
LM4
3
9
0
0.019227
3.091731315407742
e
-
010
LM5
3
9
0
0.020131
1.788683664960544
e
-
006
f
1
:
x
0
=
2
NM
5
10
0
0.057231
5.020497351182485
e
-
010
WF
4
12
0
0.039384
4.440892098500626
e
-
016
PP
4
12
0
0.044239
7.949196856316121
e
-
014
CHM
4
12
0
0.042782
2.065014825802791
e
-
014
HM
3
9
0
0.017733
3.107350415199051
e
-
006
LM1
4
12
0
0.032638
2.706235235905297
e
-
011
LM2
4
12
0
0.037830
5.551115123125783
e
-
015
LM3
4
12
0
0.034087
3.108624468950438
e
-
015
LM4
4
12
0
0.038399
2.220446049250313
e
-
016
LM5
4
12
0
0.035809
5.551115123125783
e
-
015
f
2
:
x
0
=
0
NM
4
8
0
0.035702
2.665312415217613
e
-
012
WF
3
9
0
0.019699
7.801814749797131
e
-
012
PP
3
9
0
0.021596
1.219191414492116
e
-
012
CHM
3
9
0
0.023312
8.906764215055318
e
-
013
HM
3
9
0
0.017431
7.374600929921371
e
-
012
LM1
3
9
0
0.016351
9.004856806882344
e
-
007
LM2
3
9
0
0.017242
8.067152316715287
e
-
007
LM3
3
9
0
0.018978
8.334935165388302
e
-
007
LM4
3
9
0
0.016178
7.334822825222354
e
-
007
LM5
3
9
0
0.016795
8.067152316715287
e
-
007
f
2
:
x
0
=
0.5
NM
4
8
0
0.044392
1.791899961745003
e
-
013
WF
3
9
0
0.023132
6.424749621203318
e
-
012
PP
3
9
0
0.024293
4.607425552194400
e
-
014
CHM
3
9
0
0.026689
3.087480271446452
e
-
011
HM
3
9
0
0.022199
4.208039472430869
e
-
011
LM1
3
9
0
0.024120
2.850281847210923
e
-
007
LM2
3
9
4.440892098500626
e
-
016
0.021035
2.744941123289379
e
-
007
LM3
3
9
0
0.021310
2.778229745148408
e
-
007
LM4
3
9
0
0.020886
2.655207446689012
e
-
007
LM5
3
9
4.440892098500626
e
-
016
0.021592
2.744941123289379
e
-
007
f
3
:
x
0
=
1
NM
7
14
3.537126081266182
e
-
024
0.085178
1.085848323840232
e
-
012
WF
4
12
2.621304391538411
e
-
016
0.030341
4.330310691887267
e
-
006
PP
5
15
8.196910187379942
e
-
033
0.063971
8.806888499109001
e
-
012
CHM
5
15
6.352230116524407
e
-
022
0.053064
2.520356663650445
e
-
011
HM
5
15
7.257520328033309
e
-
029
0.055315
8.459855063117184
e
-
015
LM1
5
15
5.355132874954819
e
-
017
0.041370
4.825251676864564
e
-
007
LM2
4
12
6.348042447121620
e
-
018
0.026722
5.276208195236892
e
-
013
LM3
3
9
1.559423551656343
e
-
017
0.019476
5.545966966175517
e
-
011
LM4
4
12
6.485252546346663
e
-
018
0.030513
1.317913979899399
e
-
009
LM5
4
12
6.348042447121620
e
-
018
0.031575
5.276208195236892
e
-
013
f
3
:
x
0
=
0.5
NM
6
12
5.905159674954809
e
-
020
0.071738
1.402992074360412
e
-
010
WF
4
12
1.764824578467612
e
-
017
0.039948
4.200981459101664
e
-
009
PP
4
12
1.694834561079519
e
-
019
0.056046
2.767209186089879
e
-
007
CHM
4
12
8.798671206634291
e
-
017
0.036458
3.059650585008672
e
-
007
HM
4
12
3.044907255908736
e
-
017
0.072238
5.518067735074518
e
-
009
LM1
4
12
1.293991275513200
e
-
017
0.033751
1.322689611777333
e
-
013
LM2
4
12
2.828488650436813
e
-
016
0.032997
4.626219030783813
e
-
006
LM3
4
12
6.789745028898642
e
-
018
0.036425
1.912594303370619
e
-
014
LM4
4
12
1.506225012219799
e
-
017
0.034137
2.365541379047147
e
-
011
LM5
4
12
2.828488650436813
e
-
016
0.035087
4.626219030783813
e
-
006
f
4
:
x
0
=
2.5
NM
6
12
0
0.071353
1.154631945610163
e
-
014
WF
4
12
0
0.034140
7.314593375440381
e
-
012
PP
4
12
0
0.042598
4.221685223626537
e
-
010
CHM
4
12
0
0.036365
9.853584614916144
e
-
011
HM
4
12
0
0.038080
4.662936703425658
e
-
014
LM1
3
9
0
0.018372
9.336336148635382
e
-
010
LM2
4
12
0
0.035562
1.610311883837312
e
-
010
LM3
4
12
0
0.036985
2.312174895990893
e
-
010
LM4
4
12
0
0.035463
5.127454016928823
e
-
012
LM5
4
12
0
0.033362
1.610311883837312
e
-
010
f
4
:
x
0
=
3.5
NM
7
14
0
2.151288
2.877564853065451
e
-
011
WF
5
15
0
0.054803
6.550315845288424
e
-
013
PP
5
15
0
0.063380
4.512221707386743
e
-
010
CHM
5
15
0
0.061753
4.188738245147761
e
-
011
HM
4
12
0
0.046608
4.485352507632712
e
-
006
LM1
5
15
6.661338147750939
e
-
016
0.050909
2.781577124189028
e
-
008
LM2
5
15
0
0.048688
9.426681657487279
e
-
011
LM3
10
30
0
0.111640
1.081885248055414
e
-
009
LM4
4
12
0
0.040109
2.154232348061669
e
-
011
LM5
5
15
0
0.043482
9.426681657487279
e
-
011
f
5
:
x
0
=
0
NM
5
10
0
0.050123
1.701233598438989
e
-
010
WF
3
9
0
0.023397
7.792236328407753
e
-
007
PP
4
12
0
0.038373
1.500558566291943
e
-
010
CHM
4
12
0
0.037498
5.327979279989847
e
-
009
HM
4
12
0
0.035701
1.121325254871408
e
-
014
LM1
4
12
3.330669073875470
e
-
016
0.028362
4.678117765610779
e
-
006
LM2
4
12
0
0.028666
7.181037887660224
e
-
007
LM3
6
18
0
0.053924
6.676881270095691
e
-
013
LM4
4
12
0
0.028247
2.163133006050089
e
-
010
LM5
4
12
0
0.029150
7.181037887660224
e
-
007
f
5
:
x
0
=
1
NM
4
8
0
0.039328
1.701233598438989
e
-
010
WF
2
6
4.440892098500626
e
-
016
0.003092
2.674277017133964
e
-
005
PP
3
9
0
0.020461
9.809075773858922
e
-
011
CHM
3
9
0
0.019378
1.600380383770528
e
-
009
HM
3
9
0
0.018186
6.624212289807474
e
-
010
LM1
3
9
0
0.018199
7.500783549829748
e
-
008
LM2
3
9
1.110223024625157
e
-
016
0.017732
5.330569952111119
e
-
008
LM3
3
9
0
0.013523
8.804162354714151
e
-
008
LM4
3
9
0
0.013782
3.531516445942629
e
-
008
LM5
3
9
1.110223024625157
e
-
016
0.016581
5.330569952111119
e
-
008
f
6
:
x
0
=
1
NM
6
12
3.330669073875470
e
-
016
0.065520
3.059774655866931
e
-
013
WF
4
12
4.440892098500626
e
-
016
0.030929
1.793023507445923
e
-
010
PP
16
48
4.440892098500626
e
-
016
0.266014
1.531728257564424
e
-
007
CHM
5
15
4.440892098500626
e
-
016
0.061006
6.883094094689568
e
-
010
HM
4
12
4.440892098500626
e
-
016
0.043655
2.686739719592879
e
-
013
LM1
4
12
3.330669073875470
e
-
016
0.030877
1.659126835917846
e
-
008
LM2
4
12
3.330669073875470
e
-
016
0.031156
9.103828801926284
e
-
014
LM3
4
12
3.330669073875470
e
-
016
0.033076
1.261090962767497
e
-
006
LM4
5
15
4.440892098500626
e
-
016
0.041102
8.215650382226158
e
-
014
LM5
4
12
3.330669073875470
e
-
016
0.030923
9.103828801926284
e
-
014
f
6
:
x
0
=
2.5
NM
6
12
3.330669073875470
e
-
016
0.064979
1.404654170755748
e
-
012
WF
4
12
3.330669073875470
e
-
016
0.038404
4.229505634611996
e
-
012
PP
4
12
3.330669073875470
e
-
016
0.042164
1.030850205196998
e
-
008
CHM
4
12
3.330669073875470
e
-
016
0.043437
1.475204565171140
e
-
007
HM
4
12
4.440892098500626
e
-
016
0.043896
9.462626682221753
e
-
009
LM1
5
15
3.330669073875470
e
-
016
0.053132
6.201483770951199
e
-
012
LM2
4
12
3.330669073875470
e
-
016
0.036191
4.450250368215336
e
-
008
LM3
6
18
3.330669073875470
e
-
016
0.064485
8.881784197001252
e
-
016
LM4
4
12
3.330669073875470
e
-
016
0.035393
5.568381311604753
e
-
009
LM5
4
12
3.330669073875470
e
-
016
0.036946
4.450250368215336
e
-
008
f
7
:
x
0
=
3.25
NM
8
16
0
0.114189
9.720393379097914
e
-
010
WF
6
18
0
0.065645
1.691979889528739
e
-
013
PP
6
18
0
0.073469
1.131490456884876
e
-
010
CHM
6
18
0
0.082677
2.398081733190338
e
-
014
HM
5
15
0
0.055630
3.082423205569285
e
-
012
LM1
4
12
0
0.036191
2.613820271335499
e
-
011
LM2
5
15
0
0.049599
4.440892098500626
e
-
016
LM3
4
12
0
0.038172
1.546353226355990
e
-
006
LM4
6
18
0
0.051577
4.440892098500626
e
-
016
LM5
5
15
0
0.048340
4.440892098500626
e
-
016
f
7
:
x
0
=
3.45
NM
11
22
0
0.147308
4.008793297316515
e
-
011
WF
8
24
0
0.106276
7.105427357601002
e
-
015
PP
8
24
0
0.116999
2.160227552394645
e
-
010
CHM
7
21
0
0.100725
1.268533917908599
e
-
006
HM
6
18
0
0.062876
1.694565332499565
e
-
008
LM1
6
18
0
0.060920
4.440892098500626
e
-
016
LM2
6
18
0
0.060951
1.753264200488047
e
-
011
LM3
6
18
0
0.060627
1.518252190635394
e
-
011
LM4
7
21
0
0.087873
9.144596191390519
e
-
011
LM5
6
18
0
0.060607
1.753264200488047
e
-
011We used the following test functions and display the computed approximate zerox* [16]:
(61)f1(x)=x3+4x2-10,x*=1.3652300134140969,f2(x)=x2-ex-3x+2,x*=0.25753028543986076,f3(x)=sin(x)ex+ln(1+x2),x*=0,f4(x)=(x-1)3-1,x*=2,f5(x)=cosx-x,x*=0.73908513321516067,f6(x)=sin2x-x2+1,x*=1.4044916482153411,f7(x)=ex2+7x-30-1,x*=3.


## 5. Conclusions

In this paper, we presented two new families of iterative methods for solving nonlinear equations. One is developed by fitting the modelm(x)=epx(Ax2+Bx+C)to the functionf(x)and its derivativef′(x),f′′(x)at a pointxn. The other family of iterative methods was constructed by approximating the equationf(x)=0around the point(xn,f(xn))with the quadratic equation to avoid the calculation of the second derivatives. Analysis of convergence shows that the new methods have third-order or higher convergence: if m≥2, a2=a3=p=0, then the methods defined by (35) are shown to converge the order four. We observed from numerical examples that the proposed methods are efficient and demonstrate equal or better performance as compared with other well-known methods.



---

*Source: 547438-2013-12-18.xml*