Next Article in Journal
Solving Integral Equations via Fixed Point Results Involving Rational-Type Inequalities
Next Article in Special Issue
Geometry of Harmonic Nearly Trans-Sasakian Manifolds
Previous Article in Journal
New Results on the Sequence Spaces Inclusion Equations of the Form FE+Fx Where F, F′ ∈ {w0, w, w}
Previous Article in Special Issue
Mathematical Analysis of the Vibratory Pile Driving Rate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Analysis of New Hybrid Algorithms for Solving Nonlinear Equations

by
Miguel Vivas-Cortez
1,
Naseem Zulfiqar Ali
2,
Awais Gul Khan
2 and
Muhammad Uzair Awan
2,*
1
Escuela de Ciencias Fsicas y Matemticas, Facultad de Ciencias Exactas y Naturales, Pontificia Universidad Catlica del Ecuador, Av. 12 de Octubre 1076, Apartado, Quito 17-01-2184, Ecuador
2
Department of Mathematics, Government College University, Allama Iqbal Road, Faisalabad 38000, Pakistan
*
Author to whom correspondence should be addressed.
Submission received: 16 May 2023 / Revised: 12 June 2023 / Accepted: 19 June 2023 / Published: 12 July 2023

Abstract

:
In this paper, we propose two new hybrid methods for solving nonlinear equations, utilizing the advantages of classical methods (bisection, trisection, and modified false position), i.e., bisection-modified false position (Bi-MFP) and trisection-modified false position (Tri-MFP). We implemented the proposed algorithms for several benchmark problems. We discuss the graphical analysis of these problems with respect to the number of iterations and the average CPU time.

1. Introduction

It is widely acknowledged that one of the most challenging problems in the mathematical sciences, particularly in numerical analysis, is solving the nonlinear equation Ψ ( ξ ) = 0 . In multiple improvements, numerous researchers have suggested, examined, and developed a variety of numerical approaches, employing various strategies, including Taylor’s series, modified homotopy perturbation methods, decomposition methods, variational iterative methods, and quadrature formulae. There is a wealth of literature available that highlight various approaches to solving nonlinear equations; for example, see [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15]. There are various ways to determine the roots of nonlinear equations, which include bracketing methods (bisection, trisection, false position, and modified false position) and open methods (Newton–Raphson, secant, Steffenson) hybrid methods (Badr et al. [16] and Sabharwal [17]), and metaheuristic algorithms, etc. There are benefits and drawbacks to each of the strategies mentioned above (there is no specific way to solve nonlinear equations effectively). Open approaches are quick but do not converge, whereas closed approaches are known for being slow but close. This research was motivated by the desire to propose fast and convergent methods, as opposed to the conventional methods, because hybrid methods are characterized by combining speed and convergence.
In fact, the idea of developing a hybrid technique by fusing two traditional methodologies is not new but rather has a long history. In 1995, Novak et al. [18] introduced a hybrid secant–bisection method. A novel hybrid bracketing strategy (bisection-false position) was proposed by Sabhrwal [17]. On the other hand, Badr et al. [16] created a hybrid method called Tri-FP that combines two closed algorithms. On fifteen benchmark linear and non-linear equations, they evaluated the method’s performance. Sabhrwal’s [17] approach (Bi-FP) was outperformed by the hybrid algorithm (Tri-FP) proposed by Badr et al. [16]. Recently, two new blended algorithms that make use of the advantages of open approaches and bracketing methods were suggested by Badr et al. [19].
The main objective of this research was to combine bracketing methods to propose two hybrid algorithms. It is well known that the modified false position method is better than the classical false position method. The first algorithm is a combination of bisection and modified false position methods (Bi-MFP). The second algorithm combines trisection and modified false position methods (Tri-MFP). We evaluated the proposed algorithms using fifteen benchmark functions against bisection, trisection, and modified false position methods. The numerical results showed that the proposed algorithms outperformed the classical bisection, trisection, and modified false position methods in terms of the number of iterations.

2. Classical Methods

In this section, we discuss four classical methods for finding the roots of non-linear equations. These methods are the bisection method, the trisection method, the false position method, and the modified false position method. This section includes details of these classical methods.

2.1. Bisection Method

