Generalization of Liu–Zhou Method for Multiple Roots of Applied Science Problems

Abstract

Some optimal and non-optimal iterative approaches for computing multiple zeros of nonlinear functions have recently been published in the literature when the multiplicity $\theta$ of the root is known. Here, we present a new family of iterative algorithms for multiple zeros that are distinct from the existing approaches. Some special cases of the new family are presented and it is found that existing Liu-Zhou methods are the special cases of the new family. To check the consistency and stability of the new methods, we consider the continuous stirred tank reactor problem, isentropic supersonic flow problem, eigenvalue problem, complex root problem, and standard test problem in the numerical section and we find that the new methods are more competitive with other existing fourth-order methods. In the numerical section, the error of the new methods confirms their robust character.

1. Introduction

Numerous physical and technical applications [1,2,3] have demonstrated the significance of solving nonlinear equations in the numerical field with the rapid growth. Such problems exist in a variety of domains within the natural and physical sciences, such as those involving heat and fluid movement, initial and boundary value issues, and problems with global positioning systems. Except for a few nonlinear equations, finding the solution via an analytical approach is practically impossible. Iterative methods thus offer a desirable substitute for solving problems of this type.

To find the multiple roots of a nonlinear equation of the form $f\left(t\right)=0$, where $f\left(t\right)$ is a real function defined in a domain $D\subseteq \mathbb{R}$, the modified Newton method [4,5,6] is used very commonly. The modified Newton method is given by

(1)$t_{i+1}=t_i-\theta {\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}},i=0,1,2,\dots$

with the given multiplicity $\theta \ge 1$ and $\theta \in \mathbb{N}$; the scheme provided by (1) is optimal following the Kung–Traub conjecture [7] and can be used to find the desired multiple roots with quadratic convergence.

Many scientists have been working on iterative methods for finding multiple roots with higher-order convergence and efficiency in recent decades. Numerous higher-order optimal and non-optimal methods have been developed in the literature by Behl et al. [8], Behl et al. [9], Behl et al. [10], Dong [11], Geum et al. [12], Hansen [13], Hueso et al. [14], Kansal et al. [15], Li et al. [16], Li et al. [17], Liu and Zhou [18], Neta [19,20], Osada [21], Sharma and Kumar [22,23], Sharma and Sharma [24], Sharifi et al. [25], Soleymani et al. [26], Soleymani and Babajee [27], Thukral [28], Victory and Neta [29], and Zhou et al. [30,31]. These methods are two-step and three-step methods with convergence orders of three, four, six, and eight. Thus, motivated by this, in this work, our aim is to develop optimal multi-point iterative methods of a higher convergence order that may use computations that are small in number, as needed.

Taking into account the aforementioned concerns, we offer a fourth-order family that, according to Kung–Traub hypothesis [7], has optimal fourth-order convergence and requires three additional functions of information per iteration. The proposed approach is made up of two steps, the first of which uses Newton’s iteration (1) and the second of which uses Newton-type iteration. An iterative scheme is distinct in that each iteration calls for one function and two derivatives. The main benefit of the new family of methods is that the Liu–Zhou scheme [18] is a special case of the family.

The information contained in the rest of this article is outlined as follows. In Section 2, simple definitions are given. Section 3 develops the fourth-order generalized iterative scheme and examines its convergence. Section 4 examines a few practical science problems to examine the methodological stability and validate the theoretical findings. This part also includes a comparison with existing methods and some graphs displaying the calculated outcomes. Concluding remarks are reported in Section 5.

2. Basic Definitions

2.1. Multiple Root

A root $t^{\ast }$ of $f\left(t\right)=0$ is a zero of multiplicity $\theta$ if $t\ne t^{\ast }$—we can write $f\left(t\right)={(t-t^{\ast })}^{\theta }g\left(t\right),$ where $g\left(t^{\ast }\right)\ne 0$. The function $f\in {\mathbb{C}}^{\theta }[a,b]$ has a zero of multiplicity $\theta$ at $t^{\ast }$ in $(a,b)$ if

$f\left(t^{\ast }\right)=f^{\prime }\left(t^{\ast }\right)=f^{\prime \prime }\left(t^{\ast }\right)=\dots =f^{(\theta -1)}\left(t^{\ast }\right)=0,$

but $f^{\left(\theta \right)}\left(t^{\ast }\right)\ne 0$. If $\theta =1$, the root is called a simple zero.

2.2. Order of Convergence

Let ${\left\{t_i\right\}}_{i\ge 0}$ be a sequence of iterative points which converges to $t^{\ast }$. Then, the convergence is said to be of order p, $p>1,$ if there exists M, $M>0,$ and $k_0$ such that

$\parallel t_{i+1}-t^{\ast }\parallel \le M\parallel t_i-t^{\ast }{\parallel }^p\mathrm{for}\mathrm{all}i\ge i_0$

or

$\parallel e_{i+1}\parallel \le M\parallel e_i{\parallel }^p\mathrm{for}\mathrm{all}i\ge i_0,$

where $e_i=t_i-t^{\ast }$. The convergence is linear if $p=1$, quadratic if $p=2$, and there exists M such that $0<M<1$.

2.3. Error Equation

Let $e_i=t_i-t^{\ast }$ be the error in the ith iteration; we designate the relation

$e_{i+1}=L{e}_i^p+O\left(e_i^{p+1}\right)$

as the error equation. Here, L is an asymptotic error constant, p is the order of convergence, and $O\left(e_i^{p+1}\right)$ denotes the higher power of $e_i^p.$

2.4. Computational Order of Convergence

Assume that $t_{i+2}$, $t_{i+1}$, and $t_i$ are three consecutive iterations that are near to $t^{\ast }$, with $t^{\ast }$ being the zero of the function f. Then, the following formula is used to approximate the computational order of convergence (COC) (see [32]):

(2)$\mathrm{COC}={\displaystyle \frac{ln|(t_{i+2}-t^{\ast })/(t_{i+1}-t^{\ast })|}{ln|(t_{i+1}-t^{\ast })/(t_i-t^{\ast })|}},i=1,2,\dots$

2.5. Kung–Traub Hypothesis

The Kung–Traub hypothesis [7] states that iterative methods have convergence order $2^p$ if they require $p+1$ function evaluations per step. Such methods are called optimal methods.

