Table of Contents Author Guidelines Submit a Manuscript
International Journal of Mathematics and Mathematical Sciences
VolumeΒ 2012Β (2012), Article IDΒ 932420, 18 pages
Research Article

Optimized Steffensen-Type Methods with Eighth-Order Convergence and High Efficiency Index

Department of Applied Mathematics, School of Mathematical Sciences, Ferdowsi University of Mashhad, Mashhad, Iran

Received 21 March 2012; Revised 23 May 2012; Accepted 6 June 2012

Academic Editor: V. R.Β Khalilov

Copyright Β© 2012 F. Soleymani. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Steffensen-type methods are practical in solving nonlinear equations. Since, such schemes do not need derivative evaluation per iteration. Hence, this work contributes two new multistep classes of Steffensen-type methods for finding the solution of the nonlinear equation 𝑓(π‘₯)=0. New techniques can be taken into account as the generalizations of the one-step method of Steffensen. Theoretical proofs of the main theorems are furnished to reveal the eighth-order convergence. Per computing step, the derived methods require only four function evaluations. Experimental results are also given to add more supports on the underlying theory of this paper as well as lead us to draw a conclusion on the efficiency of the developed classes.

1. Introduction

Finding rapidly and accurately the zeros of nonlinear functions is a common problem in various areas of numerical analysis. This problem has fascinated mathematicians for centuries, and the literature is full of ingenious methods, and discussions of their merits and shortcomings [1–13].

Over the last decades, there exist a large number of different methods for finding the solution of nonlinear equations either iteratively or simultaneously. Thus far, some better modified methods for finding roots of nonlinear equations cover mainly the pioneering work of Kung and Traub [14] or the efficient ways to build high-order methods through the procedures discussed in [15–18].

In this paper, we consider the problem of finding a numerical procedure to solve a simple root 𝛼 of the nonlinear equation: 𝑓(π‘₯)=0,(1.1) by derivative-free high-order iterations without memory. We here note that higher order of convergence is only possible through employing multistep methods. In what follows, we first give two definitions concerning iterative processes.

Definition 1.1. Let the sequence {π‘₯𝑛} tends to 𝛼 such that limπ‘›β†’βˆžξ€·π‘₯𝑛+1ξ€Έβˆ’π›Όξ€·π‘₯π‘›ξ€Έβˆ’π›Όπ‘=𝐢≠0,𝑛β‰₯1.(1.2) Therefore, the order of convergence of the sequence {π‘₯𝑛} is 𝑝, and 𝐢 is known as the asymptotic error constant. If 𝑝=1, 𝑝=2, or 𝑝=3, the sequence is said to converge linearly, quadratically, or cubically, respectively.

Definition 1.2. The efficiency of a method [2] is measured by the concept of efficiency index and is defined by EI=𝑝1/𝛽,(1.3) where 𝑝 is the order of the method and 𝛽 is the whole number of (functional) evaluations per iteration. Note that in (1.3) we consider that all function and derivative evaluations have the same computational cost. Moreover, we should remind that by Kung-Traub conjecture [14]: an iterative scheme without memory for solving nonlinear equations has the optimal efficiency index 2(π›½βˆ’1)/𝛽 and optimal rate (speed) of convergence 2(π›½βˆ’1).

The so-called method of Newton is basically taken into account to solve (1.1). This scheme is an optimal one-step one-point iteration without memory. The method is in fact the essence of all the improvements of root solvers. We should remind that iterations are themselves divided into two main categories of derivative-involved methods (in which at least one derivative evaluation is needed to proceed; see, e.g., [19, 20]) and derivative-free methods (do not have a direct derivative evaluation per iteration [21–23]) which are more economic in terms of derivative evaluation. Clearly, the arising real problems in science and engineering are not normally of the type to let the users to calculate their first or second derivatives. Since, they include hard structures, which make the procedure of finding the derivatives so difficult or time consuming. Due to this, derivative-free methods now comes to attention to solve the related problems as easily as possible. For further reading, one may refer to [24–28].

This paper tries to overcome on the matter of derivative evaluation for nonlinear solvers by giving two general classes of three-step eighth-order convergence methods, which have the optimal efficiency index, optimal order, accurate performance in solving numerical examples as well as are totally free from derivative calculation per full cycle to proceed. Hence, after this short introduction in this section, we organize the rest of the paper as comes next. Section 2 shortly provides one of the most applicable usages of root solvers in Mechanical Engineering in order to manifest their applicability. This section is followed by Section 3, where some of the available derivative-free schemes in root-finding are discussed and presented. Section 4 gives the heart of this paper by contributing two novel classes of three-step derivative-free without memory iterations. The theoretical proofs of the main theorems are given therein to show that each member of our classes without memory attain as much as possible of the efficiency index by using as small as possible of the number of evaluations. In Section 5, we give large number of numerical experiments to reveal the accuracy of the obtained methods from the proposed classes. And finally, Section 6 contains a short conclusion of this study.

2. Describing an Application of Nonlinear Equations Solvers

In Mechanical Engineering, a trunnion (a cylindrical protrusion used as a mounting and/or pivoting point; in a cannon, the trunnions are two projections cast just forward of the center of mass of the cannon are fixed to a two-wheeled movable gun carriage) has to be cooled before it is shrink fitted into a steel hub.

The equation that gives the temperature 𝑇𝑓 to which the trunnion has to be cooled to attain the desired contraction is given by 𝑓𝑇𝑓=βˆ’0.50598Γ—10βˆ’10𝑇3𝑓+0.38292Γ—10βˆ’7𝑇2𝑓+0.74363Γ—10βˆ’4𝑇𝑓+0.88318Γ—10βˆ’2=0.(2.1)

Clearly (2.1) could be solved using nonlinear equation solvers. This was one application of the matter of nonlinear equation solving by iterative processes in the scalar case. At this time, consider that the obtained nonlinear scalar equations in other problems are complicated; therefore their first and second derivatives are not at hand and moreover, better accuracy by low elapsed time is needed. Consequently, a close look at the derivative-free high-order methods should be paid by considering multi step iterations. Also note that such iterative schemes can be extended to solve systems of nonlinear equations which have a lot of applications in engineering problems, see for more [29].

3. Available Derivative-Free Methods in the Literature

For the first time, Steffensen in [30] coined the following derivative-free form of Newton's method: π‘₯𝑛+1=π‘₯π‘›βˆ’π‘“ξ€·π‘₯𝑛𝑓π‘₯𝑛,𝑀𝑛,𝑀𝑛=π‘₯𝑛π‘₯+𝑓𝑛,(3.1) which possesses the same rate of convergence and efficiency index as Newton's.