The intermediate value theorem (IVT) says that Ψ ( ξ ) has at least a root or zero of Ψ ( ξ ) = 0 inside the interval [ ω 1 , ω 2 ] if it is continuous and defined on the interval [ ω 1 , ω 2 ] , fulfilling the relation Ψ ( ω 1 ) Ψ ( ω 2 ) < 0 with Ψ ( ω 1 ) and Ψ ( ω 2 ) are of opposite signs. Finding two real numbers ω 1 and ω 2 that lie in the interval [ ω 1 , ω 2 ] and the root in ω 1 ξ ω 2 with each step’s length being half of the interval’s initial length is the function of this method.
This process is continued until the interval imparts the desired precision as a result of obtaining the necessary root. The algorithms only work when the multiplicity of a root is odd. Additionally, the root converges very slowly and linearly. However, we can stop iterating when Ψ ( c ) is very close to zero or is very small. Consequently, the basic formula is
c = ω 1 + ω 2 2 .
The size of the interval is reduced by half at each iteration. Therefore, the tolerance ( e p s ) is determined from the following formula:
e p s = ω 2 ω 1 2 ϖ ,
where ϖ is the number of iterations. From (2), the number of iterations is found by
ϖ = log 2 ω 2 ω 1 e p s .
As a bracketing approach, the bisection method brackets the root in the range [ ω 1 , ω 2 ] , halving the range’s size after each iteration. As a result, this decreases the error between an iteration’s approximation root and precise root. However, the bisection method only works quickly if the approximate root is far from the interval’s endpoint; otherwise, more iterations are required to reach the root [3].

2.2. Trisection Method

The trisection method [7] is like the bisection method, except that it divides the interval [ ω 1 , ω 2 ] into three subintervals, while the bisection method divides the interval into two partial periods. This algorithm divides the interval [ ω 1 , ω 2 ] into three equal subintervals and searches for the root in the subinterval that contains different signs of the function values at the endpoints of these subintervals. If the condition of termination is true, then the iteration has finished its task; otherwise, the algorithm repeats the calculations.
To divide the distance [ ω 1 , ω 2 ] by ξ 1 into thirds and ξ 2 , one needs to know the location of ξ 1 and ξ 2 , as follows:
ξ 1 ω 1 = ω 2 ξ 2
ξ 2 ξ 1 = ξ 1 ω 1 .
By solving Equations (4) and (5), we obtain
ξ 1 = 2 ω 1 + ω 2 3
and
ξ 2 = ω 1 + 2 ω 2 3 .
The size of the interval [ ω 1 , ω 2 ] decreases to a third with each repetition. Therefore, the value e p s is determined from the following formula:
e p s = ω 2 ω 1 3 ϖ ,
where ϖ is the number of iterations. From (6) the number of iterations is found by
ϖ = log 3 ω 2 ω 1 e p s .
When we compare Equations (2) and (7), we can conclude that the iterations number of the trisection algorithm is less than the iteration number of the bisection algorithm. We might think that the trisection algorithm is better than the bisection algorithm, since it requires fewer iterations. However, it might be the case that one iteration of the trisection algorithm has an execution time longer than the execution time of one iteration of the bisection algorithm. Therefore, we will consider both execution time and the number of iterations to evaluate the different algorithms, see [2,3,5,6,20] and references therein.

2.3. False Position (Regula Falsi) Method

In this technique, one uses results that are known to be false to converge to the true root. This method chooses initial approximations ξ 0 and ξ 1 , such that Ψ ( ξ 0 ) Ψ ( ξ 1 ) < 0 . The new approximation value is be then obtained using the following relation:
ξ ϖ = ξ ϖ 1 Ψ ξ ϖ 1 ξ ϖ 1 ξ ϖ 2 Ψ ξ ϖ 1 Ψ ξ ϖ 2 ,
for ϖ 2 .
After that, and in order to decide which secant line to use, the product of Ψ ( ξ 2 ) and Ψ ( ξ 1 ) should be taken. If Ψ ( ξ 1 ) Ψ ( ξ 2 ) < 0 then choose ξ 3 as a line joining ( ξ 1 , Ψ ( ξ 1 ) ) and ( ξ 2 , Ψ ( ξ 2 ) ) , and if Ψ ( ξ 1 ) Ψ ( ξ 2 ) > 0 then choose ξ 3 as a line joining ( ξ 0 , Ψ ( ξ 0 ) ) and ( ξ 2 , Ψ ( ξ 2 ) ) . The process continues until the stopping criteria are satisfied.

