On-line version ISSN 1807-0302
Comput. Appl. Math. vol.29 no.1 São Carlos 2010
Mehdi Dehghan*; Masoud Hajarian
Department of Applied Mathematics, Faculty of Mathematics and Computer Science, Amirkabir University of Technology, No. 424, Hafez Avenue, Tehran 15914, Iran E-mails: firstname.lastname@example.org, email@example.com / firstname.lastname@example.org, email@example.com
Finding the zeros of a nonlinear equation (x) = 0, is a classical problem which has nice applications in various branches of science and engineering. In this paper, we introduce four iterative methods which is based on the central-difference and forward-difference approximations to derivatives. It is proved that three of the four methods have cubic convergence and another method has quadratic convergence. The best property of these methods are that do not need to calculate any derivative. In order to demonstrate convergence properties of the introduced methods, several numerical examples are given.
Mathematical subject classification: 65H05; 41A25.
Keywords: Newton’s method,cubic convergence, quadratic convergence, nonlinear equations, iterative method.
A large number of papers have been written about iterative methods for the solution of the the nonlinear equations [3, 7, 8, 9, 10, 12, 13]. In this paper, we consider the problem of finding a simple root x* of a function : D ⊂ R → R i.e., (x*) = 0 and ' (x*) ≠ 0. The famous Newton's method for finding x* uses the iterative method
starting from some initial value x0. The Newton's method is an important and basic method where converges quadratically in some neighborhood of simple root x*. Chun  obtained the iterative method with convergence cubically given by
Also there is an modification of the Newton's method with third-order convergence as 
It is well-known that the forward-difference approximation for ' (x) at x is
If the derivative ' (xn) is replaced by the forward-difference approximation with h = (xn) i.e.
the Newton's method becomes
which is the famous Steffensen's method . The Steffensen's method is based on forward-difference approximation to derivative. This method is a tough competitor of Newton's method. Both the methods are of quadratic convergence, both require two functions evaluation per iteration but in contrast to Newton's method, Steffensen's method is derivative free. Chen  studied a particular class of these methods which contain the Steffensen's method as a special case. In , a modified Steffensen's method for the numerical solution of the system of nonlinear equations is studied. Amat et al.  considered a class of the generalized Steffensen iterations procedures for solving nonlinear equations on Banach space without any derivative.
Sometimes the applications of the iterative methods depending on derivatives are restricted in engineering. In this paper we introduce some methods which are based on the approximations to the derivative ' (xn) in each iteration. These methods are based on the central-difference and forward-difference approximations to the derivatives, respectively. The central-difference approximation for ' (x) at x is
We know that the leading errors in forward and central-difference formulae are (h) and (h2), respectively. It follows that the central-difference approximation is a more efficient than the forward-difference approximation to the derivative ' (x). If in (1.1) and (1.2) we replace derivatives ' (xn) by
we obtain two free derivative methods as follows:
Now we use the forward-difference approximation. If in (1.1) and (1.2), we replace derivatives ' (xn) by
we get two free derivative methods by the following
In the next section, we derive the convergence results of the iterative methods given by (1.3)-(1.6) and (1.8)-(1.11).
2 Main results
In this section we give the main results of this paper. We will give here the mathematical proof for the order of convergence of the methods given by (1.3)-(1.6) and (1.8)-(1.11).
Theorem 2.1. Let x* ∈ D be a simple zero of sufficiently differentiable function : D → R for an open interval D. If x0 is sufficiently close to x*, then the method defined by (1.3) and (1.4) has cubic convergence, and it satisfies the error equation
where en = xn– x*.
Proof. Let en and be the errors in xn and zn+1 respectively, i.e.
By using Taylor' theorem, we can get
It follows that
Now by substituting (2.3) and (2.6) into (1.4), we have
Again using Taylor's theorem we can write
Thus we get
Now it is not difficult to obtain
Now by using all the previous expressions, we obtain
Therefore we get the error equation
|The proof is finished.|
Theorem 2.2. Under the assumptions of Theorem 2.1, the method given by (1.5) and (1.6) has cubic convergence, and it verifies the error equation
Proof. Let = yn+1 – x*, by using the obtained equations in the proof of Theorem 2.1, we get
It is not difficult to get
By some calculations, we can show that
Hence we obtain the error equation
|The proof is completed.|
In the next theorems, we show that if the derivatives of the methods (1.1) and (1.2) are replaced by the forward-difference approximations, only the rate of convergence (1.8) and (1.9) is decreased and the obtained method has quadratic convergence. But the method defined by (1.10), (1.11) has cubic convergence.
Theorem 2.3. Under the assumptions of Theorem 2.1, the method defined by (1.8) and (1.9) has quadratic convergence, and it satisfies the error equation
Proof. Let en = xn– x* and = zn+1– x*, by substituting (2.3) and (2.4) into (1.9) and some calculations, we obtain
By using Taylor's theorem, we have
Also we obtain
Now by some calculations, it is not difficult to obtain
|The proof is finished.|
Theorem 2.4. Under the assumptions of Theorem 2.1, the method given by (1.10) and (1.11) has cubic convergence, and it verifies the error equation
Proof. Let en = xn– x* and = zn+1– x*, by substituting (2.3) and (2.4) into (1.11) and some calculations, we have
Hence we obtain
It follows from the above equations that
|The proof is completed.|
3 Numerical examples
In this section, in order to compare the methods introduced in this paper with Newton's method and Steffensen's method, we present some numerical examples. For this purpose we take two examples from the literature. Consider two nonlinear equations as
Tables 1 and 2 demonstrate the comparison of these methods for f (x) and g (x), respectively. All the tests are performed using MATLAB 7 which has a machine precision of around 10-16 on a Pentium IV. The numerical results indicate that the proposed iterative methods may be very efficient.
The problem of locating roots of nonlinear equations (or zeros of functions) occurs frequently in scientific work. In this paper, we have introduced some techniques for solving nonlinear equations. The techniques were based on the central-difference and forward-difference approximations to derivatives. We have shown that that three of the four methods have cubic convergence and another method has quadratic convergence. The introduced methods can be used for solving nonlinear equations without computing derivatives. Meanwhile, the methods introduced in this paper can be used to more class of nonlinear equations. The numerical examples shown in this paper illustrated the the efficiency of the new methods. We used the well known software MATLAB 7 to calculate the numerical results obtained from the proposed techniques.
 V. Alarcón, S. Amat, S. Busquier and D.J. López, A Steffensen's type method in Banach spaces with applications on boundary-value problems. J. Comput. Appl. Math., 216 (2008), 243-250. [ Links ]
 S. Amat, S. Busquier and J. M. Gutiérrez, Geometric constructions of iterative functions to solve nonlinear equations. J. Comput. Appl. Math., 157 (2003), 197-205. [ Links ]
 K.E. Atkinson, An Introduction to Numerical Analysis, second ed., John Wiley & Sons, New York (1989). [ Links ]
 D. Chen, On the convergence of a class of generalized Steffensen's iterative procedures and error analysis. Int. J. Comput. Math., 31 (1989), 195-203. [ Links ]
 C. Chun, A geometric construction of iterative functions of order three to solve nonlinear equations. Comput. Math. Appl., 53 (2007), 972-976. [ Links ]
 S.D. Conte and C. de Boor, Elementary Numerical Analysis: An Algorithmic Approach, 3rd edition, McGraw-Hill, Auckland (1986). [ Links ]
 M. Dehghan and M. Hajarian, On some cubic convergence iterative formulae without derivatives for solving nonlinear equations. Communications in Numerical Methods in Engineering, in press. [ Links ]
 M. Dehghan and M. Hajarian, New iterative method for solving non-linear equations with fourth-order convergence. Int. J. Comput. Math., in press. [ Links ]
 M. Frontini and E. Sormani, Modified Newton's method with third-order convergence and multiple roots. J. Comput. Appl. Math., 156 (2003), 345-354. [ Links ]
 H.H.H. Homeier, On Newton-type methods with cubic convergence. J. Comput. Appl. Math., 176 (2005), 425-432. [ Links ]
 D. Kincaid and W. Cheney, Numerical Analysis, second ed., Brooks/Cole, Pacific Grove, CA (1996). [ Links ]
 J.M. Ortega and W.C. Rheinboldt, Iterative Solution of Nonlinear Equations in Several Variables. Academic Press (1975). [ Links ]
 A.Y. Özban, Some new variants of Newton's method. Appl. Math. Lett., 17 (2004), 677-682. [ Links ]
 F.A. Potra and V. Pták, Nondiscrete induction and iterative processes. Research Notes in Mathematics, vol. 103, Pitman, Boston (1984). [ Links ]
* Corresponding author.