We here remind the well-written family of derivative-free three-step methods, which was given by Kung and Traub [14] as comes next 𝑧𝑛=𝑦𝑛𝑓π‘₯βˆ’π›½π‘›ξ€Έπ‘“ξ€·π‘¦π‘›ξ€Έπ‘“ξ€·π‘¦π‘›ξ€Έξ€·π‘₯βˆ’π‘“π‘›ξ€Έ,𝑦𝑛=π‘₯𝑛π‘₯+𝛽𝑓𝑛,𝑀𝑛=π‘§π‘›βˆ’π‘“ξ€·π‘₯𝑛𝑓𝑦𝑛𝑓𝑧𝑛π‘₯βˆ’π‘“π‘›ξ€Έξƒ¬1𝑓𝑦𝑛,π‘₯π‘›ξ€»βˆ’1𝑓𝑧𝑛,𝑦𝑛,π‘₯𝑛+1=π‘€π‘›βˆ’π‘“ξ€·π‘₯𝑛𝑓𝑦𝑛𝑓𝑧𝑛𝑓𝑀𝑛π‘₯βˆ’π‘“π‘›ξ€ΈΓ—ξƒ¬1π‘“ξ€·π‘€π‘›ξ€Έξ€·π‘¦βˆ’π‘“π‘›ξ€Έξƒ―1𝑓𝑀𝑛,π‘§π‘›ξ€»βˆ’1𝑓𝑧𝑛,π‘¦π‘›ξ€»ξƒ°βˆ’1𝑓𝑧𝑛π‘₯βˆ’π‘“π‘›ξ€Έξƒ―1𝑓𝑧𝑛,π‘¦π‘›ξ€»βˆ’1𝑓𝑦𝑛,π‘₯𝑛,ξƒ°ξƒ­(3.2) where π›½βˆˆβ„β§΅{0}. This family of one-parameter methods possesses the eighth-order convergence utilizing four pieces of information, namely, 𝑓(π‘₯𝑛), 𝑓(𝑦𝑛), 𝑓(𝑧𝑛) and 𝑓(𝑀𝑛). Therefore, its classical efficiency index is 1.682. Note that the first two steps of (3.2), are an optimal fourth-order derivative-free uniparametric family, which could also be given as comes next 𝑦𝑛=π‘₯π‘›βˆ’π‘“ξ€·π‘₯𝑛𝑓π‘₯𝑛,𝑀𝑛,𝑀𝑛=π‘₯𝑛π‘₯+𝛽𝑓𝑛π‘₯,π›½βˆˆβ„β§΅{0},𝑛+1=π‘¦π‘›βˆ’π‘“ξ€·π‘¦π‘›ξ€Έπ‘“ξ€·π‘€π‘›ξ€Έξ€·π‘“ξ€·π‘€π‘›ξ€Έξ€·π‘¦βˆ’π‘“π‘›π‘“ξ€Ίπ‘₯𝑛,𝑦𝑛.(3.3)

Another uniparametric three-step derivative-free iteration was presented recently by Zheng et al. in [21] as follows: 𝑦𝑛=π‘₯π‘›βˆ’π‘“ξ€·π‘₯𝑛𝑓π‘₯𝑛,𝑀𝑛,𝑀𝑛=π‘₯𝑛π‘₯+𝛽𝑓𝑛,𝑧𝑛=π‘¦π‘›βˆ’π‘“ξ€·π‘¦π‘›ξ€Έπ‘“ξ€Ίπ‘₯𝑛,𝑦𝑛𝑦+𝑓𝑛,𝑀𝑛π‘₯βˆ’π‘“π‘›,𝑀𝑛,π‘₯𝑛+1=π‘§π‘›βˆ’π‘“ξ€·π‘§π‘›ξ€Έπ‘“ξ€Ίπ‘§π‘›,𝑦𝑛𝑧+𝑓𝑛,𝑦𝑛,π‘₯π‘›π‘§ξ€»ξ€·π‘›βˆ’π‘¦π‘›ξ€Έξ€Ίπ‘§+𝑓𝑛,𝑦𝑛,π‘₯𝑛,π‘€π‘›π‘§ξ€»ξ€·π‘›βˆ’π‘¦π‘›π‘§ξ€Έξ€·π‘›βˆ’π‘₯𝑛,(3.4) where 𝑓[π‘₯𝑛,π‘₯π‘›βˆ’1,…,π‘₯π‘›βˆ’π‘–] is the divided differences of 𝑓(π‘₯). And we recall that they can be defined recursively via 𝑓π‘₯𝑖π‘₯=𝑓𝑖π‘₯;𝑓𝑖,π‘₯𝑗=𝑓π‘₯𝑖π‘₯βˆ’π‘“π‘—ξ€»π‘₯π‘–βˆ’π‘₯𝑗,π‘₯𝑖≠π‘₯𝑗(3.5) and, for π‘š>𝑖+1, via 𝑓π‘₯𝑖,π‘₯𝑖+1,…,π‘₯π‘šξ€»=𝑓π‘₯𝑖,π‘₯𝑖+1,…,π‘₯π‘šβˆ’1ξ€»ξ€Ίπ‘₯βˆ’π‘“π‘–+1,π‘₯𝑖+2,…,π‘₯π‘šξ€»π‘₯π‘–βˆ’π‘₯π‘š,π‘₯𝑖≠π‘₯π‘š.(3.6)

4. Main Contributions

As usual to build high-order iterations, we must consider a multipoint cycle. Now in order to contribute, we pay heed to the follow up three-step scheme in which the first step is Steffensen, while the second and the third steps are Newton's iterations. This procedure is completely inefficient. Since it possesses 81/4β‰ˆ1.414 as its efficiency index, which is the same to Steffensen's and Newton's: 𝑦𝑛=π‘₯π‘›βˆ’π‘“ξ€·π‘₯𝑛𝑓π‘₯𝑛,𝑀𝑛,𝑧𝑛=π‘¦π‘›βˆ’π‘“ξ€·π‘¦π‘›ξ€Έξ€·π‘¦π‘“β€²π‘›ξ€Έ,π‘₯𝑛+1=π‘§π‘›βˆ’π‘“ξ€·π‘§π‘›ξ€Έξ€·π‘§π‘“β€²π‘›ξ€Έ.(4.1)