2.4. Modified False Position Method

In this method, an improvement over the false position method is obtained by replacing the secant with straight lines with an even smaller slope, until ξ falls to the other side of the zero of Ψ ( ξ ) . To begin, we take the interval [ ω 1 , ω 2 ] in which the root lies and apply the false position to determine the value of ξ 1 , as given by
ξ 1 = ω 2 Ψ ω 2 ω 2 ω 1 Ψ ω 2 Ψ ω 1 .
Then, the products of Ψ ( ξ 1 ) , Ψ ( ω 1 ) and Ψ ( ω 2 ) decide the root interval. If Ψ ( ξ 1 ) Ψ ( ω 1 ) > 0 , then this implies that root lies in interval [ ξ 1 , ω 2 ] and we take ω 1 = ξ 1 in interval [ ω 1 , ω 2 ] with a fixed value of ω 2 in a complete iteration. Then, to find the next approximation ξ 2 , impose the half function value of ω 2 using the following relation:
ξ 2 = ω 2 Ψ ω 2 2 ω 2 ω 1 Ψ ω 2 2 Ψ ω 1 .
If Ψ ( ξ 1 ) Ψ ( ω 2 ) > 0 then we take ω 2 = ξ 1 and the fixed value of ω 1 in whole iteration. The next value of approximation ξ 2 can be found by imposing half function value of ω 1 , as given below
ξ 2 = ω 2 Ψ ω 2 ω 2 ω 1 Ψ ω 2 Ψ ω 1 2 .
This process continues until the stopping criteria satisfied.

3. Hybrid Algorithms

Instead of the classical methods (bisection, trisection, false position, modified false position), we propose two new hybrid algorithms: a bisection-modified false position method (Bi-MFP) and a trisection-modified false position method (Tri-MFP). These methods have the benefits of bisection, trisection, and modified false position methods. Badr et al. [16,17] proposed blended algorithms containing bisection, trisection, and false position methods. It is shown in this paper that our proposed algorithms performed better than their component algorithms, in the sense of the number of iterations.

3.1. Bisection-Modified False Position Method

This proposed algorithm (bisection-modified false position method) contains the advantages of both the bisection and modified false position methods. It is well known that the modified false position method outperforms the classical false position method, so we took advantage of the superiority of the modified false position method in our proposed algorithm. A flowchart and pseudocode for the Bi-MFP method are given in Figure 1 and Algorithm 1, respectively.
Algorithm 1: Hybrid Bi-MRF ( Ψ , ω 1 , ω 2 , e p s ) .
  • INPUT: The function Ψ ( ξ ) . The interval [ ω 1 , ω 2 ] where the root lies in, The absolute  error ( e p s )
  • OUTPUT: The root ( ξ ) . The value of Ψ ( ξ ) . Number of iterations ( ϖ ) . The interval [ ω 1 , ω 2 ] where the root lies in.
  • ϖ : = 0
  • ω 1 , a : = ω 1
  • ω 1 , b : = ω 1
  • ω 2 , a : = ω 2
  • ω 2 , b : = ω 2
  • while true do
  •     ϖ : = ϖ + 1
  •     ξ B : = ( ω 1 + ω 2 ) / 2
  •     ξ F : = ω 1 ( Ψ ( ω 1 ) ( ω 2 ω 1 ) ) / ( Ψ ( ω 2 ) Ψ ( ω 1 ) )
  •    if  | Ψ ( ξ B ) | < | Ψ ( ξ F ) |  then
  •       ξ : = ξ B
  •    else
  •       ξ : = ξ F
  •    end if
  •    if  | Ψ ( ξ ) | < = e p s  then
  •      return  ξ , Ψ ( ξ ) , ϖ , ω 1 , ω 2
  •    end if
  •    if  Ψ ( ω 1 ) Ψ ( ξ B ) < 0  then
  •       ω 2 , a : = ξ B
  •    else
  •       ω 1 , a : = ξ B
  •    end if
  •    if  Ψ ( ω 1 ) Ψ ( ξ F ) > 0  then
  •       ω 1 , b : = ξ F
  •       Ψ ( ω 1 ) : = Ψ ( ξ F )
  •       Ψ ( ω 2 ) : = Ψ ( ω 2 ) / 2
  •    else
  •      if  Ψ ( ω 2 ) Ψ ( ξ F ) > 0  then
  •          ω 2 , b : = ξ F
  •          Ψ ( ω 2 ) : = Ψ ( ξ F )
  •      end if
  •    end if
  •     ω 1 : = m a x ( ω 1 , a , ω 1 , b )
  •     ω 2 : = m i n ( ω 2 , a , ω 2 , b )
  •    end while
  • end while