3. Formulation of Scheme

The design and convergence analysis of the suggested scheme, which is the primary contribution of this paper, are covered in this section. To find multiple zeros with multiplicity $\theta >1$, we take into account the following ideal fourth-order family:

(3)$\begin{array}{cc}{u}_i\hfill & =t_i-\theta {\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}},\\ t_{i+1}\hfill & =u_i-\theta G\left(h_i\right){\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}},\end{array}$

where $h_i={\displaystyle \frac{x_i}{a_1+a_2{x}_i}}$, $a_1+a_2{x}_i\ne 0$, $x_i={\left({\displaystyle \frac{f^{\prime }\left(u_i\right)}{f^{\prime }\left(t_i\right)}}\right)}^{{\displaystyle \frac{1}{\theta -1}}}$, $a_1$ and $a_2$ are not simultaneously zero, and the function $G\left(h\right)$ is analytic in the neighborhood of 0. Note that the second factor $G\left(h\right)$ is multiplied, so the factor is called the weight function.

In the followings section, we will examine certain circumstances in which Scheme (3) achieves the highest feasible order of convergence with the minimum number of functions. The software Mathematica (v. 12.0.0.0 has been used) and other computer algebra systems were used to manage lengthy calculations.

Let $e_i=t_i-t^{\ast }$ be the error at i-th iteration. Taylor’s development about $\alpha$ yields

(4)$f\left(t_i\right)={\displaystyle \frac{f^{\left(\theta \right)}\left(t^{\ast }\right)}{\theta !}}{e}_i^{\theta }\left(1+D_1{e}_i+D_2{e}_i^2+D_3{e}_i^3+D_4{e}_i^4+\cdots \right)$

and

(5)$f^{\prime }\left(t_i\right)={\displaystyle \frac{f^{\left(\theta \right)}\left(t^{\ast }\right)}{\theta !}}{e}_i^{\theta -1}\left(\theta +(\theta +1)D_1{e}_i+(\theta +2)D_2{e}_i^2+(\theta +3)D_3{e}_i^3+(\theta +4)D_4{e}_i^4+\cdots \right),$

where $D_n={\displaystyle \frac{\theta !}{(\theta +n)!}}{\displaystyle \frac{f^{(\theta +n)}\left(t^{\ast }\right)}{f^{\left(\theta \right)}\left(t^{\ast }\right)}}$ for $n\in \mathbb{N}$.

Using (4) and (5), we have

(6)$\begin{array}{cc}\hfill e_{u_i}=& u_i-\alpha \\ \hfill =& {\displaystyle \frac{D_1}{\theta }}{e}_i^2+{\displaystyle \frac{2\theta D_2-(1+\theta )D_1^2}{{\theta }^2}}{e}_i^3\\ & +{\displaystyle \frac{1}{{\theta }^3}}\left({(1+\theta )}^2{D}_1^3-\theta (4+3\theta )D_1{D}_2+3{\theta }^2{D}_3\right)e_i^4+O\left(e_i^5\right).\end{array}$

Expanding $f^{\prime }\left(u_i\right)$ about $t^{\ast }$ gives

(7)$f^{\prime }\left(u_i\right)={\displaystyle \frac{f^{\left(\theta \right)}\left(t^{\ast }\right)}{\theta !}}{e}_{u_i}^{\theta -1}\left(\theta +(\theta +1)D_1{e}_{u_i}+(\theta +2)D_2{e}_{u_i}^2+(\theta +3)D_3{e}_{u_i}^3+(\theta +4)D_4{e}_{u_i}^4+\cdots \right).$

Then, we obtain that

(8)$\begin{array}{c}\hfill x_i={\left({\displaystyle \frac{f^{\prime }\left(u_i\right)}{f^{\prime }\left(t_i\right)}}\right)}^{{\displaystyle \frac{1}{\theta -1}}}={\displaystyle \frac{D_1}{\theta }}{e}_i+{\displaystyle \frac{2(\theta -1)D_2-(\theta +1)D_1^2}{\theta (\theta -1)}}{e}_i^2+{\eta }_1{e}_i^3+{\eta }_2{e}_i^4+O\left(e_i^5\right),\end{array}$

where ${\eta }_1={\eta }_1(\theta ,D_1,D_2,D_3)$ and ${\eta }_2={\eta }_2(\theta ,D_1,D_2,D_3,D_4)$.

Expanding function $G\left(h_i\right)$ in the neighborhood of the origin, we have that

(9)$G\left(h_i\right)\approx G\left(0\right)+h_i{G}^{\prime }\left(0\right)+{\displaystyle \frac{1}{2}}{h}_i^2{G}^{\prime \prime }\left(0\right)+{\displaystyle \frac{1}{6}}{h}_i^3{G}^{\prime \prime \prime }\left(0\right)+O\left(h_i^4\right),$

where $h_i={\displaystyle \frac{x_i}{a_1+a_2{x}_i}}$.

By using Equations (4)–(6), (8) and (9) in Scheme (3), we have

(10)$\begin{array}{cc}\hfill e_{i+1}=& -G\left(0\right)e_i+{\displaystyle \frac{(a_1+a_{1}G\left(0\right)-G^{\prime }\left(0\right))}{a_1\theta }}{D}_1{e}_i^2+{\displaystyle \frac{1}{2a_1^2(\theta -1){\theta }^2}}(\left(\right(2a_2{G}^{\prime }\left(0\right)\hfill \\ \hfill & -G^{\prime \prime }\left(0\right))(\theta -1)-2a_1^2(1+G\left(0\right))({\theta }^2-1)+2a_1{G}^{\prime }\left(0\right)({\theta }^2+2\theta -1))D_1^2\hfill \\ \hfill & +4a_1(a_1+a_{1}G\left(0\right)-G^{\prime }\left(0\right))(\theta -1)\theta D_2)e_i^3+\gamma e_i^4+O\left(e_i^5\right),\hfill \end{array}$

where

