An Efficient Class of Weighted-Newton Multiple Root Solvers with Seventh Order Convergence

Abstract

In this work, we construct a family of seventh order iterative methods for finding multiple roots of a nonlinear function. The scheme consists of three steps, of which the first is Newton’s step and last two are the weighted-Newton steps. Hence, the name of the scheme is ‘weighted-Newton methods’. Theoretical results are studied exhaustively along with the main theorem describing convergence analysis. Stability and convergence domain of the proposed class are also demonstrated by means of using a graphical technique, namely, basins of attraction. Boundaries of these basins are fractal like shapes through which basins are symmetric. Efficacy is demonstrated through numerical experimentation on variety of different functions that illustrates good convergence behavior. Moreover, the theoretical result concerning computational efficiency is verified by computing the elapsed CPU time. The overall comparison of numerical results including accuracy and CPU-time shows that the new methods are strong competitors for the existing methods.

1. Introduction

Finding numerically a root of an equation is an interesting and challenging problem. It is also very important in many diverse areas such as Mathematical Biology, Physics, Chemistry, Economics and Engineering, to name a few [1,2,3,4]. This is due to the fact that many problems from these disciplines are ultimately reduced to finding the root of an equation. Researchers are using iterative methods for approximating root since closed form solutions cannot be obtained in general. In particular, here we consider the problem of computing multiple roots of equation $f\left(x\right)=0$ by iterative methods. A root (say, $\alpha$) of $f\left(x\right)=0$ is called multiple root with multiplicity m, if $f^{\left(j\right)}\left(\alpha \right)=0,j=0,1,2,\dots ,m-1$ and $f^{\left(m\right)}\left(\alpha \right)\ne 0$.

A basic and widely used iterative method is the well-known modified Newton’s method

$$x_{n+1}=x_n-m\frac{f\left(x_n\right)}{f^{\prime }\left(x_n\right)}\forall n=0,1,2,\dots .$$

(1)

This method efficiently locates the required multiple root with quadratic order of convergence provided that the initial value $x_0$ is sufficiently close to root [5]. In terms of Traub’s classification (see [1]), Newton’s method (1) is called one-point method. Some other important methods that belong to this class have been developed in [6,7,8,9].

Recently, numerous higher order methods, either independent or based on the modified Newton’s method (1), have been proposed and analyzed in the literature, see e.g., [10,11,12,13,14,15,16,17,18,19,20,21,22,23] and references cited therein. Such methods belong to the category of multipoint methods [1]. Multipoint iterative methods compute new approximations to root $\alpha$ by sampling the function $f\left(x\right)$, and its derivatives at several points of the independent variable, per each step. These methods have the strategy similar to Runge–Kutta methods for solving differential equations and Gaussian quadrature integration rules in the sense that they possess free parameters which can be used to ensure that the convergence speed is of a certain order, and that the sampling is done at some suitable points.

In particular, Geum et al. in [22,23] have proposed two- and three-point Newton-like methods with convergence order six for finding multiple roots. The two-point method [22], applicable for $m>1$, is given as

$$\begin{array}{c}{y}_n=x_n-m\frac{f\left(x_n\right)}{f^{\prime }\left(x_n\right)}\hfill \end{array}$$

$$\begin{array}{c}\hfill x_{n+1}=y_n-Q(u,s)\frac{f\left(y_n\right)}{f^{\prime }\left(y_n\right)}\end{array}$$

(2)

where $u={\left(\frac{f\left(y_n\right)}{f\left(x_n\right)}\right)}^{\frac{1}{m}}$ and $s={\left(\frac{f^{\prime }\left(y_n\right)}{f^{\prime }\left(x_n\right)}\right)}^{\frac{1}{m-1}}$ and $Q:{\mathbb{C}}^2\to \mathbb{C}$ is a holomorphic function in some neighborhood of origin $(0,0)$. The three-point method [23] for $m\ge 1$ is given as

$$\begin{array}{c}{y}_n=x_n-m\frac{f\left(x_n\right)}{f^{\prime }\left(x_n\right)}\hfill \\ z_n=x_n-m{Q}_f\left(u\right)\frac{f\left(x_n\right)}{f^{\prime }\left(x_n\right)}\hfill \\ \hfill x_{n+1}=x_n-m{K}_f(u,v)\frac{f\left(x_n\right)}{f^{\prime }\left(x_n\right)}\end{array}$$

(3)

wherein $u={\left(\frac{f\left(y_n\right)}{f\left(x_n\right)}\right)}^{\frac{1}{m}}$ and $v={\left(\frac{f\left(z_n\right)}{f\left(x_n\right)}\right)}^{\frac{1}{m}}$. The function $Q_f:\mathbb{C}\to \mathbb{C}$ is analytic in a neighborhood of 0 and $K_f:{\mathbb{C}}^2\to \mathbb{C}$ is holomorphic in a neighborhood of $(0,0)$. Both schemes (2) and (3) require four function evaluations to obtain sixth order convergence with the efficiency index (see [24]), $6^{1/4}\approx 1.565$.

The goal and motivation in constructing iterative methods is to attain convergence of order as high as possible by using function evaluations as small as possible. With these considerations, here we propose a family of three-point methods that attain seventh order of convergence for locating multiple roots. The methodology is based on Newton’s and weighted-Newton iterations. The algorithm requires four evaluations of function per iteration and, therefore, possesses the efficiency index $7^{1/4}\approx 1.627$. This shows that the proposed methods have better efficiency $\left(1.627\right)$ than the efficiency $\left(1.565\right)$ of existing methods (2) and (3). Theoretical results concerning convergence order and computational efficiency are verified by performing numerical tests. In the comparison of numerical results with existing techniques, the proposed methods are observed computationally more efficient since they require less computing time (CPU-time) to achieve the solution of required accuracy.

Contents of the article are summarized as follows. In Section 2, we describe the approach to develop new methods and prove their seventh order convergence. In Section 3, stability of the methods is checked by means of using a graphical technique called basins of attraction. In Section 4, some numerical tests are performed to verify the theoretical results by implementing the methods on some examples. Concluding remarks are reported in Section 5.

2. Formulation of Method

Let $m\ge 1$ be the multiplicity of a root of the equation $f\left(x\right)=0$. To compute the root let us consider the following three-step iterative scheme:

$$\left\{\begin{array}{c}{y}_n=x_n-m\frac{f\left(x_n\right)}{f^{\prime }\left(x_n\right)}\hfill \\ z_n=y_n-muH\left(u\right)\frac{f\left(x_n\right)}{f^{\prime }\left(x_n\right)}\hfill \\ x_{n+1}=z_n-mvG(u,w)\frac{f\left(x_n\right)}{f^{\prime }\left(x_n\right)}\hfill \end{array}\right.$$

(4)

where $u={\left(\frac{f\left(y_n\right)}{f\left(x_n\right)}\right)}^{\frac{1}{m}}$, $v={\left(\frac{f\left(z_n\right)}{f\left(x_n\right)}\right)}^{\frac{1}{m}}$, $w={\left(\frac{f\left(z_n\right)}{f\left(y_n\right)}\right)}^{\frac{1}{m}}$, and the function $H:\mathbb{C}\to \mathbb{C}$ is analytic in some neighborhood of 0 and $G:{\mathbb{C}}^2\to \mathbb{C}$ is holomorphic in a neighborhood of $(0,0)$. Notice that the first step is Newton iteration (1) whereas second and third steps are weighted by employing the factors $H\left(u\right)$ and $G(u,w)$, and so we call the algorithm (4) by the name weighted-Newton method. Factors H and G are called weight factors or more appropriately weight functions.

In the sequel we shall find conditions under which the algorithm (4) achieves high convergence order. Thus, the following theorem is stated and proved:

Theorem 1.

Assume that $f:\mathbb{C}\to \mathbb{C}$ is an analytic function in a domain enclosing a root α with multiplicity m. Suppose that initial point $x_0$ is closer enough to the root α, then the iterative formula defined by (4) has seventh order of convergence, if the functions $H\left(u\right)$ and $G(u,w)$ verify the conditions: $H\left(0\right)=1,H^{\prime }\left(0\right)=2$, $H^{″}\left(0\right)=-2$, $G(0,0)=1,G_{10}(0,0)=2$, $G_{01}(0,0)=1,G_{20}(0,0)=0,|H^{‴}\left(0\right)|<\infty \mathrm{and}|G_{11}(0,0)|<\infty$, where $G_{ij}(0,0)=\frac{{\partial }^{i+j}}{\partial u^i\partial w^j}{G(u,w)|}_{(0,0)}$.

Proof. 

Let $e_n=x_n-\alpha$ be the error at n-th iteration. Taking into account that $f^{\left(j\right)}\left(\alpha \right)=0$, $j=0,1,2,\dots ,m-1$, we have by the Taylor’s expansion of $f\left(x_n\right)$ about $\alpha$

$$\begin{array}{cc}\hfill f\left(x_n\right)=& {\displaystyle \frac{f^{\left(m\right)}\left(\alpha \right)}{m!}}{e}_n^m+{\displaystyle \frac{f^{(m+1)}\left(\alpha \right)}{(m+1)!}}{e}_n^{m+1}+{\displaystyle \frac{f^{(m+2)}\left(\alpha \right)}{(m+2)!}}{e}_n^{m+2}+{\displaystyle \frac{f^{(m+3)}\left(\alpha \right)}{(m+3)!}}{e}_n^{m+3}+{\displaystyle \frac{f^{(m+4)}\left(\alpha \right)}{(m+4)!}}{e}_n^{m+4}\hfill \\ & +{\displaystyle \frac{f^{(m+5)}\left(\alpha \right)}{(m+5)!}}{e}_n^{m+5}+{\displaystyle \frac{f^{(m+6)}\left(\alpha \right)}{(m+6)!}}{e}_n^{m+6}+{\displaystyle \frac{f^{(m+7)}\left(\alpha \right)}{(m+7)!}}{e}_n^{m+7}+O\left(e_n^{m+8}\right)\hfill \end{array}$$

or:

$$f\left(x_n\right)={\displaystyle \frac{f^{\left(m\right)}\left(\alpha \right)}{m!}}{e}_n^m\left(1+C_1{e}_n+C_2{e}_n^2+C_3{e}_n^3+C_4{e}_n^4+C_5{e}_n^5+C_6{e}_n^6+C_7{e}_n^7+O\left(e_n^8\right)\right)$$

(5)

where $C_k={\displaystyle \frac{m!}{(m+k)!}}{\displaystyle \frac{f^{(m+k)}\left(\alpha \right)}{f^{\left(m\right)}\left(\alpha \right)}}$ for $k\in \mathbb{N}$.

also

$$\begin{array}{cc}\hfill f^{\prime }\left(x_n\right)=& \frac{f^{\left(m\right)}\left(\alpha \right)}{m!}{e}_n^{m-1}(m+C_1(m+1)e_n\hfill \\ & +C_2(m+2)e_n^2+C_3(m+3)e_n^3+C_4(m+4)e_n^4\hfill \\ & +C_5(m+5)e_n^5+C_6(m+6)e_n^6+C_7(m+7)e_n^7+O\left(e_n^8\right)),\hfill \end{array}$$

(6)

where $C_k=\frac{m!}{(m+k)!}\frac{f^{(m+k)}\left(\alpha \right)}{f^{\left(m\right)}\left(\alpha \right)}$ for $k\in \mathbb{N}$.

Using (5) and (6) in first step of (4), it follows that

$$\begin{array}{cc}\hfill y_n-\alpha =& \frac{C_1}{m}{e}_n^2+\sum _{i=1}^5{\omega }_i{e}_n^{i+2}+O\left(e_n^8\right),\hfill \end{array}$$

(7)

where ${\omega }_i={\omega }_i(m,C_1,C_2,\dots ,C_7)$ are given in terms of $m,C_1,C_2,\dots ,C_7$ with explicitly written two coefficients ${\omega }_1=\frac{2m{C}_2-(m+1)C_1^2}{m^2}$, ${\omega }_2=\frac{1}{m^3}\left(3m^2{C}_3+{(m+1)}^2{C}_1^3-m(4+3m)C_1{C}_2\right)$. Here, rest of the expressions of ${\omega }_i$ are not being produced explicitly since they are very lengthy.

Expansion of $f\left(y_n\right)$ about $\alpha$ yields

$$\begin{array}{cc}\hfill f\left(y_n\right)=& \frac{f^{\left(m\right)}\left(\alpha \right)}{m!}{\left(\frac{C_1}{m}\right)}^m{e}_n^{2m}(1+\frac{2C_{2}m-C_1^2(m+1)}{C_1}{e}_n+\hfill \\ & \frac{1}{2m{C}_1^2}((3+3m+3m^2+m^3)C_1^4-2m(2+3m+2m^2)C_1^2{C}_2\hfill \\ & +4(-1+m)m^2{C}_2^2+6m^2{C}_1{C}_3)e_n^2+\sum _{i=1}^4{\overline{\omega }}_i{e}_n^{i+2}+O\left(e_n^8\right)),\hfill \end{array}$$

(8)

where ${\overline{\omega }}_i={\overline{\omega }}_i(m,C_1,C_2,\dots ,C_7).$

Using (5) and (8) in the expression of u, it follows that

$$\begin{array}{c}\hfill u=\frac{C_1}{m}{e}_n+\frac{2C_{2}m-C_1^2(m+2)}{m^2}{e}_n^2+\sum _{i=1}^5{\eta }_i{e}_n^{i+2}+O\left(e_n^8\right),\end{array}$$

(9)

where ${\eta }_i={\eta }_i(m,C_1,C_2,\dots ,C_7)$ with explicitly written one coefficient ${\eta }_1=\frac{1}{2m^3}\left(C_1^3(2m^2+7m+7)+6C_3{m}^2-2C_2{C}_{1}m(3m+7)\right)$.

Developing weight function $H\left(u\right)$ in neighborhood 0,

$$H\left(u\right)\approx H\left(0\right)+u{H}^{\prime }\left(0\right)+\frac{1}{2!}{u}^2{H}^{″}\left(0\right)+\frac{1}{2!}{u}^3{H}^{‴}\left(0\right).$$

(10)

Inserting Equations (5), (8) and (10) in the second step of (4), after some simplifications we have that

$$\begin{array}{cc}\hfill z_n-\alpha =& -\frac{A}{m}{C}_1{e}_n^2+\frac{1}{m^2}\left(-2mA{C}_2+C_1^2(-1+mA+3H\left(0\right)-H^{\prime }\left(0\right))\right)e_n^3\hfill \\ & +\frac{1}{2m^3}(-6A{m}^2{C}_3+2m{C}_1{C}_2(-4+3Am+11H\left(0\right)-4H^{\prime }\left(0\right))\hfill \\ & +C_1^3(2-2A{m}^2-13H\left(0\right)+10H^{\prime }\left(0\right)\hfill \\ & +m(4-11H\left(0\right)+4H^{\prime }\left(0\right))-H^{″}\left(0\right)\left)\right)e_n^4\hfill \\ & +\sum _{i=1}^3{\gamma }_i{e}_n^{i+4}+O\left(e_n^8\right),\hfill \end{array}$$

(11)

where $A=-1+H\left(0\right)$ and ${\gamma }_i={\gamma }_i(m,H\left(0\right),H^{\prime }\left(0\right),H^{″}\left(0\right),H^{‴}\left(0\right),C_1,C_2,\dots ,C_7)$.

In order to accelerate convergence, the coefficients of $e_n^2$ and $e_n^3$ should be equal to zero. That is possible only if we have

$$H\left(0\right)=1,H^{\prime }\left(0\right)=2.$$

(12)

By using the above values in (11), we obtain that

$$\begin{array}{c}\hfill z_n-\alpha =\frac{-2m{C}_1{C}_2+C_1^3(9+m-H^{″}\left(0\right))}{2m^3}{e}_n^4+\sum _{i=1}^3{\gamma }_i{e}_n^{i+4}+O\left(e_n^8\right).\end{array}$$

(13)

Expansion of $f\left(z_n\right)$ about $\alpha$ yields

$$\begin{array}{c}\hfill f\left(z_n\right)=\frac{f^{\left(m\right)}\left(\alpha \right)}{m!}{(z_n-\alpha )}^m\left(1+C_1(z_n-\alpha )+C_2{(z_n-\alpha )}^2+O\left({(z_n-\alpha )}^3\right)\right).\end{array}$$

(14)

From (5), (8) and (14), we obtain forms of v and w as

$$\begin{array}{c}\hfill v=\frac{(9+m)C_1^3-2m{C}_1{C}_2}{2m^3}{e}_n^3+\sum _{i=1}^4{\tau }_i{e}_n^{i+3}+O\left(e_n^8\right),\end{array}$$

(15)

where ${\tau }_i={\tau }_i(m,H^{″}\left(0\right),H^{‴}\left(0\right),C_1,C_2,\dots ,C_7)$ and

$$\begin{array}{c}\hfill w=\frac{(9+m-H^{″}\left(0\right))C_1^2-2m{C}_2}{2m^3}{e}_n^2+\sum _{i=1}^5{\varsigma }_i{e}_n^{i+2}+O\left(e_n^8\right),\end{array}$$

(16)

where ${\varsigma }_i={\varsigma }_i(m,H^{″}\left(0\right),H^{‴}\left(0\right),C_1,C_2,\dots ,C_7)$.

Expanding $G(u,w)$ in neighborhood of origin $(0,0)$ by Taylor series, it follows that

$$G(u,w)\approx G_{00}(0,0)+u{G}_{10}(0,0)+\frac{1}{2}{u}^2{G}_{20}(0,0)+w(G_{01}(0,0)+u{G}_{11}(0,0)),$$

(17)

where $G_{ij}(0,0)=\frac{{\partial }^{i+j}}{\partial u^i\partial w^j}{G(u,w)|}_{(0,0)}$.

Then by substituting (5), (6), (15)–(17) into the last step of scheme (4), we obtain that

$$\begin{array}{cc}\hfill e_{n+1}=& \frac{1}{2m^3}\left((-1+G_{00}(0,0))C_1\left(2m{C}_1-(9+m-H^{″}\left(0\right))C_1^2\right)\right)e_n^4\hfill \\ & +\sum _{i=1}^3{\xi }_i{e}_n^{i+4}+O\left(e_n^8\right),\hfill \end{array}$$

(18)

where ${\xi }_i={\xi }_i(m,H^{″}\left(0\right),H^{‴}\left(0\right),G_{00}(0,0),G_{10}(0,0),G_{20}(0,0),G_{01}(0,0),G_{11}(0,0),C_1,C_2,\dots ,C_7)$.

From Equation (18) it is clear that we can obtain at least fifth order convergence when $G_{00}(0,0)=1$. In addition, using this value in ${\xi }_1=0$, we will obtain that

$$G_{10}(0,0)=2.$$

(19)

By using $G_{00}=1$ and (19) in ${\xi }_2=0$, the following equation is obtained

$$\begin{array}{c}{C}_1\left(2m{C}_2-C_1^2(9+m-H^{″}\left(0\right))\right)\hfill \\ (-2m{C}_2(-1+G_{01}(0,0))+C_1^2(-11+m(-1+G_{01}(0,0))\hfill \\ -(-9+H^{″}\left(0\right))G_{01}(0,0)+G_{20}(0,0)\left)\right)\hfill \\ =0,\hfill \end{array}$$

(20)

which further yields

$$G_{01}(0,0)=1,G_{20}(0,0)=0\mathrm{and}{H}^{″}\left(0\right)=-2.$$

(21)

Using the above values in (18), we obtain the error equation

$$\begin{array}{cc}\hfill e_{n+1}=& \frac{1}{360m^6}(360m^3\left((39+5m)C_2^3-6m{C}_3^2-10m{C}_2{C}_4\right)+120m^3{C}_1((515+78m)C_2{C}_3\hfill \\ & -12m{C}_5)-60m^2{C}_1^3{C}_3\left(1383+845m+78m^2+12H^{‴}\left(0\right)\right)+10m{C}_1^4{C}_2(21571\hfill \\ & +8183m^2+558m^3+515H^{‴}\left(0\right)+324G_{11}(0,0)+36m\left(667+6H^{‴}\left(0\right)+G_{11}(0,0)\right))\hfill \\ & -60m^2{C}_1^2\left(-6m(55+9m)C_4+C_2^2\left(2619+1546m+135m^2+24H^{‴}\left(0\right)+6G_{11}(0,0)\right)\right)\hfill \\ & -C_1^6(55017+17005m+978m^4+2775H^{‴}\left(0\right)+7290G_{11}(0,0)+15m^2(4463+40H^{‴}\left(0\right)\hfill \\ & +6G_{11}(0,0))+5m(21571+515H^{‴}\left(0\right)+324G_{11}(0,0))\left)\right)e_n^7+O\left(e_n^8\right).\hfill \end{array}$$

(22)

Thus, the seventh order convergence is established. □

Based on the conditions on $H\left(u\right)$ and $G(u,w)$ as shown in Theorem 1, we can generate numerous methods of the family (4). However, we restrict to the following simple forms:

2.1. Some Concrete Forms of $H\left(u\right)$

Case 1.Considering $H\left(u\right)$ a polynomial function, i.e.,

$$H\left(u\right)=A_0+A_{1}u+A_2{u}^2.$$

Using the conditions of Theorem 1, we get $A_0=1$, $A_1=2$ and $A_2=-1$. Then

$$H\left(u\right)=1+2u-u^2.$$

Case 2.When $H\left(u\right)$ is a rational function, i.e.,

$$H\left(u\right)=\frac{1+A_{0}u}{A_1+A_{2}u}.$$

Using the conditions of Theorem 1, we get that $A_0=\frac{5}{2}$, $A_1=1$ and $A_2=\frac{1}{2}$. So

$$H\left(u\right)=\frac{2+5u}{2+u}.$$

Case 3.Consider $H\left(u\right)$ as another rational weight function, e.g.,

$$H\left(u\right)=\frac{1+A_{0}u+A_1{u}^2}{1+A_{2}u}.$$

Using the conditions of Theorem 1, we obtain $A_0=3$, $A_1=2$ and $A_2=1$. Then $H\left(u\right)$ becomes

$$H\left(u\right)=\frac{1+3u+u^2}{1+u}.$$

Case 4.When $H\left(u\right)$ is a yet another rational function of the form

$$H\left(u\right)=\frac{1+A_{0}u}{1+A_{1}u+A_2{u}^2}.$$

Using the conditions of Theorem 1, we have $A_0=1$, $A_1=-1$ and $A_2=1$. Then

$$H\left(u\right)=\frac{1+u}{1-u+3u^2}.$$

2.2. Some Concrete Forms of $G(u,w)$

Case 5.Considering $G(u,w)$ a polynomial function, e.g.,

$$G(u,w)=A_0+A_{1}u+A_2{u}^2+A_{3}w.$$

From the conditions of Theorem 1, we get $A_0=1$, $A_1=2$, $A_2=0$ and $A_3=1$. So

$$G(u,w)=1+2u+w.$$

Case 6.Considering $G(u,w)$ a sum of two rational functions, that is

$$G(u,w)=\frac{A_0+2u}{1+A_{1}u}+\frac{B_0}{1+B_{1}w}.$$

By using the conditions of Theorem 1, we find that $A_0=0$, $A_1=0$, $B_0=1$ and $B_1=-1$. $G(u,w)$ becomes

$$G(u,w)=2u+\frac{1}{1-w}.$$

Case 7.When $G(u,w)$ is a product of two rational functions, that is

$$G(u,w)=\frac{1+A_{0}u}{1+A_{1}u}\times \frac{B_0}{1+B_{1}w}.$$

Then the conditions of Theorem 1 yield $A_0=2$, $A_1=0$, $B_0=1$ and $B_1=-1$. So

$$G(u,w)=\frac{1+2u}{1-w}.$$

3. Complex Dynamics of Methods

Here we analyze the complex dynamics of new methods based on a graphical tool ‘the basins of attraction’ of the roots of a polynomial $p\left(z\right)$ in Argand plane. Analysis of the basins gives an important information about the stability and convergence region of iterative methods. Wider is the convergence region (i.e., basin), better is the stability. The idea of complex dynamics was introduced initially by Vrscay and Gilbert [25]. In recent times, many authors have used this concept in their work, see, for example [26,27] and references therein. We consider some of the cases corresponding to the previously obtained forms of $H\left(u\right)$ and $G(u,w)$) of family (4) to assess the basins of attraction. Let us select the combinations: cases 1 and 2 of $H\left(u\right)$ with the cases 5, 6 and 7 of $G(u,w)$) in the scheme (4), and denote the corresponding methods by NM-i(j), i = I, II and j = a, b, c.