3.2. Trisection-Modified False Position Method

Badr et al. [16] proposed a new blended algorithm called the trisection-false position algorithm, which is better than Sabharwal’s bisection-false position algorithm [17] in terms of iteration number and average CPU time. We proposed a new method having supremacy over the modified false position method, called the trisection-modified false position method (Tri-MFP). A flowchart and pseudocode of the proposed methods are given below in Figure 2 & Algorithm 2 respectively:

3.3. Empirical Study

MAPLE version 18.0, 64-bit, windows 7, Core(TM)2, T5600, 1.83 GHz, and 3.00 GB of memory were the elements that made up the software environment. In this paper, we used an absolute error of ( e p s = 10 10 ) to terminate all the algorithms. The iteration number and CPU time are effective tools for the comparison of algorithms, so we ran every algorithm ten times and calculated the average of the running time to obtain an accurate running time and avoid the problem of operating systems.
We evaluate the suggested technique for numerous benchmark problems as shown in Table 1, since it was not accurate enough to draw conclusions from only one problem. The abbreviations AppRoot, Error, LowerB, and UpperB in Table 2, Table 3, Table 4, Table 5 and Table 6 represent the approximate root, the difference between two consecutive roots, and the lower and upper bounds, respectively.
Figure 1 and Figure 2 shows flow charts of the proposed algorithms Bi-MFP and Tri-MFP. Meanwhile, Figure 3 and Figure 4 display the analysis of the benchmark problems shown in Table 1, according to the number of iterations and average CPU time, respectively.
Algorithm 2: Hybrid Tri-MRF ( Ψ , ω 1 , ω 2 , e p s ) .
  • INPUT: The function Ψ ( ξ ) . The interval [ ω 1 , ω 2 ] where the root lies in, The absolute error ( e p s )
  • OUTPUT: The root ( ξ ) . The value of Ψ ( ξ ) . Number of iterations ( ϖ ) . The interval [ ω 1 , ω 2 ] where the root lies in.
  • ϖ : = 0 , ω 1 , a : = ω 1 , ω 1 , b : = ω 1 , ω 2 , a : = ω 2 , ω 2 , b : = ω 2
  • while true do
  •     ϖ : = ϖ + 1
  •     ξ T 1 : = ( 2 ω 1 + ω 2 ) / 3
  •     ξ T 2 : = ( ω 1 + 2 ω 2 ) / 3
  •     ξ F : = ω 1 ( Ψ ( ω 1 ) ( ω 2 ω 1 ) ) / ( Ψ ( ω 2 ) Ψ ( ω 1 ) )
  •     ξ : = ξ T 1
  •     Ψ ξ : = Ψ ( ξ T 1 )
  •    if  | Ψ ( ξ T 2 ) | < | Ψ ( ξ ) |  then
  •       ξ : = ξ T 2
  •    end if
  •    if  | Ψ ( ξ F ) | < | Ψ ( ξ ) |  then
  •       ξ : = ξ F
  •    end if
  •    if  | Ψ ( ξ ) | < = e p s  then
  •      return  ξ , Ψ ( ξ ) , ϖ , ω 1 , ω 2
  •    end if
  •    if  Ψ ( ω 1 ) Ψ ( ξ T 1 ) < 0  then
  •       ω 2 , a : = ξ T 1
  •    else
  •      if  Ψ ( ξ T 1 ) Ψ ( ξ T 2 ) < 0  then
  •          ω 1 , a : = ξ T 1
  •          ω 2 , a : = ξ T 2
  •      else
  •          ω 1 , a : = ξ T 2
  •      end if
  •    end if
  •    if  Ψ ( ω 1 ) Ψ ( ξ F ) > 0  then
  •       ω 1 , b : = ξ F
  •       Ψ ( ω 1 ) : = Ψ ( ξ F )
  •       Ψ ( ω 2 ) : = Ψ ( ω 2 ) / 2
  •    else
  •      if  Ψ ( ω 2 ) Ψ ( ξ F ) > 0  then
  •          ω 2 , b : = ξ F
  •          Ψ ( ω 2 ) : = Ψ ( ξ F )
  •      end if
  •    end if
  •     ω 1 : = m a x ( ω 1 , a , ω 1 , b )
  •     ω 2 : = m i n ( ω 2 , a , ω 2 , b )
  •    end while
  • end while

