Mathematical Problems in EngineeringVolume 2011 (2011), Article ID 469512, 10 pagesdoi:10.1155/2011/469512

Research Article

A General Three-Step Class of Optimal Iterations for Nonlinear Equations

Fazlollah Soleymani,¹ Solat Karimi Vanani,² and Abtin Afghani²

¹Young Researchers Club, Islamic Azad University, Zahedan Branch, Zahedan 98168, Iran²Department of Mathematics, Islamic Azad University, Zahedan Branch, Zahedan 98168, Iran

Received 7 August 2011; Accepted 16 September 2011

Academic Editor: Hung Nguyen-Xuan

Copyright © 2011 Fazlollah Soleymani et al. 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.

Abstract.

Many of the engineering problems are reduced to solve a nonlinear equation numerically, and as a result, an especial attention to suggest efficient and accurate root solvers is given in literature. Inspired and motivated by the research going on in this area, this paper establishes an efficient general class of root solvers, where per computing step, three evaluations of the function and one evaluation of the first-order derivative are used to achieve the optimal order of convergence eight. The without-memory methods from the developed class possess the optimal efficiency index 1.682. In order to show the applicability and validity of the class, some numerical examples are discussed.

1. Introduction

Numerical solution of nonlinear scalar equations plays a crucial role in many optimization and engineering problems. For example, many engineering systems can be modeled as neutral delay differential equations (NDDEs) that involve a time delay in the derivative of the highest order, which are different from retarded delay differential equations (RDDEs) that do not involve a time delay in the derivative of the highest order. To illustrate more, a system, which consists of a mass mounted on a linear spring to which a pendulum is attached via a hinged massless rod, is used to predict the dynamic response of structures to external forces using a set of actuators, and it is modeled as an NDDE if the delay in actuators is taken into consideration [1]. On the other hand, the stability of a delay differential equation can be investigated on the basis of the root location of the characteristic function. This simple example shows the importance of numerical root solvers in engineering problems.

There are numerical methods, which find one root at a time, such as Newton’s iteration or its variant, and the schemes, which find all the roots at a time, namely, simultaneous methods, such as Weierstrass method. Recently many journals such as Numerical Algorithms, Mathematical Problems in Engineering, Applied Mathematics and Computation, etc., have published new findings; see, for example, [2–5] and the references therein in this active topic of study. To shortly provide some of the newest findings in this field, we mention the following.

Noor et al. in [3] developed the follow-up quartically iterative scheme consisting of three steps and eight numbers of evaluation per full iteration as comes next

In 2010, an eighth-order method is provided in [6] using Ostrowski's method in the first two steps of a three-step cycle as follows: wherein are the divided differences of the function .

Soleymani and Mousavi in [7] suggested an iteration without memory scheme including three steps and only four functional evaluations per iteration in what follows: where denotes , and denotes .

For further reading, one may consult [8], where a complete review of the methods given in literature from 2000 to 2010 was furnished, and also [9] for obtaining a background on the application use of such root solvers. We here remark that the efficiency of different methods could be assessed by the measure of efficiency index, which could be defined as , wherein is the order of convergence and is the total number of evaluations per iteration. Now, we should remark that Kung and Traub in [10] conjectured that an iterative scheme without memory by using evaluations per cycle can arrive at the maximum order of convergence . Any without memory iteration, which satisfies this hypothesis, is named as an optimal method in literature.

After providing a short background of this research in this section, we give the main contribution in Section 2. The convergence study of our general three-step class is also furnished therein. We will also produce different optimal three-step iterations from the contributed class. Section 3 discusses some numerical comparisons with the existing methods in literature, and finally Section 4 draws a conclusion of this research paper.

2. New Class of Iteration Methods