To annihilate the derivative evaluations of the structure (4.1) and also keep the order at eight we must reconstruct our structure in which the first two steps provide the optimal order four using three evaluations and subsequently to the eighth-order convergence by consuming four function evaluations. Toward this end, we make use of weight function approach as comes next (and also by replacing π‘“ξ…ž(𝑦𝑛)β‰ˆπ‘“[π‘₯𝑛,𝑀𝑛] and π‘“ξ…ž(𝑧𝑛)β‰ˆπ‘“[π‘₯𝑛,𝑀𝑛]): 𝑦𝑛=π‘₯π‘›βˆ’π‘“ξ€·π‘₯𝑛𝑓π‘₯𝑛,𝑀𝑛,𝑀𝑛=π‘₯𝑛π‘₯+𝛽𝑓𝑛,𝑧𝑛=π‘¦π‘›βˆ’π‘“ξ€·π‘¦π‘›ξ€Έπ‘“ξ€Ίπ‘₯𝑛,𝑀𝑛π‘₯{𝐺(A)×𝐻(B)},𝑛+1=π‘§π‘›βˆ’π‘“ξ€·π‘§π‘›ξ€Έπ‘“ξ€Ίπ‘₯𝑛,𝑀𝑛{𝐾(Ξ“)×𝐿(Ξ”)×𝑃(E)×𝑄(B)×𝐽(A)},(4.2) wherein π›½βˆˆβ„β§΅{0},  𝐺(A),  𝐻(B),  𝐾(Ξ“),  𝐿(Ξ”),  𝑃(E),  𝑄(B), and 𝐽(A) are seven real-valued weight functions when A=(𝑓(𝑦)/𝑓(π‘₯)), B=(𝑓(𝑦)/𝑓(𝑀)), Ξ“=(𝑓(𝑧)/𝑓(π‘₯)), Ξ”=(𝑓(𝑧)/𝑓(𝑀)), and E=(𝑓(𝑧)/𝑓(𝑦)), without the index 𝑛, should be chosen such that the order of convergence arrives at the optimal level eight. Theorem 4.1 indicates the way of selecting the weight functions in order to reach the optimal efficiency index by using the smallest possible number of function evaluations.

Theorem 4.1. Let us consider 𝛼 as a simple root of the nonlinear equation 𝑓(π‘₯)=0 in the domain 𝐷. And assume that 𝑓(π‘₯) is sufficiently smooth in the neighborhood of the root. Then the derivative-free iterative class without memory defined by (4.2) is of optimal order eight, when the weight functions satisfy 𝐺(0)=πΊξ…ž(0)=1,πΊξ…žξ…ž(0)=𝐺(3)||𝐺(0)=0,(4)||(0)<∞,𝐻(0)=π»ξ…ž(0)=1,π»ξ…žξ…žπ»(0)=0,(3)ξ€Ίπ‘₯(0)=18+6𝛽𝑓𝑛,𝑀𝑛π‘₯ξ€»ξ€·5+𝛽𝑓𝑛,𝑀𝑛π‘₯ξ€»ξ€·4+𝛽𝑓𝑛,𝑀𝑛,||𝐻(4)||𝐾(0)<∞,(0)=πΎξ…ž(0)=1,𝐿(0)=πΏξ…ž(0)=1,𝑃(0)=π‘ƒξ…ž||𝑃(0)=1,ξ…žξ…ž||𝑄(0)<∞,(0)=π‘„ξ…ž(0)=1,π‘„ξ…žξ…žξ€Ίπ‘₯(0)=2+2𝛽𝑓𝑛,𝑀𝑛,𝑄(3)||𝑄(0)=0,(4)||(0)<∞,𝐽(0)=π½ξ…ž(0)=1,π½ξ…žξ…ž(0)=𝐽(3)||𝐽(0)=0,(4)||(0)<∞.(4.3)

Proof. Using Taylor's series and symbolic computations, we can determine the asymptotic error constant of the three-step uniparametric class (4.2)-(4.3). Furthermore, assume 𝑒𝑛=π‘₯π‘›βˆ’π›Ό be the error in the 𝑛th iterate and take into account 𝑓(𝛼)=0, π‘π‘˜=𝑓(π‘˜)(𝛼)/π‘˜!,  for all π‘˜=1,2,3,…. Now, we expand 𝑓(π‘₯𝑛) around the simple zero 𝛼. Hence, we have 𝑓π‘₯𝑛=𝑐1𝑒𝑛+𝑐2𝑒2𝑛+𝑐3𝑒3𝑛+𝑐4𝑒4𝑛+𝑐5𝑒5𝑛+𝑐5𝑒5𝑛+𝑐6𝑒6𝑛+𝑐7𝑒7𝑛+𝑐8𝑒8𝑛𝑒+𝑂9𝑛.(4.4) By considering (4.4) and the first step of (4.2), we attain π‘¦π‘›βˆ’π›Ό=𝑐2ξ‚΅1𝑐1𝑒+𝛽2𝑛𝑒+β‹―+𝑂9𝑛.(4.5) In the same vein, by considering (4.3) and (4.5), we obtain for the second step that π‘§π‘›π‘βˆ’π›Ό=2ξ€·1+𝑐1π›½ξ€Έξ€·βˆ’π‘1𝑐3ξ€·1+𝑐1𝛽+𝑐22(2+𝑐1𝛽)2𝑐31𝑒4𝑛𝑒+β‹―+𝑂9𝑛,(4.6) and 𝑓(𝑧𝑛)=((𝑐2(1+𝑐1𝛽)(βˆ’π‘1𝑐3(1+𝑐1𝛽)+𝑐22(2+𝑐1𝛽)2))/𝑐21)𝑒4𝑛+β‹―+𝑂(𝑒9𝑛). We also have (𝑓(𝑧𝑛)/𝑓[π‘₯𝑛,𝑀𝑛])=(1/𝑐31)𝑐2(1+𝑐1𝛽)(βˆ’π‘1𝑐3(1+𝑐1𝛽)+𝑐22(2+𝑐1𝛽)2)𝑒4𝑛+β‹―+𝑂(𝑒9𝑛). Using symbolic computation in the last step of (4.2) and considering (4.3) and (4.6), we have π‘§π‘›βˆ’π‘“ξ€·π‘§π‘›ξ€Έπ‘“ξ€Ίπ‘₯𝑛,π‘€π‘›ξ€»π‘βˆ’π›Ό=22ξ€·1+𝑐1𝛽2+𝑐1π›½ξ€Έξ‚€βˆ’π‘1𝑐3ξ€·1+𝑐1𝛽+𝑐22ξ€·2+𝑐1𝛽2𝑐41𝑒5𝑛𝑒+β‹―+𝑂9𝑛.(4.7) Furthermore using (4.3), we have 𝑓𝑧𝑛𝑓π‘₯𝑛,𝑀𝑛=𝑐{𝐾(Ξ“)×𝐿(Ξ”)×𝑃(E)×𝑄(B)×𝐽(A)}2ξ€·1+𝑐1π›½ξ€Έξ‚€βˆ’π‘1𝑐3ξ€·1+𝑐1𝛽+𝑐22ξ€·2+𝑐1𝛽2𝑒4𝑛𝑐31βˆ’1𝑐41𝑐21𝑐23ξ€·1+𝑐1𝛽2ξ€·2+𝑐1𝛽+𝑐21𝑐2𝑐4ξ€·1+𝑐1𝛽2ξ€·2+𝑐1π›½ξ€Έβˆ’π‘1𝑐22𝑐3ξ€·1+𝑐1𝛽26+𝑐1𝛽37+𝑐1𝛽19+3𝑐1𝛽+𝑐42ξ€·29+𝑐1𝛽65+𝑐1𝛽57+4𝑐1𝛽6+𝑐1𝛽𝑒5𝑛𝑒+β‹―+𝑂9𝑛.(4.8) And finally, Taylor series expansion around the simple root in the last iterate by using (4.3) and the above relation (values of higher order derivatives of 𝐾(Ξ“) and 𝐿(Ξ”); not explicitly given in (4.9), can be arbitrary at the point 0), will result in 𝑒𝑛+11=βˆ’24𝑐71𝑐2ξ€·1+𝑐1π›½ξ€Έξ‚€βˆ’π‘1𝑐3ξ€·1+𝑐1𝛽+𝑐22ξ€·2+𝑐1𝛽2ξ‚Γ—ξ‚€βˆ’24𝑐21𝑐2𝑐4ξ€·1+𝑐1𝛽2+12𝑐21𝑐23ξ€·1+𝑐1𝛽2ξ€·βˆ’2+π‘ƒξ…žξ…žξ€Έ(0)+24𝑐1𝑐22𝑐3ξ€·1+𝑐1𝛽2ξ€·6+𝑐1𝛽6+𝑐1π›½βˆ’ξ€·ξ€Έξ€Έ2+𝑐1𝛽2π‘ƒξ…žξ…žξ‚(0)+𝑐42ξ€·βˆ’504+192π‘ƒξ…žξ…ž(0)βˆ’πΊ(4)(0)βˆ’π»(4)(0)+𝐽(4)(0)+𝑐1𝛽4ξ€·βˆ’270+96π‘ƒξ…žξ…ž(0)βˆ’πΊ(4)(0)+𝐽(4)ξ€Έ(0)+𝑐1𝛽6ξ€·βˆ’144+48π‘ƒξ…žξ…ž(0)βˆ’πΊ(4)(0)+𝐽(4)ξ€Έ(0)+𝑐1𝛽𝑐1π›½ξ€·βˆ’48+12π‘ƒξ…žξ…ž(0)βˆ’πΊ(4)(0)+𝐽(4)ξ€Έξ€·(0)+4βˆ’78+24π‘ƒξ…žξ…ž(0)βˆ’πΊ(4)(0)+𝐽(4)(0)ξ€Έξ€Έξ€Έξ€Έ+𝑄(4)𝑒(0)ξ€Έξ€Έξ€Έ8𝑛𝑒+𝑂9𝑛,(4.9) which shows that the three-step without memory class (4.2)-(4.3) reaches the order of convergence eight by using only four pieces of information. This completes that proof.

