 Research note
 Open Access
 Published:
Quadratically convergent algorithm for computing real root of nonlinear transcendental equations
BMC Research Notes volume 11, Article number: 909 (2018)
Abstract
Objectives
The present paper describes a new algorithm to find a root of nonlinear transcendental equations. It is found that RegulaFalsi method always gives guaranteed result but slow convergence. However, Newton–Raphson method does not give guaranteed result but faster than RegulaFalsi method. Therefore, the present paper used these two ideas and developed a new algorithm which has better convergence than RegulaFalsi and guaranteed result. One of the major issue in Newton–Raphson method is, it fails when first derivative is zero or approximately zero.
Results
The proposed method implemented the failure condition of Newton–Raphson method with better convergence. Error calculation has been discussed for certain real life examples using Bisection, RegulaFalsi, Newton–Raphson method and new proposed method. The computed results show that the new proposed quadratically convergent method provides better convergence than other methods.
Introduction
Most of the real lifeproblems are nonlinear in nature therefore it is a challenging task for the mathematician and engineer to find the exact solution of such problems [1, 2]. In this reference, a number of methods have been proposed/implemented in the last two decades [1, 3,4,5,6,7,8]. Analytical solutions of such nonlinear equations are very difficult, therefore only numerical method based iterative techniques are the way to find approximate solution. In the literature, there are some numerical methods such as Bisection, Secant, RegulaFalsi, Newton–Raphson, Mullers methods, etc., to calculate an approximate root of the nonlinear transcendental equations. It is well known [1, 3,4,5,6,7,8,9,10,11, 14] that all the iterative methods require one or more initial guesses for the initial approximations.
In RegulaFalsi method, two initial guesses are taken in such a way that the corresponding function values have opposite signs. Then these two points are connected through the straight line and next approximation is the point where this line intersect the xaxis. This method gives guaranteed result but slow convergence therefore several researchers have improved this standard RegulaFalsi method into different hybrid models to speed up the convergence [1, 3,4,5, 7, 10, 11, 15, 16]. Thus previously published works have revised/implemented RegulaFalsi method in several ways to obtain better convergence. However, it is found that modified form of RegualFalsi method becomes more complicated from computational point of view. Therefore, in the present work RegualFalsi method has been used as its standard form with Newton–Raphson method and found better convergence. Newton–Raphson method is generally used to improve the result obtained by one of the above methods. This method uses the concept of tangent at the initial approximation point. The next approximate root is taken those value where the tangent intersect the xaxis. So this method fails where tangent is parallel to xaxis, i.e. the derivative of the function is zero or approximately zero. The order of convergence of Newton–Raphson method is two, therefore it converges very rapidly than other methods (Bisection, RegulaFalsi, etc.). However it does not always give guaranteed root. Many scientists and engineers have been proposed different hybrid models on Newton–Raphson method [8, 9, 12,13,14, 17,18,19,20,21,22].
It is clear from the survey [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22], that the most of new algorithms are either based on three classical methods namely Bisection, RegulaFalsi and Newton–Raphson or created by hybrid processes. In the present work, the proposed new algorithm is based on standard RegulaFalsi and Newton–Raphson methods, which provides guaranteed results and higher order convergence over RegulaFalsi method. The new proposed algorithm will work even the first derivative equals to zero where Newton–Raphson method fails.
Main text
Consider a continuous function f(x) between a and b such that f(a) and f(b) having opposite signs, consequently \(f(a) \cdot f(b) < 0\). Without loss of generality, assume that f(a) is negative and f(b) is positive, and \(f(a) < f(b)\), hence at least one root lies between a and b. From RegulaFalsi method, the first approximate root can be calculated by using the formula
and, the first approximate root by using Newton–Raphson method is
where \(f'(a)\) indicates the derivative of f(x) at \(x=a\).
Now, in the present proposed algorithm, we take the average of the iterations in Eqs. (1) and (2) as our first approximate root \({\widehat{x}}\) and follow the conditions given below for further iterations:

Choose two values a and b where the root exists as in RegulaFalsi method.

Select the value such that the corresponding function is closer to zero as a and the other one as b, i.e. \(f(a) < f(b)\).

If first derivative at a is zero (i.e., \(f'(a) =0\)) then interchange the values of a and b i.e., interchange (a, b) to (b, a).
The generalization of this process is described in the following section.
Formulation of proposed algorithm
Recall the Eqs. (1) and (2) in terms of iteration formulae by replacing a, b, x by \(x_{n1}, x_{n+1}, x_n\) respectively, as follows
and
where n is the iteration number and \(f(x_{n1}) < f(x_{n+1})\). Now the average of Eqs. (3) and (4) is
After simplification of (5), we get
The value \(x_n\) in Eq. (6) gives the iterative formula with \(f(x_{n1}) < f(x_{n+1})\). If \(f'(x_{n1}) = 0\), then Eq. (6) gives undefined value, then we have to interchange the values \(x_{n1}\) and \(x_{n+1}\).
The following Theorem gives the generalization of above formulation.
Lemma 1
Let f(x) be a continuous function and (a, b) be a sufficiently small interval such that \(f(a)f(b) < 0\), and \(f'(x)\) exists on [a, b]. Then the approximation of a root of f(x) can be find using the iterative formula given in Eq. (6).
Steps for calculating a root

I.
Select two initial approximations \(x_{n1}\) and \(x_{n+1}\) such that product of the corresponding function values must be negative, i.e. \(f(x_{n1})f(x_{n+1})<0\).

II.
Now calculate \(x_n\) using the formula given in Eq. (6). Check \(f(x_n) = 0\), if so, then \(x_n\) is required root and process stop. Otherwise we check the following possible conditions.

(i)
For \(f(x_n)f(x_{n1})<0\), suppose \(f(x_{n1}) < f(x_{n})\) then \(x_n\) replace by \(x_{n1}\) and \(x_{n1}\) replace by \(x_{n}\).

(ii)
For \(f(x_n)f(x_{n+1})<0\), suppose \(f(x_{n}) < f(x_{n+1})\) then \(x_{n+1}\) replace by \(x_{n}\).

(i)

III.
If \(f'(x_{n1}) \approx 0\) then interchange \(x_{n1}\) and \(x_{n+1}\).

IV.
Repeat steps I, II and III until we get required approximate solution.
The implementation of the proposed algorithm in Matlab is also discussed (See, Additional file 1). This algorithm would help to implement the manual calculations in commercial packages such as Maple, Mathematica, SCILab, Singular,etc.
Order of convergence
The order of converges of any iterative method is defined as
where p is the order of convergence and c is a positive finite constant. The following Theorem shows the order of convergence of the proposed algorithm is quadratic.
Theorem 2
Let \(\beta\) be a exact root of a continuous function f(x) and [a, b] be a sufficiently small neighbourhood of \(\beta\). Then the sequence \(\{x_n\}\) generated by the iterative formula (6) is at least quadratically convergent.
Proof
If \(x_n\) is an approximate value of \(\beta\) and \(e_n\) is the error of the \(x_n\) then
In the similar way, \(x_{n+1} = \beta + e_{n+1}\) and \(x_{n+2} = \beta + e_{n+2}\). By proposed algorithm,
putting values of \(x_{n},x_{n+1}\) and \(x_{n+2}\) in above equation, we get
After simplification of above equation using Taylor’s series, we get
Putting \(\frac{f''(\beta )}{2f'(\beta )} = A\) (constant), then
We have, \(e_{n+2} = c e_{n+1}^p\), \(c>0\); \(e_{n+1} = c e_{n}^p\); and \(e_{n} = c^{1/p} e_{n+1}^{1/p}\). Substituting in Eq. (8) and after simplification, we get
If \(p=1\), then
where \(c^*=c_1^*+c_2^*\). From Eq. (7), it shows that the iterative formula (6) has quadratic convergent. \(\square\)
Results
This section provides three examples to discuss the algorithm presented in “Main text” section and comparisons are taken into account to conform that the algorithm is more efficient than other existing methods. Moreover, it is also observed that the proposed method takes less time in comparison of Regula Falsi method but takes more convergence time in comparison of Newton–Raphson method.
Example 3
Consider a transcendental equation of the form
We compute a root of Eq. (9) in the interval (0,1) using Bisection, RegulaFalsi, Newton–Raphson and proposed algorithm.
Table 1 shows that the comparison between Bisection, RegulaFalsi, Newton–Raphson and proposed method. The errors given in table are indicating the difference between two consecutive iterations. It is clear that the proposed method rapid convergence towards exact root than bisection and Regulafasi method. It is not as speed convergent as Newton–Raphson method, but provide guaranteed result.
It is found that bisection and RegulaFalsi methods converged after 22 and 14 iterations respectively (Table 1), while the proposed algorithm converged in 7th iteration. Thus proposed method is efficient over bisection and RegulaFalsi methods. It is also clear that both of the methods (proposed and Newton–Raphson) are converged in 7th iteration. But one of the main advantage of proposed method is that it gives guaranteed result over the Newton–Raphson method. Therefore, the proposed method is not only reduce the computational affords but also provide the guaranteed result for solving the real life problem.
The error estimation after the 3rd iteration (Table 1), show that the proposed method having \(5.4\%\) error in comparison to Bisection (\(20\%\)), RegulaFalsi (\(9.57\%\)) and Newton–Raphson (\(22.91\%\)) methods. Thus, the proposed method is also efficient for error estimation.
Most of the real life problems take too much computational time for convergence because of the complex flow physics and higher degree polynomial equations. Therefore, the proposed method is useful also for solving such the real life problem.
Example 4
Consider a transcendental equation of the form
We compute a root of above equation in the interval (1,3) using Bisection, RegulaFalsi, Newton–Raphson and proposed algorithm.
Table 2 shows that bisection method converged in 21st iteration however the remaining three methods (RegulaFalsi, Newton–Raphson, proposed) converged after 3rd iteration. Thus, as similar to the Example 3, proposed method is efficient to solve this logarithmic problem.
Example 5
Consider the real root of \(f(x) = 1x^2\) in the interval (0, 2).
In this example, Newton–Raphson method unable to find the real root because of \(f'(x)\) is zero at initial approximation \(x=0\) as in Eq. (2). However, the proposed method give approximate root with any range of initial approximation as given in Table 3. Thus, the proposed method is also applicable to such equation where \(f'(x)\) is zero with any one of the initial approximations. Since \(f'(x)=0\) at \(a=0\), therefore according to the proposed method, the initial approximations automatically interchange i.e., \(a=2\) and \(b=0\) as shown in Table 3.
Limitations
The order of the presented method is quadratic. Even though there are methods with a higher order of convergence, the proposed method is simple and efficient quadratic convergence method.
Abbreviations
 BM:

Bisection method
 RF:

RegulaFalsi
 N–R:

Newton–Raphson
 PM:

proposed method
References
 1.
Sharma JR, Goyal RK. Fourth order derivative methods for solving nonlinear equations. Int J Comput Math. 2006;83(1):101–6.
 2.
Gottlieb RG, Thompson BF. Bisected direct quadratic RegulaFalsi. Appl Math Sci. 2010;4(15):709–18.
 3.
Wu W, Wu H. On a class of quadratic convergence iteration formula without derivatives. Appl Math Comput. 2000;107:77–80.
 4.
Wu X, Shen Z, Jianlin X. An improved Regula Falsi method with quadratic convergence of both diameter and point for enclosing simple zeros of nonlinear equations. Appl Math Comput. 2003;144:381–8.
 5.
Mamta VK, Kukreja VK, Singh S. On some thirdorder iterative methods for solving nonlinear equations. Appl Math Comput. 2005;171:272–80.
 6.
Noor MA, Ahmad F. Numerical comparison of iterative methods for solving nonlinear equations. Appl Math Comput. 2006;180:167–72.
 7.
Noor MA, Noor KI, Khan WA, Ahmad F. On iterative methods for nonlinear equations. Appl Math Comput. 2006;183:128–33.
 8.
Thota S, Srivastav VK. Interpolation based hybrid algorithm for computing real root of nonlinear transcendental functions. Int J Math Comput Res. 2014;2(11):729–35.
 9.
Hussain S, Srivastav VK, Thota S. Assessment of interpolation methods for solving the real life problem. Int J Math Sci Appl. 2015;5(1):91–5.
 10.
Chen J, Li W. An exponential Regula Falsi method for solving nonlinear equations. Numer Algorithms. 2006;41:327–38.
 11.
Chen J, Li W. An improved exponential Regula Falsi methods with quadratic convergence of both diameter and point for solving nonlinear equations. Appl Numer Math. 2007;57:80–8.
 12.
Johan, Ronald. The Newton–Raphson method. Int J Math Educ Sci Technol. 1995;26(2):177–93.
 13.
Wu XY, Fu D. New highorder convergence iteration methods without employing derivatives for solving nonlinear equations. Comput Math Appl. 2001;41:489–95.
 14.
Frontini M, Sormani E. Modified Newton’s method with thirdorder convergence and multiple roots. J Comput Appl Math. 2003;156:345–54.
 15.
Dowell M, Jarratt P. A modified RegulaFalsi method for computing the real root of an equation. BIT Numer Math. 1971;11:168–74.
 16.
Gottlieb RG, Thompson BF. Bisected direct quadratic RegulaFalsi. Appl Math Sci. 2010;4(15):709–18.
 17.
Wu XY, Xia JL, Shao R. Quadratically convergent multiple roots finding method without derivatives. Comput Math Appl. 2001;42:115–9.
 18.
Amata S, Busquiera S, Gutierrezb JM. Geometric constructions of iterative functions to solve nonlinear equations. J Comput Appl Math. 2003;157:197–205.
 19.
Kou J, Li Y, Wang X. A variant of superHalley method with accelerated fourthorder convergence. Appl Math Comput. 2007;186:535–9.
 20.
Mamta VK, Kukreja VK, Singh S. On a class of quadratically convergent iteration formulae. Appl Math Comput. 2005;166:633–7.
 21.
Sagraloff M. Computing real roots of real polynomials. 2015. arXiv:1308.4088v2 [cs.SC].
 22.
Abbott J. Quadratic interval refinement for real roots. ACM Commun Comput Algebra. 2014;48(1):3–12.
 23.
Chapra Steven C. Solutions Manual to accompany. Applied numerical methods With MATLAB for engineers and scientists; Pennsylvania: Mc Graw Hill Higher education. ISBN 9780073401065. 2010.
Authors' contributions
ST is involved in creating the proposed algorithm and implementation of the algorithm in Matlab. VKS is participated in modifications, corrections and writing of the manuscript. Both authors read and approved the final manuscript.
Acknowledgements
The author is thankful to the editor for providing valuable inputs to improve the present format of manuscript.
Competing interests
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated and analyzed during the current study are available from the corresponding author on reasonable request.
Consent to publish
Not applicable.
Ethics approval and consent to participate
Not applicable.
Funding
Not applicable.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Affiliations
Corresponding author
Additional file
Additional file 1.
Implementation of the proposed method in Matlab. In the Additional file, we provide the implementation of the proposed method inMatlab code similar to RegulaFalsi method in [23] by creating a data type NewAlgorithm (f, a, b, esp, n),as given in Additional file, where f is a given transcendental equation, a, b are the initial approximationof the root, esp is the relative error and n is the number of iterations required.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Thota, S., Srivastav, V.K. Quadratically convergent algorithm for computing real root of nonlinear transcendental equations. BMC Res Notes 11, 909 (2018). https://0doiorg.brum.beds.ac.uk/10.1186/s131040184008z
Received:
Accepted:
Published:
Keywords
 Root of transcendental equations
 RegulaFalsi method
 Newton–Raphson method
 Quadratic convergence
Mathematics Subject Classification
 65Hxx
 65H04