To start with we take the initial point $z_0$ in a rectangular region $R\in \mathbb{C}$ that contains all the roots of a polynomial $p\left(z\right).$ The iterative method when starts from point $z_0$ in a rectangle either converges to the root $P\left(z\right)$ or eventually diverges. The stopping criterion for convergence is considered to be ${10}^{-3}$ up to a maximum of 25 iterations. If the required accuracy is not achieved in 25 iterations, we conclude that the method with initial point $z_0$ does not converge to any root. The strategy adopted is as follows: A color is allocated to each initial point $z_0$ lying in the basin of attraction of a root. If the iteration initiating at $z_0$ converges, it represents the attraction basin painted with assigned color to it, otherwise, the non-convergent cases are painted by the black color.

To view the geometry in complex plane, we characterize attraction basins associated with the methods NM-I(a–c) and NM-II(a–c) considering the following four polynomials:

Problem 1.

Consider the polynomial $p_1\left(z\right)={(z^2-1)}^3$, which has roots $\{-1,1\}$ with multiplicity three. We use a grid of $400\times 400$ points in a rectangle $R\in \mathbb{C}$ of size $[-3,3]\times [-3,3]$ and assign red color to each initial point in the attraction basin of root $-1$ and green color to each point in the attraction basin of root 1. The basins so plotted for NM-I(a–c) and NM-II(a–c) are displayed in Figure 1. Looking at these graphics, we conclude that the method NM-II(c) possesses better stability followed by NM-I(c) and NM-II(b). Black zones in the figures show the divergent nature of a method when it starts assuming initial point from such zones.