In order to contribute and give a general class of methods consistent with the optimality conjecture of Kung-Traub, an iteration eighth-order scheme without memory in this section should be constructed such that four evaluations per computing step are used. Such schemes are also known as predictor-corrector methods in which the first step is (Newton's step) predictor, while the other two steps correct the obtained solution. To achieve our goal, we consider the following three-step scheme on which the first two steps are the King's fourth-order family with one free parameter in real numbers, ,

Clearly in (2.1), should be annihilated as the order of convergence remains at the highest level by the smallest use of number of evaluations per iteration. Toward this end, we approximate it by a polynomial of degree two that fits in , , and . Therefore, we take into account where . Subsequently, by considering , , and , we attain ,

Solving the system of two linear equations with two unknowns, (2.2) gives us and . Using the obtained relations for the unknowns in the approximation and simplifying, we have

Considering (2.3) in (2.1) and using weight function approach, we have the following general class of three-step without-memory iteration: wherein , , and are three real-valued weight functions with , , and (without the index ) which should be chosen such that the order of convergence arrives at the optimal level eight. We summarize this in the following theorem.

Theorem 2.1. Let be a simple zero of a sufficiently differentiable function   for an open interval , which contains as the initial approximation of , then the three-step iteration (2.4), which includes four evaluations per full cycle, has the optimal convergence rate eight, when and satisfies the error equation below

Proof. By defining as the error of the iterative scheme in the th iterate, applying the Taylor's series expansion for (2.4), and taking into account , we have where . Furthermore, we have Dividing (2.7) by (2.8) gives us . Again by substituting this relation in the first step of (2.4) and writing the Taylor's series expansion for , we obtain, respectively, Furthermore, we find Similarly, we have We moreover by using (2.11) and (2.5) attain that . Considering this new relation, (2.12) and (2.5) in the last step of (2.4) will end in This concludes the proof. And it shows that our suggested general class of three-step without-memory methods (2.4)-(2.5) possesses the eighth order of convergence.

Remark 2.2. The class of three-step methods (2.4)-(2.5) requires four evaluations and has the order of convergence eight. Therefore, this class is of optimal order and supports the Kung-Traub conjecture [10]. Hence, the efficiency index of the eighth-order derivative-involved methods from the class is .Some efficient methods from the contributed optimal three-step class are given below. Per computing step, these methods are free from second or higher order derivative computations. The new contributed methods are where , where is its error equation. We also here mention some of the typical forms of the weight functions , , and in iteration (2.4), which satisfy (2.5) for making the order optimal. These forms are listed in Table 1.Other than the very efficient methods (2.14) and (2.15) of optimal order eight, many more three-step without-memory iterations can be constructed using Table 1, that is, (2.5) in (2.4) and also with different values for the free parameter . Thus, now, in order to save the space and also giving some of the other such optimal eighth-order methods according to (2.4) and (2.5), we list the interesting ones in Table 2. Note that we consider first that the weight functions satisfy (2.5), and then we try to make new error equations based on the available data in (2.13).Future researches in this field of study can now be turned to finding optimal sixteenth-order four-step without-memory iterations based on the general class (2.4)-(2.5). Furthermore, producing with-memory iterations according to this class can also be of researcher's interest for future studies.

Table 1: Typical forms of the weight functions satisfying (2.5).

Weight function
Form 1
Form 2
Form 3

,

, , (without the index n), and .

Table 2: Interesting choices of , , and in (2.13), which provide very efficient optimal root solvers.

Method

Error equation

1

−219

2

0
−120

3

−1
−54

4

0
0

5
0
−2
−96

3. Computational Examples

The contribution given in Section 2 is supported here through numerical works. We check the effectiveness of the novel methods (2.14) and (2.15) from our class of methods. For this reason, we have compared our new methods with Newton's method (NM), (1.1), (1.2), and (1.3). The nonlinear test functions are furnished in Table 3. The results of comparisons are given in Table 4 in terms of the number significant digits for each test function after some specified iterations.

Table 3: Test function, their roots, and the starting points.

Test functions
Roots
Guess

0.2588298273352688443917065956960
0.4

2.55419595283704303782966617379187
2

−1.25644837202179610636347166071080
−1.3

1.37318816571803393434603174397956
1.5

−0.94316510784112785413792698729205
−0.8

−1.41582823891434741754690479611753
−1.4

−1.33663345183152516454641803558884
−1.31

0.931908190324955905662501231761206
1.5

0.671280744928817559650708950129257
1

−0.785398163397448309615660845819875
−0.6

0.714180223545189655580041854721814
1

1.004983564135382106859677582107315
1

Table 4: Comparison of different methods for finding the simple roots of test functions.

Test functions
Methods
NM
(1.1)
(1.2)
(1.3)
(2.14)
(2.15)

IN
8
4
3
3
3
3

TNE
16
32
12
12
12
12

IN
9
4
3
3
3
3

TNE
18
32
12
12
12
12

IN
9
4
3
3
3
3

TNE
18
32
12
12
12
12

IN
8
4
3
3
3
3

TNI
16
32
12
12
12
12

IN
8
5
3
3
3
3

TNE
16
40
12
12
12
12

IN
9
4
3
3
3
3

TNE
18
32
12
12
12
12

IN
9
4
3
3
3
3

TNE
18
32
12
12
12
12

IN
8
4
3
3
3
3

TNE
16
32
12
12
12
12

IN
9
4
3
3
3
3

TNE
18
32
12
12
12
12

IN
8
4
3
3
3
3

TNE
16
32
12
12
12
12

IN
8
4
3
3
3
3

TNE
16
32
12
12
12
12

IN
9
4
3
3
3
3

TNE
18
32
12
12
12
12

All computations in this paper were performed in MATLAB 7.6 using variable precision arithmetic (VPA) to increase the number of significant digits. We have considered the following stopping criterion . In Table 4, shows that the absolute value of the given nonlinear function after three iterations is zero up to 448 decimal places. In Table 4, IN and TNE stand for iteration number and total number of evaluation. As shown in Table 4, the proposed method (2.14) is preferable to Newton’s method and some methods with fourth- and eighth-order of convergences. It is evident that (2.14) is more robust than the other competent from various orders. We also recall an important concern in using multi point iterations, which indicates that the high-order root solvers are very sensitive for initial guesses far from the root. And they are so powerful for starting points in the vicinity of the sought zero and so close.

Remark 3.1. If we need to solve a lot of equations from a large system of boundary-value problems, then the cost of function evaluations becomes important. Therefore, the proposed class (2.4)-(2.5) is valuable for solving such problems.

4. Concluding Remarks

In recent years, numerous works have been focusing on the development of more advanced and efficient methods for nonlinear scalar equations. Many methods have been developed, which improve the convergence rate of the Newton’s method. One practical drawback of so many methods is their slow rate of convergence. This paper has developed and established a rapid class of eighth-order iteration methods. Per iteration, the methods from our class require three evaluations of the function and one of its first derivatives; and therefore, the efficiency of the methods is equal to , which is better than that of the classical Newton’s method. Kung and Traub [10] conjectured that a multipoint iteration without memory based on evaluations of or its derivatives could achieve optimal convergence order . Newton’s method is an example, which agrees with Kung-Traub’s conjecture for , and the class of methods (2.4)-(2.5) is another example, which agrees with Kung-Traub’s hypothesis for . Thus, the suggested class (2.4)-(2.5) is effective and attracts the attention of researchers.

References

Z. H. Wang, “Numerical stability test of neutral delay differential equations,” Mathematical Problems in Engineering, vol. 2008, Article ID 698043, 10 pages, 2008.

N. A. Mir, R. Muneer, and I. Jabeen, “Some families of two-step simultaneous methods for determining zeros of nonlinear equations,” ISRN Applied Mathematics, Article ID 817174, 11 pages, 2011.

M. A. Noor, K. I. Noor, E. Al-Said, and M. Waseem, “Some new iterative methods for nonlinear equations,” Mathematical Problems in Engineering, Article ID 198943, 12 pages, 2010.

P. Sargolzaei and F. Soleymani, “Accurate fourteenth-order methods for solving nonlinear equations,” Numerical Algorithms, vol. 58, pp. 513–527, 2011.

F. Soleymani, M. Sharifi, and B. S. Mousavi, “An improvement of Ostrowski's and King's techniques with optimal convergence order eight,” Journal of Optimization Theory and Applications. In press

J. R. Sharma and R. Sharma, “A new family of modified Ostrowski's methods with accelerated eighth order convergence,” Numerical Algorithms, vol. 54, no. 4, pp. 445–458, 2010.

F. Soleymani and B. S. Mousavi, “A novel computational technique for finding simple roots of nonlinear equations,” International Journal of Mathematical Analysis, vol. 5, pp. 1813–1819, 2011.

A. Iliev and N. Kyurkchiev, Nontrivial Methods in Numerical Analysis (Selected Topics in Numerical Analysis), Lambeth Academy, 2010.

G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, Kluwer Academic, Boston, Mass, USA, 1994.

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.