A simple computational example from the class (4.2)-(4.3) can be 𝑦𝑛=π‘₯π‘›βˆ’π‘“ξ€·π‘₯𝑛𝑓π‘₯𝑛,𝑀𝑛,𝑀𝑛=π‘₯𝑛π‘₯+𝑓𝑛,𝑧𝑛=π‘¦π‘›βˆ’π‘“ξ€·π‘¦π‘›ξ€Έπ‘“ξ€Ίπ‘₯𝑛,π‘€π‘›ξ€»Γ—βŽ§βŽͺ⎨βŽͺβŽ©βŽ›βŽœβŽœβŽπ‘“ξ€·π‘¦1+𝑛𝑓π‘₯𝑛+𝑓𝑦𝑛𝑓π‘₯𝑛ξƒͺ5βŽžβŽŸβŽŸβŽ Γ—βŽ›βŽœβŽœβŽπ‘“ξ€·π‘¦1+𝑛𝑓𝑀𝑛+ξ€·ξ€Ίπ‘₯3+𝑓𝑛,𝑀𝑛π‘₯ξ€»ξ€·5+𝑓𝑛,𝑀𝑛π‘₯ξ€»ξ€·4+𝑓𝑛,𝑀𝑛𝑓𝑦𝑛𝑓𝑀𝑛ξƒͺ3⎞⎟⎟⎠⎫βŽͺ⎬βŽͺ⎭,π‘₯𝑛+1=π‘§π‘›βˆ’π‘“ξ€·π‘§π‘›ξ€Έπ‘“ξ€Ίπ‘₯𝑛,π‘€π‘›ξ€»ξ€½π‘Š1ξ€Ύ,(4.10) where π‘Š1=βŽ›βŽœβŽœβŽπ‘“ξ€·π‘§1+𝑛𝑓π‘₯𝑛+𝑓𝑧𝑛𝑓π‘₯𝑛ξƒͺ3βŽžβŽŸβŽŸβŽ βŽ›βŽœβŽœβŽπ‘“ξ€·π‘§1+𝑛𝑓𝑀𝑛+𝑓𝑧𝑛𝑓𝑀𝑛ξƒͺ3βŽžβŽŸβŽŸβŽ βŽ›βŽœβŽœβŽπ‘“ξ€·π‘§1+𝑛𝑓𝑦𝑛+𝑓𝑧𝑛𝑓𝑦𝑛ξƒͺ2βŽžβŽŸβŽŸβŽ Γ—ξƒ©π‘“ξ€·π‘¦1+𝑛𝑓𝑀𝑛+ξ€·ξ€Ίπ‘₯1+𝑓𝑛,𝑀𝑛𝑓(𝑦𝑛)𝑓(𝑀𝑛)ξ‚Ά2𝑓𝑦ξƒͺ1+𝑛𝑓π‘₯𝑛+𝑓(𝑦𝑛)𝑓(π‘₯𝑛)ξ‚Ά5ξƒͺ,(4.11) and its error equation is as comes next 𝑒𝑛+1=ξ€·1+𝑐1ξ€Έ3𝑐22ξ‚€ξ€·2+𝑐1ξ€Έ2𝑐22βˆ’π‘1ξ€·1+𝑐1𝑐35+𝑐1ξ€·3+𝑐1𝑐32βˆ’4𝑐1𝑐2𝑐3+𝑐21𝑐4𝑐71𝑒8𝑛𝑒+𝑂9𝑛.(4.12)

Remark 4.2. Although the structure (4.10) is hard, it could be coded easily because 𝑓[π‘₯𝑛,𝑀𝑛] and the factors 𝑓(𝑦𝑛)/𝑓(π‘₯𝑛), 𝑓(𝑦𝑛)/𝑓(𝑀𝑛), 𝑓(𝑧𝑛)/𝑓(π‘₯𝑛), 𝑓(𝑧𝑛)/𝑓(𝑀𝑛), and 𝑓(𝑧𝑛)/𝑓(𝑦𝑛) should be computed only once per computing step and their values will be used throughout and we thus have some operational calculations in implementing (4.10).