Problem 2.

Let $p_2\left(z\right)={(z^3-1)}^2$ that has three roots $\{-0.5\pm 0.866025i,1\}$ each with multiplicity two. To plot the graphics, we use a grid of $400\times 400$ points in a rectangle $R\in \mathbb{C}$ of size $[-3,3]\times [-3,3]$ and assign the colors blue, green and red corresponding to each point in the basins of attraction of 1, $-0.5+0.866025i$ and $-0.5-0.866025i$. Basins drawn for the methods NM-I(a–c) and NM-II(a–c) are shown in Figure 2. As can be observed from the pictures, the method NM-I(c) and NM-II(c) possess a small number of divergent points and therefore have better convergence than the remaining methods.

Problem 3.

Let $p_3\left(z\right)={(z^6-1)}^3$ with six roots $\{\pm 1,-0.5\pm 0.866025i,0.5\pm 0.866025i\}$ each with multiplicity $m=3$. Basins obtained for the considered methods are presented in Figure 3. To draw the pictures, the red, blue, green, pink, cyan and magenta colors have been assigned to the attraction basins of the six roots. We observe from the graphics that the method NM-I(c) and NM-II(c) have better convergence behavior since they have lesser number of divergent points. On the other hand NM-I(a) and NM-II(a) contain large black regions followed by NM-I(b) and NM-II(b) indicating that the methods do not converge in 25 iterations starting at those points.

Problem 4.

Consider the polynomial $p_4\left(z\right)=z^4-6z^2+8$ that has four simple roots $\{\pm 2,\pm 1.414\dots \}$. In this case also, we use a grid of $400\times 400$ points in a rectangle $R\in \mathbb{C}$ of size $[-3,3]\times [-3,3]$ and allocate the red, blue, green and yellow colors to the basins of attraction of these four roots. Basins obtained for the methods are shown in Figure 4. Observing the basins, we conclude that the method NM-II(c) possesses better stability followed by NM-I(c). Remaining methods show chaotic nature along the boundaries of the attraction basins.

Looking at the graphics, one can easily judge the stable behavior and so the better convergence of any method. We reach to a root, if we start the iteration choosing $z_0$ anywhere in the basin of that root. However, if we choose an initial guess $z_0$ in a region wherein different basins of attraction meet each other, it is difficult to predict which root is going to be attained by the iterative method that starts from $z_0$. So, the choice of $z_0$ in such a region is not a good one. Both black regions and the regions with different colors are not suitable to assume the initial guess as $z_0$ when we are required to achieve a particular root. The most intricate geometry is between the basins of attraction, which corresponds to the cases where the method is more demanding with respect to the initial point. From the basins, one can conclude that the method NM-II(c) possesses better stability followed by NM-I(c) than the remaining methods.

4. Numerical Tests

In this section, we apply the special cases NM-i(j), i = I, II and j = a, b, c of the scheme (4), corresponding to the combinations of $H\left(u\right)$: cases 1 and 2 with that of $G(u,w)$): cases 5, 6 and 7, to solve some nonlinear equations for validation of the theoretical results that we have derived. The theoretical seventh order of convergence is verified by calculating the computational order of convergence (COC)