$\begin{array}{cc}\hfill \gamma =& {\displaystyle \frac{1}{6a_1^3{(\theta -1)}^2{\theta }^3}}\left(\right(-(6a_2^2{G}^{\prime }\left(0\right)-6a_2{G}^{\prime \prime }\left(0\right)+G^{\prime \prime \prime }\left(0\right)){(\theta -1)}^2+6a_1^3(1+G\left(0\right)){({\theta }^2-1)}^2\hfill \\ \hfill & -3a_1(2a_2{G}^{\prime }\left(0\right)-G^{\prime \prime }\left(0\right))(1-4\theta +{\theta }^2+2{\theta }^3)-3a_1^2{G}^{\prime }\left(0\right)\theta (-5+7{\theta }^2+2{\theta }^3))D_1^3\hfill \\ \hfill & -6a_1(\theta -1)\theta (-2(2a_2{G}^{\prime }\left(0\right)-G^{\prime \prime }\left(0\right))(\theta -1)+a_1{G}^{\prime }\left(0\right)(4-8\theta -3{\theta }^2)\hfill \\ \hfill & +a_1^2(1+G\left(0\right))(-4+\theta +3{\theta }^2))D_1{D}_2+18a_1^2(a_1+a_{1}G\left(0\right)-G^{\prime }\left(0\right)){(\theta -1)}^2{\theta }^2{D}_3).\hfill \end{array}$

The vanishing of the coefficients of $e_i$, $e_i^2$, and $e_i^3$ in Equation (10) is clear. We thus have the optimal fourth-order convergence. So, after simple calculations, we have the following:

(11)$G\left(0\right)=0,G^{\prime }\left(0\right)=a_1\mathrm{and}{G}^{\prime \prime }\left(0\right)={\displaystyle \frac{2(2a_1^2\theta +a_1{a}_2\theta -a_1{a}_2)}{\theta -1}},\theta \ne 1.$

Thus, the error Equation (10) is given by

(12)$\begin{array}{cc}\hfill e_{i+1}=& {\displaystyle \frac{D_1}{6a_1^3{(\theta -1)}^2{\theta }^3}}\left(\right(6a_1{a}_2^2{(\theta -1)}^2-G^{\prime \prime \prime }\left(0\right){(\theta -1)}^2+24a_1^2{a}_2(\theta -1)\theta \\ & +3a_1^3(2+\theta +8{\theta }^2+{\theta }^3))D_1^2-6a_1^3(\theta -1){\theta }^2{D}_2)e_i^4+O\left(e_i^5\right).\end{array}$

The following theorem states the above results:

Theorem 1.

Let $f:\mathbb{C}\to \mathbb{C}$ represent an analytical function in the neighborhood of a multiple-zero $t^{\ast }$ with multiplicity $\theta >1$. If the initial guess $t_0$ is sufficiently near to $t^{\ast }$, then Scheme (3) has a local order of convergence that is at least four, if $G\left(0\right)=0$, $G^{\prime }\left(0\right)=a_1$, $G^{″}\left(0\right)={\displaystyle \frac{2(2a_1^2\theta +a_1{a}_2\theta -a_1{a}_2)}{\theta -1}}$, and $|G^{‴}\left(0\right)|<\infty$.

Special Members of Scheme (3)

By distributing various values of weight functions $G\left(h_i\right)$ that satisfy (11), we discuss a few particular examples of our suggested Scheme (3) in this section. Thus, here we have specified various members of the suggested family in this regard. The corresponding simple forms of $G\left(h_i\right)$ are given by

(13)$\begin{array}{cc}\hfill G\left(h_i\right)=& {\displaystyle \frac{a_1{h}_i(a_2{h}_i(\theta -1)+\theta +2a_1{h}_i\theta -1)}{\theta -1}},\end{array}$

(14)$\begin{array}{cc}\hfill G\left(h_i\right)=& {\displaystyle \frac{a_1{h}_i(\theta -1)}{\theta +h_i(a_2-2a_1\theta -a_2\theta )-1}}.\end{array}$

Based on the values of parameters $a_1$ and $a_2$, we present the following special members of the family in (3):

• $\left(1\right)$

  Combining $a_1=1$, $a_2=0$, and (13) in Expression (3), we have

  (15)$\begin{array}{cc}\hfill t_{i+1}=& u_i-\theta \left(h_i+{\displaystyle \frac{2\theta h_i^2}{\theta -1}}\right){\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}}.\end{array}$

  It is very important to remember that method (15) above is a Liu–Zhou Method [18]. This demonstrates that the Liu–Zhou method [18] is a special case of our family, which is given in (3).

• $\left(2\right)$

  Combining $a_1=1$, $a_2=0$, and (14) in Expression (3), we have

  (16)$\begin{array}{cc}\hfill t_{i+1}=& u_i-\theta {\displaystyle \frac{h_i(\theta -1)}{\theta -2h_i\theta -1}}{\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}}.\end{array}$

  Again, it is important to notice that method (16) above is a Liu–Zhou method [18]. This demonstrates that Liu–Zhou methods [18] are a special case of our family (3).

• $\left(3\right)$

  Using $a_1=1$, $a_2=1$, and (13) in (3), we have

  (17)$\begin{array}{cc}\hfill t_{i+1}=& u_i-\theta \left(h_i+{\displaystyle \frac{(3\theta -1)h_i^2}{\theta -1}}\right){\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}}.\end{array}$

• $\left(4\right)$

  Use $a_1=-1$, $a_2=1$, and (13) in (3), we obtain

  (18)$\begin{array}{cc}\hfill t_{i+1}=& u_i-\theta \left(-h_i+{\displaystyle \frac{(1+\theta )h_i^2}{\theta -1}}\right){\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}}.\end{array}$

• $\left(5\right)$

  Let $a_1=1$, $a_2=-2$, and (13) in (3), we have

  (19)$\begin{array}{cc}\hfill t_{i+1}=& u_i-\theta \left(h_i+{\displaystyle \frac{2h_i^2}{\theta -1}}\right){\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}}.\end{array}$

• $\left(6\right)$

  Let $a_1=1$, $a_2=-3$, and (13) in (3), we have

  (20)$\begin{array}{cc}\hfill t_{i+1}=& u_i-\theta \left(h_i-{\displaystyle \frac{(\theta -3)h_i^2}{\theta -1}}\right){\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}}.\end{array}$

In each of the above cases, $u_i=t_i-\theta {\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}}$. For numerical work in the following, the proposed methods (15)–(20) are denoted by LZ-1, LZ-2, M-1, M-2, M-3, and M-4, respectively.

Remark 1.