The contributed class (4.2)-(4.3) uses the forward finite difference approximation in its first step. If we use backward finite difference estimation, then another novel three-step eighth-order derivative-free iteration without memory can be constructed. For this reason, our second general derivative-free class could be given as 𝑦𝑛=π‘₯π‘›βˆ’π‘“ξ€·π‘₯𝑛𝑓π‘₯𝑛,𝑀𝑛,𝑀𝑛=π‘₯𝑛π‘₯βˆ’π›½π‘“π‘›ξ€Έ,𝑧𝑛=π‘¦π‘›βˆ’π‘“ξ€·π‘¦π‘›ξ€Έπ‘“ξ€Ίπ‘₯𝑛,𝑀𝑛π‘₯{𝐴(𝑑)×𝐡(𝑖)},𝑛+1=π‘§π‘›βˆ’π‘“ξ€·π‘§π‘›ξ€Έπ‘“ξ€Ίπ‘₯𝑛,𝑀𝑛{𝑃(π‘Ÿ)+𝑄(𝑖)+𝐽(𝑑)+𝐿(π‘˜)},(4.13) where π›½βˆˆβ„β§΅{0},𝐴(𝑑),𝐡(𝑖),𝑃(π‘Ÿ),𝑄(𝑖),𝐽(𝑑), and 𝐿(π‘˜)are six real-valued weight functions when 𝑑=𝑓(𝑦)/𝑓(π‘₯), 𝑖=𝑓(𝑦)/𝑓(𝑀), π‘Ÿ=𝑓(𝑧)/𝑓(𝑦), π‘˜=𝑓(𝑧)/𝑓(𝑀), without the index 𝑛, and they must read 𝐴(0)=π΄ξ…ž(0)=1,π΄ξ…žξ…ž(0)=2,𝐴(3)||𝐴(0)=0,(4)||(0)<∞,𝐡(0)=π΅ξ…ž(0)=1,π΅ξ…žξ…žξ€Ίπ‘₯(0)=6βˆ’4𝛽𝑓𝑛,𝑀𝑛,𝐡(3)||𝐡(0)=0,(4)||||𝑃(0)<∞,𝑃(0)=𝑃′(0)=1,ξ…žξ…ž||(0)<∞,𝑄(0)=0,π‘„ξ…ž(0)=1,π‘„ξ…žξ…žξ€Ίπ‘₯(0)=10βˆ’8𝛽𝑓𝑛,𝑀𝑛,𝑄(3)ξ€Ίπ‘₯(0)=60+12𝛽𝑓𝑛,𝑀𝑛π‘₯ξ€»ξ€·βˆ’8+3𝛽𝑓𝑛,𝑀𝑛,||𝑄(4)||(0)<∞,𝐽(0)=0,π½ξ…ž(0)=1,π½ξ…žξ…ž(0)=2,𝐽(3)(||𝐽0)=0,(4)(||0)<∞,𝐿(0)=0,πΏξ…žξ€Ίπ‘₯(0)=4βˆ’2𝛽𝑓𝑛,𝑀𝑛.(4.14)

The scheme (4.13)-(4.14) defines a new family of multipoint methods. To obtain the solution of (1.1) by the new derivative-free class, we must set a particular guess π‘₯0, ideally close to the simple root. In numerical analysis, it is very useful and essential to know the behavior of an approximate method. Therefore, we will prove the order of convergence of the new eighth-order class.

Theorem 4.3. Let us consider 𝛼 as a simple root of the nonlinear equation 𝑓(π‘₯)=0 in the domain 𝐷. And assume that 𝑓(π‘₯) is sufficiently smooth in the neighborhood of the root. Then the derivative-free iterative class defined by (4.13)-(4.14) is of optimal order eight.

Proof. Applying Taylor series and symbolic computations, we can determine the asymptotic error constant of the three-step uni-parametric family (4.13)-(4.14). Assume 𝑒𝑛=π‘₯π‘›βˆ’π›Ό be the error in the nth iterate and take into account 𝑓(𝛼)=0, π‘π‘˜=𝑓(π‘˜)(𝛼)/π‘˜!,  for all π‘˜=1,2,3,…. Then, the procedure of this proof is similar to the proof of Theorem 4.1. Thus, we below give the final error equation of (4.13)-(4.14): 𝑒𝑛+11=βˆ’24𝑐61𝑐2𝑐3ξ€·βˆ’1+𝑐1𝛽2Γ—ξ‚€24𝑐21𝑐2𝑐4ξ€·βˆ’1+𝑐1𝛽2+48𝑐1𝑐22𝑐3ξ€·βˆ’1+𝑐1𝛽7+𝑐1π›½ξ€·βˆ’7+𝑐1π›½ξ€Έξ€Έβˆ’12𝑐21𝑐23ξ€·βˆ’1+𝑐1𝛽2ξ€·βˆ’2+π‘ƒξ…žξ…žξ€Έ(0)+𝑐42ξ€·264+𝐴(4)(0)+𝐡(4)(0)βˆ’π½(4)(0)+𝑐1π›½ξ€·ξ€·βˆ’4156+𝐴(4)(0)βˆ’π½(4)(ξ€Έ0)+𝑐1𝛽6ξ€·60+𝐴(4)(0)βˆ’π½(4)ξ€Έ(0)+𝑐1π›½ξ€·βˆ’4𝐴(4)(0)+𝑐1π›½ξ€·βˆ’24+𝐴(4)(0)βˆ’π½(4)ξ€Έ(0)+4𝐽(4)(0)ξ€Έξ€Έξ€Έβˆ’π‘„(4)𝑒(0)ξ€Έξ€Έξ€Έ8𝑛𝑒+𝑂9𝑛.(4.15) This shows that the contributed class (4.13)-(4.14) achieves the optimal order eight by using only four pieces of information. The proof is complete.