3.4. Some New Numerical Experiments

We now consider some nonlinear polynomial equations from applied biomedical engineering and practical sciences.
Example 1 
(Blood rheology and fractional nonlinear equations Model [21]). Blood is modeled as a “Casson Fluid” as it is a non-Newtonian fluid. According to the Casson fluid model, basic fluids such as blood and water flow through a tube in such a way that the velocity gradient near the wall and the center core of the fluids travel as plugs, with minimal deformation. The following nonlinear fractional equation, which measures the drop in flow rate, was used to elaborate the plug flow of Casson fluids:
Ψ ( ξ ) = 1 16 7 ξ + 4 3 ξ 1 21 ξ 4 G ,
or where the reduction in flow rate is measured by G = 0.40 . Table 7 and Figure 5 and Figure 6 show numerical and graphical comparisons of the classical and proposed methods.
Example 2 
(Fluid Permeability in Biogels [21]). The following nonlinear equation comes from the specific hydraulic permeability relationship between the pressure gradient and fluid velocity in porous media (agarose gel or extracellular fiber matrix):
Ψ ( ξ ) = e ξ 3 20 κ ( 1 ξ ) 2 ,
Comparison of the proposed methods with the classical methods are numerically and graphically displayed in Table 8 and Figure 7 & Figure 8, respectively.
Example 3 
(beam position model [22]). Consider a mechanical issue with beam positioning that results in the following nonlinear equation:
Ψ ( ξ ) = ξ 4 + 4 ξ 3 24 ξ 2 + 16 ξ + 16 .
Comparison of the proposed methods with the classical methods are numerically and graphically displayed in Table 9 and Figure 9 & Figure 10, respectively.

4. Conclusions

This research study introduced a pair of innovative hybrid algorithms that enhanced the efficiency of bisection, trisection, and modified false position methods. These new hybrid algorithms—bisection-modified false position (Bi-MFP) and trisection-modified false position (Tri-MFP)—outperformed the methods previously proposed by [16,17]. Moreover, we tested some examples and arranged the results in a numerical and graphical manner. Moving forward, in the future, we will focus on some open and bracketing methods for solving non-linear equations.

Author Contributions

Conceptualization, N.Z.A. and A.G.K.; methodology, M.V.-C., N.Z.A., A.G.K. and M.U.A.; software, N.Z.A., A.G.K. and M.U.A.; validation, M.V.-C., N.Z.A., A.G.K. and M.U.A.; formal analysis, N.Z.A. and A.G.K.; investigation, M.V.-C., N.Z.A., A.G.K. and M.U.A.; writing—original draft preparation, M.V.-C., N.Z.A., A.G.K. and M.U.A.; writing—review and editing, A.G.K. and M.U.A.; visualization, M.V.-C., N.Z.A., A.G.K. and M.U.A.; supervision, A.G.K. and M.U.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by Pontificia Universidad Católica del Ecuador Proyect Título: “Algunos resultados Cualitativos sobre Ecuaciones diferenciales fraccionales y desigualdades integrales” Cod UIO2022.

Acknowledgments