The new family (3) only requires three functional evaluations (viz., $f\left(t_i\right)$, $f^{\prime }\left(t_i\right)$, and $f^{\prime }\left(u_i\right)$) per iteration to reach fourth-order convergence. According to the Kung and Traub [7] hypothesis, the approaches have optimal fourth-order convergence.

4. Numerical Simulation

This section applies the new methods LZ-1, LZ-2, M-1, M-2, M-3, and M-4 to a few basic science problems and illustrates their convergence behavior and computational effectiveness. Their performance is also contrasted with current approaches. For instance, we chose Li et al. [16,17], Sharma–Sharma [24], Zhou et al. [30], Soleymani et al. [26], and Kansal et al. [15]. Now, these methods are expressed as follows:

Li–Liao–Cheng method (LLC):

$\begin{array}{cc}\hfill u_i=& t_i-{\displaystyle \frac{2\theta }{\theta +2}}{\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}},\\ \hfill t_{i+1}=& t_i-{\displaystyle \frac{\theta (\theta -2){\left(\frac{\theta }{\theta +2}\right)}^{-\theta }{f}^{\prime }\left(u_i\right)-{\theta }^2{f}^{\prime }\left(t_i\right)}{f^{\prime }\left(t_i\right)-{\left(\frac{\theta }{\theta +2}\right)}^{-\theta }{f}^{\prime }\left(u_i\right)}}{\displaystyle \frac{f\left(t_i\right)}{2f^{\prime }\left(t_i\right)}}.\end{array}$

Li–Cheng–Neta method (LCN):

$\begin{array}{cc}\hfill u_i=& t_i-{\displaystyle \frac{2\theta }{\theta +2}}{\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}},\\ \hfill t_{i+1}=& t_i-{\alpha }_1{\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(u_i\right)}}-{\displaystyle \frac{f\left(t_i\right)}{{\alpha }_2{f}^{\prime }\left(t_i\right)+{\alpha }_3{f}^{\prime }\left(u_i\right)}},\end{array}$

where

$\begin{array}{cc}\hfill {\alpha }_1=& -{\displaystyle \frac{1}{2}}{\displaystyle \frac{{\left(\frac{\theta }{\theta +2}\right)}^{\theta }\theta ({\theta }^4+4{\theta }^3-16\theta -16)}{{\theta }^3-4\theta +8}},\\ \hfill {\alpha }_2=& -{\displaystyle \frac{{({\theta }^3-4\theta +8)}^2}{\theta ({\theta }^4+4{\theta }^3-4{\theta }^2-16\theta +16)({\theta }^2+2\theta -4)}},\\ \hfill {\alpha }_3=& {\displaystyle \frac{{\theta }^2({\theta }^3-4\theta +8)}{{\left(\frac{\theta }{\theta +2}\right)}^{\theta }({\theta }^4+4{\theta }^3-4{\theta }^2-16\theta +16)({\theta }^2+2\theta -4)}}.\end{array}$

Sharma–Sharma method (SM):

$\begin{array}{cc}\hfill u_i=& t_i-{\displaystyle \frac{2\theta }{\theta +2}}{\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}},\\ \hfill t_{i+1}=& t_i-{\displaystyle \frac{\theta }{8}}[({\theta }^3-4\theta +8)-{(\theta +2)}^2{\left({\displaystyle \frac{\theta }{\theta +2}}\right)}^{\theta }{\displaystyle \frac{f^{\prime }\left(t_i\right)}{f^{\prime }\left(u_i\right)}}(2(\theta -1)\\ & -(\theta +2){\left({\displaystyle \frac{\theta }{\theta +2}}\right)}^{\theta }{\displaystyle \frac{f^{\prime }\left(t_i\right)}{f^{\prime }\left(u_i\right)}}\left)\right]{\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}}.\end{array}$

Zhou–Chen–Song method (ZCS):

$\begin{array}{cc}\hfill u_i=& t_i-{\displaystyle \frac{2\theta }{\theta +2}}{\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}},\\ \hfill t_{i+1}=& t_i-{\displaystyle \frac{\theta }{8}}[{\theta }^3{\left({\displaystyle \frac{\theta +2}{\theta }}\right)}^{2\theta }{\left({\displaystyle \frac{f^{\prime }\left(u_i\right)}{f^{\prime }\left(t_i\right)}}\right)}^2-2{\theta }^2(\theta +3){\left({\displaystyle \frac{\theta +2}{\theta }}\right)}^{\theta }{\displaystyle \frac{f^{\prime }\left(u_i\right)}{f^{\prime }\left(t_i\right)}}\\ & +({\theta }^3+6{\theta }^2+8\theta +8)]{\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}}.\end{array}$

Soleymani–Babajee–Lotfi method (SBL):

$\begin{array}{cc}\hfill u_i=& t_i-{\displaystyle \frac{2\theta }{\theta +2}}{\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}},\\ \hfill t_{i+1}=& t_i-{\displaystyle \frac{f^{\prime }\left(u_i\right)f\left(t_i\right)}{q_1{\left(f^{\prime }\left(u_i\right)\right)}^2+q_2{f}^{\prime }\left(u_i\right)f^{\prime }\left(t_i\right)+q_3{\left(f^{\prime }\left(t_i\right)\right)}^2}},\end{array}$

where

$\begin{array}{cc}\hfill q_1=& {\displaystyle \frac{1}{16}}{\theta }^{3-\theta }{(\theta +2)}^{\theta },\\ \hfill q_2=& {\displaystyle \frac{8-\theta (\theta +2)({\theta }^2-2)}{8\theta }},\\ \hfill q_3=& {\displaystyle \frac{1}{16}}(\theta -2){\theta }^{\theta -1}{(\theta +2)}^{3-\theta }.\end{array}$

Kansal–Kanwar–Bhatia method (KKB):

$\begin{array}{cc}\hfill u_i=& t_i-{\displaystyle \frac{2\theta }{\theta +2}}{\displaystyle \frac{f\left(t_i\right)}{f^{\prime }\left(t_i\right)}},\\ \hfill u_{i+1}=& t_i-{\displaystyle \frac{\theta }{4}}f\left(t_i\right)\left(1+{\displaystyle \frac{{\theta }^4{p}^{-2\theta }{\left(p^{\theta -1}-\frac{f^{\prime }\left(u_i\right)}{f^{\prime }\left(t_i\right)}\right)}^2(p^{\theta }-1)}{8(2p^{\theta }+\theta (p^{\theta }-1))}}\right)\\ & \times \left({\displaystyle \frac{4-2\theta +{\theta }^2(p^{-\theta }-1)}{f^{\prime }\left(t_i\right)}}-{\displaystyle \frac{p^{-\theta }{(2p^{\theta }+\theta (p^{\theta }-1))}^2}{f^{\prime }\left(t_i\right)-f^{\prime }\left(u_i\right)}}\right),\end{array}$