A very efficient example from our novel class (4.13)-(4.14) can be the following iteration without memory 𝑦𝑛=π‘₯π‘›βˆ’π‘“ξ€·π‘₯𝑛𝑓π‘₯𝑛,𝑀𝑛,𝑀𝑛=π‘₯𝑛π‘₯βˆ’π‘“π‘›ξ€Έ,𝑧𝑛=π‘¦π‘›βˆ’π‘“ξ€·π‘¦π‘›ξ€Έπ‘“ξ€Ίπ‘₯𝑛,π‘€π‘›ξ€»Γ—βŽ§βŽͺ⎨βŽͺβŽ©βŽ›βŽœβŽœβŽπ‘“ξ€·π‘¦1+𝑛𝑓π‘₯𝑛+𝑓𝑦𝑛𝑓π‘₯𝑛ξƒͺ2βŽžβŽŸβŽŸβŽ βŽ›βŽœβŽœβŽπ‘“ξ€·π‘¦1+𝑛𝑓𝑀𝑛+ξ€·ξ€Ίπ‘₯3βˆ’2𝑓𝑛,𝑀𝑛𝑓𝑦𝑛𝑓𝑀𝑛ξƒͺ2⎞⎟⎟⎠⎫βŽͺ⎬βŽͺ⎭,π‘₯𝑛+1=π‘§π‘›βˆ’π‘“ξ€·π‘§π‘›ξ€Έπ‘“ξ€Ίπ‘₯𝑛,π‘€π‘›ξ€»ξ€½π‘Š2ξ€Ύ,(4.16) where π‘Š2𝑓𝑧=1+𝑛𝑓𝑦𝑛+𝑓𝑧𝑛𝑓𝑦𝑛ξƒͺ2+𝑓𝑦𝑛𝑓𝑀𝑛+ξ€·ξ€Ίπ‘₯5βˆ’4𝑓𝑛,𝑀𝑛𝑓𝑦𝑛𝑓𝑀𝑛ξƒͺ2+ξ€·ξ€Ίπ‘₯10+2𝑓𝑛,𝑀𝑛π‘₯ξ€»ξ€·βˆ’8+3𝑓𝑛,𝑀𝑛𝑓𝑦𝑛𝑓𝑀𝑛ξƒͺ3+ξ‚€ξ€Ίπ‘₯11βˆ’π‘“π‘›,𝑀𝑛π‘₯26βˆ’15𝑓𝑛,𝑀𝑛+𝑓π‘₯𝑛,𝑀𝑛3𝑓𝑦𝑛𝑓𝑀𝑛ξƒͺ4+𝑓𝑦𝑛𝑓π‘₯𝑛+𝑓𝑦𝑛𝑓π‘₯𝑛ξƒͺ2+ξ€·ξ€Ίπ‘₯4βˆ’2𝑓𝑛,𝑀𝑛𝑓𝑧𝑛𝑓𝑀𝑛,(4.17) and its very simple error equation comes next as 𝑒𝑛+1=βˆ’(βˆ’1+𝑐1)3𝑐22𝑐3ξ€·2ξ€·ξ€·7+βˆ’7+𝑐1𝑐1𝑐2𝑐3+ξ€·βˆ’1+𝑐1𝑐1𝑐4𝑐51𝑒8𝑛𝑒+𝑂9𝑛.(4.18) Mostly, and according to the assumptions made in the proof of Theorem 4.1, 𝑐71 should be in the denominator of the asymptotic error constant of the optimal eighth-order methods, and if one obtains some forms like 𝑐61 or 𝑐51, then the derived methods will be mostly finer than the other existing forms of optimal eighth-order methods.

Remark 4.4. Each method from the proposed derivative-free classes in this paper reaches the efficiency index 4√8β‰ˆ1.682, which is greater than 3√4β‰ˆ1.587 of the optimal fourth-order techniques' and 2√2β‰ˆ1.414 of optimal one-point methods' without memory.

Some other examples from the class (4.13)-(4.14) can be easily constructed by changing the involved weight function in (4.14). For example, we can have 𝑦𝑛=π‘₯π‘›βˆ’π‘“ξ€·π‘₯𝑛𝑓π‘₯𝑛,𝑀𝑛,𝑀𝑛=π‘₯𝑛π‘₯βˆ’π‘“π‘›ξ€Έ,𝑧𝑛=π‘¦π‘›βˆ’π‘“ξ€·π‘¦π‘›ξ€Έπ‘“ξ€Ίπ‘₯𝑛,π‘€π‘›ξ€»Γ—βŽ§βŽͺ⎨βŽͺβŽ©βŽ›βŽœβŽœβŽπ‘“ξ€·π‘¦1+𝑛𝑓π‘₯𝑛+𝑓𝑦𝑛𝑓π‘₯𝑛ξƒͺ2βŽžβŽŸβŽŸβŽ βŽ›βŽœβŽœβŽπ‘“ξ€·π‘¦1+𝑛𝑓𝑀𝑛+ξ€·ξ€Ίπ‘₯3βˆ’2𝑓𝑛,𝑀𝑛𝑓𝑦𝑛𝑓𝑀𝑛ξƒͺ2⎞⎟⎟⎠⎫βŽͺ⎬βŽͺ⎭,π‘₯𝑛+1=π‘§π‘›βˆ’π‘“ξ€·π‘§π‘›ξ€Έπ‘“ξ€Ίπ‘₯𝑛,π‘€π‘›ξ€»ξ€½π‘Š3ξ€Ύ,(4.19) with π‘Š3𝑓𝑧=1+𝑛𝑓𝑦𝑛+𝑓𝑦𝑛𝑓𝑀𝑛+ξ€·ξ€Ίπ‘₯5βˆ’4𝑓𝑛,𝑀𝑛𝑓𝑦𝑛𝑓𝑀𝑛ξƒͺ2+ξ€·ξ€Ίπ‘₯10+2𝑓𝑛,𝑀𝑛π‘₯ξ€»ξ€·βˆ’8+3𝑓𝑛,𝑀𝑛𝑓𝑦𝑛𝑓𝑀𝑛ξƒͺ3+ξ‚€ξ€Ίπ‘₯11βˆ’π‘“π‘›,𝑀𝑛π‘₯26βˆ’15𝑓𝑛,𝑀𝑛+𝑓π‘₯𝑛,𝑀𝑛3×𝑓𝑦𝑛𝑓𝑀𝑛ξƒͺ4+𝑓𝑦𝑛𝑓π‘₯𝑛+𝑓𝑦𝑛𝑓π‘₯𝑛ξƒͺ2+ξ€·ξ€Ίπ‘₯4βˆ’2𝑓𝑛,𝑀𝑛𝑓𝑧𝑛𝑓𝑀𝑛,(4.20) where its error equation comes next as 𝑒𝑛+1=βˆ’(βˆ’1+𝑐1)3𝑐2𝑐3ξ€·2ξ€·ξ€·7+βˆ’7+𝑐1𝑐1𝑐22𝑐3+ξ€·βˆ’1+𝑐1𝑐1𝑐23+ξ€·βˆ’1+𝑐1𝑐1𝑐2𝑐4𝑐51𝑒8𝑛𝑒+𝑂9𝑛.(4.21) And also we could construct 𝑦𝑛=π‘₯π‘›βˆ’π‘“ξ€·π‘₯𝑛𝑓π‘₯𝑛,𝑀𝑛,𝑀𝑛=π‘₯𝑛π‘₯βˆ’π‘“π‘›ξ€Έ,𝑧𝑛=π‘¦π‘›βˆ’π‘“ξ€·π‘¦π‘›ξ€Έπ‘“ξ€Ίπ‘₯𝑛,π‘€π‘›ξ€»Γ—βŽ§βŽͺ⎨βŽͺβŽ©βŽ›βŽœβŽœβŽπ‘“ξ€·π‘¦1+𝑛𝑓π‘₯𝑛+𝑓𝑦𝑛𝑓π‘₯𝑛ξƒͺ2βŽžβŽŸβŽŸβŽ βŽ›βŽœβŽœβŽπ‘“ξ€·π‘¦1+𝑛𝑓𝑀𝑛+ξ€·ξ€Ίπ‘₯3βˆ’2𝑓𝑛,𝑀𝑛𝑓𝑦𝑛𝑓𝑀𝑛ξƒͺ2⎞⎟⎟⎠⎫βŽͺ⎬βŽͺ⎭,π‘₯𝑛+1=π‘§π‘›βˆ’π‘“ξ€·π‘§π‘›ξ€Έπ‘“ξ€Ίπ‘₯𝑛,π‘€π‘›ξ€»ξ€½π‘Š4ξ€Ύ,(4.22) where π‘Š4𝑓𝑧=1+𝑛𝑓𝑦𝑛+𝑓𝑧𝑛𝑓𝑦𝑛ξƒͺ2+𝑓𝑦𝑛𝑓𝑀𝑛+ξ€·ξ€Ίπ‘₯5βˆ’4𝑓𝑛,𝑀𝑛𝑓𝑦𝑛𝑓𝑀𝑛ξƒͺ2+ξ€·ξ€Ίπ‘₯10+2𝑓𝑛,𝑀𝑛π‘₯ξ€»ξ€·βˆ’8+3𝑓𝑛,𝑀𝑛𝑓𝑦𝑛𝑓𝑀𝑛ξƒͺ3+ξ‚€ξ€Ίπ‘₯11βˆ’π‘“π‘›,𝑀𝑛π‘₯26βˆ’15𝑓𝑛,𝑀𝑛+𝑓π‘₯𝑛,𝑀𝑛3𝑓𝑦𝑛𝑓𝑀𝑛ξƒͺ4+𝑓𝑦𝑛𝑓π‘₯𝑛+𝑓𝑦𝑛𝑓π‘₯𝑛ξƒͺ2+𝑓𝑦𝑛𝑓π‘₯𝑛ξƒͺ4+ξ€·ξ€Ίπ‘₯4βˆ’2𝑓𝑛,𝑀𝑛𝑓𝑧𝑛𝑓𝑀𝑛,(4.23) and its error equation comes next as 𝑒𝑛+1=(βˆ’1+𝑐1)3𝑐22𝑐3ξ€·(βˆ’1+𝑐1)3𝑐32βˆ’2𝑐1ξ€·ξ€·7+βˆ’7+𝑐1𝑐1𝑐2𝑐3βˆ’ξ€·βˆ’1+𝑐1𝑐21𝑐4𝑐61𝑒8𝑛𝑒+𝑂9𝑛.(4.24)

