Achieving Optimal Order in a Novel Family of Numerical Methods: Insights from Convergence and Dynamical Analysis Results

Abstract

In this manuscript, we introduce a novel parametric family of multistep iterative methods designed to solve nonlinear equations. This family is derived from a damped Newton’s scheme but includes an additional Newton step with a weight function and a “frozen” derivative, that is, the same derivative than in the previous step. Initially, we develop a quad-parametric class with a first-order convergence rate. Subsequently, by restricting one of its parameters, we accelerate the convergence to achieve a third-order uni-parametric family. We thoroughly investigate the convergence properties of this final class of iterative methods, assess its stability through dynamical tools, and evaluate its performance on a set of test problems. We conclude that there exists one optimal fourth-order member of this class, in the sense of Kung–Traub’s conjecture. Our analysis includes stability surfaces and dynamical planes, revealing the intricate nature of this family. Notably, our exploration of stability surfaces enables the identification of specific family members suitable for scalar functions with a challenging convergence behavior, as they may exhibit periodical orbits and fixed points with attracting behavior in their corresponding dynamical planes. Furthermore, our dynamical study finds members of the family of iterative methods with exceptional stability. This property allows us to converge to the solution of practical problem-solving applications even from initial estimations very far from the solution. We confirm our findings with various numerical tests, demonstrating the efficiency and reliability of the presented family of iterative methods.

1. Introduction

A multitude of challenges in Computational Sciences and other fields in Science and Technology can be effectively represented as nonlinear equations through mathematical modeling, see for example [1,2,3]. Finding solutions $\xi$ to nonlinear equations of the form $f\left(x\right)=0$ stands as a classical yet formidable problem in the realms of Numerical Analysis, Applied Mathematics, and Engineering. Here, the function $f:I\subset \mathbb{R}\to \mathbb{R}$ is assumed to be differentiable enough within the open interval I. Extensive overviews of iterative methods for solving nonlinear equations published in recent years can be found in [4,5,6], and their associated references.

In recent years, many iterative methods have been developed to solve nonlinear equations. The essence of these methods is as follows: if one knows a sufficiently small domain that contains only one root $\xi$ of the equation $f\left(x\right)=0$, and we select a sufficiently close initial estimate of the root $x_0$, we generate a sequence of iterates $x_1,x_2,\dots ,x_k,\dots$, by means of a fixed point function $g\left(x\right)$, which under certain conditions converges to $\xi$. The convergence of the sequence is guaranteed, among other elements, by the appropriate choice of the function g and the initial approximation $x_0$.

The method described by the iteration function $g:I\subseteq \mathbb{R}\to \mathbb{R}$ such that

$$x_{k+1}=g\left(x_k\right),k=0,1,\dots ,$$

(1)

starting from a given initial estimate $x_0$, includes a large number of iterative schemes. These differ from each other by the way the iteration function g is defined.

Among these methods, Newton’s scheme is widely acknowledged as the most renowned approach for locating a solution $\xi \in I$. This scheme is defined by the iterative formula:

$$\begin{array}{cc}\hfill x_{k+1}& =x_k-{\displaystyle \frac{f\left(x_k\right)}{f^{\prime }\left(x_k\right)}},\hfill \end{array}$$

where $k=0,1,2,\dots$, and $f^{\prime }\left(x_k\right)$ denotes the derivative of the function f evaluated in the kth iteration.

A very important concept of iterative methods is their order of convergence, which provides a measure of the speed of convergence of the iterates. Let ${\left\{x_k\right\}}_{k_{\ge }0}$ be a sequence of real numbers such that ${lim}_{k\to \infty }{x}_k=\xi$. The convergence is called (see [7]):

(a)

Linear, if there exist C, $0<C<1$ and $k_0\in \mathbb{N}$ such that

$$\frac{|x_k-\xi |}{|x_{k-1}-\xi |}\le C, \mathrm{for} \mathrm{all} k>k_0;$$

(b)

Is of order p, if there exist $C>0$ and $k_0\in \mathbb{N}$ such that

$$\frac{|x_k-\xi |}{|x_{k-1}{-\xi |}^p}\le C, \mathrm{for} \mathrm{all} k>k_0.$$

We denote by $e_k=x_k-\xi$ the error of the k-th iteration. Moreover, equation $e_{k+1}=C{e}_k^p+O\left(e_k^{p+1}\right)$ is called the error equation of the iterative method, where p is its order of convergence and C is called the asymptotic error constant.

It is known (see, for example, [4]) that if $x_{k+1}=g\left(x_k\right)$ is an iterative point-to-point method with d functional evaluations per step, then the order of convergence of the method is, at most, $p=d$. On the other hand, Traub proves in [4] that to design a point-to-point method of order p, the iterative expression must contain derivatives of the nonlinear function whose zero we are looking for, at least of order $p-1$. This is why point-to-point methods are not efficient if we seek to simultaneously increase the order of convergence and computational efficiency.

These restrictions of point-to-point methods are the starting point for the growing interest of researchers in multipoint methods, see for example [4,5,6]. In such schemes, also called predictor–corrector, the $(k+1)$-th iterate is obtained by using functional evaluations of the k-th iterate and also of other intermediate points. For example, a two-step multipoint method has the expression

$$\begin{array}{ccc}\hfill y_k& =& \Psi \left(x_k\right),\hfill \\ \hfill x_{k+1}& =& \Phi (x_k,y_k),k=0,1,2,\dots \hfill \end{array}$$

Thus, the main motivation for designing new iterative schemes is to increase the order of convergence without adding many functional evaluations. The first multipoint schemes were designed by Traub in [4]. At that time the concept of optimality had not yet been defined and the fact of designing multipoint schemes with the same order as classical schemes such as Halley or Chebyshev, but with a much simpler iterative expression and without using second derivatives, was of great importance. The techniques used then have been the seed of those that allowed the appearance of higher-order methods.

In recent years, different authors have developed a large number of optimal schemes for solving nonlinear equations [6,8]. A common way to increase the convergence order of an iterative scheme is to use the composition of methods, based on the following result (see [4]).

Theorem 1.

Let $g_1\left(x\right)$ and $g_2\left(x\right)$ be the fixed-point functions of orders $p_1$ and $p_2$, respectively. Then, the iterative method resulting from composing them, $x_{k+1}=g_1\left(g_2\left(x_k\right)\right)$, $k=0,1,2,\dots$, has an order of convergence $p_1{p}_2$.

However, this composition necessarily increases the number of functional evaluations. So, to preserve optimality, the number of evaluations must be reduced. There are many techniques used for this purpose by different authors, such as approximating some of the evaluations that have appeared with the composition by means of interpolation polynomials, Padé approximants, inverse interpolation, Adomian polynomials, etc. (see, for example, [6,9,10]). If after the reduction of functional evaluations the resulting method is not optimal, the weight function technique, introduced by Chun in [11], can be used to increase its order of convergence.

There are also other ways in the literature to compare different iterative methods with each other. Traub in [4] defined the information efficiency of an iterative method as

$$I\left(M\right)=\frac{p}{d},$$

where p is the order of convergence and d is the number of functional evaluations per iteration. On the other hand, Ostrowski in [12] introduced the so-called efficiency index,

$$EI\left(M\right)=p^{1/d},$$

which, in turn, gives rise to the concept of optimality of an iterative method.

Regarding the order of convergence, Kung and Traub in their conjecture (see [13]) establish what is the highest order that a multipoint iterative scheme without memory can reach. Schemes that attain this limit are called optimal methods. Such a conjecture states that the order of convergence of any memoryless multistep method cannot exceed $2^{d-1}$ (called optimal order), where d is the number of functional evaluations per iteration, with efficiency index $2^{(d-1)/d}$ (called optimal index). In this sense, Newton is an optimal method.

Furthermore, in order to numerically test the behavior of the different iterative methods, Weerakoon and Fernando in [14] introduced the so-called computational order of convergence (COC),

$$p\approx COC=\frac{ln\frac{|x_{k+1}-\xi |}{|x_k-\xi |}}{ln\frac{|x_k-\xi |}{|x_{k-1}-\xi |}},k=1,2,\dots ,$$

where $x_{k+1}$, $x_k$ and $x_{k-1}$ are three consecutive approximations of the root of the nonlinear equation, obtained in the iterative process. However, the value of the zero $\xi$ is not known in practice, which motivated the definition in [15] of the approximate computational convergence order ACOC,

$$p\approx ACOC=\frac{ln\frac{|x_{k+1}-x_k|}{|x_k-x_{k-1}|}}{ln\frac{|x_k-x_{k-1}|}{|x_{k-1}-x_{k-2}|}},k=2,3,\dots .$$

(2)

On the other hand, the dynamical analysis of rational operators derived from iterative schemes, particularly when applied to low-degree nonlinear polynomial equations, has emerged as a valuable tool for assessing the stability and reliability of these numerical methods. This approach is detailed, for instance, in Refs. [16,17,18,19,20] and their associated references.

Using the tools of complex discrete dynamics, it is possible to compare different algorithms in terms of their basins of attraction, the dynamical behavior of the rational functions associated with the iterative method on low-degree polynomials, etc. Varona [21], Amat et al. [22], Neta et al. [23], Cordero et al. [24], Magreñán [25], Geum et al. [26], among others, have analyzed many schemes and parametric families of methods under this point of view, obtaining interesting results about their stability and reliability.

The dynamical analysis of an iterative method focuses on the study of the asymptotic behavior of the fixed points (roots, or not, of the equation) of the operator, as well as on the basins of attraction associated with them. In the case of parametric families of iterative methods, the analysis of the free critical points (points where the derivative of the operator cancels out that are not roots of the nonlinear function) and stability functions of the fixed points allows us to select the most stable members of these families. Some of the existing works in the literature related to this approach are Refs. [27,28], among others.

In this paper, we introduce a novel parametric family of multistep iterative methods tailored for solving nonlinear equations. This family is constructed by enhancing the traditional Newton’s scheme, incorporating an additional Newton step with a weight function and a frozen derivative. As a result, the family is characterized by a two-step iterative expression that relies on four arbitrary parameters.

Our approach yields a third-order uni-parametric family and a fourth-order member. However, in the course of developing these iterative schemes, we initially start with a first-order quad-parametric family. By selectively setting just one parameter, we manage to accelerate its convergence to a third-order scheme, and for a specific value of this parameter, we achieve an optimal member. To substantiate these claims, we conduct a comprehensive convergence analysis for all classes.

The stability of this newly introduced family is rigorously examined using dynamical tools. We construct stability surfaces and dynamical planes to illustrate the intricate behavior of this class. These stability surfaces help us to identify specific family members with exceptional behavior, making them well-suited for practical problem-solving applications. To further demonstrate the efficiency and reliability of these iterative schemes, we conduct several numerical tests.

The rest of the paper is organized as follows. In Section 2, we present the proposed class of iterative methods depending on several parameters, which is step-by-step modified in order to achieve the highest order of convergence. Section 3 is devoted to the dynamical study of the uni-parametric family; by means of this analysis, we find the most stable members, less dependent from their initial estimation. In Section 4, the previous theoretical results are checked by means of numerical tests on several nonlinear problems, using a wide variety of initial guesses and parameter values. Finally, some conclusions are presented.

2. Convergence Analysis of the Family

In this section, we conduct a convergence analysis of the newly introduced quad-parametric iterative family, with the following iterative expression:

$$\left\{\begin{array}{cc}\hfill y_k& =x_k-\alpha {\displaystyle \frac{f\left(x_k\right)}{f^{\prime }\left(x_k\right)}},\hfill \\ \hfill x_{k+1}& =y_k-\left(\beta +\gamma {\displaystyle \frac{f\left(y_k\right)}{f\left(x_k\right)}}+\delta {\left({\displaystyle \frac{f\left(y_k\right)}{f\left(x_k\right)}}\right)}^2\right){\displaystyle \frac{f\left(x_k\right)}{f^{\prime }\left(x_k\right)}},\hfill \end{array}\right.$$

(3)

where $\alpha$, $\beta$, $\gamma$, $\delta$ are arbitrary parameters and $k=0,1,2,\dots$.

Additionally, we present a strategy for simplifying it into a uni-parametric class to enhance convergence speed. Consequently, even though the quad-parametric family has a first-order convergence rate, we employ higher-order Taylor expansions in our proof, as they are instrumental in establishing the convergence rate of the uni-parametric subfamily. In Appendix A, the Mathematica code used for checking it is available.

Theorem 2 (quad-parametric family).

Let $f:I\subseteq \mathbb{R}\to \mathbb{R}$ be a sufficiently differentiable function in an open interval I and $\xi \in I$ a simple root of the nonlinear equation $f\left(x\right)=0$. Let us suppose that $f^{\prime }\left(x\right)$ is continuous and nonsingular at ξ, and $x_0$ is an initial estimate close enough to ξ. Then, the sequence ${\left\{x_k\right\}}_{k\ge 0}$ obtained by using the expression (3) converges to ξ with order one, being its error equation

$$\begin{array}{cc}\hfill e_{k+1}& =\left(-{\alpha }^2\delta +\alpha (\gamma +2\delta -1)-\beta -\gamma -\delta +1\right)e_k+𝒪\left(e_k^2\right),\hfill \end{array}$$

where $e_k=x_k-\xi$, and α, β, γ, and δ are free parameters.

Proof. 

Let us consider $\xi$ as the simple root of nonlinear function $f\left(x\right)$ and $x_k=\xi +e_k$. We calculate the Taylor expansion of $f\left(x_k\right)$ and $f^{\prime }\left(x_k\right)$ around the root $\xi$, we get

$$\begin{array}{cc}\hfill f\left(x_k\right)& =f\left(\xi \right)+f^{\prime }\left(\xi \right)e_k+\frac{1}{2!}{f}^{″}\left(\xi \right)e_k^2+\frac{1}{3!}{f}^{‴}\left(\xi \right)e_k^3+\frac{1}{4!}{f}^{\left(iv\right)}\left(\xi \right)e_k^4+𝒪\left(e_k^5\right)\hfill \\ \hfill & =f^{\prime }\left(\xi \right)\left[e_k+\frac{1}{2!}\frac{f^{″}\left(\xi \right)}{f^{\prime }\left(\xi \right)}{e}_k^2+\frac{1}{3!}\frac{f^{‴}\left(\xi \right)}{f^{\prime }\left(\xi \right)}{e}_k^3+\frac{1}{4!}\frac{f^{\left(iv\right)}\left(\xi \right)}{f^{\prime }\left(\xi \right)}{e}_k^4\right]+𝒪\left(e_k^5\right)\hfill \\ \hfill & =f^{\prime }\left(\xi \right)\left[e_k+C_2{e}_k^2+C_3{e}_k^3+C_4{e}_k^4\right]+𝒪\left(e_k^5\right),\hfill \end{array}$$

(4)

and

$$\begin{array}{cc}\hfill f^{\prime }\left(x_k\right)& =f^{\prime }\left(\xi \right)+f^{″}\left(\xi \right)e_k+\frac{1}{2!}{f}^{‴}\left(\xi \right)e_k^2+\frac{1}{3!}{f}^{\left(iv\right)}\left(\xi \right)e_k^3+𝒪\left(e_k^4\right)\hfill \\ \hfill & =f^{\prime }\left(\xi \right)\left[1+\frac{f^{″}\left(\xi \right)}{f^{\prime }\left(\xi \right)}{e}_k+\frac{1}{2!}\frac{f^{‴}\left(\xi \right)}{f^{\prime }\left(\xi \right)}{e}_k^2+\frac{1}{3!}\frac{f^{\left(iv\right)}\left(\xi \right)}{f^{\prime }\left(\xi \right)}{e}_k^3\right]+𝒪\left(e_k^4\right)\hfill \\ \hfill & =f^{\prime }\left(\xi \right)\left[1+2C_2{e}_k+3C_3{e}_k^2+4C_4{e}_k^3\right]+𝒪\left(e_k^4\right),\hfill \end{array}$$

(5)

where $C_p={\displaystyle \frac{1}{p!}}{\displaystyle \frac{f^{\left(p\right)}\left(\xi \right)}{f^{\prime }\left(\xi \right)}}$, $p=2,3,...$

By a direct division of (4) and (5),

$$\begin{array}{cc}\hfill \frac{f\left(x_k\right)}{f^{\prime }\left(x_k\right)}& =e_k-C_2{e}_k^2+2\left(C_2^2-C_3\right)e_k^3-\left(4C_2^3-7C_2{C}_3+3C_4\right)e_k^4+𝒪\left(e_k^5\right).\hfill \end{array}$$

(6)

Replacing (6) in (3), we have

$$\begin{array}{cc}\hfill y_k& =\xi +(1-\alpha )e_k+\alpha C_2{e}_k^2-2\alpha \left(C_2^2-C_3\right)e_k^3+\alpha \left(4C_2^3-7C_2{C}_3+3C_4\right)e_k^4+𝒪\left(e_k^5\right).\hfill \end{array}$$

(7)

Again a Taylor expansion of $f\left(y_k\right)$ around $\xi$ allows us to get

$$\begin{array}{cc}\hfill f\left(y_k\right)& =f^{\prime }\left(\xi \right)\left[(1-\alpha )e_k+\left({\alpha }^2-\alpha +1\right)C_2{e}_k^2+\left(-2{\alpha }^2{C}_2^2-\left({\alpha }^3-3{\alpha }^2+\alpha -1\right)C_3\right)e_k^3\right.\hfill \\ \hfill & \left.+\left(5{\alpha }^2{C}_2^3+{\alpha }^2(3\alpha -10)C_2{C}_3+\left({\alpha }^4-4{\alpha }^3+6{\alpha }^2-\alpha +1\right)C_4\right)e_k^4\right]+𝒪\left(e_k^5\right).\hfill \end{array}$$

(8)

Dividing (8) by (4), we obtain

$$\begin{array}{cc}\hfill \frac{f\left(y_k\right)}{f\left(x_k\right)}& =(1-\alpha )+{\alpha }^2{C}_2{e}_k-{\alpha }^2\left((\alpha -3)C_3+3C_2^2\right)e_k^2\hfill \\ \hfill & +{\alpha }^2\left(\left({\alpha }^2-4\alpha +6\right)C_4+2(2\alpha -7)C_2{C}_3+8C_2^3\right)e_k^3+𝒪\left(e_k^4\right).\hfill \end{array}$$

(9)

Finally, substituting (6), (7) and (9), in the second step of family (3), we have

$$\begin{array}{cc}\hfill x_{k+1}& =\xi +A_1{e}_k+A_2{e}_k^2+A_3{e}_k^3+A_4{e}_k^4+𝒪\left(e_k^5\right),\hfill \end{array}$$

(10)

where

$$\begin{array}{cc}\hfill A_1& =-{\alpha }^2\delta +\alpha (\gamma +2\delta -1)-\beta -\gamma -\delta +1,\hfill \\ \hfill A_2& =\left(2{\alpha }^3\delta -{\alpha }^2(\gamma +\delta )-\alpha (\gamma +2\delta -1)+\beta +\gamma +\delta \right)C_2,\hfill \\ \hfill A_3& =\left(-2{\alpha }^4\delta +{\alpha }^3(\gamma +8\delta )-{\alpha }^2(3\gamma +4\delta )-2\alpha (\gamma +2\delta -1)+2(\beta +\gamma +\delta )\right)C_3\hfill \\ \hfill & -\left({\alpha }^4\delta +8{\alpha }^3\delta -2{\alpha }^2(2\gamma +3\delta )-2\alpha (\gamma +2\delta -1)+2(\beta +\gamma +\delta )\right)C_2^2,\hfill \\ \hfill A_4& =\left(7{\alpha }^4\delta +26{\alpha }^3\delta -{\alpha }^2(13\gamma +22\delta )-4\alpha (\gamma +2\delta -1)+4(\beta +\gamma +\delta )\right)C_2^3\hfill \\ \hfill & +\left(2{\alpha }^5\delta +4{\alpha }^4\delta -{\alpha }^3(5\gamma +48\delta )+{\alpha }^2(19\gamma +31\delta )+7\alpha (\gamma +2\delta -1)-7(\beta +\gamma +\delta )\right)C_2{C}_3\hfill \\ \hfill & +\left(2{\alpha }^5\delta -{\alpha }^4(\gamma +10\delta )+4{\alpha }^3(\gamma +5\delta )-3{\alpha }^2(2\gamma +3\delta )-3\alpha (\gamma +2\delta -1)+3(\beta +\gamma +\delta )\right)C_4,\hfill \end{array}$$

(11)

being the error equation

$$\begin{array}{cc}\hfill e_{k+1}& =A_1{e}_k+A_2{e}_k^2+A_3{e}_k^3+A_4{e}_k^4+𝒪\left(e_k^5\right)\hfill \\ \hfill & =\left(-{\alpha }^2\delta +\alpha (\gamma +2\delta -1)-\beta -\gamma -\delta +1\right)e_k+𝒪\left(e_k^2\right),\hfill \end{array}$$

(12)

and the proof is finished.    □

From Theorem 2, it is evident that the newly introduced quad-parametric family exhibits a convergence order of one, irrespective of the values assigned to $\alpha$, $\beta$, $\gamma$, and $\delta$. Nevertheless, we can expedite convergence by holding only two parameters constant, effectively reducing the family to a bi-parametric iterative scheme. In Appendix B, the Mathematica code used for checking it is available.

Theorem 3 (bi-parametric family).

Let $f:I\subseteq \mathbb{R}\to \mathbb{R}$ be a sufficiently differentiable function in an open interval I and $\xi \in I$ a simple root of the nonlinear equation $f\left(x\right)=0$. Let us suppose that $f^{\prime }\left(x\right)$ is continuous and nonsingular at ξ, and $x_0$ is an initial estimate close enough to ξ. Then, the sequence ${\left\{x_k\right\}}_{k\ge 0}$ obtained by using the expression (3) converges to ξ with order three, provided that ${\displaystyle \beta =\frac{{(\alpha -1)}^2\left({\alpha }^2\delta -\alpha -1\right)}{{\alpha }^2}}$ and ${\displaystyle \gamma =\frac{2{\alpha }^3\delta -2{\alpha }^2\delta +1}{{\alpha }^2}}$, being its error equation

$$\begin{array}{cc}\hfill e_{k+1}& =\left(-\left({\alpha }^4\delta -2\right)C_2^2+(\alpha -1)C_3\right)e_k^3+𝒪\left(e_k^4\right),\hfill \end{array}$$

where $e_k=x_k-\xi$, $C_q={\displaystyle \frac{1}{q!}}{\displaystyle \frac{f^{\left(q\right)}\left(\xi \right)}{f^{\prime }\left(\xi \right)}}$, $q=2,3,...$, and α, δ are arbitrary parameters.

Proof. 

Using the results of Theorem 2 to cancel $A_1$ and $A_2$ accompanying $e_k$ and $e_k^2$ in (12), respectively, it must be satisfied that

$$\left\{\begin{array}{cc}\hfill -{\alpha }^2\delta +\alpha (\gamma +2\delta -1)-\beta -\gamma -\delta +1& =0,\hfill \\ \hfill 2{\alpha }^3\delta -{\alpha }^2(\gamma +\delta )-\alpha (\gamma +2\delta -1)+\beta +\gamma +\delta & =0.\hfill \end{array}\right.$$

(13)

It is clear that system (13) has infinite solutions for

$$\beta =\frac{{(\alpha -1)}^2\left({\alpha }^2\delta -\alpha -1\right)}{{\alpha }^2}\mathrm{and}\gamma =\frac{2{\alpha }^3\delta -2{\alpha }^2\delta +1}{{\alpha }^2},$$

(14)

where $\alpha$ and $\delta$ are free parameters. Therefore, replacing (14) in (11), we obtain that

$$\begin{array}{cc}\hfill A_1& =0,\hfill \\ \hfill A_2& =0,\hfill \\ \hfill A_3& =-\left({\alpha }^4\delta -2\right)C_2^2+(\alpha -1)C_3,\hfill \\ \hfill A_4& =\left(7{\alpha }^4\delta -9\right)C_2^3+\left(2(\alpha -3){\alpha }^4\delta -5\alpha +12\right)C_2{C}_3-(\alpha -3)(\alpha -1)C_4,\hfill \end{array}$$

(15)

being the error equation

$$\begin{array}{cc}\hfill e_{k+1}& =A_3{e}_k^3+𝒪\left(e_k^4\right)\hfill \\ \hfill & =\left(-\left({\alpha }^4\delta -2\right)C_2^2+(\alpha -1)C_3\right)e_k^3+𝒪\left(e_k^4\right),\hfill \end{array}$$

(16)

and the proof is finished.    □

According to the findings in Theorem 3, it is evident that the newly introduced bi-parametric family

$$\left\{\begin{array}{cc}\hfill y_k& =x_k-\alpha {\displaystyle \frac{f\left(x_k\right)}{f^{\prime }\left(x_k\right)}},\hfill \\ \hfill x_{k+1}& =y_k-\left(\beta +\gamma {\displaystyle \frac{f\left(y_k\right)}{f\left(x_k\right)}}+\delta {\left({\displaystyle \frac{f\left(y_k\right)}{f\left(x_k\right)}}\right)}^2\right)\left({\displaystyle \frac{f\left(x_k\right)}{f^{\prime }\left(x_k\right)}}\right),\hfill \end{array}\right.$$

(17)

where $k=0,1,2,...$, ${\displaystyle \beta =\frac{{(\alpha -1)}^2\left({\alpha }^2\delta -\alpha -1\right)}{{\alpha }^2}}$ and ${\displaystyle \gamma =\frac{2{\alpha }^3\delta -2{\alpha }^2\delta +1}{{\alpha }^2}}$ consistently exhibits a third-order convergence across all values of $\alpha$ and $\delta$. Nevertheless, it is noteworthy that by restricting one of the parameters while transitioning to a uni-parametric iterative scheme, not only can we sustain convergence, but we can also enhance performance. This improvement arises from the reduction in the error equation complexity, resulting in more efficient computations.

Corollary 1 (uni-parametric family).

Let $f:I\subseteq \mathbb{R}\to \mathbb{R}$ be a sufficiently differentiable function in an open interval I and $\xi \in I$ a simple root of the nonlinear equation $f\left(x\right)=0$. Let us suppose that $f^{\prime }\left(x\right)$ is continuous and nonsingular at ξ and $x_0$ is an initial estimate close enough to ξ. Then, the sequence ${\left\{x_k\right\}}_{k\ge 0}$ obtained by using the expression (17) converges to ξ with order three, provided that $ϵ={\alpha }^4\delta =2$, being its error equation

$$\begin{array}{cc}\hfill e_{k+1}& =(\alpha -1)C_3{e}_k^3+𝒪\left(e_k^4\right),\hfill \end{array}$$

where $e_k=x_k-\xi$, $C_q={\displaystyle \frac{1}{q!}}{\displaystyle \frac{f^{\left(q\right)}\left(\xi \right)}{f^{\prime }\left(\xi \right)}}$, $q=2,3,...$, and α is an arbitrary parameter. Indeed, $\alpha =1$ and, therefore, $\delta =ϵ=2$ provides the only member of the family of the optimal fourth-order of convergence.

Proof. 

Using the results of Theorem 3 to reduce the expression of $A_3$ accompanying $e_k^3$ in (15), it must be satisfied that ${\alpha }^4\delta -2=0$ and/or $\alpha -1=0$. It is easy to show that the first equation has infinite solutions for

$$ϵ={\alpha }^4\delta =2.$$

(18)

Therefore, replacing (18) in (15), we obtain that

$$\begin{array}{cc}\hfill e_{k+1}& =A_3{e}_k^3+𝒪\left(e_k^4\right)\hfill \\ \hfill & =(\alpha -1)C_3{e}_k^3+𝒪\left(e_k^4\right),\hfill \end{array}$$

(19)

and the proof is finished.    □

Based on the outcomes derived from Corollary 1, it becomes apparent that the recently introduced uni-parametric family, which we will call MCCTU($\alpha$),

$$\left\{\begin{array}{cc}\hfill y_k& =x_k-\alpha {\displaystyle \frac{f\left(x_k\right)}{f^{\prime }\left(x_k\right)}},\hfill \\ \hfill x_{k+1}& =y_k-\left(\beta +\gamma {\displaystyle \frac{f\left(y_k\right)}{f\left(x_k\right)}}+\delta {\left({\displaystyle \frac{f\left(y_k\right)}{f\left(x_k\right)}}\right)}^2\right)\left({\displaystyle \frac{f\left(x_k\right)}{f^{\prime }\left(x_k\right)}}\right),\hfill \end{array}\right.$$

(20)

where $k=0,1,2,...$, ${\displaystyle \beta =\frac{{(\alpha -1)}^2\left({\alpha }^2\delta -\alpha -1\right)}{{\alpha }^2}}$, ${\displaystyle \gamma =\frac{2{\alpha }^3\delta -2{\alpha }^2\delta +1}{{\alpha }^2}}$ and ${\displaystyle \delta =\frac{2}{{\alpha }^4}}$ consistently exhibits a convergence order of three, regardless of the chosen value for $\alpha$. Nevertheless, a remarkable observation emerges when $\alpha =1$: in such a case, a member of this family attains an optimal convergence order of four.

Due to the previous results, we have chosen to concentrate our efforts solely on the MCCTU($\alpha$) class of iterative schemes moving forward. To pinpoint the most effective members within this family, we will utilize dynamical techniques outlined in Section 3.

3. Stability Analysis

This section delves into the examination of the dynamical characteristics of the rational operator linked to the iterative schemes within the MCCTU($\alpha$) family. This exploration provides crucial insights into the stability and dependence of the members of the family with respect to the initial estimations used. To shed light on the performance, we create rational operators and visualize their dynamical planes. These visualizations enable us to discern the behavior of specific methods in terms of the attraction basins of periodic orbits, fixed points, and other relevant dynamics.

Now, we introduce the basic concepts of complex dynamics used in the dynamical analysis of iterative methods. The texts [29,30], among others, provide extensive and detailed information on this topic.

Given a rational function $R:\widehat{\mathbb{C}}\to \widehat{\mathbb{C}}$, where $\widehat{\mathbb{C}}$ is the Riemann sphere, the orbit of a point $z_0\in \widehat{\mathbb{C}}$ is defined as:

$$\{z_0,R^1\left(z_0\right),R^2\left(z_0\right),...,R^n\left(z_0\right),\dots \}.$$

We are interested in the study of the asymptotic behavior of the orbits depending on the initial estimate $z_0,$ analyzed in the dynamical plane of the rational function R defined by the different iterative methods.

To obtain these dynamical planes, we must first classify the fixed or periodic points of the rational operator R. A point $z_0\in \widehat{\mathbb{C}}$ is called a fixed point if it satisfies $R\left(z_0\right)=z_0$. If the fixed point is not a solution of the equation, it is called a strange fixed point. $z_0$ is said to be a periodic point of period $p>1$ if $R^p\left(z_0\right)=z_0$ and $R^k\left(z_0\right)\ne z_0$, $k<p$. A critical point $z_C$ is a point where $R^{\prime }\left(z_C\right)=0$.

On the other hand, a fixed point $z_0$ is called attracting if $|R^{\prime }\left(z_0\right)|<1$, superattracting if $|R^{\prime }\left(z_0\right)|=0$, repulsive if $|R^{\prime }\left(z_0\right)|>1$, and parabolic if $|R^{\prime }\left(z_0\right)|=1$.

The basin of attraction of an attractor $\overline{z}$ is defined as the set of pre-images of any order:

$$\mathcal{A}\left(\overline{z}\right)=\{z_0\in \widehat{\mathbb{C}}:R^n\left(z_0\right)\to \overline{z},n\to \infty \}.$$

The Fatou set consists of the points whose orbits have an attractor (fixed point, periodic orbit or infinity). Its complementary in $\widehat{\mathbb{C}}$ is the Julia set, $\mathcal{J}$. Therefore, the Julia set includes all the repulsive fixed points and periodic orbits, and also their pre-images. So, the basin of attraction of any fixed point belongs to the Fatou set. Conversely, the boundaries of the basins of attraction compose the Julia set.

The following classical result, which is due to Fatou [31] and Julia [32], includes both periodic points (of any period) and fixed points, considered as periodic points of the unit period.

Theorem 4. ([31,32]).

Let R be a rational function. The immediate basins of attraction of each attracting periodic point contain at least one critical point.

By means of this key result, all the attracting behavior can be found using the critical points as a seed.

3.1. Rational Operator

While the fixed-point operator can be formulated for any nonlinear function, our focus here lies on constructing this operator for low-degree nonlinear polynomial equations, in order to get a rational function. This choice stems from the fact that the stability or instability criteria applied to methods on these equations can often be extended to other cases. Therefore, we introduce the following nonlinear equation represented by $f\left(x\right)$:

$$f\left(x\right)=(x-a)(x-b)=0,$$

(21)

where $a,b\in \mathbb{R}$ are the roots of the polynomial.

Let us remark that when MCCTU($\alpha$) is directly applied to $f\left(x\right)$, parameter $\alpha$ disappears in the resulting rational expression; so, no dynamical analysis can be made. However, if we use parameter $ϵ={\alpha }^4\delta$ appearing in Corollary 1 the same class of iterative methods can be expressed as MCCTU($ϵ$) and the dynamical analyisis can be made depending on $ϵ$.

Proposition 1 (rational operator${\mathit{R}}_{\mathit{f}}$).

Let the polynomial equation $f\left(x\right)$ given in (21), for $a,b\in \mathbb{C}$. Rational operator $R_f$ related to the MCCTU(ϵ) family given in (20) on $f\left(x\right)$ is

$$\begin{array}{cc}\hfill R_f(x,ϵ)& =\frac{x^3\left(ϵ-x^3-4x^2-5x-2\right)}{x^3\left(ϵ-2\right)-5x^2-4x-1},\hfill \end{array}$$

(22)

with $ϵ\in \mathbb{C}$ being an arbitrary parameter.

Proof. 

Let $f\left(x\right)$ be a generic quadratic polynomial function with roots $a,b\in \mathbb{C}$. We apply the iterative scheme MCCTU($ϵ$) given in (20) on $f\left(x\right)$ and obtain a rational function $A_f(x,ϵ)$ that depends on the roots $a,b\in \mathbb{C}$ and the parameters $ϵ\in \mathbb{C}$. Then, by using a Möbius transformation (see [22,33,34]) on $A_f(x,ϵ)$ with

$$\begin{array}{c}\hfill h\left(w\right)=\frac{w-a}{w-b},\end{array}$$

satisfying $h(\infty )=1,h\left(a\right)=0$ and $h\left(b\right)=\infty$, we get

$$\begin{array}{cc}\hfill R_f(x,ϵ)& =\left(h\circ A_f(x,ϵ)\circ h^{-1}\right)\left(x\right)=\frac{x^3\left(ϵ-x^3-4x^2-5x-2\right)}{x^3\left(ϵ-2\right)-5x^2-4x-1},\hfill \end{array}$$

(23)

which depends on two arbitrary parameters $ϵ\in \mathbb{C}$, thus completing the proof.    □

From Proposition 1, if we set $ϵ-2=0$, we obtain

$$\begin{array}{cc}\hfill \delta & ={\displaystyle \frac{2}{{\alpha }^4}},\hfill \end{array}$$

(24)

and, then, it is easy to show that the rational operator $R_f(x,ϵ)$ simplifies to the expression

$$\begin{array}{cc}\hfill R_f\left(x\right)& ={\displaystyle \frac{x^4\left(x^2+4x+5\right)}{5x^2+4x+1}},\hfill \end{array}$$

(25)

which does not depend on any free parameters.

3.2. Fixed Points

Now, we calculate all the fixed points of $R_f(x,ϵ)$ given by (22), to afterwards analyze their character (attracting, repulsive, or neutral or parabolic).

Proposition 2.

The fixed points of $R_f(x,ϵ)$ are $x=0$, $x=\infty$, and also five strange fixed points:

• $e{x}_1=1$,
• $e{x}_{2,3}\left(ϵ\right)=-{\displaystyle \frac{5}{4}}-{\displaystyle \frac{1}{4}}\sqrt{1-4ϵ}\pm {\displaystyle \frac{1}{2}}\sqrt{{\displaystyle \frac{5}{2}}-ϵ-{\displaystyle \frac{20\left(ϵ+8\right)-165}{2\sqrt{1-4ϵ}}}}$, and 
• $e{x}_{4,5}\left(ϵ\right)=-{\displaystyle \frac{5}{4}}+{\displaystyle \frac{1}{4}}\sqrt{1-4ϵ}\pm {\displaystyle \frac{1}{2}}\sqrt{{\displaystyle \frac{5}{2}}-ϵ+{\displaystyle \frac{20\left(ϵ+8\right)-165}{2\sqrt{1-4ϵ}}}}$.

By using Equation (24), the strange fixed points $e{x}_{2,3}\left(ϵ\right)$ and $e{x}_{4,5}\left(ϵ\right)$ do not depend on any free parameter,

• $e{x}_{2,3}\left(2\right)=-2.1943\pm 1.5370i$, and 
• $e{x}_{4,5}\left(2\right)=-0.3057\pm 0.2142i$.

Morover, strange fixed points depending on ϵ are conjugated, $e{x}_{2,3}\left(ϵ\right)$ and $e{x}_{4,5}\left(ϵ\right)$. If $ϵ=\frac{1}{4}$, $e{x}_1\left(ϵ\right)=e{x}_3\left(ϵ\right)$ and $e{x}_2\left(ϵ\right)=e{x}_2\left(ϵ\right)$, so the amount of strange fixed points is three. Indeed, $e{x}_3(-20)=e{x}_4(-20)=1$ and $e{x}_3\left(0\right)=e{x}_4\left(0\right)=-1$.

From Proposition 2, we establish that there are seven fixed points. Among these, 0 and ∞ come from the roots a and b of $f\left(x\right)$. $e{x}_1=1$ comes from the divergence of the original scheme, previous to the Möbius transformation.

Proposition 3.

The strange fixed point $e{x}_1=1$, $\forall ϵ\in \mathbb{C}$, has the following character:

(i) 

If $\left|ϵ-12\right|>32$, then $e{x}_1$ is an attractor.

(ii) 

If $\left|ϵ-12\right|<32$, then $e{x}_1$ is a repulsor.

(iii) 

If $\left|ϵ-12\right|=32$, then $e{x}_1$ is parabolic.

Moreover, $e{x}_1$ can be attracting but not superattracting. The superattracting fixed points of $R_f$ are $x=0$, $x=\infty$, and the strange fixed points $e{x}_{4,5}\left(ϵ\right)$ for $ϵ=\frac{1}{9}\left(-5\sqrt{97}-47\right)$ and $ϵ=\frac{1}{9}\left(5\sqrt{97}-47\right)$.

In the particular case of $ϵ=2$ (using the Equation (24)), all the strange fixed points are repulsive.

Proof. 

We prove this result by analyzing the stability of the fixed points found in Proposition 2. It must be done by evaluating $\left|R_f^{\prime }(x,ϵ)\right|$ at each fixed point and, if it is lower, equal, or greater than one it is called attracting, neutral, or repulsive, respectively.

The cases of $x=0$ and ∞ are straightforward from the expression of $R_f(x,ϵ)$. When $e{x}_1\left(ϵ\right)$ is studied, then

$$\left|R_f^{\prime }(1,ϵ)\right|=\left|{\displaystyle \frac{32}{12-ϵ}}\right|,$$

so it is attracting, repelling or neutral if $\left|ϵ-12\right|$ is greater, lower, or equal to 32. It can be graphically viewed in Figure 1.

By a graphical and numerical study of $\left|R_f^{\prime }(e{x}_i\left(ϵ\right),ϵ)\right|$, $i=1,2,3,4$, it can be deduced that $e{x}_{2,3}\left(ϵ\right)$ are repulsive for all $ϵ$, meanwhile $e{x}_{4,5}\left(ϵ\right)$ are superattracting for $ϵ=\frac{1}{9}\left(-5\sqrt{97}-47\right)\approx -10.6938$ or $ϵ=\frac{1}{9}\left(5\sqrt{97}-47\right)\approx 0.249365$. Their stability function is presented in Figure 2a,b. Moreover, $e{x}_1$ can not be a superattractor as $\left|R_f^{\prime }(1,ϵ)\right|\ne 0$.    □

It is clear that 0 and ∞ are always superattracting fixed points, but the stability of the remaining fixed points depends on the values of $ϵ$. According to Proposition 3, two strange fixed points can become superattractors. This implies that there would exist basins of attraction for them, potentially causing the method to fail to converge to the solution. However, even when they are only attracting (that can be the case of $e{x}_1$), these basins of attraction exist.

As we have stated previously, Figure 1 represents the stability function of the strange fixed point $e{x}_1$. In this figure, the zones of attraction are the yellow area and the repulsion zone corresponds to the grey area. For values of $ϵ$ within the disk, $e{x}_1$ is repulsive; whereas for values of $ϵ$ outside the grey disk, $e{x}_1$ becomes attracting. So, it is natural to select values within the grey disk, as a repulsive divergence improves the performance of the iterative scheme.

Similar conclusions can be stated from the stability region of strange fixed points $e{x}_{4,5}\left(ϵ\right)$, appearing in Figure 2. When a value of parameter $ϵ$ is taken in the yellow area of Figure 2, both points are simultaneously attracting, so there are at least four different basins of attraction.

However, the basins of attraction also appear when there exist attracting periodic orbits of any period. To detect this kind of behavior, the role of critical points is crucial.

3.3. Critical Points

Now, we obtain the critical points of $R_f(x,ϵ)$.

Proposition 4.

The critical points of $R_f(x,ϵ)$ are $x=0$, $x=\infty$ and also:

• $c{r}_1=-1$, and 
• $c{r}_{2,3}\left(ϵ\right)={\displaystyle \frac{2ϵ+6\pm \sqrt{5}\sqrt{12ϵ-{ϵ}^2}}{3\left(ϵ-2\right)}}$.

Morover, if $ϵ=2$, critical points are not free $c{r}_{2,3}\left(2\right)=0$. In any other case, $c{r}_{2,3}\left(ϵ\right)$ are conjugated free critical points.

From Proposition 4, we establish that, in general, there are five critical points. The free critical point $c{r}_1=-1$ is a pre-image of the strange fixed point $e{x}_1=1$. Therefore, the stability of $c{r}_1$ corresponds to the stability of $e{x}_1$ (see Section 3.2). Note that if the Equation (24) is satisfied, the only remaining free critical point is $c{r}_1$. Since $c{r}_1$ is the pre-image of $e{x}_1$, it would be a repulsor.

Then, we use the only independent free critical point $c{r}_2\left(ϵ\right)$ (conversely, $c{r}_3\left(ϵ\right)$, as they are conjugates) to generate the parameter plane. This a graphical representation of the global stability performance of the member of the class of iterative methods. In a definite area of the complex plane, a mesh of $500\times 500$ points is generated. Each one of these points is used as a value of parameter $ϵ$, i.e., we get a particular element of the family. For each one of these values, we get as our initial guess the critical point $c{r}_2\left(ϵ\right)$ and calculate its orbit. If it converges to $x=0$ or $x=\infty$, then the point corresponding to this value of $ϵ$ is represented using a red color. In other case, it is left in black. So, convergent schemes to the original roots of the quadratic equations appear in the red stable area and the black area corresponds to schemes of the classes that are not able to converge to them, by reason of an attracting strange fixed point or periodic orbit. This performance can be seen in Figure 3, representing the domain $D_1=[-30,50]\times [-40,40]$, where a wide area of stable performance can be found around the origin, $D_2=[-5,15]\times [-10,10]$ (Figure 3b).

3.4. Dynamical Planes

A dynamical plane is defined as a mesh in a limited domain of the complex plane, where each point corresponds to a different initial estimate $x_0$. The graphical representation shows the method’s convergence starting from $x_0$ within a maximum of 80 iterations and ${10}^{-3}$ as the tolerance. Fixed points appear as a white circle ‘○’, critical points are ‘□’, and  a white asterisk ‘∗’ symbolizes an attracting point. Additionally, the basins of attraction are depicted in different colors. To generate this graph, we use MATLAB R2020b with a resolution of $400\times 400$ pixels.

Here, we analyze the stability of various MCCTU($ϵ$) methods using dynamical planes. We consider methods with $ϵ$ values both inside and outside the stability surface of $e{x}_1$, specifically, in the red and black areas of the parameter plane represented in Figure 3a.

Firstly, examples of methods within the stability region are provided for $ϵ\in \{1,2,10,5+5i\}$. Their dynamical planes, along with their respective basins of attraction, are shown in Figure 4. Let us remark that all selected values of $ϵ$ lie in the red area of the parameter plane and have only two basins of attraction, corresponding to $x=0$ (in orange color in the figures) and $x=\infty$ (blue in the figures).

Secondly, some schemes outside the stability region (in black in the parameter plane) are provided for $ϵ\in \{100,15,-15,30\}$. Their dynamical planes are shown in Figure 5. Each of these members have specific characteristics: in Figure 5a, the widest basin of attraction (in green color) corresponds to $e{x}_1=1$, which is attracting for this value of $ϵ$, the basin of $x=0$ is a very narrow area around the point; for $ϵ=15$, we observe in Figure 5b three different basins of attraction, the third of the two being attracting periodic orbits of period 2 (one of them is plotted in yellow in the figure); Figure 5c corresponds to $ϵ=-15$, inside the stability area of $e{x}_{4,5}\left(ϵ\right)$ (see Figure 2), where both are simultaneously attracting; finally, for $ϵ=30$, the widest basin of attraction corresponds to an attracting periodic orbit of period 2, see Figure 5d.

4. Numerical Results

In this section, we conduct several numerical tests to validate the theoretical convergence and stability results of the MCCTU($\alpha$) family obtained in previous sections. We use both stable and unstable methods from (20) and apply them to ten nonlinear test equations, with their expressions and corresponding roots provided in

Table 1. Nonlinear test equations and corresponding roots.

Nonlinear Test Equations	Roots
$f_1\left(x\right)=sin\left(x\right)-x^2+1=0$	$\xi \approx -0.63673$
$f_2\left(x\right)=x^2-e^x-3x+2=0$	$\xi \approx 0.25753$
$f_3\left(x\right)=cos\left(x\right)-x{e}^x+x^2=0$	$\xi \approx 0.63915$
$f_4\left(x\right)=e^x-1.5-arctan\left(x\right)=0$	$\xi \approx -14.10127$
$f_5\left(x\right)=x^3+4x^2-10=0$	$\xi \approx 1.36523$
$f_6\left(x\right)=8x-cos\left(x\right)-2x^2=0$	$\xi \approx 0.12808$
$f_7\left(x\right)=x{e}^{x^2}-{sin}^2\left(x\right)+3cos\left(x\right)+5=0$	$\xi \approx -1.20765$
$f_8\left(x\right)=\sqrt{x^2+2x+5}-2sin\left(x\right)-x^2+3=0$	$\xi \approx 2.33197$
$f_9\left(x\right)=x^4+sin\left({\displaystyle \frac{\pi }{x^2}}\right)-5=0$	$\xi \approx -1.41421$
$f_{10}\left(x\right)=\sqrt{x^4}+sin\left({\displaystyle \frac{\pi }{x^2}}\right)-{\displaystyle \frac{3}{16}}=0$	$\xi \approx -0.90599$

.

We aim to demonstrate the theoretical results by testing the MCCTU($\alpha$) family. Specifically, we evaluate three representative members of the family with $\delta =\frac{2}{{\alpha }^4}$ and $\alpha =1$, $\alpha =2$, and $\alpha =100$. Therefore, in all cases, $ϵ=2$.

We conduct two experiments. In the first experiment, we analyze the stability of the MCCTU($\alpha$) family using two of its methods, chosen based on stable and unstable values of the parameter $\alpha$. In the second experiment, we perform an efficiency analysis of the MCCTU($\alpha$) family through a comparative study between its optimal stable member and fifteen different fourth-order methods from the literature: Ostrowski (OS) in [12,35], King (KI) in [35,36], Jarratt (JA) in [35,37], Özban and Kaya (OK1, OK2, OK3) in [8], Chun (CH) in [38], Maheshwari (MA) in [39], Behl et al. (BMM) in [40], Chun et al. (CLND1, CLND2) in [41], Artidiello et al. (ACCT1, ACCT2) in [42], Ghanbari (GH) in [43], and Kou et al. (KLW) in [44].

While performing these numerical tests, we start the iterations with different initial estimates: close ($x_0\approx \xi$), far ($x_0\approx 3\xi$), and very far ($x_0\approx 10\xi$) from the root $\xi$. This approach allows us to evaluate how sensitive the methods are to the initial estimation when finding a solution.

The calculations are performed using the MATLAB R2020b programming package with variable precision arithmetic set to 200 digits of mantissa (in Appendix C, an example with double-precision arithmetics is included). For each method, we analyze the number of iterations (iter) required to converge to the solution, with stopping criteria defined as $|x_{k+1}-x_k|<{10}^{-100}$ or $|f\left(x_{k+1}\right)|<{10}^{-100}$. Here, $|x_{k+1}-x_k|$ represents the error estimation between two consecutive iterations, and $\left|f\right(x_{k+1}\left)\right|$ is the residual error of the nonlinear test function.

To check the theoretical order of convergence (p), we calculate the approximate computational order of convergence (ACOC) as described by Cordero and Torregrosa in [15]. In the numerical results, if the ACOC values do not stabilize throughout the iterative process, it is marked as ‘-’; and if any method fails to converge within a maximum of 50 iterations, it is marked as ‘nc’.

4.1. First Experiment: Stability Analysis of MCCTU($\alpha$) Family

In this experiment, we conducted a stability analysis of the MCCTU($\alpha$) family by considering values of $\alpha$ both within the stability regions ($\alpha =2$) and outside of them ($\alpha =100$), setting $\delta =\frac{2}{{\alpha }^4}$. The methods analyzed are of order 3, consistent with the theoretical convergence results. A special case occurs when $\alpha =0$, where the associated method never converges to the solution because the denominator in the relation $\delta =2/{\alpha }^4$ becomes zero, causing $\delta$ to grow indefinitely.

The numerical performance of the iterative methods MCCTU(2) and MCCTU(100) is presented in

Table 2. Numerical performance of MCCTU(2) method on nonlinear equations (“nc” means non-convergence).

Function	${\mathit{x}}_0$	$|{\mathit{x}}_{\mathit{k}+1}-{\mathit{x}}_{\mathit{k}}|$	$\left|\mathit{f}\right({\mathit{x}}_{\mathit{k}+1}\left)\right|$	Iter	ACOC
	Close to $\xi$
$f_1$	−0.6	2.2252 $\times {10}^{-54}$	1.4765 $\times {10}^{-162}$	4	3
$f_2$	0.2	1.8447 $\times {10}^{-50}$	1.3536 $\times {10}^{-150}$	4	3
$f_3$	0.6	2.3846 $\times {10}^{-44}$	1.4235 $\times {10}^{-131}$	4	3
$f_4$	−14.1	5.1414 $\times {10}^{-36}$	3.3633 $\times {10}^{-111}$	3	3
$f_5$	1.3	1.6295 $\times {10}^{-53}$	4.3267 $\times {10}^{-159}$	4	3
$f_6$	0.1	4.6096 $\times {10}^{-78}$	2.4334 $\times {10}^{-208}$	4	3
$f_7$	−1.2	3.6237 $\times {10}^{-54}$	1.9349 $\times {10}^{-159}$	4	3
$f_8$	2.3	3.0861 $\times {10}^{-54}$	6.9791 $\times {10}^{-162}$	4	3
$f_9$	−1.4	7.0858 $\times {10}^{-51}$	3.8746 $\times {10}^{-150}$	4	3
$f_{10}$	−0.9	8.9456 $\times {10}^{-45}$	9.7874 $\times {10}^{-131}$	4	3
	Far from $\xi$
$f_1$	−1.8	1.5223  $\times {10}^{-92}$	0	5	3
$f_2$	0.6	6.6012 $\times {10}^{-87}$	0	5	3
$f_3$	1.8	3.8851 $\times {10}^{-45}$	6.1565 $\times {10}^{-134}$	6	3
$f_4$	−42.3	nc	nc	nc	nc
$f_5$	3.9	1.0792 $\times {10}^{-59}$	1.2569 $\times {10}^{-177}$	6	3
$f_6$	0.3	1.0805 $\times {10}^{-48}$	2.6855 $\times {10}^{-146}$	4	3
$f_7$	−3.6	2.2394 $\times {10}^{-55}$	4.5662 $\times {10}^{-163}$	14	3
$f_8$	6.9	1.1722 $\times {10}^{-41}$	3.8248 $\times {10}^{-124}$	6	3
$f_9$	−4.2	1.3408 $\times {10}^{-101}$	0	8	3
$f_{10}$	−2.7	4.3149 $\times {10}^{-78}$	3.1147 $\times {10}^{-207}$	8	3
	Very far from $\xi$
$f_1$	−6.0	1.5491 $\times {10}^{-52}$	4.9812 $\times {10}^{-157}$	6	3
$f_2$	2.0	1.6192 $\times {10}^{-89}$	0	6	3
$f_3$	6.0	7.1447 $\times {10}^{-57}$	3.8290 $\times {10}^{-169}$	10	3
$f_4$	−141.0	nc	nc	nc	nc
$f_5$	13.0	1.6531 $\times {10}^{-82}$	0	8	3
$f_6$	1.0	1.6423 $\times {10}^{-56}$	9.4291 $\times {10}^{-170}$	5	3
$f_7$	−12.0	nc	nc	nc	nc
$f_8$	23.0	1.2648 $\times {10}^{-44}$	4.8043 $\times {10}^{-133}$	7	3
$f_9$	−14.0	2.3358 $\times {10}^{-43}$	1.3880 $\times {10}^{-127}$	10	3
$f_{10}$	−9.0	3.0298 $\times {10}^{-44}$	1.2080 $\times {10}^{-128}$	6	3

and

Table 3. Numerical performance of MCCTU(100) method on nonlinear equations (“nc” means non-convergence).

Function	${\mathit{x}}_0$	$|{\mathit{x}}_{\mathit{k}+1}-{\mathit{x}}_{\mathit{k}}|$	$\left|\mathit{f}\right({\mathit{x}}_{\mathit{k}+1}\left)\right|$	Iter	ACOC
	Close to $\xi$
$f_1$	−0.6	6.1808 $\times {10}^{-99}$	7.6768 $\times {10}^{-113}$	9	-
$f_2$	0.2	2.1827 $\times {10}^{-88}$	4.9309 $\times {10}^{-102}$	9	-
$f_3$	0.6	6.0791 $\times {10}^{-94}$	8.8104 $\times {10}^{-108}$	9	-
$f_4$	−14.1	4.5379 $\times {10}^{-95}$	1.3573 $\times {10}^{-111}$	8	-
$f_5$	1.3	4.9631 $\times {10}^{-94}$	4.8998 $\times {10}^{-107}$	9	-
$f_6$	0.1	3.0953 $\times {10}^{-100}$	1.4092 $\times {10}^{-113}$	9	-
$f_7$	−1.2	8.7126 $\times {10}^{-95}$	1.0578 $\times {10}^{-107}$	9	-
$f_8$	2.3	2.1622 $\times {10}^{-95}$	3.1373 $\times {10}^{-109}$	9	-
$f_9$	−1.4	4.0458 $\times {10}^{-95}$	2.7366 $\times {10}^{-108}$	9	-
$f_{10}$	−0.9	6.2830 $\times {10}^{-95}$	3.1368 $\times {10}^{-108}$	9	-
	Far from $\xi$
$f_1$	−1.8	2.7746 $\times {10}^{-92}$	3.4462 $\times {10}^{-106}$	10	-
$f_2$	0.6	6.8191 $\times {10}^{-99}$	1.5405 $\times {10}^{-112}$	10	-
$f_3$	1.8	8.0835 $\times {10}^{-90}$	1.1715 $\times {10}^{-103}$	12	-
$f_4$	−42.3	nc	nc	nc	nc
$f_5$	3.9	nc	nc	nc	nc
$f_6$	0.3	4.0669 $\times {10}^{-95}$	1.8516 $\times {10}^{-108}$	9	-
$f_7$	−3.6	nc	nc	nc	nc
$f_8$	6.9	1.5980 $\times {10}^{-88}$	2.3186 $\times {10}^{-102}$	11	-
$f_9$	−4.2	nc	nc	nc	nc
$f_{10}$	−2.7	1.5127 $\times {10}^{-97}$	3.0929 $\times {10}^{-110}$	11	-
	Very far from $\xi$
$f_1$	−6.0	1.2947 $\times {10}^{-94}$	1.6081 $\times {10}^{-108}$	11	-
$f_2$	2.0	3.5429 $\times {10}^{-94}$	8.0036 $\times {10}^{-108}$	11	-
$f_3$	6.0	4.5426 $\times {10}^{-97}$	6.5836 $\times {10}^{-111}$	18	-
$f_4$	−141.0	nc	nc	nc	nc
$f_5$	13.0	nc	nc	nc	nc
$f_6$	1.0	1.4843 $\times {10}^{-94}$	6.7580 $\times {10}^{-108}$	10	-
$f_7$	−12.0	nc	nc	nc	nc
$f_8$	23.0	7.4725 $\times {10}^{-92}$	1.0842 $\times {10}^{-105}$	12	-
$f_9$	−14.0	nc	nc	nc	nc
$f_{10}$	−9.0	6.5629e $\times {10}^{-95}$	3.2765 $\times {10}^{-108}$	12	-

, using initial estimates that are close, far, and very far from the root. This approach enables us to assess the stability and reliability of the methods under various initial conditions.

From the analysis of the first experiment, it is evident that the MCCTU(2) method exhibits robust performance. For initial estimates close to the root ($x_0\approx \xi$), the method consistently converges to the solution with very low errors, achieving convergence in three or four iterations, and the ACOC value stabilizes at 3. For initial estimates that are far ($x_0\approx 3\xi$), the number of iterations increases, but the method still converges to the solution in nine out of ten cases. For initial estimates that are very far ($x_0\approx 10\xi$), the method holds a similar performance, converging to the solution in eight out of ten cases. It is notable that as the initial condition moves further away, the method shows a slight difficulty in finding the solution. This slight dependence is understandable given the complexity of the nonlinear functions $f_4$ and $f_7$. Nonetheless, the method is shown to be stable and robust, with a convergence order of 3, verifying the theoretical results.

On the other hand, MCCTU(100) method encounters significant difficulties in finding the solution. As the initial conditions move further away, the number of iterations increases. Despite lacking good stability characteristics, the method converges to the solution for initial estimates close to the root. However, for initial estimates that are far and very far from the root, it fails to converge in four out of ten cases. Additionally, the method never stabilizes the ACOC value in any case. These results confirm the theoretical instability of the method, as $\alpha =100$ lies outside the stability surface studied in Section 3.

4.2. Second Experiment: Efficiency Analysis of MCCTU($\alpha$) Family

In this experiment, we conduct a comparative study between an optimal method of the MCCTU($\alpha$) family and the fifteen fourth-order methods mentioned in the introduction of Section 4, to contrast their numerical performances on nonlinear equations. We consider the method associated with $\alpha =1$ and $\delta =2$, denoted as MCCTU(1), as the optimal stable member of the MCCTU($\alpha$) family with fourth-order of convergence.

Thus, in 

Table 4. Numerical performance of iterative methods on nonlinear equations for $x_0$ close to $\xi$ (1/4).

Function	Method	$|{\mathit{x}}_{\mathit{k}+1}-{\mathit{x}}_{\mathit{k}}|$	$\left|\mathit{f}\right({\mathit{x}}_{\mathit{k}+1}\left)\right|$	Iter	ACOC
$f_1$	MCCTU(1)	8.4069 $\times {10}^{-27}$	2.2344 $\times {10}^{-105}$	3	4.0111
$x_0=$ −0.6	OS	1.2193 $\times {10}^{-29}$	2.7787 $\times {10}^{-117}$	3	4.0062
	KI	3.9435 $\times {10}^{-29}$	3.8183 $\times {10}^{-115}$	3	4.0070
	JA	1.3498 $\times {10}^{-29}$	4.2651 $\times {10}^{-117}$	3	4.0061
	OK1	5.0547 $\times {10}^{-32}$	2.9443 $\times {10}^{-127}$	3	3.9991
	OK2	4.0266 $\times {10}^{-30}$	2.6729 $\times {10}^{-119}$	3	4.0052
	OK3	2.5735 $\times {10}^{-30}$	4.5908 $\times {10}^{-120}$	3	3.9937
	CH	1.6691 $\times {10}^{-28}$	1.6213 $\times {10}^{-112}$	3	4.0081
	MA	3.0371 $\times {10}^{-27}$	3.1217 $\times {10}^{-107}$	3	4.0103
	BMM	1.2299 $\times {10}^{-28}$	4.4824 $\times {10}^{-113}$	3	4.0084
	CLND1	8.643 $\times {10}^{-27}$	2.5116 $\times {10}^{-105}$	3	4.0110
	CLND2	1.6691 $\times {10}^{-28}$	1.6213 $\times {10}^{-112}$	3	4.0081
	ACCT1	8.4069 $\times {10}^{-27}$	2.2344 $\times {10}^{-105}$	3	4.0111
	ACCT2	7.4417 $\times {10}^{-32}$	1.0756 $\times {10}^{-126}$	3	4.0294
	GH	1.9739 $\times {10}^{-26}$	8.0112 $\times {10}^{-104}$	3	4.0119
	KLW	8.4441 $\times {10}^{-28}$	1.4567 $\times {10}^{-109}$	3	4.0092
$f_2$	MCCTU(1)	4.0916 $\times {10}^{-36}$	1.3257 $\times {10}^{-144}$	3	3.9624
$x_0=$ 0.2	OS	2.6718 $\times {10}^{-32}$	8.6963 $\times {10}^{-129}$	3	3.9998
	KI	1.7333 $\times {10}^{-32}$	1.4291 $\times {10}^{-129}$	3	3.9987
	JA	1.1553 $\times {10}^{-31}$	4.1074 $\times {10}^{-126}$	3	3.9990
	OK1	2.4295 $\times {10}^{-31}$	9.1464 $\times {10}^{-125}$	3	4.0008
	OK2	1.4863 $\times {10}^{-31}$	1.1754 $\times {10}^{-125}$	3	3.9997
	OK3	1.3844 $\times {10}^{-31}$	8.8054 $\times {10}^{-126}$	3	3.9988
	CH	5.002 $\times {10}^{-32}$	1.2502 $\times {10}^{-127}$	3	3.9969
	MA	2.1425 $\times {10}^{-34}$	1.6464 $\times {10}^{-137}$	3	3.9844
	BMM	5.5838 $\times {10}^{-31}$	2.8585 $\times {10}^{-123}$	3	4.0057
	CLND1	9.7338 $\times {10}^{-34}$	9.6229 $\times {10}^{-135}$	3	3.9830
	CLND2	5.002 $\times {10}^{-32}$	1.2502 $\times {10}^{-127}$	3	3.9969
	ACCT1	4.0916 $\times {10}^{-36}$	1.3257 $\times {10}^{-144}$	3	3.9624
	ACCT2	1.4832 $\times {10}^{-31}$	1.1243 $\times {10}^{-125}$	3	4.0029
	GH	2.5248 $\times {10}^{-38}$	6.6868 $\times {10}^{-154}$	3	3.9675
	KLW	1.8553 $\times {10}^{-33}$	1.2914 $\times {10}^{-133}$	3	3.9925
$f_3$	MCCTU(1)	2.2096 $\times {10}^{-83}$	0	4	4
$x_0=$ 0.6	OS	1.7622 $\times {10}^{-27}$	3.3439 $\times {10}^{-108}$	3	3.9992
	KI	1.6297 $\times {10}^{-100}$	0	4	4
	JA	2.9708 $\times {10}^{-27}$	2.989 $\times {10}^{-107}$	3	3.9996
	OK1	6.1743 $\times {10}^{-100}$	1.9467 $\times {10}^{-208}$	4	4
	OK2	1.0137 $\times {10}^{-33}$	6.7493 $\times {10}^{-135}$	3	4.0975
	OK3	1.0148 $\times {10}^{-27}$	3.9188 $\times {10}^{-110}$	3	4.2357
	CH	2.1262 $\times {10}^{-94}$	0	4	4
	MA	6.9765 $\times {10}^{-86}$	0	4	4
	BMM	1.6076 $\times {10}^{-85}$	1.9467 $\times {10}^{-208}$	4	4
	CLND1	2.5512 $\times {10}^{-83}$	0	4	4
	CLND2	2.1262 $\times {10}^{-94}$	0	4	4
	ACCT1	2.2096 $\times {10}^{-83}$	6.8135 $\times {10}^{-208}$	4	4
	ACCT2	2.5202 $\times {10}^{-91}$	0	4	4
	GH	2.2217 $\times {10}^{-81}$	0	4	4
	KLW	2.4336 $\times {10}^{-89}$	0	4	4

,

Table 5. Numerical performance of iterative methods on nonlinear equations for $x_0$ close to $\xi$ (2/4).

Function	Method	$|{\mathit{x}}_{\mathit{k}+1}-{\mathit{x}}_{\mathit{k}}|$	$\left|\mathit{f}\right({\mathit{x}}_{\mathit{k}+1}\left)\right|$	Iter	ACOC
$f_4$	MCCTU(1)	2.4812 $\times {10}^{-61}$	0	3	4
$x_0=$ −14.1	OS	5.7494 $\times {10}^{-76}$	0	3	4
	KI	2.6178 $\times {10}^{-66}$	0	3	4
	JA	4.5662 $\times {10}^{-69}$	3.8934 $\times {10}^{-208}$	3	4
	OK1	1.6181 $\times {10}^{-64}$	0	3	4
	OK2	1.2341 $\times {10}^{-67}$	0	3	4
	OK3	4.782 $\times {10}^{-68}$	3.8934 $\times {10}^{-208}$	3	3.9998
	CH	4.1273 $\times {10}^{-64}$	0	3	4
	MA	5.9003 $\times {10}^{-62}$	0	3	4
	BMM	2.4555 $\times {10}^{-61}$	3.8934 $\times {10}^{-208}$	3	4
	CLND1	2.8374 $\times {10}^{-61}$	0	3	4
	CLND2	4.1273 $\times {10}^{-64}$	0	3	4
	ACCT1	2.4812 $\times {10}^{-61}$	0	3	4
	ACCT2	7.6144 $\times {10}^{-63}$	0	3	4
	GH	7.562 $\times {10}^{-61}$	0	3	4
	KLW	7.8025 $\times {10}^{-63}$	0	3	4
$f_5$	MCCTU(1)	1.5146 $\times {10}^{-80}$	0	4	4
$x_0=$ 1.3	OS	4.0399 $\times {10}^{-98}$	0	4	4
	KI	4.6142 $\times {10}^{-94}$	0	4	4
	JA	4.0399 $\times {10}^{-98}$	0	4	4
	OK1	1.6263 $\times {10}^{-26}$	3.9339 $\times {10}^{-104}$	3	4.0265
	OK2	3.7251 $\times {10}^{-26}$	1.5538 $\times {10}^{-102}$	3	4.0049
	OK3	2.6244 $\times {10}^{-29}$	4.1697 $\times {10}^{-115}$	3	3.8563
	CH	5.0966 $\times {10}^{-90}$	0	4	4
	MA	8.4188 $\times {10}^{-83}$	0	4	4
	BMM	8.6757 $\times {10}^{-85}$	0	4	4
	CLND1	1.5146 $\times {10}^{-80}$	0	4	4
	CLND2	5.0966 $\times {10}^{-90}$	0	4	4
	ACCT1	1.5146 $\times {10}^{-80}$	0	4	4
	ACCT2	1.0557 $\times {10}^{-91}$	0	4	4
	GH	1.0682 $\times {10}^{-78}$	0	4	4
	KLW	8.3547 $\times {10}^{-86}$	0	4	4
$f_6$	MCCTU(1)	1.1439 $\times {10}^{-32}$	5.0948 $\times {10}^{-129}$	3	3.9969
$x_0=$ 0.1	OS	5.058 $\times {10}^{-36}$	4.1154 $\times {10}^{-143}$	3	3.9980
	KI	2.4554 $\times {10}^{-35}$	3.1386 $\times {10}^{-140}$	3	3.9979
	JA	7.2379 $\times {10}^{-36}$	1.8516 $\times {10}^{-142}$	3	3.9981
	OK1	8.6178 $\times {10}^{-41}$	3.6529 $\times {10}^{-163}$	3	4.0021
	OK2	1.3158 $\times {10}^{-36}$	1.4361 $\times {10}^{-145}$	3	3.9982
	OK3	2.2299 $\times {10}^{-36}$	1.2384 $\times {10}^{-144}$	3	4.0031
	CH	1.6189 $\times {10}^{-34}$	8.6635 $\times {10}^{-137}$	3	3.9977
	MA	3.8478 $\times {10}^{-33}$	5.2367 $\times {10}^{-131}$	3	3.9971
	BMM	7.9902 $\times {10}^{-34}$	7.003 $\times {10}^{-134}$	3	3.9985
	CLND1	1.2375 $\times {10}^{-32}$	7.0855 $\times {10}^{-129}$	3	3.9970
	CLND2	1.6189 $\times {10}^{-34}$	8.6635 $\times {10}^{-137}$	3	3.9977
	ACCT1	1.1439 $\times {10}^{-32}$	5.0948 $\times {10}^{-129}$	3	3.9969
	ACCT2	2.0595 $\times {10}^{-36}$	9.7992 $\times {10}^{-145}$	3	3.9952
	GH	2.7869 $\times {10}^{-32}$	2.1488 $\times {10}^{-127}$	3	3.9967
	KLW	9.5356 $\times {10}^{-34}$	1.49 $\times {10}^{-133}$	3	3.9974

,

Table 6. Numerical performance of iterative methods on nonlinear equations for $x_0$ close to $\xi$ (3/4).

Function	Method	$|{\mathit{x}}_{\mathit{k}+1}-{\mathit{x}}_{\mathit{k}}|$	$\left|\mathit{f}\right({\mathit{x}}_{\mathit{k}+1}\left)\right|$	Iter	ACOC
$f_7$	MCCTU(1)	8.109 $\times {10}^{-29}$	1.224 $\times {10}^{-110}$	3	4.0025
$x_0=$ −1.2	OS	1.0259 $\times {10}^{-36}$	8.588 $\times {10}^{-144}$	3	3.9987
	KI	2.2 $\times {10}^{-33}$	8.266 $\times {10}^{-130}$	3	4.0003
	JA	1.3275 $\times {10}^{-35}$	3.987 $\times {10}^{-139}$	3	3.9993
	OK1	2.8995 $\times {10}^{-32}$	4.1375 $\times {10}^{-125}$	3	4.0011
	OK2	4.5559 $\times {10}^{-36}$	4.3528 $\times {10}^{-141}$	3	4.0014
	OK3	3.3899 $\times {10}^{-35}$	9.9763 $\times {10}^{-138}$	3	4.0602
	CH	1.5282 $\times {10}^{-31}$	4.4541 $\times {10}^{-122}$	3	4.0011
	MA	1.9806 $\times {10}^{-29}$	3.2969 $\times {10}^{-113}$	3	4.0021
	BMM	5.13 $\times {10}^{-29}$	1.8531 $\times {10}^{-111}$	3	3.9988
	CLND1	8.8475 $\times {10}^{-29}$	1.7657 $\times {10}^{-110}$	3	4.0024
	CLND2	1.5282 $\times {10}^{-31}$	4.4541 $\times {10}^{-122}$	3	4.0011
	ACCT1	8.109 $\times {10}^{-29}$	1.224 $\times {10}^{-110}$	3	4.0025
	ACCT2	1.7685 $\times {10}^{-30}$	1.2708 $\times {10}^{-117}$	3	4.0037
	GH	2.4542 $\times {10}^{-28}$	1.2766 $\times {10}^{-108}$	3	4.0029
	KLW	2.7616 $\times {10}^{-30}$	8.4579 $\times {10}^{-117}$	3	4.0014
$f_8$	MCCTU(1)	3.2362 $\times {10}^{-36}$	1.278 $\times {10}^{-144}$	3	4.0010
$x_0=$ 2.3	OS	4.7781 $\times {10}^{-35}$	1.1082 $\times {10}^{-139}$	3	3.9959
	KI	3.867 $\times {10}^{-35}$	4.5395 $\times {10}^{-140}$	3	3.9962
	JA	6.4886 $\times {10}^{-36}$	2.6103 $\times {10}^{-143}$	3	3.9934
	OK1	1.3631 $\times {10}^{-35}$	5.9439 $\times {10}^{-142}$	3	3.9927
	OK2	8.3354 $\times {10}^{-36}$	7.4954 $\times {10}^{-143}$	3	3.9931
	OK3	8.2958 $\times {10}^{-36}$	7.3117 $\times {10}^{-143}$	3	3.9935
	CH	2.8017 $\times {10}^{-36}$	7.5934 $\times {10}^{-145}$	3	3.9943
	MA	7.3689 $\times {10}^{-36}$	4.1437 $\times {10}^{-143}$	3	3.9992
	BMM	2.6822 $\times {10}^{-34}$	1.5979 $\times {10}^{-136}$	3	3.9934
	CLND1	5.1248 $\times {10}^{-38}$	3.528 $\times {10}^{-152}$	3	4.0005
	CLND2	2.8017 $\times {10}^{-36}$	7.5934 $\times {10}^{-145}$	3	3.9943
	ACCT1	3.2362 $\times {10}^{-36}$	1.278 $\times {10}^{-144}$	3	4.0010
	ACCT2	1.2316 $\times {10}^{-34}$	5.9984 $\times {10}^{-138}$	3	3.9946
	GH	1.2035 $\times {10}^{-36}$	1.9404 $\times {10}^{-146}$	3	4.0036
	KLW	1.4928 $\times {10}^{-35}$	8.1719 $\times {10}^{-142}$	3	3.9978
$f_9$	MCCTU(1)	1.2504 $\times {10}^{-28}$	6.0286 $\times {10}^{-111}$	3	3.9982
$x_0=$ −1.4	OS	2.2297 $\times {10}^{-33}$	5.9539 $\times {10}^{-131}$	3	4.0107
	KI	3.571 $\times {10}^{-39}$	4.8453 $\times {10}^{-155}$	3	3.9663
	JA	6.6365 $\times {10}^{-33}$	5.9006 $\times {10}^{-129}$	3	4.0095
	OK1	1.7043 $\times {10}^{-30}$	8.4881 $\times {10}^{-119}$	3	4.0019
	OK2	8.1242 $\times {10}^{-32}$	2.3078 $\times {10}^{-124}$	3	4.0049
	OK3	1.3061 $\times {10}^{-31}$	1.4689 $\times {10}^{-123}$	3	4.0184
	CH	5.6961 $\times {10}^{-33}$	3.922 $\times {10}^{-129}$	3	3.9887
	MA	2.4063 $\times {10}^{-29}$	5.9988 $\times {10}^{-114}$	3	3.9973
	BMM	2.911 $\times {10}^{-28}$	2.1169 $\times {10}^{-109}$	3	3.9971
	CLND1	1.0887 $\times {10}^{-28}$	3.3751 $\times {10}^{-111}$	3	3.9980
	CLND2	5.6961 $\times {10}^{-33}$	3.922 $\times {10}^{-129}$	3	3.9887
	ACCT1	1.2504 $\times {10}^{-28}$	6.0286 $\times {10}^{-111}$	3	3.9982
	ACCT2	1.7434 $\times {10}^{-29}$	1.473 $\times {10}^{-114}$	3	4.0025
	GH	4.3546 $\times {10}^{-28}$	1.1301 $\times {10}^{-108}$	3	3.9989
	KLW	2.0248 $\times {10}^{-30}$	1.8702 $\times {10}^{-118}$	3	3.9955

,

Table 7. Numerical performance of iterative methods on nonlinear equations for $x_0$ close to $\xi$ (4/4).

Function	Method	$|{\mathit{x}}_{\mathit{k}+1}-{\mathit{x}}_{\mathit{k}}|$	1c$\left|\mathit{f}\right({\mathit{x}}_{\mathit{k}+1}\left)\right|$	Iter	ACOC
$f_{10}$	MCCTU(1)	1.3096 $\times {10}^{-27}$	1.0557 $\times {10}^{-105}$	3	4.0263
$x_0=$ −0.9	OS	1.2157 $\times {10}^{-28}$	5.2236 $\times {10}^{-110}$	3	4.0178
	KI	1.6268 $\times {10}^{-28}$	1.7588 $\times {10}^{-109}$	3	4.0189
	JA	2.5808 $\times {10}^{-28}$	1.2592 $\times {10}^{-108}$	3	4.0158
	OK1	1.2566 $\times {10}^{-28}$	6.3023 $\times {10}^{-110}$	3	4.0126
	OK2	2.0733 $\times {10}^{-28}$	5.0608 $\times {10}^{-109}$	3	4.0149
	OK3	1.9638 $\times {10}^{-28}$	4.0898 $\times {10}^{-109}$	3	4.0133
	CH	4.7545 $\times {10}^{-28}$	1.6033 $\times {10}^{-107}$	3	4.0184
	MA	7.9356 $\times {10}^{-28}$	1.3045 $\times {10}^{-106}$	3	4.0246
	BMM	1.3934 $\times {10}^{-30}$	4.5027 $\times {10}^{-118}$	3	3.9969
	CLND1	2.1208 $\times {10}^{-27}$	8.164 $\times {10}^{-105}$	3	4.0242
	CLND2	4.7545 $\times {10}^{-28}$	1.6033 $\times {10}^{-107}$	3	4.0184
	ACCT1	1.3096 $\times {10}^{-27}$	1.0557 $\times {10}^{-105}$	3	4.0263
	ACCT2	1.9256 $\times {10}^{-29}$	2.4651 $\times {10}^{-113}$	3	4.0090
	GH	2.0723 $\times {10}^{-27}$	7.171 $\times {10}^{-105}$	3	4.0278
	KLW	4.546 $\times {10}^{-28}$	1.2771 $\times {10}^{-107}$	3	4.0226

,

Table 8. Numerical performance of iterative methods on nonlinear equations for $x_0$ far from $\xi$ (“nc” means non-convergence) (1/4).

Function	Method	$|{\mathit{x}}_{\mathit{k}+1}-{\mathit{x}}_{\mathit{k}}|$	$\left|\mathit{f}\right({\mathit{x}}_{\mathit{k}+1}\left)\right|$	Iter	ACOC
$f_1$	MCCTU(1)	3.15 $\times {10}^{-28}$	4.4044 $\times {10}^{-111}$	4	3.9913
$x_0=$ −1.8	OS	5.8375 $\times {10}^{-36}$	1.46 $\times {10}^{-142}$	4	3.9979
	KI	1.5765 $\times {10}^{-34}$	9.7538 $\times {10}^{-137}$	4	3.9972
	JA	5.3832 $\times {10}^{-35}$	1.0789 $\times {10}^{-138}$	4	3.9976
	OK1	9.9392 $\times {10}^{-40}$	4.4016 $\times {10}^{-158}$	4	4.0001
	OK2	2.4525 $\times {10}^{-36}$	3.6785 $\times {10}^{-144}$	4	3.9982
	OK3	1.1878 $\times {10}^{-33}$	2.0829 $\times {10}^{-133}$	4	4.0017
	CH	2.7866 $\times {10}^{-32}$	1.2595 $\times {10}^{-127}$	4	3.9958
	MA	2.0068 $\times {10}^{-29}$	5.9514 $\times {10}^{-116}$	4	3.9929
	BMM	1.3126 $\times {10}^{-31}$	5.814 $\times {10}^{-125}$	4	4.0050
	CLND1	7.4349 $\times {10}^{-28}$	1.3753 $\times {10}^{-109}$	4	3.9907
	CLND2	2.7866 $\times {10}^{-32}$	1.2595 $\times {10}^{-127}$	4	3.9958
	ACCT1	3.15 $\times {10}^{-28}$	4.4044 $\times {10}^{-111}$	4	3.9913
	ACCT2	6.7276 $\times {10}^{-44}$	7.1842 $\times {10}^{-175}$	4	3.9954
	GH	2.7189 $\times {10}^{-27}$	2.8837 $\times {10}^{-107}$	4	3.9896
	KLW	7.9363 $\times {10}^{-31}$	1.1367 $\times {10}^{-121}$	4	3.9946
$f_2$	MCCTU(1)	6.8509 $\times {10}^{-86}$	0	4	4
$x_0=$ 0.6	OS	7.8707 $\times {10}^{-82}$	0	4	4
	KI	4.1628 $\times {10}^{-82}$	0	4	4
	JA	5.9451 $\times {10}^{-78}$	7.7869 $\times {10}^{-208}$	4	4
	OK1	1.7391 $\times {10}^{-77}$	7.7869 $\times {10}^{-208}$	4	4
	OK2	8.4717 $\times {10}^{-78}$	0	4	4
	OK3	8.9827 $\times {10}^{-78}$	0	4	4
	CH	1.8951 $\times {10}^{-78}$	0	4	4
	MA	1.8733 $\times {10}^{-84}$	0	4	4
	BMM	1.2206 $\times {10}^{-79}$	0	4	4
	CLND1	2.1212 $\times {10}^{-80}$	0	4	4
	CLND2	1.8951 $\times {10}^{-78}$	0	4	4
	ACCT1	6.8509 $\times {10}^{-86}$	0	4	4
	ACCT2	1.366 $\times {10}^{-80}$	0	4	4
	GH	1.4879 $\times {10}^{-88}$	0	4	4
	KLW	2.1154 $\times {10}^{-83}$	0	4	4
$f_3$	MCCTU(1)	6.0868 $\times {10}^{-31}$	6.9016 $\times {10}^{-121}$	5	3.9978
$x_0=$ 1.8	OS	7.2812 $\times {10}^{-73}$	0	5	4
	KI	8.2846 $\times {10}^{-59}$	0	5	4
	JA	6.0259 $\times {10}^{-71}$	0	5	4
	OK1	1.1879 $\times {10}^{-82}$	0	5	4
	OK2	6.2111 $\times {10}^{-27}$	9.5138 $\times {10}^{-108}$	4	4.1522
	OK3	7.5783 $\times {10}^{-53}$	0	5	4.0205
	CH	9.6923 $\times {10}^{-49}$	1.3715 $\times {10}^{-192}$	5	3.9999
	MA	1.3275 $\times {10}^{-34}$	1.1979 $\times {10}^{-135}$	5	3.9989
	BMM	nc	nc	nc	nc
	CLND1	9.5034 $\times {10}^{-31}$	4.1315 $\times {10}^{-120}$	5	3.9978
	CLND2	9.6923 $\times {10}^{-49}$	1.3715 $\times {10}^{-192}$	5	3.9999
	ACCT1	6.0868 $\times {10}^{-31}$	6.9016 $\times {10}^{-121}$	5	3.9978
	ACCT2	6.1271 $\times {10}^{-31}$	2.8102 $\times {10}^{-121}$	4	3.9953
	GH	1.4039 $\times {10}^{-28}$	2.4077 $\times {10}^{-111}$	5	3.9965
	KLW	4.5965 $\times {10}^{-39}$	1.1996 $\times {10}^{-153}$	5	3.9996

,

Table 9. Numerical performance of iterative methods on nonlinear equations for $x_0$ far from $\xi$ (“nc” means non-convergence) (2/4).

Function	Method	$|{\mathit{x}}_{\mathit{k}+1}-{\mathit{x}}_{\mathit{k}}|$	$\left|\mathit{f}\right({\mathit{x}}_{\mathit{k}+1}\left)\right|$	Iter	ACOC
$f_4$	MCCTU(1)	nc	nc	nc	nc
$x_0=$ −42.3	OS	2.602 $\times {10}^{-54}$	0	6	4.0004
	KI	nc	nc	nc	nc
	JA	1.0645 $\times {10}^{-51}$	0	6	4
	OK1	nc	nc	nc	nc
	OK2	nc	nc	nc	nc
	OK3	nc	nc	nc	nc
	CH	nc	nc	nc	nc
	MA	nc	nc	nc	nc
	BMM	nc	nc	nc	nc
	CLND1	nc	nc	nc	nc
	CLND2	nc	nc	nc	nc
	ACCT1	nc	nc	nc	nc
	ACCT2	nc	nc	nc	nc
	GH	nc	nc	nc	nc
	KLW	nc	nc	nc	nc
$f_5$	MCCTU(1)	5.1192 $\times {10}^{-33}$	6.3445 $\times {10}^{-129}$	5	3.9976
$x_0=$ 3.9	OS	1.8922 $\times {10}^{-60}$	0	5	4
	KI	1.6925 $\times {10}^{-53}$	0	5	3.9999
	JA	1.8922 $\times {10}^{-60}$	0	5	4
	OK1	4.3746 $\times {10}^{-85}$	0	5	4
	OK2	8.1261 $\times {10}^{-70}$	0	5	4
	OK3	8.7491 $\times {10}^{-49}$	5.1503 $\times {10}^{-193}$	5	4.0015
	CH	1.68 $\times {10}^{-47}$	2.7094 $\times {10}^{-187}$	5	3.9998
	MA	3.351 $\times {10}^{-36}$	9.1961 $\times {10}^{-142}$	5	3.9986
	BMM	nc	nc	nc	nc
	CLND1	5.1192 $\times {10}^{-33}$	6.3445 $\times {10}^{-129}$	5	3.9976
	CLND2	1.68 $\times {10}^{-47}$	2.7094 $\times {10}^{-187}$	5	3.9998
	ACCT1	5.1192 $\times {10}^{-33}$	6.3445 $\times {10}^{-129}$	5	3.9976
	ACCT2	1.7037 $\times {10}^{-75}$	0	5	4
	GH	8.0066 $\times {10}^{-31}$	4.5963 $\times {10}^{-120}$	5	3.9964
	KLW	5.3477 $\times {10}^{-40}$	4.373 $\times {10}^{-157}$	5	3.9993
$f_6$	MCCTU(1)	2.8249 $\times {10}^{-77}$	1.2167 $\times {10}^{-208}$	4	4
$x_0=$ 0.3	OS	4.615 $\times {10}^{-92}$	1.2167 $\times {10}^{-208}$	4	4
	KI	4.375 $\times {10}^{-89}$	1.2167 $\times {10}^{-208}$	4	4
	JA	1.7544 $\times {10}^{-91}$	1.2167 $\times {10}^{-208}$	4	4
	OK1	3.5822 $\times {10}^{-29}$	1.0907 $\times {10}^{-116}$	3	3.9593
	OK2	1.0602 $\times {10}^{-94}$	1.2167 $\times {10}^{-208}$	4	4
	OK3	1.3907 $\times {10}^{-101}$	1.2167 $\times {10}^{-208}$	4	4
	CH	1.6778 $\times {10}^{-85}$	1.2167 $\times {10}^{-208}$	4	4
	MA	2.2127 $\times {10}^{-79}$	1.2167 $\times {10}^{-208}$	4	4
	BMM	3.3893 $\times {10}^{-83}$	1.2167 $\times {10}^{-208}$	4	4
	CLND1	3.6933 $\times {10}^{-77}$	1.2167 $\times {10}^{-208}$	4	4
	CLND2	1.6778 $\times {10}^{-85}$	1.2167 $\times {10}^{-208}$	4	4
	ACCT1	2.8249 $\times {10}^{-77}$	2.4334 $\times {10}^{-208}$	4	4
	ACCT2	3.1138 $\times {10}^{-91}$	1.2167 $\times {10}^{-208}$	4	4
	GH	1.6004 $\times {10}^{-75}$	1.2167 $\times {10}^{-208}$	4	4
	KLW	3.9889 $\times {10}^{-82}$	1.2167 $\times {10}^{-208}$	4	4

,

Table 10. Numerical performance of iterative methods in nonlinear equations for $x_0$ far from $\xi$ (“nc” means non-convergence) (3/4).

Function	Method	$|{\mathit{x}}_{\mathit{k}+1}-{\mathit{x}}_{\mathit{k}}|$	$\left|\mathit{f}\right({\mathit{x}}_{\mathit{k}+1}\left)\right|$	Iter	ACOC
$f_7$	MCCTU(1)	2.1695 $\times {10}^{-40}$	6.2709 $\times {10}^{-157}$	12	3.9997
$x_0=$ −3.6	OS	8.7445 $\times {10}^{-42}$	4.5328 $\times {10}^{-164}$	9	4.0005
	KI	3.5832 $\times {10}^{-59}$	0	10	4
	JA	1.445 $\times {10}^{-33}$	5.5984 $\times {10}^{-131}$	9	4.0010
	OK1	4.5904 $\times {10}^{-56}$	0	9	4
	OK2	2.2993 $\times {10}^{-100}$	0	9	4
	OK3	4.4822 $\times {10}^{-55}$	0	11	3.9955
	CH	2.5759 $\times {10}^{-93}$	0	11	4
	MA	1.0752 $\times {10}^{-70}$	0	12	4
	BMM	nc	nc	nc	nc
	CLND1	5.1735 $\times {10}^{-38}$	2.0643 $\times {10}^{-147}$	12	3.9995
	CLND2	2.5759 $\times {10}^{-93}$	0	11	4
	ACCT1	2.1695 $\times {10}^{-40}$	6.2709 $\times {10}^{-157}$	12	3.9997
	ACCT2	7.4341 $\times {10}^{-57}$	0	8	4
	GH	2.2938 $\times {10}^{-28}$	9.7412 $\times {10}^{-109}$	12	3.9971
	KLW	6.4489 $\times {10}^{-31}$	2.515 $\times {10}^{-119}$	11	3.9988
$f_8$	MCCTU(1)	8.9717 $\times {10}^{-34}$	7.5485 $\times {10}^{-135}$	5	3.9964
$x_0=$ 6.9	OS	3.5465 $\times {10}^{-42}$	3.3638 $\times {10}^{-168}$	5	3.9988
	KI	3.0788 $\times {10}^{-46}$	1.8241 $\times {10}^{-184}$	5	3.9994
	JA	3.8134 $\times {10}^{-44}$	3.1142 $\times {10}^{-176}$	5	3.9984
	OK1	4.9365 $\times {10}^{-41}$	1.0225 $\times {10}^{-163}$	5	3.9972
	OK2	1.6379 $\times {10}^{-42}$	1.1175 $\times {10}^{-169}$	5	3.9979
	OK3	9.2522 $\times {10}^{-52}$	1.0123 $\times {10}^{-206}$	5	4.0004
	CH	6.7803 $\times {10}^{-74}$	1.5574 $\times {10}^{-207}$	5	4
	MA	4.1752 $\times {10}^{-36}$	4.2707 $\times {10}^{-144}$	5	4.0002
	BMM	1.4528 $\times {10}^{-98}$	1.5574 $\times {10}^{-207}$	5	4
	CLND1	2.6243 $\times {10}^{-36}$	2.4259 $\times {10}^{-145}$	5	3.9960
	CLND2	6.7803 $\times {10}^{-74}$	1.5574 $\times {10}^{-207}$	5	4
	ACCT1	8.9717 $\times {10}^{-34}$	7.5485 $\times {10}^{-135}$	5	3.9964
	ACCT2	7.0924 $\times {10}^{-38}$	6.5962 $\times {10}^{-151}$	5	3.9970
	GH	9.3051 $\times {10}^{-33}$	6.9333 $\times {10}^{-131}$	5	3.9867
	KLW	1.1619 $\times {10}^{-39}$	2.9996 $\times {10}^{-158}$	5	4.0009
$f_9$	MCCTU(1)	1.9461 $\times {10}^{-27}$	3.5371 $\times {10}^{-106}$	6	4.0014
$x_0=$ −4.2	OS	4.9716 $\times {10}^{-100}$	0	6	4
	KI	2.3907 $\times {10}^{-73}$	0	6	4.0003
	JA	1.5785 $\times {10}^{-91}$	0	6	4
	OK1	7.3314 $\times {10}^{-101}$	0	6	4
	OK2	5.2967 $\times {10}^{-94}$	0	6	4
	OK3	3.3512 $\times {10}^{-48}$	6.3655 $\times {10}^{-190}$	6	3.9987
	CH	1.4014 $\times {10}^{-54}$	0	6	4.0003
	MA	2.5559 $\times {10}^{-32}$	7.6354 $\times {10}^{-126}$	6	4.0012
	BMM	nc	nc	nc	nc
	CLND1	1.866 $\times {10}^{-27}$	2.9131 $\times {10}^{-106}$	6	4.0016
	CLND2	1.4014 $\times {10}^{-54}$	0	6	4.0003
	ACCT1	1.9461 $\times {10}^{-27}$	3.5371 $\times {10}^{-106}$	6	4.0014
	ACCT2	7.4372 $\times {10}^{-40}$	4.8779 $\times {10}^{-156}$	5	3.9995
	GH	3.5153 $\times {10}^{-95}$	0	7	4
	KLW	1.7919 $\times {10}^{-38}$	1.147 $\times {10}^{-150}$	6	4.0009

,

Table 11. Numerical performance of iterative methods on nonlinear equations for $x_0$ far from $\xi$ (4/4).

Function	Method	$|{\mathit{x}}_{\mathit{k}+1}-{\mathit{x}}_{\mathit{k}}|$	$\left|\mathit{f}\right({\mathit{x}}_{\mathit{k}+1}\left)\right|$	Iter	ACOC
$f_{10}$	MCCTU(1)	1.4573 $\times {10}^{-79}$	1.0707 $\times {10}^{-207}$	6	3.9998
$x_0=$ −2.7	OS	5.9731 $\times {10}^{-31}$	3.7644 $\times {10}^{-111}$	10	4.0223
	KI	3.5019 $\times {10}^{-84}$	3.3094 $\times {10}^{-207}$	6	4.0006
	JA	2.7724 $\times {10}^{-41}$	1.3193 $\times {10}^{-154}$	5	3.9916
	OK1	4.672 $\times {10}^{-37}$	3.3675 $\times {10}^{-141}$	5	4.0067
	OK2	6.9502 $\times {10}^{-38}$	2.3984 $\times {10}^{-144}$	5	3.9545
	OK3	4.1808 $\times {10}^{-48}$	2.2923 $\times {10}^{-179}$	5	4.0014
	CH	1.8136 $\times {10}^{-97}$	3.8389 $\times {10}^{-205}$	6	4
	MA	2.2689 $\times {10}^{-93}$	3.8934 $\times {10}^{-208}$	5	4
	BMM	2.6823 $\times {10}^{-41}$	6.1829 $\times {10}^{-161}$	6	3.9999
	CLND1	8.2625 $\times {10}^{-101}$	2.7254 $\times {10}^{-207}$	5	4
	CLND2	1.8136 $\times {10}^{-97}$	3.8389 $\times {10}^{-205}$	6	4
	ACCT1	1.4573 $\times {10}^{-79}$	1.0707 $\times {10}^{-207}$	6	3.9998
	ACCT2	3.3917 $\times {10}^{-58}$	2.3908 $\times {10}^{-204}$	6	4.0003
	GH	2.2337 $\times {10}^{-72}$	1.0707 $\times {10}^{-207}$	5	3.9994
	KLW	2.4673 $\times {10}^{-32}$	6.302 $\times {10}^{-122}$	4	3.8791

,

Table 12. Numerical performance of iterative methods on nonlinear equations for $x_0$ very far from $\xi$ (“nc” means non-convergence) (1/4).

Function	Method	$|{\mathit{x}}_{\mathit{k}+1}-{\mathit{x}}_{\mathit{k}}|$	$\left|\mathit{f}\right({\mathit{x}}_{\mathit{k}+1}\left)\right|$	Iter	ACOC
$f_1$	MCCTU(1)	3.9494 $\times {10}^{-95}$	0	6	4
$x_0=$ −6.0	OS	7.5454 $\times {10}^{-40}$	4.0753 $\times {10}^{-158}$	5	3.9989
	KI	2.2846 $\times {10}^{-36}$	4.3008 $\times {10}^{-144}$	5	3.9980
	JA	4.2789 $\times {10}^{-41}$	4.3067 $\times {10}^{-163}$	5	3.9992
	OK1	2.8437 $\times {10}^{-53}$	0	5	4
	OK2	2.0962 $\times {10}^{-46}$	1.9632 $\times {10}^{-184}$	5	3.9997
	OK3	1.9553 $\times {10}^{-33}$	1.5296 $\times {10}^{-132}$	5	4.0018
	CH	1.6161 $\times {10}^{-33}$	1.4248 $\times {10}^{-132}$	5	3.9966
	MA	2.632 $\times {10}^{-26}$	1.7609 $\times {10}^{-103}$	5	3.9875
	BMM	nc	nc	nc	nc
	CLND1	5.8255 $\times {10}^{-95}$	0	6	4
	CLND2	1.6161 $\times {10}^{-33}$	1.4248 $\times {10}^{-132}$	5	3.9966
	ACCT1	3.9494 $\times {10}^{-95}$	0	6	4
	ACCT2	5.5395 $\times {10}^{-58}$	0	5	3.9996
	GH	6.2374 $\times {10}^{-89}$	0	6	4
	KLW	1.0186 $\times {10}^{-28}$	3.0849 $\times {10}^{-113}$	5	3.9921
$f_2$	MCCTU(1)	1.0368 $\times {10}^{-34}$	5.4646 $\times {10}^{-139}$	4	4.0222
$x_0=$ 2.0	OS	2.3862 $\times {10}^{-93}$	0	5	4
	KI	1.2873 $\times {10}^{-25}$	4.3485 $\times {10}^{-102}$	4	3.9933
	JA	4.1797 $\times {10}^{-95}$	0	5	4
	OK1	8.2892 $\times {10}^{-82}$	0	5	4
	OK2	5.051 $\times {10}^{-87}$	0	5	4
	OK3	4.2138 $\times {10}^{-33}$	7.557 $\times {10}^{-132}$	4	3.9991
	CH	1.639 $\times {10}^{-29}$	1.4412 $\times {10}^{-117}$	4	3.9949
	MA	9.1807 $\times {10}^{-43}$	5.5512 $\times {10}^{-171}$	4	4.0028
	BMM	2.9675 $\times {10}^{-52}$	0	6	3.9998
	CLND1	9.0974 $\times {10}^{-32}$	7.3425 $\times {10}^{-127}$	4	4.0155
	CLND2	1.639 $\times {10}^{-29}$	1.4412 $\times {10}^{-117}$	4	3.9949
	ACCT1	1.0368 $\times {10}^{-34}$	5.4646 $\times {10}^{-139}$	4	4.0222
	ACCT2	1.9358 $\times {10}^{-73}$	0	5	4
	GH	6.4242 $\times {10}^{-33}$	2.803 $\times {10}^{-132}$	4	4.0791
	KLW	1.6753 $\times {10}^{-39}$	8.5848 $\times {10}^{-158}$	4	3.9976
$f_3$	MCCTU(1)	1.0037 $\times {10}^{-74}$	0	9	4
$x_0=$ 6.0	OS	1.4888 $\times {10}^{-45}$	1.7033 $\times {10}^{-180}$	7	4
	KI	3.4193 $\times {10}^{-27}$	1.1139 $\times {10}^{-106}$	7	3.9978
	JA	2.4788 $\times {10}^{-42}$	1.4488 $\times {10}^{-167}$	7	4
	OK1	7.5531 $\times {10}^{-64}$	0	7	4
	OK2	5.2399 $\times {10}^{-91}$	0	7	4
	OK3	2.024 $\times {10}^{-56}$	0	8	4.0130
	CH	3.2307 $\times {10}^{-66}$	6.8135 $\times {10}^{-208}$	8	4
	MA	2.7511 $\times {10}^{-26}$	2.21 $\times {10}^{-102}$	8	3.9951
	BMM	nc	nc	nc	nc
	CLND1	5.7924 $\times {10}^{-73}$	0	9	4
	CLND2	3.2307 $\times {10}^{-66}$	6.8135 $\times {10}^{-208}$	8	4
	ACCT1	1.0037 $\times {10}^{-74}$	0	9	4
	ACCT2	1.643 $\times {10}^{-31}$	1.4531 $\times {10}^{-123}$	6	4.0042
	GH	1.8548 $\times {10}^{-60}$	0	9	4
	KLW	1.1156 $\times {10}^{-35}$	4.1633 $\times {10}^{-140}$	8	3.9992

,

Table 13. Numerical performance of iterative methods on nonlinear equations for $x_0$ very far from $\xi$ (“nc” means non-convergence) (2/4).

Function	Method	$|{\mathit{x}}_{\mathit{k}+1}-{\mathit{x}}_{\mathit{k}}|$	$\left|\mathit{f}\right({\mathit{x}}_{\mathit{k}+1}\left)\right|$	Iter	ACOC
$f_4$	MCCTU(1)	nc	nc	nc	nc
$x_0=$ −141.0	OS	nc	nc	nc	nc
	KI	nc	nc	nc	nc
	JA	nc	nc	nc	nc
	OK1	nc	nc	nc	nc
	OK2	nc	nc	nc	nc
	OK3	nc	nc	nc	nc
	CH	nc	nc	nc	nc
	MA	nc	nc	nc	nc
	BMM	nc	nc	nc	nc
	CLND1	nc	nc	nc	nc
	CLND2	nc	nc	nc	nc
	ACCT1	nc	nc	nc	nc
	ACCT2	nc	nc	nc	nc
	GH	nc	nc	nc	nc
	KLW	nc	nc	nc	nc
$f_5$	MCCTU(1)	1.2254 $\times {10}^{-58}$	0	7	4
$x_0=$ 13.0	OS	4.3174 $\times {10}^{-43}$	5.0572 $\times {10}^{-170}$	6	3.9996
	KI	5.4113 $\times {10}^{-35}$	1.9154 $\times {10}^{-137}$	6	3.9985
	JA	4.3174 $\times {10}^{-43}$	5.0572 $\times {10}^{-170}$	6	3.9996
	OK1	1.2884 $\times {10}^{-87}$	0	6	4
	OK2	1.1488 $\times {10}^{-54}$	0	6	4
	OK3	1.3547 $\times {10}^{-27}$	2.9602 $\times {10}^{-108}$	6	4.0338
	CH	1.1569 $\times {10}^{-28}$	6.0929 $\times {10}^{-112}$	6	3.9948
	MA	2.1036 $\times {10}^{-69}$	6.2295 $\times {10}^{-207}$	7	4
	BMM	nc	nc	nc	nc
	CLND1	1.2254 $\times {10}^{-58}$	0	7	4
	CLND2	1.1569 $\times {10}^{-28}$	6.0929 $\times {10}^{-112}$	6	3.9948
	ACCT1	1.2254 $\times {10}^{-58}$	0	7	4
	ACCT2	7.7193 $\times {10}^{-68}$	0	6	4
	GH	6.6848 $\times {10}^{-52}$	2.2302 $\times {10}^{-204}$	7	3.9999
	KLW	1.4904 $\times {10}^{-82}$	6.2295 $\times {10}^{-207}$	7	4
$f_6$	MCCTU(1)	1.4106 $\times {10}^{-84}$	1.2167 $\times {10}^{-208}$	5	4
$x_0=$ 1.0	OS	1.0011 $\times {10}^{-40}$	6.3155 $\times {10}^{-162}$	4	3.9991
	KI	6.4033 $\times {10}^{-37}$	1.4516 $\times {10}^{-146}$	4	3.9984
	JA	2.3288 $\times {10}^{-40}$	1.9843 $\times {10}^{-160}$	4	3.9991
	OK1	1.8287 $\times {10}^{-53}$	1.2167 $\times {10}^{-208}$	4	3.9997
	OK2	1.2544 $\times {10}^{-44}$	1.1864 $\times {10}^{-177}$	4	3.9995
	OK3	7.2892 $\times {10}^{-32}$	1.4139 $\times {10}^{-126}$	4	3.9897
	CH	9.2308 $\times {10}^{-32}$	9.1584 $\times {10}^{-126}$	4	3.9962
	MA	1.5346 $\times {10}^{-95}$	1.2167 $\times {10}^{-208}$	5	4
	BMM	2.958 $\times {10}^{-31}$	1.3155 $\times {10}^{-123}$	4	4.0024
	CLND1	1.5451 $\times {10}^{-84}$	1.2167 $\times {10}^{-208}$	5	4
	CLND2	9.2308 $\times {10}^{-32}$	9.1584 $\times {10}^{-126}$	4	3.9962
	ACCT1	1.4106 $\times {10}^{-84}$	1.2167 $\times {10}^{-208}$	5	4
	ACCT2	5.8091 $\times {10}^{-73}$	1.2167 $\times {10}^{-208}$	5	4
	GH	4.0743 $\times {10}^{-77}$	1.2167 $\times {10}^{-208}$	5	4
	KLW	3.7745 $\times {10}^{-28}$	3.658 $\times {10}^{-111}$	4	3.9931

and

Table 14. Numerical performance of iterative methods on nonlinear equations for $x_0$ very far from $\xi$ (“nc” means non-convergence) (3/4).

Function	Method	$|{\mathit{x}}_{\mathit{k}+1}-{\mathit{x}}_{\mathit{k}}|$	$\left|\mathit{f}\right({\mathit{x}}_{\mathit{k}+1}\left)\right|$	iter	ACOC
$f_7$	MCCTU(1)	nc	nc	nc	nc
$x_0=$ −12.0	OS	nc	nc	nc	nc
	KI	nc	nc	nc	nc
	JA	nc	nc	nc	nc
	OK1	nc	nc	nc	nc
	OK2	nc	nc	nc	nc
	OK3	nc	nc	nc	nc
	CH	nc	nc	nc	nc
	MA	nc	nc	nc	nc
	BMM	nc	nc	nc	nc
	CLND1	nc	nc	nc	nc
	CLND2	nc	nc	nc	nc
	ACCT1	nc	nc	nc	nc
	ACCT2	1.2624 $\times {10}^{-39}$	3.2997 $\times {10}^{-154}$	50	4.0007
	GH	nc	nc	nc	nc
	KLW	nc	nc	nc	nc
$f_8$	MCCTU(1)	7.1071 $\times {10}^{-32}$	2.9726 $\times {10}^{-127}$	6	3.9944
$x_0=$ 23.0	OS	3.9961 $\times {10}^{-44}$	5.4218 $\times {10}^{-176}$	6	3.9991
	KI	1.9961 $\times {10}^{-46}$	3.223 $\times {10}^{-185}$	6	3.9995
	JA	2.8208 $\times {10}^{-93}$	1.5574 $\times {10}^{-207}$	6	4
	OK1	1.2245 $\times {10}^{-30}$	3.8716 $\times {10}^{-122}$	5	3.9812
	OK2	1.3604 $\times {10}^{-44}$	5.3174 $\times {10}^{-178}$	5	3.9985
	OK3	1.477 $\times {10}^{-56}$	1.5574 $\times {10}^{-207}$	6	3.9998
	CH	3.6894 $\times {10}^{-61}$	1.5574 $\times {10}^{-207}$	6	3.9999
	MA	5.2575 $\times {10}^{-34}$	1.0738 $\times {10}^{-135}$	6	4.0001
	BMM	3.7259 $\times {10}^{-45}$	3.1549 $\times {10}^{-179}$	7	4.0007
	CLND1	1.0076 $\times {10}^{-36}$	5.271 $\times {10}^{-147}$	6	3.9964
	CLND2	3.6894 $\times {10}^{-61}$	1.5574 $\times {10}^{-207}$	6	3.9999
	ACCT1	7.1071 $\times {10}^{-32}$	2.9726 $\times {10}^{-127}$	6	3.9944
	ACCT2	1.1468 $\times {10}^{-52}$	1.5574 $\times {10}^{-207}$	6	3.9997
	GH	3.8246 $\times {10}^{-31}$	1.9787 $\times {10}^{-124}$	6	3.9800
	KLW	1.5819 $\times {10}^{-37}$	1.0306 $\times {10}^{-149}$	6	4.0012
$f_9$	MCCTU(1)	7.6712 $\times {10}^{-52}$	8.5422 $\times {10}^{-204}$	9	4
$x_0=$ −14.0	OS	2.3748 $\times {10}^{-99}$	0	8	4
	KI	7.9344 $\times {10}^{-69}$	0	8	4.0006
	JA	1.0139 $\times {10}^{-94}$	0	8	4
	OK1	3.2135 $\times {10}^{-50}$	1.0728 $\times {10}^{-197}$	7	3.9999
	OK2	1.5049 $\times {10}^{-32}$	2.7173 $\times {10}^{-127}$	7	3.9952
	OK3	7.1691 $\times {10}^{-34}$	1.3332 $\times {10}^{-132}$	8	3.9822
	CH	8.1205 $\times {10}^{-44}$	1.62 $\times {10}^{-172}$	8	4.0015
	MA	5.7604 $\times {10}^{-70}$	0	9	4
	BMM	nc	nc	nc	nc
	CLND1	9.8978 $\times {10}^{-52}$	2.306 $\times {10}^{-203}$	9	4
	CLND2	8.1205 $\times {10}^{-44}$	1.62 $\times {10}^{-172}$	8	4.0015
	ACCT1	7.6712 $\times {10}^{-52}$	8.5422 $\times {10}^{-204}$	9	4
	ACCT2	7.9723 $\times {10}^{-43}$	6.4405 $\times {10}^{-168}$	7	4.0003
	GH	3.8751 $\times {10}^{-42}$	7.0871 $\times {10}^{-165}$	9	4.0001
	KLW	1.1099 $\times {10}^{-94}$	0	9	4

, we present the numerical results for the sixteen known methods, considering initial estimates that are close, far, and very far from the root, as well as the ten test equations.

In Table 4, Table 5, Table 6 and Table 7, we observe that MCCTU(1) consistently converges to the solution for initial estimates close to the root ($x_0\approx \xi$), with a similar number of iterations as other methods across all equations. The theoretical convergence order is confirmed by the ACOC, which is close to 4. However, what about the dependence of MCCTU(1) on initial estimates? To answer this, we analyze the method for initial estimates far and very far from the solution, specifically for $x_0\approx 3\xi$ and $x_0\approx 10\xi$, respectively. The results are shown in Table 8, Table 9, Table 10, Table 11, Table 12, Table 13, Table 14 and

Table 15. Numerical performance of iterative methods in nonlinear equations for $x_0$ very far from $\xi$ (“nc” means non-convergence) (4/4).

Function	Method	$|{\mathit{x}}_{\mathit{k}+1}-{\mathit{x}}_{\mathit{k}}|$	$\left|\mathit{f}\right({\mathit{x}}_{\mathit{k}+1}\left)\right|$	Iter	ACOC
$f_{10}$	MCCTU(1)	1.2776 $\times {10}^{-70}$	1.0707 $\times {10}^{-207}$	6	4.0008
$x_0=$ −9.0	OS	2.9225 $\times {10}^{-29}$	1.7446 $\times {10}^{-112}$	5	3.9821
	KI	7.9476 $\times {10}^{-94}$	1.5574 $\times {10}^{-207}$	8	4
	JA	6.1519 $\times {10}^{-29}$	1.4803 $\times {10}^{-108}$	8	4.0703
	OK1	2.3282 $\times {10}^{-63}$	9.7336 $\times {10}^{-208}$	6	4
	OK2	2.4314 $\times {10}^{-75}$	1.9467 $\times {10}^{-208}$	7	4
	OK3	1.9369 $\times {10}^{-93}$	4.9203 $\times {10}^{-206}$	6	4
	CH	1.0491 $\times {10}^{-37}$	1.9596 $\times {10}^{-135}$	8	4.0092
	MA	9.5564 $\times {10}^{-89}$	3.8934 $\times {10}^{-208}$	6	4
	BMM	nc	nc	nc	nc
	CLND1	1.5067 $\times {10}^{-28}$	4.9584 $\times {10}^{-107}$	12	3.7072
	CLND2	1.0491 $\times {10}^{-37}$	1.9596 $\times {10}^{-135}$	8	4.0092
	ACCT1	1.2776 $\times {10}^{-70}$	1.0707 $\times {10}^{-207}$	6	4.0008
	ACCT2	1.463 $\times {10}^{-49}$	5.2801 $\times {10}^{-191}$	5	4.0594
	GH	2.4513 $\times {10}^{-45}$	6.3226 $\times {10}^{-174}$	7	4.0044
	KLW	5.7956 $\times {10}^{-39}$	1.4646 $\times {10}^{-150}$	6	4.0153

.

The results presented in Table 8, Table 9, Table 10 and Table 11 are promising. MCCTU(1) converges to the solution in nine out of the ten nonlinear equations, even when the initial estimate is far from the root ($x_0\approx 3\xi$). In these cases, the ACOC consistently stabilizes and approaches 4. Only in one instance, for the function $f_4$, does MCCTU(1) fail to converge, similar to the other thirteen methods. For this particular equation, only two methods successfully approximate the root. In the remaining equations, MCCTU(1) converges to the solution with a comparable number of iterations to other methods and even requires fewer iterations than Ostrowski’s method, as seen with function $f_{10}$. Therefore, we confirm that this method is robust, consistent with the stability results shown in previous sections.

The results presented in Table 12, Table 13, Table 14 and Table 15 confirm the exceptional robustness of the MCCTU(1) method for initial estimates that are very far from the root ($x_0\approx 10\xi$), as the method converges in eight out of ten cases. A slight dependence on the initial estimate is observed for functions $f_4$ and $f_7$, where the method does not converge; however, in these two cases, the other methods also fail to approximate the solution, except for the ACCT2 method, which converges to the root of function $f_7$ with 50 iterations. The complexity of the nonlinear equations plays a significant role in finding their solutions. Moreover, in the cases where the MCCTU(1) method converges to the roots, it does so with a comparable number of iterations to other methods and often with fewer iterations, as seen in function $f_2$. Additionally, for these cases, the ACOC consistently stabilizes at values close to 4.

Therefore, based on the results of the second experiment, we conclude that the MCCTU($\alpha$) family demonstrates impressive numerical performance when using the optimal stable member with $\alpha =1$ as a representative, highlighting its robustness and efficiency even with challenging initial conditions. Overall, the selected MCCTU(1) method exhibits low errors and requires a similar or fewer number of iterations compared to other methods. In certain cases, as the complexity of the nonlinear equation increases, the MCCTU(1) method outperforms Ostrowski’s method and others. The theoretical convergence order is also confirmed by the ACOC, which is always close to 4.

5. Conclusions

The development of the parametric family of multistep iterative schemes MCCTU($\alpha$) based on the damped Newton scheme has proven to be an effective strategy for solving nonlinear equations. The inclusion of an additional Newton step with a weight function and a “frozen” derivative significantly improved the convergence speed from a first-order class to a uniparametric third-order family.

The numerical results confirm the robustness of the MCCTU(2) method for initial estimates close to the root ($x_0\approx \xi$), with very low errors and convergence in three or four iterations. As the initial estimates move further away ($x_0\approx 3\xi$) and ($x_0\approx 10\xi$), the method continues to show solid performance, converging in most cases and confirming its theoretical stability and robustness.

Through the analysis of stability surfaces and dynamical planes, specific members of the MCCTU($\alpha$) family with exceptional stability were identified. These members are particularly suitable for scalar functions with challenging convergence behavior, exhibiting attractive periodic orbits and strange fixed points in their corresponding dynamical planes. The MCCTU(1) member stood out for its optimal and stable performance.

In the comparative analysis, the MCCTU(1) method demonstrated superior numerical performance in many cases, requiring a similar or fewer number of iterations compared to well-established fourth-order methods such as Ostrowski’s method. This superior performance is especially notable in more complex nonlinear equations, where MCCTU(1) outperforms several alternative methods.

The theoretical convergence order of the MCCTU($\alpha$) family was confirmed by calculating the approximate computational order of convergence (ACOC). In most cases, the ACOC value stabilized close to 3, validating the effectiveness and accuracy of the proposed methods both theoretically and practically. Additionally, it was confirmed that the convergence order of the method associated with $\alpha =1$ is optimal, achieving a fourth-order convergence.

Finally, the analysis revealed that certain members of the MCCTU($\alpha$) family, particularly those with $\alpha$ values outside the stability surface, exhibited significant instability. These methods struggled to converge to the solution, especially when initial estimates were far or very far from the root. For instance, the method with $\alpha =100$ failed to stabilize and did not meet the convergence criteria in four out of ten cases. Additionally, the ACOC values for this method did not stabilize, confirming its theoretical instability. This highlights the importance of selecting appropriate parameter values within the stability regions to ensure reliable performance.