where $p={\displaystyle \frac{\theta }{\theta +2}}.$

The various problems considered for numerical testing are shown in Table 1. Computations were compiled in the programming package of the software Mathematica using multiple-precision arithmetic. Numerical results displayed in Table 2, Table 3, Table 4, Table 5 and Table 6 include the following: (i) number of iterations $\left(i\right)$ required to obtain the desired solution using the stopping criterion $|t_{i+1}-t_i|+|f\left(t_i\right)|<{10}^{-350}$, (ii) estimated errors $|t_{i+1}-t_i|$ for the first three iterations, and (iii) computational order of convergence (COC). Table 7 displays the CPU time utilized in the execution of a programm which is computed by the Mathematica command “TimeUsed[ ]”. The computational order of convergence (COC) is calculated using Formula (2).

We note that the increased accuracy of the proposed approaches exhibits increasing precision in the successive approximations based on the numerical data shown in Table 2, Table 3, Table 4, Table 5 and Table 6. This explains the method’s excellent convergence nature. The theoretical fourth-order convergence of the new methods is strongly supported by the computational order of convergence shown in the penultimate columns of the tables. In addition, the CPU time taken by the techniques as displayed in Table 7 demonstrates the computationally efficient nature of the new technique as compared to the CPU time of the considered existing techniques of the same order. We also show the time required by each method in bar graphs. Figure 1a–e show the graphical representation of data in Table 7. Similar numerical testing carried out for many other problems has confirmed the above conclusions to a large extent.

5. Conclusions

We have proposed a fourth-order family of iterative algorithms that are computationally effective for identifying multiple roots in applied science problems. According to the Kung–Traub conjecture, the approaches converge to the required root with fourth-order convergence and with three function evaluations per iteration, so the new family is optimal. The procedure is unique in the sense that there is no such algorithm available in the literature. The main benefit of the new family is that the current Liu–Zhou approach is a special case of the new family. Analysis of the convergence was carried out, which proves fourth-order convergence under standard assumptions of the function whose zeros we are looking for. Numerical testing was checked to evaluate performance. Additionally, the new methods were also applied to many other problems, further supporting the effectiveness of the new methods.

Author Contributions

Conceptualization, methodology, S.K. and A.K.; software, writing, M.K. and M.V.; draft preparation, P.D.; formal analysis, validation, resources, L.J. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

The original contributions presented in this study are included in the article. Further inquiries can be directed to the corresponding author.

Acknowledgments

We would like to express our gratitude to the anonymous reviewers for their help with the publication of this paper.

Conflicts of Interest

The authors declare no conflicts of interest.

Footnotes

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.

Figure and Tables

Figure 1. (a) Bar diagram for function f1(t). (b) Bar diagram for function f2(t). (c) Bar diagram for function f3(t). (d) Bar diagram for function f4(t). (e) Bar diagram for function f5(t).

Table 1

Test functions.

Functions	Root (${\mathit{t}}^{\ast }$)	$\mathit{\theta }$	${\mathit{t}}_0$
Continuous stirred tank reactor problem [33]
$f_1\left(t\right)=t^4+11.50t^3+47.49t^2+83.06325t+51.23266875$	−2.85	2	−3
Standard test problem
$f_2\left(t\right)=-{\displaystyle \frac{t^4}{12}}+{\displaystyle \frac{t^2}{2}}+t+e^t(t-3)+sin\left(t\right)+3$	0	3	0.6
Isentropic supersonic flow problem [34]
$f_3\left(t\right)=[{tan}^{-1}\left({\displaystyle \frac{\sqrt{5}}{2}}\right)-{tan}^{-1}\left(\sqrt{t^2-1}\right)+\sqrt{6}({tan}^{-1}\left(\sqrt{{\displaystyle \frac{t^2-1}{6}}}\right)$
$-{tan}^{-1}\left({\displaystyle \frac{1}{2}}\sqrt{{\displaystyle \frac{5}{6}}}\right))-{\displaystyle \frac{11}{63}}{]}^3$	1.8411294068…	3	1.5
Eigen value problem [10]
$f_4\left(t\right)=(t-8)(t-5)(t-4){(t-3)}^4(t^2-1)$	3	4	2.3
Complex root problem [35]
$f_5\left(t\right)=t(t^2+1)(2e^{t^2+1}+t^2-1){cosh}^4\left({\displaystyle \frac{\pi t}{2}}\right)$	i	6	1.1 i

Table 2

Performance of methods for function $f_1\left(t\right)$.

Methods	i	$|{\mathit{t}}_2-{\mathit{t}}_1|$	$|{\mathit{t}}_3-{\mathit{t}}_2|$	$|{\mathit{t}}_4-{\mathit{t}}_3|$	COC
LLC	4	$1.18\times {10}^{-5}$	$2.22\times {10}^{-22}$	$2.75\times {10}^{-89}$	4.000
LCN	4	$1.18\times {10}^{-5}$	$2.22\times {10}^{-22}$	$2.75\times {10}^{-89}$	4.000
SM	4	$1.18\times {10}^{-5}$	$2.23\times {10}^{-22}$	$2.80\times {10}^{-89}$	4.000
ZCS	4	$1.19\times {10}^{-5}$	$2.26\times {10}^{-22}$	$2.98\times {10}^{-89}$	4.000
SBL	4	$1.18\times {10}^{-5}$	$2.22\times {10}^{-22}$	$2.75\times {10}^{-89}$	4.000
KKB	4	$1.18\times {10}^{-5}$	$2.18\times {10}^{-22}$	$2.59\times {10}^{-89}$	4.000
LZ-1	4	$1.26\times {10}^{-5}$	$2.97\times {10}^{-22}$	$9.04\times {10}^{-89}$	4.000
LZ-2	4	$1.17\times {10}^{-5}$	$2.16\times {10}^{-22}$	$2.48\times {10}^{-89}$	4.000
M-1	4	$1.32\times {10}^{-5}$	$3.54\times {10}^{-22}$	$1.85\times {10}^{-89}$	4.000
M-2	4	$1.22\times {10}^{-5}$	$2.58\times {10}^{-22}$	$5.14\times {10}^{-89}$	4.000
M-3	4	$1.19\times {10}^{-5}$	$2.34\times {10}^{-22}$	$3.43\times {10}^{-89}$	4.000
M-4	4	$1.18\times {10}^{-5}$	$2.20\times {10}^{-22}$	$2.69\times {10}^{-89}$	4.000