The authors are thankful to the editor and anonymous reviewers for their valuable comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Burden, R.L.; Faires, J.D. Numerical Analysis, 9th ed.; Brooks/Cole, Cengage Learning: Boston, MA, USA, 2011. [Google Scholar]
  2. Chapra, S.C.; Canale, R.P. Numerical Methods for Engineers, 7th ed.; McGraw-Hill: Boston, MA, USA, 2015. [Google Scholar]
  3. Ehiwario, J.C.; Aghamie, S.O. Comparative Study of Bisection, Newton-Raphson and Secant Methods of Root-Finding Problems. IOSR J. Eng. 2014, 4, 1–7. [Google Scholar]
  4. Ernst, T. A method for q-calculus. J. Nonlinear Math. Phys. 2003, 10, 487–525. [Google Scholar] [CrossRef] [Green Version]
  5. Esfandiari, R.S. Numerical Methods for Engineers and Scientists Using MATLAB; CRC Press: Boca Raton, FL, USA, 2013. [Google Scholar]
  6. Harder, D.W.; Khoury, R. Numerical Methods and Modelling for Engineering; Springer International Publishing: Berlin/Heidelberg, Germany, 2016. [Google Scholar]
  7. Joe, D.H. Numerical Methods for Engineers and Scientists, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2001. [Google Scholar]
  8. Kac, V.; Cheung, P. Quantum Calculus; Springer: New York, NY, USA, 2002. [Google Scholar]
  9. Khan, A.G.; Ameen, F.; Awan, M.U.; Nonlaopon, K. Some new numerical schemes for finding the solutions to nonlinear equations. AIMS Math. 2022, 7, 18616–18631. [Google Scholar] [CrossRef]
  10. Mathews, J.H.; Fink, K.D. Numerical Methods Using Matlab, 4th ed.; Prentice-Hall Inc.: Upper Saddle River, NJ, USA, 2004; ISBN 0-13-065248-2. [Google Scholar]
  11. Nonlaopon, K.; Khan, A.G.; Ameen, F.; Awan, M.U.; Cesarano, C. Multi-Step Quantum Numerical Techniques for Finding the Solutions of Nonlinear Equations. Mathematics 2022, 10, 2595. [Google Scholar] [CrossRef]
  12. Nonlaopon, K.; Khan, A.G.; Ameen, F.; Awan, M.U.; Cesarano, C. Some New Quantum Numerical Techniques for Solving Nonlinear Equations. Symmetry 2022, 14, 1829. [Google Scholar] [CrossRef]
  13. Noor, M.A. Fifth-order convergent iterative method for solving nonlinear equations using quadrature formula. J. Math. Control Sci. Appl. 2018, 4, 95–104. [Google Scholar]
  14. Sana, G.; Mohammed, P.O.; Shin, D.Y.; Noor, M.A.; Oudat, M.S. On iterative methods for solving nonlinear equations in quantum calculus. Fractal Fract. 2021, 5, 60. [Google Scholar] [CrossRef]
  15. Srivastava, R.B.; Srivastava, S. Comparison of numerical rate of convergence of bisection, Newton and secant methods. J. Chem. Biol. Phys. Sci. 2011, 2, 472–479. [Google Scholar]
  16. Badr, E.; Almotairi, S.; Ghamry, A.E. A Comparative Study among New Hybrid Root Finding Algorithms and Traditional Methods. Mathematics 2021, 9, 1306. [Google Scholar] [CrossRef]
  17. Sabharwal, C.L. Blended root finding algorithm outperforms bisection and Regula Falsi Algorithms. Mathematics 2019, 7, 1118. [Google Scholar] [CrossRef] [Green Version]
  18. Novak, E.; Ritter, K.; Wozniakowski, H. Average-case ompitmality of a hybrid secant-bisection method. Math. Comput. 1995, 64, 1517–1539. [Google Scholar] [CrossRef]
  19. Badr, H.; Attiya, S.; Ghamry, A.E. Novel hybrid algorithms for root determining using advantages of open methods and bracketing methods. Alex. Eng. J. 2022, 61, 11579–11588. [Google Scholar] [CrossRef]
  20. Hasan, A.; Ahmad, N. Compartive study of a new iterative method with that Newtons Method for solving algebraic and transcesental equations. Int. J. Comput. Math. Sci. 2015, 4, 32–37. [Google Scholar]
  21. Fournier, R.L. Basic Transport Phenomena in Biomedical Engineering; Taylor and Francis: New York, NY, USA, 2007. [Google Scholar]
  22. Zafar, F.; Cordero, A.; Torregrosa, J.R. An efficient family of optimal eighth-order multiple root finders. Mathematics 2018, 6, 310. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flow chart of bisection-modified false position.