5. Numerical Examples

We check the effectiveness of the novel derivative-free classes of iterative methods (4.2)-(4.3) and (4.13)-(4.14) in this section. In order to do this, we choose (4.10) as the representative from the optimal class (4.2)-(4.3) and (4.16) as the representative of the novel three-step class (4.13)-(4.14). We have compared (4.10) and (4.16) with Steffensen's method (3.1), the fourth-order family of Kung and Traub (3.3) with 𝛽=0.01, the eighth-order technique of Kung and Traub (3.2) with 𝛽=1, and the optimal eighth-order family of Zheng et al. (3.4) with 𝛽=1, using the examples listed in Table 1.

Table 1: The examples considered in this study.

The results of comparisons are given in Table 2 in terms of the number significant digits for each test function after the specified number of iterations, that is, for example, 0.1π‘’βˆ’1190 shows that the absolute value of the given nonlinear function 𝑓1, after three iterations, is zero up to 1190 decimal places. For numerical comparisons, the stopping criterion is |𝑓(π‘₯𝑛)|<1.πΈβˆ’1200. In Table 2, IT and TNE stand for the number of iterations and the total number of (function) evaluations. We used four different initial guesses to analyze the behaviors of the methods totally. In Table 2, F stands for failure, for instance, when the iteration for the particular initial guess becomes divergence, or finds another root, needs more numerical computing steps to find an acceptable solution of the nonlinear equations.

Table 2: Results of convergence under fair circumstances for different derivative-free methods.

It can be observed from Table 2 that almost in most cases our contributed methods from the classes of derivative free without memory iterations are superior in solving nonlinear equations. Numerical computations have been carried out using variable precision arithmetic in MATLAB 7.6.

We here remark that the eighth-order iterative methods such as (4.10) and (4.16) improve the number of correct digits in the convergence phase for the simple roots of (1.1) by a factor of eight and in order to show this and also the asymptotic error constant, we have applied 1200 digits floating point.

Note that experimental results show that whatever the value of 𝛽≠0 is small, then the output results of solving nonlinear equations will be more reliable. As a matter of fact, experimental results for the contributed methods from the three-step derivative-free classes (4.2)-(4.3) and (4.13)-(4.14) can give better feedbacks in most cases, for instance they provide better accuracy than those illustrated in Table 2, by choosing very small values for 𝛽. By doing this for 𝛽, the error equation will be narrowed. Also note that, if we approximate 𝛽 by an iteration through the data of the first step per cycle, then with memory iterations from the suggested classes will be attained, which is the topic of the forthcoming papers in this field.

A simple glance at Table 2 reveals that (4.16) is mostly better than its competitors. The reason indeed is that whatever the error equation is finer, the better numerical results will be attained. The error equations correspond to (4.16), for instance (4.18) is very small, in fact in its denominator we have 𝑐51 which clearly shows these refinements. In general, and according to the assumptions made in the proof of Theorem 4.1, 𝑐71 should be in the denominator of the optimal eighth-order methods, and if one obtains forms like 𝑐61 or 𝑐51, then the derived methods will be mostly finer than the other existing forms of the optimal eighth-order methods.

6. Conclusions

Multipoint methods without memory are methods that use new information at a number of points. Much literature on the multipoint Newton-like methods for function of one variable and their convergence analysis can be found in Traub [2].

In this paper, two novel classes of iterations without memory were discussed fully. We have shown that each member of our contributions reach the optimal order of convergence eight by consuming only four function evaluations per full cycle. Thus, our classes support the optimality conjecture of Kung and Traub for building optimal without memory iterations. Our classes can be taken into account as the generalizations of the well-cited derivative-free method of Steffensen. We have given a lot of numerical examples in Section 5 to support the underlying theory developed in this paper. Numerical results were completely in harmony with the theory developed in this paper, and accordingly, the contributions hit the target.