Table 3

Performance of methods for function $f_2\left(t\right)$.

Methods	i	$|{\mathit{t}}_2-{\mathit{t}}_1|$	$|{\mathit{t}}_3-{\mathit{t}}_2|$	$|{\mathit{t}}_4-{\mathit{t}}_3|$	COC
LLC	5	$2.02\times {10}^{-4}$	$2.11\times {10}^{-17}$	$2.51\times {10}^{-69}$	4.000
LCN	5	$2.02\times {10}^{-4}$	$2.12\times {10}^{-17}$	$2.54\times {10}^{-69}$	4.000
SM	5	$2.02\times {10}^{-4}$	$2.12\times {10}^{-17}$	$2.60\times {10}^{-69}$	4.000
ZCS	5	$2.02\times {10}^{-4}$	$2.15\times {10}^{-17}$	$2.75\times {10}^{-69}$	4.000
SBL	5	$2.02\times {10}^{-4}$	$2.13\times {10}^{-17}$	$2.62\times {10}^{-69}$	4.000
KKB	5	$2.02\times {10}^{-4}$	$2.08\times {10}^{-17}$	$2.31\times {10}^{-69}$	4.000
LZ-1	5	$1.39\times {10}^{-4}$	$5.16\times {10}^{-18}$	$9.75\times {10}^{-72}$	4.000
LZ-2	5	$1.37\times {10}^{-4}$	$3.05\times {10}^{-18}$	$7.40\times {10}^{-73}$	4.000
M-1	5	$1.40\times {10}^{-4}$	$6.93\times {10}^{-18}$	$4.12\times {10}^{-71}$	4.000
M-2	5	$1.38\times {10}^{-4}$	$3.96\times {10}^{-18}$	$2.68\times {10}^{-72}$	4.000
M-3	5	$1.38\times {10}^{-4}$	$3.27\times {10}^{-18}$	$1.05\times {10}^{-72}$	4.000
M-4	5	$1.37\times {10}^{-4}$	$3.05\times {10}^{-18}$	$7.40\times {10}^{-73}$	4.000

Table 4

Performance of methods for function $f_3\left(t\right)$.

Methods	i	$|{\mathit{t}}_2-{\mathit{t}}_1|$	$|{\mathit{t}}_3-{\mathit{t}}_2|$	$|{\mathit{t}}_4-{\mathit{t}}_3|$	COC
LLC	5	$8.19\times {10}^{-4}$	$2.64\times {10}^{-15}$	$2.87\times {10}^{-61}$	4.000
LCN	5	$8.19\times {10}^{-4}$	$2.61\times {10}^{-15}$	$2.73\times {10}^{-61}$	4.000
SM	5	$8.19\times {10}^{-4}$	$2.55\times {10}^{-15}$	$2.44\times {10}^{-61}$	4.000
ZCS	5	$8.19\times {10}^{-4}$	$2.40\times {10}^{-15}$	$1.79\times {10}^{-61}$	4.000
SBL	5	$8.19\times {10}^{-4}$	$2.53\times {10}^{-15}$	$2.34\times {10}^{-61}$	4.000
KKB	5	$8.19\times {10}^{-4}$	$2.85\times {10}^{-15}$	$4.22\times {10}^{-61}$	4.000
LZ-1	5	$3.99\times {10}^{-5}$	$6.65\times {10}^{-20}$	$5.13\times {10}^{-79}$	4.000
LZ-2	5	$3.97\times {10}^{-5}$	$3.86\times {10}^{-20}$	$3.44\times {10}^{-80}$	4.000
M-1	5	$4.01\times {10}^{-5}$	$8.91\times {10}^{-20}$	$2.17\times {10}^{-78}$	4.000
M-2	5	$3.98\times {10}^{-5}$	$5.09\times {10}^{-20}$	$1.36\times {10}^{-79}$	4.000
M-3	5	$3.97\times {10}^{-5}$	$4.16\times {10}^{-20}$	$5.02\times {10}^{-80}$	4.000
M-4	5	$3.97\times {10}^{-5}$	$3.86\times {10}^{-20}$	$3.44\times {10}^{-80}$	4.000

Table 5

Performance of methods for function $f_4\left(t\right)$.

Methods	i	$|{\mathit{t}}_2-{\mathit{t}}_1|$	$|{\mathit{t}}_3-{\mathit{t}}_2|$	$|{\mathit{t}}_4-{\mathit{t}}_3|$	COC
LLC	5	$1.42\times {10}^{-3}$	$6.70\times {10}^{-13}$	$3.33\times {10}^{-50}$	4.000
LCN	5	$1.42\times {10}^{-3}$	$6.72\times {10}^{-13}$	$3.37\times {10}^{-50}$	4.000
SM	5	$1.42\times {10}^{-3}$	$6.76\times {10}^{-13}$	$3.48\times {10}^{-50}$	4.000
ZCS	5	$1.42\times {10}^{-3}$	$6.83\times {10}^{-13}$	$3.67\times {10}^{-50}$	4.000
SBL	5	$1.42\times {10}^{-3}$	$6.87\times {10}^{-13}$	$3.78\times {10}^{-50}$	4.000
KKB	5	$1.42\times {10}^{-3}$	$6.53\times {10}^{-13}$	$2.93\times {10}^{-50}$	4.000
LZ-1	5	$7.07\times {10}^{-4}$	$4.35\times {10}^{-14}$	$6.26\times {10}^{-55}$	4.000
LZ-2	5	$7.02\times {10}^{-4}$	$1.94\times {10}^{-14}$	$1.12\times {10}^{-56}$	4.000
M-1	5	$7.10\times {10}^{-4}$	$6.59\times {10}^{-14}$	$4.90\times {10}^{-54}$	4.000
M-2	5	$7.04\times {10}^{-4}$	$2.86\times {10}^{-14}$	$7.88\times {10}^{-56}$	4.000
M-3	5	$7.02\times {10}^{-4}$	$2.08\times {10}^{-14}$	$1.62\times {10}^{-56}$	4.000
M-4	5	$7.02\times {10}^{-4}$	$1.97\times {10}^{-14}$	$1.23\times {10}^{-56}$	4.000