$$\mathrm{COC}=\frac{ln|(x_{n+1}-\alpha )/(x_n-\alpha )|}{ln|(x_n-\alpha )/(x_{n-1}-\alpha )|},$$

which is given in (see [28]). Comparison of performance is also done with some existing methods such as the sixth order methods by Geum et al. [22,23], which are already expressed by (2) and (3). To represent $Q_f(u,s)$, we choose the following four special cases in the formula (2) and denote the respective methods by GKN-I(j), j = a, b, c, d:

(a)

$Q_f(u,s)=m(1+2(m-1)(u-s)-4us+s^2).$

(b)

$Q_f(u,s)=m(1+2(m-1)(u-s)-u^2-2us).$

(c)

$Q_f(u,s)=\frac{m+au}{1+bu+cs+dus},$ where $a=\frac{2m}{m-1}$, $b=2-2m$, $c=\frac{2(2-2m+m^2)}{m-1}$, $d=-2m(m-1)$.

(d)

$Q_f(u,s)=\frac{m+a_{1}u}{1+b_{1}u+c_1{u}^2}\frac{1}{1+d_{1}s},$ where $a_1=\frac{2m(4m^4-16m^3+31m^2-30m+13}{(m-1)(4m^2-8m+7)}$, $b_1=\frac{4(2m^2-4m+3)}{(m-1)(4m^2-8m+7)}$,

$c_1=-\frac{4m^2-8m+3}{4m^2-8m+7}$, $d_1=2(m-1)$.

For the formula (3), considering the following four combinations of the functions $Q_f\left(u\right)$ and $K_f(u,v)$, and denoting the corresponding methods by GKN-II(j), j = a, b, c, d:

(a)

$Q_f\left(u\right)=\frac{1+u^2}{1-u}$, $K_f(u,v)=\frac{1+u^2-v}{1-u+(u-2)v}$.

(b)

$Q_f\left(u\right)=1+u+2u^2$, $K_f(u,v)=1+u+2u^2+(1+2u)v$.

(c)

$Q_f\left(u\right)=\frac{1+u^2}{1-u}$, $K_f(u,v)=1+u+2u^2+2u^3+2u^4+(2u+1)v$.

(d)

$Q_f\left(u\right)=\frac{(2u-1)(4u-1)}{1-7u+13u^2}$, $K_f(u,v)=\frac{(2u-1)(4u-1)}{1-7u+13u^2-(1-6u)v}$.

Computational work is compiled in the programming package of Mathematica software using multiple-precision arithmetic. Numerical results as displayed in

Table 1. Comparison of the numerical results for Example 1.

Methods	n	$|{\mathit{e}}_{\mathit{n}-3}|$	$|{\mathit{e}}_{\mathit{n}-2}|$	$|{\mathit{e}}_{\mathit{n}-1}|$	COC	CPU-t (s)
GKN-I(a)	4	$1.06\times {10}^{-9}$	$3.86\times {10}^{-56}$	$9.03\times {10}^{-335}$	6.0000	0.1567
GKN-I(b)	4	$1.06\times {10}^{-9}$	$3.91\times {10}^{-56}$	$9.85\times {10}^{-335}$	6.0000	0.1583
GKN-I(c)	4	$1.06\times {10}^{-9}$	$4.34\times {10}^{-56}$	$2.02\times {10}^{-334}$	6.0000	0.1525
GKN-I(d)	4	$1.07\times {10}^{-9}$	$1.17\times {10}^{-55}$	$2.02\times {10}^{-331}$	6.0000	0.1600
GKN-II(a)	4	$1.19\times {10}^{-6}$	$5.39\times {10}^{-38}$	$4.56\times {10}^{-226}$	5.9999	0.1835
GKN-II(b)	4	$1.20\times {10}^{-6}$	$1.61\times {10}^{-37}$	$9.49\times {10}^{-223}$	5.9999	0.1640
GKN-II(c)	4	$1.20\times {10}^{-6}$	$1.12\times {10}^{-37}$	$7.51\times {10}^{-224}$	5.9999	0.1718
GKN-II(d)	4	$1.20\times {10}^{-6}$	$1.87\times {10}^{-37}$	$2.76\times {10}^{-222}$	5.9999	0.1680
NM-I(a)	3	$9.83\times {10}^{-8}$	$4.34\times {10}^{-51}$	0	7.0000	0.1562
NM-I(b)	3	$1.16\times {10}^{-9}$	$1.38\times {10}^{-64}$	0	7.0000	0.1170
NM-I(c)	3	$6.30\times {10}^{-10}$	$7.75\times {10}^{-67}$	0	7.0000	0.1485
NM-II(a)	3	$9.83\times {10}^{-8}$	$4.41\times {10}^{-51}$	0	7.0000	0.1367
NM-II(b)	3	$1.16\times {10}^{-9}$	$1.40\times {10}^{-64}$	0	7.0000	0.1562
NM-II(c)	3	$6.30\times {10}^{-10}$	$8.07\times {10}^{-67}$	0	7.0000	0.1405

,

Table 2. Comparison of the numerical results for Example 2.

Methods	n	$|{\mathit{e}}_{\mathit{n}-3}|$	$|{\mathit{e}}_{\mathit{n}-2}|$	$|{\mathit{e}}_{\mathit{n}-1}|$	COC	CPU-t (s)
GKN-I(a)	4	$2.17\times {10}^{-8}$	$4.61\times {10}^{-25}$	$1.01\times {10}^{-152}$	6.0000	1.4218
GKN-I(b)	4	$2.17\times {10}^{-8}$	$4.60\times {10}^{-25}$	$2.27\times {10}^{-151}$	6.0000	1.4923
GKN-I(c)	4	$2.11\times {10}^{-8}$	$4.21\times {10}^{-25}$	$1.03\times {10}^{-150}$	6.0000	1.4532
GKN-I(d)	4	$1.77\times {10}^{-8}$	$2.48\times {10}^{-25}$	$2.68\times {10}^{-151}$	6.0000	1.4960
GKN-II(a)	4	$4.83\times {10}^{-7}$	$1.36\times {10}^{-41}$	$6.84\times {10}^{-249}$	6.0000	1.3867
GKN-II(b)	4	$4.90\times {10}^{-7}$	$2.89\times {10}^{-41}$	$1.21\times {10}^{-246}$	6.0000	1.3790
GKN-II(c)	4	$4.88\times {10}^{-7}$	$2.22\times {10}^{-41}$	$1.98\times {10}^{-247}$	6.0000	1.4110
GKN-II(d)	4	$4.89\times {10}^{-7}$	$3.22\times {10}^{-41}$	$2.62\times {10}^{-246}$	6.0000	1.3982
NM-I(a)	3	$1.65\times {10}^{-8}$	$2.82\times {10}^{-58}$	0	7.0000	1.1367
NM-I(b)	3	$7.69\times {10}^{-9}$	$1.35\times {10}^{-60}$	0	7.0000	1.1915
NM-I(c)	3	$3.65\times {10}^{-9}$	$3.19\times {10}^{-63}$	0	7.0000	1.1407
NM-II(a)	3	$1.65\times {10}^{-9}$	$2.86\times {10}^{-58}$	0	7.0000	1.1290
NM-II(b)	3	$7.69\times {10}^{-9}$	$1.36\times {10}^{-60}$	0	7.0000	1.2540
NM-II(c)	3	$3.65\times {10}^{-9}$	$3.27\times {10}^{-63}$	0	7.0000	1.1445

,

Table 3. Comparison of numerical results for Example 3.

Methods	n	$|{\mathit{e}}_{\mathit{n}-3}|$	$|{\mathit{e}}_{\mathit{n}-2}|$	$|{\mathit{e}}_{\mathit{n}-1}|$	COC	CPU-t (s)
GKN-I(a)	4	$1.29\times {10}^{-3}$	$5.18\times {10}^{-20}$	$2.19\times {10}^{-118}$	6.0000	0.0313
GKN-I(b)	4	$1.48\times {10}^{-3}$	$1.63\times {10}^{-19}$	$2.19\times {10}^{-115}$	5.9998	0.0390
GKN-I(c)	4	$1.45\times {10}^{-3}$	$1.76\times {10}^{-19}$	$5.56\times {10}^{-115}$	5.9997	0.0352
GKN-I(d)	4	$1.97\times {10}^{-3}$	$1.80\times {10}^{-18}$	$1.07\times {10}^{-108}$	5.9996	0.0428
GKN-II(a)	4	$5.67\times {10}^{-4}$	$1.20\times {10}^{-22}$	$1.06\times {10}^{-134}$	5.9999	0.0314
GKN-II(b)	4	$2.39\times {10}^{-3}$	$5.78\times {10}^{-18}$	$1.16\times {10}^{-105}$	5.9996	0.0396
GKN-II(c)	4	$1.70\times {10}^{-3}$	$4.26\times {10}^{-19}$	$1.08\times {10}^{-112}$	5.9997	0.0392
GKN-II(d)	4	$1.55\times {10}^{-2}$	$5.18\times {10}^{-13}$	$7.23\times {10}^{-76}$	6.0000	0.0354
NM-I(a)	4	$1.13\times {10}^{-4}$	$6.52\times {10}^{-23}$	$1.41\times {10}^{-157}$	6.9998	0.0275
NM-I(b)	4	$9.26\times {10}^{-4}$	$1.63\times {10}^{-23}$	$8.75\times {10}^{-162}$	6.9998	0.0313
NM-I(c)	4	$4.64\times {10}^{-4}$	$4.44\times {10}^{-26}$	$3.23\times {10}^{-180}$	6.9998	0.0275
NM-II(a)	4	$1.13\times {10}^{-4}$	$6.83\times {10}^{-23}$	$2.00\times {10}^{-157}$	6.9998	0.0316
NM-II(b)	4	$9.33\times {10}^{-4}$	$1.77\times {10}^{-23}$	$1.58\times {10}^{-161}$	6.9998	0.0275
NM-II(c)	4	$4.78\times {10}^{-4}$	$5.86\times {10}^{-26}$	$2.43\times {10}^{-179}$	6.9998	0.0354

,

Table 4. Comparison of numerical results for Example 4.

Methods	n	$|{\mathit{e}}_{\mathit{n}-3}|$	$|{\mathit{e}}_{\mathit{n}-2}|$	$|{\mathit{e}}_{\mathit{n}-1}|$	COC	CPU-t (s)
GKN-I(a)	5	$1.90\times {10}^{-5}$	$9.03\times {10}^{-22}$	$1.05\times {10}^{-119}$	6.0000	0.0471
GKN-I(b)	5	$2.31\times {10}^{-5}$	$3.69\times {10}^{-21}$	$6.14\times {10}^{-116}$	6.0000	0.0472
GKN-I(c)	5	$2.18\times {10}^{-5}$	$3.18\times {10}^{-21}$	$3.14\times {10}^{-116}$	6.0000	0.0465
GKN-I(d)	5	$3.58\times {10}^{-5}$	$1.01\times {10}^{-19}$	$5.02\times {10}^{-107}$	6.0000	0.0483
GKN-II(a)	5	$3.00\times {10}^{-6}$	$4.91\times {10}^{-27}$	$9.51\times {10}^{-152}$	6.0000	0.0474
GKN-II(b)	5	$4.78\times {10}^{-5}$	$5.42\times {10}^{-19}$	$1.17\times {10}^{-102}$	6.0000	0.0472
GKN-II(c)	5	$2.51\times {10}^{-5}$	$6.82\times {10}^{-21}$	$2.75\times {10}^{-114}$	6.0000	0.0481
GKN-II(d)	7	$3.85\times {10}^{-11}$	$1.78\times {10}^{-55}$	$1.75\times {10}^{-321}$	6.0000	0.0625
NM-I(a)	5	$1.06\times {10}^{-5}$	$4.09\times {10}^{-26}$	$5.33\times {10}^{-169}$	7.0000	0.0368
NM-I(b)	5	$5.10\times {10}^{-6}$	$2.51\times {10}^{-28}$	$1.73\times {10}^{-184}$	7.0000	0.0322
NM-I(c)	5	$1.15\times {10}^{-6}$	$2.55\times {10}^{-33}$	$6.75\times {10}^{-220}$	7.0000	0.0327
NM-II(a)	5	$1.05\times {10}^{-5}$	$4.13\times {10}^{-23}$	$5.89\times {10}^{-169}$	7.0000	0.0316
NM-II(b)	5	$5.16\times {10}^{-6}$	$2.76\times {10}^{-23}$	$3.48\times {10}^{-184}$	7.0000	0.0323
NM-II(c)	5	$1.20\times {10}^{-6}$	$3.65\times {10}^{-26}$	$9.09\times {10}^{-219}$	7.0000	0.0314

and

Table 5. Comparison of the numerical results for Example 5.

Methods	n	$|{\mathit{e}}_{\mathit{n}-3}|$	$|{\mathit{e}}_{\mathit{n}-2}|$	$|{\mathit{e}}_{\mathit{n}-1}|$	COC	CPU-t (s)
GKN-I(a)	4	$1.12\times {10}^{-4}$	$5.78\times {10}^{-24}$	$1.10\times {10}^{-139}$	6.0000	0.2772
GKN-I(b)	4	$1.55\times {10}^{-4}$	$7.30\times {10}^{-23}$	$8.07\times {10}^{-133}$	6.0000	0.2462
GKN-I(c)	4	$1.39\times {10}^{-4}$	$4.40\times {10}^{-23}$	$4.43\times {10}^{-134}$	6.0000	0.2497
GKN-I(d)	4	$2.32\times {10}^{-4}$	$1.95\times {10}^{-21}$	$6.85\times {10}^{-124}$	6.0000	0.2812
GKN-II(a)	4	$3.36\times {10}^{-5}$	$8.72\times {10}^{-28}$	$2.66\times {10}^{-163}$	6.0000	0.3397
GKN-II(b)	4	$3.39\times {10}^{-5}$	$2.19\times {10}^{-20}$	$1.57\times {10}^{-117}$	6.0000	0.2695
GKN-II(c)	4	$2.16\times {10}^{-5}$	$7.70\times {10}^{-22}$	$1.58\times {10}^{-126}$	6.0000	0.2460
GKN-II(d)	4	$3.51\times {10}^{-3}$	$3.25\times {10}^{-14}$	$2.03\times {10}^{-80}$	6.0000	0.2342
NM-I(a)	4	$1.52\times {10}^{-4}$	$8.45\times {10}^{-26}$	$1.41\times {10}^{-174}$	6.9999	0.1445
NM-I(b)	4	$1.25\times {10}^{-4}$	$2.22\times {10}^{-26}$	$1.23\times {10}^{-178}$	6.9999	0.1522
NM-I(c)	4	$5.26\times {10}^{-4}$	$1.58\times {10}^{-29}$	$3.54\times {10}^{-201}$	6.9999	0.1640
NM-II(a)	4	$1.52\times {10}^{-4}$	$9.05\times {10}^{-26}$	$2.36\times {10}^{-174}$	6.9999	0.1482
NM-II(b)	4	$1.27\times {10}^{-4}$	$2.49\times {10}^{-26}$	$2.84\times {10}^{-178}$	6.9999	0.1492
NM-II(c)	4	$5.54\times {10}^{-4}$	$2.51\times {10}^{-29}$	$9.82\times {10}^{-200}$	6.9999	0.1642

contain: (i) number of iterations $\left(n\right)$ needed to converge to desired solution, (ii) last three successive errors $e_n=|x_{n+1}-x_n|$, (iii) computational order of convergence (COC) and (iv) CPU-time (CPU-t) in seconds elapsed during the execution of a program. Required iteration $\left(n\right)$ and elapsed CPU-time are computed by selecting $|x_{n+1}-x_n|+|f\left(x_n\right)|<{10}^{-350}$ as the stopping condition.

For numerical tests we select seven problems. The first four problems are of practical interest where as last three are of academic interest. In the problems we need not to calculate the root multiplicity m and it is set a priori, before running the algorithm.

Example 1 (Eigen value problem).

Finding Eigen values of a large sparse square matrix is a challenging task in applied mathematics and engineering sciences. Calculating the roots of a characteristic equation of matrix of order larger than 4 is even a big job. We consider the following 9 × 9 matrix.

$$A=\frac{1}{8}\left[\begin{array}{ccccccccc}-12& 0& 0& 19& -19& 76& -19& 18& 437\\ -64& 24& 0& -24& 24& 64& -8& 32& 376\\ -16& 0& 24& 4& -4& 16& -4& 8& 92\\ -40& 0& 0& -10& 50& 40& 2& 20& 242\\ -4& 0& 0& -1& 41& 4& 1& 0& 25\\ -40& 0& 0& 18& -18& 104& -18& 20& 462\\ -84& 0& 0& -29& 29& 84& 21& 42& 501\\ 16& 0& 0& -4& 4& -16& 4& 16& -92\\ 0& 0& 0& 0& 0& 0& 0& 0& 24\end{array}\right].$$

We calculate the characteristic polynomial of the matrix $\left(A\right)$ as

$$\begin{array}{c}\hfill f_1\left(x\right)=x^9-29x^8+349x^7-2261x^6+8455x^5-17663x^4+15927x^3+6993x^2-24732x+12960.\end{array}$$

This function has a multiple root $\alpha =3$ with multiplicity 4. We select initial value $x_0=2.25$ and obtain the numerical results as shown in Table 1.

Example 2 (Manning equation for fluid dynamics).

Next, the problem of isentropic supersonic flow around a sharp expansion corner is chosen (see [2]). Relation among the Mach number before the corner (say $M_1$) and after the corner (say $M_2$) is given by

$$\delta =b^{1/2}\left({tan}^{-1}{\left(\frac{M_2^2-1}{b}\right)}^{1/2}-{tan}^{-1}{\left(\frac{M_1^2-1}{b}\right)}^{1/2}\right)-\left({tan}^{-1}{(M_2^2-1)}^{1/2}-{tan}^{-1}{(M_1^2-1)}^{1/2}\right),$$

where $b=\frac{\gamma +1}{\gamma -1}$ and γ is the specific heat ratio of gas.

For a specific case, the above equation is solved for for $M_2$, given that $M_1=1.5$, $\gamma =1.4$ and $\delta ={10}^0$. Then, we have that

$${tan}^{-1}\left(\frac{\sqrt{5}}{2}\right)-{tan}^{-1}\left(\sqrt{x^2-1}\right)+\sqrt{6}\left({tan}^{-1}\left(\sqrt{\frac{x^2-1}{6}}\right)-{tan}^{-1}\left(\frac{1}{2}\sqrt{\frac{5}{6}}\right)\right)-\frac{11}{63}=0$$

where $x=M_2$.

Let us consider this particular case seven times using same values of the involved parameters and then obtain the nonlinear function

$$f_2\left(x\right)={\left[{tan}^{-1}\left(\frac{\sqrt{5}}{2}\right)-{tan}^{-1}\left(\sqrt{x^2-1}\right)+\sqrt{6}\left({tan}^{-1}\left(\sqrt{\frac{x^2-1}{6}}\right)-{tan}^{-1}\left(\frac{1}{2}\sqrt{\frac{5}{6}}\right)\right)-\frac{11}{63}\right]}^7.$$

The above function has one root at $\alpha =1.8411027704\dots$ of multiplicity 7 with initial approximations $x_0$$=1.50$. Computed numerical results are shown in Table 2.

Example 3 (Beam designing model).

We consider the problem of beam positioning (see [4]) where a beam of length r unit leans against the edge of a cubical box of sides 1 unit distance each, such that one end of the beam touches the wall and the other end touches the floor, as depicted in Figure 5.

The problem is: What will be the distance alongside the floor from the base of wall to the bottom of beam? Suppose that y is distance along the beam from the floor to the edge of the box and x is the distance from the bottom of box to the bottom of beam. For a given r, we can obtain the equation

$$f_3\left(x\right)=x^4+4x^3-24x^2+16x+16=0.$$

One of the roots of this equation is the double root $x=2$. We select the initial guess $x_0=3$ to find the root. Numerical results by various methods are shown in Table 3.

Example 4 (van der Waals equation).

Consider the Van der Waals equation

$$\left(P+\frac{a_1{n}^2}{V^2}\right)(V-n{a}_2)=nRT,$$

that describes nature of a real gas by adding in the ideal gas equation two parameters, $a_1$ and $a_2$, which are specific for each gas. To find the volume V in terms of rest of the parameters one requires to solve the equation

$$P{V}^3-(n{a}_{2}P+nRT)V^2+a_1{n}^{2}V=a_1{a}_2{n}^3.$$

Given a set of values of $a_1$ and $a_2$ of a particular gas, one can find values for n, P and T, so that this equation has three real roots. Using a particular set of values (see [3]), we have the equation

$$f_4\left(x\right)=x^3-5.22x^2+9.0825x-5.2675=0,$$

where $x=V$. This equation has a multiple root $\alpha =1.75$ with multiplicity 2. The initial guess chosen to obtain the root 1.75 is $x_0=2$. Numerical results are shown in Table 4.

Example 5.

Consider now the standard nonlinear test function (see [23])

$$f_5\left(x\right)=\left(9-2x-2x^4+cos2x\right)\left(5-x-x^4-{sin}^{2}x\right).$$

The root $\alpha =1.29173329244360\dots$ of multiplicity 2 is computed with initial guess $x_0=1.5$. Numerical results are displayed in Table 5.

Example 6.

Let us assume another nonlinear test function given as (see [22])

$$f_6\left(x\right)=\left(x-\sqrt{3}{x}^{3}cos\left(\frac{\pi x}{6}\right)+\frac{1}{x^2+1}-\frac{11}{5}+4\sqrt{3}\right){(x-2)}^4.$$

The root $\alpha =2$ of this function is of multiplicity 5. This root is calculated assuming the initial approximation $x_0=1.5$. Results so obtained are shown in

Table 6. Comparison of the numerical results for Example 6.

Methods	n	$|{\mathit{e}}_{\mathit{n}-3}|$	$|{\mathit{e}}_{\mathit{n}-2}|$	$|{\mathit{e}}_{\mathit{n}-1}|$	COC	CPU-t (s)
GKN-I(a)	4	$1.20\times {10}^{-5}$	$6.82\times {10}^{-31}$	$2.31\times {10}^{-182}$	6.0000	0.6797
GKN-I(b)	4	$1.20\times {10}^{-5}$	$6.86\times {10}^{-31}$	$2.40\times {10}^{-182}$	6.0000	0.6680
GKN-I(c)	4	$1.21\times {10}^{-5}$	$7.72\times {10}^{-31}$	$5.18\times {10}^{-182}$	6.0000	0.6992
GKN-I(d)	4	$1.58\times {10}^{-5}$	$1.00\times {10}^{-29}$	$6.51\times {10}^{-175}$	6.0000	0.6720
GKN-II(a)	4	$3.17\times {10}^{-5}$	$1.64\times {10}^{-28}$	$3.21\times {10}^{-168}$	6.0000	0.8047
GKN-II(b)	4	$3.50\times {10}^{-5}$	$6.90\times {10}^{-28}$	$4.05\times {10}^{-164}$	6.0000	0.8280
GKN-II(c)	4	$3.41\times {10}^{-5}$	$4.42\times {10}^{-28}$	$2.09\times {10}^{-165}$	6.0000	0.7967
GKN-II(d)	4	$3.54\times {10}^{-5}$	$8.45\times {10}^{-28}$	$1.56\times {10}^{-163}$	6.0000	0.8242
NM-I(a)	4	$5.14\times {10}^{-6}$	$4.35\times {10}^{-38}$	$1.35\times {10}^{-262}$	7.0000	0.5625
NM-I(b)	4	$3.45\times {10}^{-6}$	$2.68\times {10}^{-39}$	$4.53\times {10}^{-271}$	7.0000	0.5782
NM-I(c)	4	$2.05\times {10}^{-6}$	$2.95\times {10}^{-41}$	$3.76\times {10}^{-285}$	7.0000	0.5277
NM-II(a)	4	$5.14\times {10}^{-6}$	$4.42\times {10}^{-38}$	$1.53\times {10}^{-262}$	7.0000	0.4805
NM-II(b)	4	$3.45\times {10}^{-6}$	$2.73\times {10}^{-39}$	$5.24\times {10}^{-271}$	7.0000	0.4725
NM-II(c)	4	$2.05\times {10}^{-6}$	$3.07\times {10}^{-41}$	$5.17\times {10}^{-285}$	7.0000	0.4610

.

Example 7.

Lastly, consider the test function

$$f_7\left(x\right)=\left(x^2+1\right)\left(2x{e}^{x^2+1}+x^3-x\right){cosh}^2\left(\frac{\pi x}{2}\right)$$

The function has multiple root $\alpha =i$ of multiplicity 4. We choose the initial approximations $x_0=1.25i$ for obtaining the root of the function. The results computed by various methods are shown in

Table 7. Comparison of the numerical results for Example 7.

Methods	n	$|{\mathit{e}}_{\mathit{n}-3}|$	$|{\mathit{e}}_{\mathit{n}-2}|$	$|{\mathit{e}}_{\mathit{n}-1}|$	COC	CPU-t (s)
GKN-I(a)	4	$2.53\times {10}^{-6}$	$3.79\times {10}^{-35}$	$4.32\times {10}^{-208}$	6.0000	1.1564
GKN-I(b)	4	$2.53\times {10}^{-6}$	$3.92\times {10}^{-35}$	$5.33\times {10}^{-208}$	6.0000	1.1577
GKN-I(c)	4	$2.68\times {10}^{-6}$	$6.07\times {10}^{-35}$	$8.23\times {10}^{-207}$	6.0000	1.1415
GKN-I(d)	4	$4.80\times {10}^{-6}$	$5.34\times {10}^{-33}$	$1.01\times {10}^{-194}$	6.0000	1.0473
GKN-II(a)	4	$5.04\times {10}^{-6}$	$1.82\times {10}^{-33}$	$4.04\times {10}^{-198}$	6.0000	1.0212
GKN-II(b)	4	$7.15\times {10}^{-6}$	$4.23\times {10}^{-32}$	$1.81\times {10}^{-189}$	6.0000	1.1215
GKN-II(c)	4	$6.39\times {10}^{-6}$	$1.51\times {10}^{-32}$	$2.64\times {10}^{-192}$	6.0000	1.2035
GKN-II(d)	4	$8.22\times {10}^{-6}$	$1.41\times {10}^{-31}$	$8.09\times {10}^{-187}$	6.0000	1.1416
NM-I(a)	4	$1.08\times {10}^{-6}$	$6.96\times {10}^{-43}$	$3.13\times {10}^{-296}$	7.0000	0.5787
NM-I(b)	4	$9.01\times {10}^{-7}$	$1.91\times {10}^{-43}$	$3.71\times {10}^{-300}$	7.0000	0.5632
NM-I(c)	4	$4.64\times {10}^{-7}$	$7.44\times {10}^{-46}$	$2.01\times {10}^{-317}$	7.0000	0.5586
NM-II(a)	4	$1.09\times {10}^{-6}$	$7.21\times {10}^{-43}$	$4.10\times {10}^{-296}$	7.0000	0.5478
NM-II(b)	4	$9.04\times {10}^{-7}$	$2.00\times {10}^{-43}$	$5.10\times {10}^{-300}$	7.0000	0.5946
NM-II(c)	4	$4.68\times {10}^{-7}$	$8.21\times {10}^{-46}$	$4.20\times {10}^{-317}$	7.0000	0.5644

.

From the numerical values of errors we observe the increasing accuracy in the values of successive approximations as the iteration proceed, which points to the stable nature of the methods. Like the existing methods, the convergence behavior of new methods is also consistent. At the stage when stopping criterion $|x_{n+1}-x_n|+|f\left(x_n\right)|<{10}^{-350}$ has been satisfied we display the value ‘0’ of $|x_{n+1}-x_n|$. From the calculation of computational order of convergence shown in the penultimate column in each table, we verify the theoretical convergence of seventh order. The entries of last column in each table show that the new methods consume less CPU-time during the execution of program than the time taken by existing methods. This confirms the computationally more efficient nature of the new methods. Among the new methods, the better performers (in terms of accuracy) are NM-I(c) and NM-II(c) since they produce approximations of the root with small error. However, this is not true when execution time is taken into account because if one method is better in some situations, then the other is better in some other situation. The main purpose of implementing the new methods for solving different type of nonlinear equations is purely to illustrate the better accuracy of the computed solution and the better computational efficiency than existing techniques. Similar numerical experimentation, performed on a variety of numerical problems of different kinds, confirmed the above remarks to a large extent.

5. Conclusions

In the present work, we have constructed a class of seventh order methods for solving nonlinear equations containing multiple roots. Analysis of the convergence has been carried out, which proves the seventh order of convergence under standard conditions of the function whose zeros we are looking for. Some particular cases of the family are presented. The stability of these cases are tested by means of visual display of the basins of attraction when the methods are applied on different polynomials. The methods are also implemented to solve nonlinear equations including those arising in practical problems. The performance is compared with existing methods in numerical testing. Superiority of proposed methods over the known techniques is endorsed by the numerical tests including the elapsed CPU-time in execution of program.