The author cheerfully acknowledges the interesting comments of the reviewer, which have led to the improvement of this paper.


  1. A. Iliev and N. Kyurkchiev, Nontrivial Methods in Numerical Analysis (Selected Topics in Numerical Analysis), Lambert Academic Publishing, 2010.
  2. J. F. Traub, Iterative Methods for the Solution of Equations, Prentice-Hall, Englewood Cliffs, NJ, USA, 1964.
  3. F. Soleymani, R. Sharma, X. Li, and E. Tohidi, β€œAn optimized derivative-free form of the Potra-Ptak method,” Mathematical and Computer Modelling, vol. 56, pp. 97–104, 2012. View at Publisher Β· View at Google Scholar
  4. F. Soleymani, β€œOptimal fourth-order iterative methods free from derivative,” Miskolc Mathematical Notes, vol. 12, pp. 255–264, 2011. View at Google Scholar
  5. F. Soleymani and B. S. Mousavi, β€œOn novel classes of iterative methods for solving nonlinear equations,” Computational Mathematics and Mathematical Physics, vol. 52, pp. 203–210, 2012. View at Publisher Β· View at Google Scholar
  6. F. Soleymani and S. K. Khattri, β€œFinding simple roots by seventh- and eighth-order derivative-free methods,” International Journal of Mathematical Models and Methods in Applied Sciences, vol. 6, pp. 45–52, 2012. View at Google Scholar
  7. F. Soleymani, β€œOn a novel optimal quartically class of methods,” Far East Journal of Mathematical Sciences, vol. 58, pp. 199–206, 2011. View at Google Scholar
  8. F. Soleymani, β€œOptimal eighth-order simple root-finders free from derivative,” WSEAS Transactions on Information Science and Applications, vol. 8, pp. 293–299, 2011. View at Google Scholar
  9. F. Soleymani and F. Soleimani, β€œNovel computational derivative-free methods for simple roots,” Fixed Point Theory, Computation and Applications, vol. 13, pp. 247–258, 2012. View at Google Scholar
  10. F. Soleymani, S. K. Vanani, and A. Afghani, β€œA general three-step class of optimal iterations for nonlinear equations,” Mathematical Problems in Engineering, vol. 2011, Article ID 469512, 10 pages, 2011. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  11. W.-X. Wang, Y.-L. Shang, W.-G. Sun, and Y. Zhang, β€œFinding the roots of system of nonlinear equations by a novel filled function method,” Abstract and Applied Analysis, vol. 2011, Article ID 209083, 9 pages, 2011. View at Publisher Β· View at Google Scholar
  12. B. Ignatova, N. Kyurkchiev, and A. Iliev, β€œMultipoint algorithms arising from optimal in the sense of Kung-Traub iterative procedures for numerical solution of nonlinear equations,” General Mathematics Notes, vol. 6, no. 2, 2011. View at Google Scholar
  13. S. Kumar, V. Kanwar, S. K. Tomar, and S. Singh, β€œGeometrically constructed families of Newton's method for unconstrained optimization and nonlinear equations,” International Journal of Mathematics and Mathematical Sciences, vol. 2011, Article ID 972537, 9 pages, 2011. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  14. H. T. Kung and J. F. Traub, β€œOptimal order of one-point and multipoint iteration,” Journal of the Association for Computing Machinery, vol. 21, pp. 643–651, 1974. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  15. P. Sargolzaei and F. Soleymani, β€œAccurate fourteenth-order methods for solving nonlinear equations,” Numerical Algorithms, vol. 58, no. 4, pp. 513–527, 2011. View at Publisher Β· View at Google Scholar
  16. M. S. Petković, J. Džunić, and B. Neta, β€œInterpolatory multipoint methods with memory for solving nonlinear equations,” Applied Mathematics and Computation, vol. 218, no. 6, pp. 2533–2541, 2011. View at Publisher Β· View at Google Scholar
  17. F. Soleymani, S. Karimi Vanani, and M. Jamali Paghaleh, β€œA class of three-step derivative-free root-solvers with optimal convergence order,” Journal of Applied Mathematics, vol. 2012, Article ID 568740, 15 pages, 2012. View at Publisher Β· View at Google Scholar
  18. M. Sharifi, D. K. R. Babajee, and F. Soleymani, β€œFinding the solution of nonlinear equations by a class of optimal methods,” Computers & Mathematics with Applications, vol. 63, pp. 764–774, 2012. View at Google Scholar
  19. F. Soleymani, S. Karimi Vanani, M. Khan, and M. Sharifi, β€œSome modifications of King's family with optimal eighth order of convergence,” Mathematical and Computer Modelling, vol. 55, pp. 1373–1380, 2012. View at Publisher Β· View at Google Scholar
  20. F. Soleymani, M. Sharifi, and B. Somayeh Mousavi, β€œAn improvement of Ostrowski's and King's techniques with optimal convergence order eight,” Journal of Optimization Theory and Applications, vol. 153, no. 1, pp. 225–236, 2012. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  21. Q. Zheng, J. Li, and F. Huang, β€œAn optimal Steffensen-type family for solving nonlinear equations,” Applied Mathematics and Computation, vol. 217, no. 23, pp. 9592–9597, 2011. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  22. F. Soleymani, β€œOn a new class of optimal eighth-order derivative-free methods,” Mathematica, vol. 3, pp. 169–180, 2011. View at Google Scholar
  23. F. Soleymani and S. Karimi Vanani, β€œOptimal Steffensen-type methods with eighth order of convergence,” Computers & Mathematics with Applications, vol. 62, no. 12, pp. 4619–4626, 2011. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  24. F. Soleymani, β€œNew optimal iterative methods in solving nonlinear equations,” International Journal of Pure and Applied Mathematics, vol. 72, pp. 195–202, 2011. View at Google Scholar
  25. F. Soleymani, β€œOn a bi-parametric class of optimal eighth-order derivative-free methods,” International Journal of Pure and Applied Mathematics, vol. 72, no. 1, pp. 27–37, 2011. View at Google Scholar
  26. F. Soleymani, S. K. Khattri, and S. Karimi Vanani, β€œTwo new classes of optimal Jarratt-type fourth-order methods,” Applied Mathematics Letters, vol. 25, pp. 847–853, 2012. View at Publisher Β· View at Google Scholar
  27. F. Soleymani, β€œNovel computational iterative methods with optimal order for nonlinear equations,” Advances in Numerical Analysis, vol. 2011, Article ID 270903, 10 pages, 2011. View at Google Scholar Β· View at Zentralblatt MATH
  28. B. I. Yun, β€œA non-iterative method for solving non-linear equations,” Applied Mathematics and Computation, vol. 198, no. 2, pp. 691–699, 2008. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  29. D. K. R. Babajee, , Analysis of higher order variants of Newton's method and their applications to differential and integral equations and in ocean acidification [Ph.D. thesis], University of Mauritius, 2010.
  30. J. F. Steffensen, β€œRemarks on iteration,” Scandinavian Actuarial Journal, vol. 16, no. 1, pp. 64–72, 1933. View at Publisher Β· View at Google Scholar