Table 6

Performance of methods for function $f_5\left(t\right)$.

Methods	i	$|{\mathit{t}}_2-{\mathit{t}}_1|$	$|{\mathit{t}}_3-{\mathit{t}}_2|$	$|{\mathit{t}}_4-{\mathit{t}}_3|$	COC
LLC	5	$1.75\times {10}^{-5}$	$3.01\times {10}^{-20}$	$2.66\times {10}^{-79}$	4.000
LCN	5	$1.75\times {10}^{-5}$	$3.02\times {10}^{-20}$	$2.68\times {10}^{-79}$	4.000
SM	5	$1.75\times {10}^{-5}$	$3.03\times {10}^{-20}$	$2.73\times {10}^{-79}$	4.000
ZCS	5	$1.75\times {10}^{-5}$	$3.05\times {10}^{-20}$	$2.79\times {10}^{-79}$	4.000
SBL	5	$1.76\times {10}^{-5}$	$3.16\times {10}^{-20}$	$3.28\times {10}^{-79}$	4.000
KKB	5	$1.74\times {10}^{-5}$	$2.94\times {10}^{-20}$	$2.40\times {10}^{-79}$	4.000
LZ-1	5	$8.39\times {10}^{-6}$	$9.14\times {10}^{-22}$	$1.29\times {10}^{-85}$	4.000
LZ-2	5	$6.91\times {10}^{-6}$	$2.76\times {10}^{-22}$	$7.01\times {10}^{-88}$	4.000
M-1	5	$9.81\times {10}^{-6}$	$2.29\times {10}^{-21}$	$6.85\times {10}^{-84}$	4.000
M-2	5	$7.43\times {10}^{-6}$	$4.34\times {10}^{-22}$	$5.06\times {10}^{-87}$	4.000
M-3	5	$6.95\times {10}^{-6}$	$2.87\times {10}^{-22}$	$8.33\times {10}^{-88}$	4.000
M-4	5	$7.01\times {10}^{-6}$	$3.02\times {10}^{-22}$	$1.04\times {10}^{-87}$	4.000

Table 7

CPU time consumed by methods.

Methods	${\mathit{f}}_1\left(\mathit{t}\right)$	${\mathit{f}}_2\left(\mathit{t}\right)$	${\mathit{f}}_3\left(\mathit{t}\right)$	${\mathit{f}}_4\left(\mathit{t}\right)$	${\mathit{f}}_5\left(\mathit{t}\right)$
LLC	0.1562	0.8491	1.8570	0.4217	1.4506
LCN	0.1557	0.8897	2.0283	0.4534	2.2164
SM	0.1884	0.8433	2.0442	0.4681	2.1681
ZCS	0.1714	0.8265	2.0129	0.4575	2.2165
SBL	0.1873	1.0342	2.3187	0.4567	2.4653
KKB	0.1723	0.8587	2.0284	0.4052	2.0287
LZ-1	0.0392	0.6867	1.5138	0.2029	1.2015
LZ-2	0.0399	0.7221	1.6533	0.2347	1.3261
M-1	0.0648	0.7024	1.6592	0.2184	1.3326
M-2	0.0637	0.6870	1.6387	0.2392	1.3427
M-3	0.0461	0.7107	1.6534	0.2501	1.3419
M-4	0.0624	0.7231	1.6239	0.2562	1.3426

References

1. Azarmanesh, M.; Dejam, M.; Azizian, P.; Yesiloz, G.; Mohamad, A.A.; Sanati-Nezhad, A. Passive microinjection within high-throughput microfluidics for controlled actuation of droplets and cells. Sci. Rep.; 2019; 9, 6723. [DOI: https://dx.doi.org/10.1038/s41598-019-43056-2]

2. Dejam, M. Advective-diffusive-reactive solute transport due to non-Newtonian fluid flows in a fracture surrounded by a tight porous medium. Int. J. Heat Mass Transf.; 2019; 128, pp. 1307-1321. [DOI: https://dx.doi.org/10.1016/j.ijheatmasstransfer.2018.09.061]

3. Nikpoor, M.H.; Dejam, M.; Chen, Z.; Clarke, M. Chemical-Gravity-Thermal Diffusion Equilibrium in Two-Phase Non-isothermal Petroleum Reservoirs. Energy Fuel; 2016; 30, pp. 2021-2034. [DOI: https://dx.doi.org/10.1021/acs.energyfuels.5b02753]

4. Ostrowski, A.M. Solution of Equations and Systems of Equations; Academic Press: New York, NY, USA, 1960.

5. Traub, J.F. Iterative Methods for the Solution of Equations; Prentice-Hall: Englewood Cliffs, NJ, USA, 1964.

6. Petkovic, M.S.; Neta, B.; Petkovic, L.D.; Dzunic, J. Multipoint Methods for Solving Nonlinear Equations; Academic Press: New York, NY, USA, 2013.

7. Kung, H.T.; Traub, J.F. Optimal order of one-point and multipoint iteration. J. Assoc. Comput. Mach.; 1974; 21, pp. 643-651. [DOI: https://dx.doi.org/10.1145/321850.321860]

8. Behl, R.; Cordero, A.; Motsa, S.S.; Torregrosa, J.R. On developing fourth-order optimal families of methods for multiple roots and their dynamics. Appl. Math. Comput.; 2015; 265, pp. 520-532. [DOI: https://dx.doi.org/10.1016/j.amc.2015.05.004]

9. Behl, R.; Cordero, A.; Motsa, S.S.; Torregrosa, J.R.; Kanwar, V. An optimal fourth-order family of methods for multiple roots and its dynamics. Numer. Algor.; 2016; 71, pp. 775-796. [DOI: https://dx.doi.org/10.1007/s11075-015-0023-5]

10. Behl, R.; Zafar, F.; Alshormani, A.S.; Junjua, M.U.D.; Yasmin, N. An optimal eighth-order scheme for multiple zeros of unvariate functions. Int. J. Comput. Meth.; 2018; 16, 1843002. [DOI: https://dx.doi.org/10.1142/S0219876218430028]

11. Dong, C. A family of multipoint iterative functions for finding multiple roots of equations. Int. J. Comput. Math.; 1987; 21, pp. 363-367. [DOI: https://dx.doi.org/10.1080/00207168708803576]

12. Geum, Y.H.; Kim, Y.I.; Neta, B. A class of two-point sixth-order multiple-zero finders of modified double-Newton type and their dynamics. Appl. Math. Comput.; 2015; 270, pp. 387-400. [DOI: https://dx.doi.org/10.1016/j.amc.2015.08.039]

13. Hansen, E.; Patrick, M. A family of root finding methods. Numer. Math.; 1977; 27, pp. 257-269. [DOI: https://dx.doi.org/10.1007/BF01396176]

14. Hueso, J.L.; Martínez, E.; Teruel, C. Determination of multiple roots of nonlinear equations and applications. J. Math. Chem.; 2015; 53, pp. 880-892. [DOI: https://dx.doi.org/10.1007/s10910-014-0460-8]

15. Kansal, M.; Kanwar, V.; Bhatia, S. On some optimal multiple root-finding methods and their dynamics. Appl. Appl. Math.; 2015; 10, pp. 349-367.

16. Li, S.; Liao, X.; Cheng, L. A new fourth-order iterative method for finding multiple roots of nonlinear equations. Appl. Math. Comput.; 2009; 215, pp. 1288-1292.

17. Li, S.G.; Cheng, L.Z.; Neta, B. Some fourth-order nonlinear solvers with closed formulae for multiple roots. Comput Math. Appl.; 2010; 59, pp. 126-135. [DOI: https://dx.doi.org/10.1016/j.camwa.2009.08.066]

18. Liu, B.; Zhou, X. A new family of fourth-order methods for multiple roots of nonlinear equations. Nonlinear Anal. Model. Control; 2013; 21, pp. 143-152. [DOI: https://dx.doi.org/10.15388/NA.18.2.14018]

19. Neta, B. New third order nonlinear solvers for multiple roots. Appl. Math. Comput.; 2008; 202, pp. 162-170. [DOI: https://dx.doi.org/10.1016/j.amc.2008.01.031]

20. Neta, B. Extension of Murakami’s high-order non-linear solver to multiple roots. Int. J. Comput. Math.; 2010; 87, pp. 1023-1031. [DOI: https://dx.doi.org/10.1080/00207160802272263]

21. Osada, N. An optimal multiple root-finding method of order three. J. Comput. Appl. Math.; 1994; 51, pp. 131-133. [DOI: https://dx.doi.org/10.1016/0377-0427(94)00044-1]

22. Sharma, J.R.; Kumar, S. An excellent numerical technique for multiple roots. Math. Comput. Simul.; 2021; 182, pp. 316-324. [DOI: https://dx.doi.org/10.1016/j.matcom.2020.11.008]

23. Sharma, J.R.; Kumar, S. A class of computationally efficient numerical algorithms for locating multiple zeros. Afr. Mat.; 2021; 32, pp. 853-864. [DOI: https://dx.doi.org/10.1007/s13370-020-00865-3]

24. Sharma, J.R.; Sharma, R. Modified Jarratt method for computing multiple roots. Appl. Math. Comput.; 2010; 217, pp. 878-881. [DOI: https://dx.doi.org/10.1016/j.amc.2010.06.031]

25. Sharifi, M.; Babajee, D.K.R.; Soleymani, F. Finding the solution of nonlinear equations by a class of optimal methods. Comput. Math. Appl.; 2012; 63, pp. 764-774. [DOI: https://dx.doi.org/10.1016/j.camwa.2011.11.040]

26. Soleymani, F.; Babajee, D.K.R.; Lotfi, T. On a numerical technique for finding multiple zeros and its dynamics. J. Egypt. Math. Soc.; 2013; 21, pp. 346-353. [DOI: https://dx.doi.org/10.1016/j.joems.2013.03.011]

27. Soleymani, F.; Babajee, D.K.R. Computing multiple zeros using a class of quartically convergent methods. Alex. Eng. J.; 2013; 52, pp. 531-541. [DOI: https://dx.doi.org/10.1016/j.aej.2013.05.001]

28. Thukral, R. A new family of fourth-order iterative methods for solving nonlinear equations with multiple roots. J. Numer. Math. Stoch.; 2014; 6, pp. 37-44.

29. Victory, H.D.; Neta, B. A higher order method for multiple zeros of nonlinear functions. Int. J. Comput. Math.; 1983; 12, pp. 329-335. [DOI: https://dx.doi.org/10.1080/00207168208803346]

30. Zhou, X.; Chen, X.; Song, Y. Constructing higher-order methods for obtaining the multiple roots of nonlinear equations. J. Comput. Appl. Math.; 2011; 235, pp. 4199-4206. [DOI: https://dx.doi.org/10.1016/j.cam.2011.03.014]

31. Zhou, X.; Chen, X.; Song, Y. Families of third and fourth order methods for multiple roots of nonlinear equations. Appl. Math. Comput.; 2013; 219, pp. 6030-6038. [DOI: https://dx.doi.org/10.1016/j.amc.2012.12.041]

32. Cordero, A.; Torregrosa, J.R. Variants of Newton’s method using fifth-order quadrature formulas. Appl. Math. Comput.; 2007; 190, pp. 686-698. [DOI: https://dx.doi.org/10.1016/j.amc.2007.01.062]

33. Douglas, J.M. Process Dynamics and Control; Prentice Hall: Englewood Cliffs, NJ, USA, 1972.

34. Hoffman, J.D. Numerical Methods for Engineers and Scientists; McGraw-Hill Book Company: New York, NY, USA, 1992.

35. Kumar, S.; Kumar, D.; Sharma, J.R.; Cesarano, C.; Agarwal, P.; Chu, Y.M. An Optimal Fourth Order Derivative-Free Numerical Algorithm for Multiple Roots. Symmetry; 2020; 12, 1038. [DOI: https://dx.doi.org/10.3390/sym12061038]