Figure 1. Flow chart of bisection-modified false position.
Axioms 12 00684 g001
Figure 2. Flow chart of trisection-modified false position.
Figure 2. Flow chart of trisection-modified false position.
Axioms 12 00684 g002
Figure 3. Comparison of six standard problems according to the number of iterations.
Figure 3. Comparison of six standard problems according to the number of iterations.
Axioms 12 00684 g003
Figure 4. Comparison of six standard problems according to CPU time.
Figure 4. Comparison of six standard problems according to CPU time.
Axioms 12 00684 g004
Figure 5. Comparison of methods according to iterations.
Figure 5. Comparison of methods according to iterations.
Axioms 12 00684 g005
Figure 6. Comparison of log of residual per iteration.
Figure 6. Comparison of log of residual per iteration.
Axioms 12 00684 g006
Figure 7. Comparison of the methods according to iterations.
Figure 7. Comparison of the methods according to iterations.
Axioms 12 00684 g007
Figure 8. Comparison of the log of residual per iteration.
Figure 8. Comparison of the log of residual per iteration.
Axioms 12 00684 g008
Figure 9. Comparison of the methods according to Iteration.
Figure 9. Comparison of the methods according to Iteration.
Axioms 12 00684 g009
Figure 10. Comparison of the log of residual per iteration.
Figure 10. Comparison of the log of residual per iteration.
Axioms 12 00684 g010
Table 1. Standard benchmark problems.
Table 1. Standard benchmark problems.
No.ProblemIntervalReferences
P1 ξ e ξ 7 [ 0 , 2 ] Calihoun [13]
P2 ξ cos ( ξ ) [ 0 , 1 ] Ehiworio [3]
P3 ξ sin ( ξ ) 1 [ 0 , 2 ] Mathews [18]
P4 sin ( ξ ) sinh ( ξ ) + 1 [ 1.5 , 4 ] Esfandiari [17]
P5 e ξ 3 ξ 2 [ 0 , 3 ] Hoffman [14]
P6 sin ( ξ ) ξ 2 [ 0.5 , 1 ] Chapra [11]
Table 2. Solutions of some problems using the bisection method.
Table 2. Solutions of some problems using the bisection method.
ProblemsBisection Method
IterAverage
CPU Time
Approx
Root
Function
Value
Lower
Bound
Upper
Bound
P1320.11181.52434520468 0.000000002 1.524345204221.52434520515
P2350.11520.739085134091 0.000000015 0.7390851340620.739085134119
P3320.11801.11415714072 0.000000000 1.114157140261.11415714119
P4330.13363.22158839905 0.000000000 3.221588398763.22158839934
P5330.13082.12539119889 0.000000000 2.125391198542.12539119923
P6330.12140.876726215414 0.000000000 0.8767262153550.876726215471
Table 3. Solutions of some problems using the trisection method.
Table 3. Solutions of some problems using the trisection method.
ProblemsTrisection Method
IterAverage
CPU Time
Approx
Root
Function
Value
Lower
Bound
Upper
Bound
P1200.12141.52434520541 0.000000003 1.524345204841.52434520656
P2220.12440.739085133325 0.000000000 0.7390851332620.739085133357
P3210.11201.11415714133 0.000000000 1.114157141141.11415714171
P4210.11803.22158839889 0.000000029 3.221588398413.22158839913
P5210.12122.12539119879 0.000000001 2.125391198232.12539119908
P6200.13380.876726215456 0.000000000 0.8767262151700.876726215600
Table 4. Solutions of some problems using the modified regula falsi method.
Table 4. Solutions of some problems using the modified regula falsi method.
ProblemsModified Regula Falsi Method
IterAverage
CPU Time
Approx
Root
Function
Value
Lower
Bound
Upper
Bound
P1200.13041.52434520493 0.0000000062 1.524345204822.0000000000
P290.13700.739085133213 0.0000000000 0.7390851331711.0000000000
P360.12441.11415714087 0.0000000000 1.114157140871.11415714304
P4310.15243.22158839849 0.0000000081 3.221588397774.0000000000
P5330.12442.12539119823 0.0000000031 2.125391197613.0000000000
P6130.13060.876726215392 0.000000000 0.8767262153721.0000000000
Table 5. Solutions of some problems using the bisection and modified regula falsi methods.
Table 5. Solutions of some problems using the bisection and modified regula falsi methods.
ProblemsBisection-Modified Regula Falsi Method
IterAverage
CPU Time
Approx
Root
Function
Value
Lower
Bound
Upper
Bound
P190.11821.52434520539 0.000000000 1.524345204451.52772146940
P270.11180.739085133226 0.000000000 0.7390851297060.745369013289
P3150.11801.11415714088 0.000000000 1.114146255001.11425614878
P4100.12143.22158839955 0.000000012 3.221588399433.22238911820
P590.10862.12539119894 0.000000001 2.125391185212.12787084220
P670.13080.876726215441 0.0000000001 0.8767262106850.877268445426
Table 6. Solutions of some problems using the trisection and modified regula falsi methods.
Table 6. Solutions of some problems using the trisection and modified regula falsi methods.
ProblemsTrisection-Modified Regula Falsi Method
IterAverage
CPU Time
Approx
Root
Function
Value
Lower
Bound
Upper
Bound
P160.11801.52434520508 0.00000255309 1.524345204131.52441127915
P260.09640.739085133236 0.0003113615 0.7390851331170.739643235290
P3120.11181.11415714132 0.000001774 1.114157140461.11415841776
P480.10563.22158839889 1.5 × 10 9 3.221588398783.22215305257
P570.12762.12539119914 1.0 × 10 9 2.125391198402.12548466697
P650.12120.876726215473 2.0 × 10 10 0.8767262151140.876727040681
Table 7. Results of the above example using Bi, Tri, MFP, Bi-MFP, and Tri-MFP methods.
Table 7. Results of the above example using Bi, Tri, MFP, Bi-MFP, and Tri-MFP methods.
MethodBlood Rheology and Fractional Nonlinear Equations Model
IterApprox
Root
Function
Value
Lower
Bound
Upper
Bound
Bisection340.104698652342 1.6 × 10 9 0.1046986522840.104698652401
Trisection230.104698652104 1.3 × 10 9 0.1046986520830.104698652115
MFP440.104698651667 2.86 × 10 10 0.0000000000000.104698651748
Bi-MFP190.104698651542 0.0000000000 0.1046981811520.104701995849
Tri-MFP130.104698651485 7.335 × 10 7 0.1046983578610.104700239538
Table 8. The results of the above example using the Bi, Tri, MFP, Bi-MFP, and Tri-MFP methods.
Table 8. The results of the above example using the Bi, Tri, MFP, Bi-MFP, and Tri-MFP methods.
MethodFluid Permeability in Biogels
IterApprox
Root
Function
Value
Lower
Bound
Upper
Bound
Bisection301.00003698747 7.51 × 10 13 1.000036986541.00003698840
Trisection201.00003698880 3.31 × 10 13 1.000036988231.00003698909
MFP Fail
Bi-MFP181.00003698850 2.09 × 10 13 1.000036988081.00003699241
Tri-MFP151.00003698885 1.41 × 10 12 1.000036988481.00003699616
Table 9. Results of the above example using the Bi, Tri, MFP, Bi-MFP, and Tri-MFP methods.
Table 9. Results of the above example using the Bi, Tri, MFP, Bi-MFP, and Tri-MFP methods.
MethodBeam Position Model
IterApprox
Root
Function
Value
Lower
Bound
Upper
Bound
Bisection34−0.535898384669 1.0 × 10 8 −0.535898384698−0.535898384639
Trisection22−0.535898384057 3.0 × 10 8 −0.535898384089−0.535898384042
MFP15−0.535898384849 6.0 × 10 10 −1.000000000000−0.535898384808
Bi-MFP15−0.535898384738 1.00 × 10 8 −0.535919189454−0.535888671875
Tri-MFP12−0.535898384898 1.06 × 10 5 −0.535899563639−0.535896741125
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Vivas-Cortez, M.; Ali, N.Z.; Khan, A.G.; Awan, M.U. Numerical Analysis of New Hybrid Algorithms for Solving Nonlinear Equations. Axioms 2023, 12, 684. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms12070684

AMA Style

Vivas-Cortez M, Ali NZ, Khan AG, Awan MU. Numerical Analysis of New Hybrid Algorithms for Solving Nonlinear Equations. Axioms. 2023; 12(7):684. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms12070684

Chicago/Turabian Style

Vivas-Cortez, Miguel, Naseem Zulfiqar Ali, Awais Gul Khan, and Muhammad Uzair Awan. 2023. "Numerical Analysis of New Hybrid Algorithms for Solving Nonlinear Equations" Axioms 12, no. 7: 684. https://0-doi-org.brum.beds.ac.uk/10.3390/axioms12070684

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop