#### Abstract

A new inverse family of the iterative method is interrogated in the present article for simultaneously estimating all distinct and multiple roots of nonlinear polynomial equations. Convergence analysis proves that the order of convergence of the newly constructed family of methods is two. The computer algebra systems CAS-Mathematica is used to determine the lower bound of convergence order, which justifies the local convergence of the newly developed method. Some nonlinear models from physics, chemistry, and engineering sciences are considered to demonstrate the performance and efficiency of the newly constructed family of inverse simultaneous methods in comparison to classical methods in the literature. The computational time in seconds and residual error graph of the inverse simultaneous methods are also presented to elaborate their convergence behavior.

#### 1. Introduction

Considering nonlinear polynomial equation of degree ,with arbitrary real or complex coefficient . Let denote all the simple or complex roots of (1) with multiplicity . Newton’s method [11] is one of the most basic and ancient methods that is used to estimate single roots of (1) at a time as below:

Iterative method (2) has local quadratic convergence. Nedzhibov et al. [13] presented corresponding inverse numerical technique of the same convergence order as

Here, we propose the following family of the optimal second-order convergence method for finding simple roots of (1) aswhere . Method (4) is optimal and the convergence order of (4) is 2 if is a simple root of (1) and . The error equation of (4) is obtained using Maple-18:or

Corresponding inverse methods of (4) is constructed as

If is an exact root of (1), , and the following is obtained:

Inverse iterative schemes (7) are second-order convergence as it is easy to prove .

Besides simple root finding methods [3–5, 13, 15, 16, 18–20] in literature, there exists another class of numerical methods which estimate all real and complex roots of (1) at a time, known as simultaneous methods. Simultaneous numerical iterative schemes are very prevalent due to their global convergence properties and its parallel execution on computers [1, 6, 8–10, 12, 14].

The most prominent method among simultaneous derivative-free iterative technique is the Weierstrass–Dochive [17] method (abbreviated as MWM1), which is defined aswhereis Weierstrass’ correction. Method (10) has local quadratic convergence. For finding all multiple roots of (1), we use the following correction [17]:where is the multiplicity of the roots.

#### 2. Construction of the Inverse Simultaneous Method

Using Weierstrass correction in (7), we get a new family of inverse modified Weierstrass method (abbreviated as MWM2):

Inverse simultaneous iterative method (13) can also be written as

Thus, we construct a new derivative-free family of inverse iterative simultaneous scheme (13), abbreviated as MWM2, for estimating all distinct roots of (1). To estimate all multiple roots of (1), we use correction (12) instead of (11) in (7).

##### 2.1. Convergence Framework

In this section, we demonstrate convergence theorem of inverse iterative scheme MWM2.

Theorem 1. *Let be single zero of (1) and for necessarily close primary distinct guess of the zero, respectively; then, MWM2 has local -order convergence.*

*Proof. *Let be the errors in and respectively. For the simplicity of the calculation, we omit the iteration index. Then,orUsing in (15), we haveThus, we obtainUsing the expression [9] in (20), we haveIf all the errors are assumed of the same order, i.e., , thenHence, it is proved.

#### 3. Lower Bound of Convergence of MWM1 and MWM2

Computer algebra system, Mathematica, has been used to find the lower bound of convergence of MWM1 and MWM2.

Considerwhere , , and are exact zeros of (23). The first component of (where ) of numerical iterative methods is for finding zeros of (23), , simultaneously. We have to express the derivatives of , i.e., the partial derivatives of with respect to are as follows:and so on.

We obtain the lower bound of convergence order till the first nonzero element of row is found. The Mathematica notebook codes are used for the following MWM1 and MWM2: Weierstrass–Dochive method, MWM1: Modified inverse family of iterative schemes, MWM2:

#### 4. Numerical Results

Some engineering problems are considered to demonstrate the performance and effectiveness of the simultaneous method, MWM2 and MWM1. For computer calculations, we use CAS-Maple-18, and the following stopping criteria for termination of computer are programmed:where signifies the absolute error. In Tables 1–5, C-Time represents computational time in second.

##### 4.1. Engineering Applications

Some engineering applications are deliberated in this section in order to show the feasibility of the present work.

*Example 1. *(see [2]). Considering a physical problem of beam positioning results in the following nonlinear polynomial equation:The exact root of (30), , is 2 with multiplicity 2 and the remaining other two roots are simple, i.e., and . We take the following initial estimates:Table 1 clearly demonstrates the superiority of MWM2 over MWM1 in terms of predicted absolute error and CPU time for guesstimating all real roots of (30) on the same number of iterations .

*Example 2. *(see [16]). In this engineering application, we consider a reactor of stirred tank. Items H and H are fed to the reactor at rates of ß and q-*ß*, respectively. Composite reaction improves in the apparatus as below:Douglas et al. [7] first examined this complex control system and obtained the following nonlinear polynomial equation:where is the gain of the proportional controller. By taking , we haveThe exact distinct roots of (34) are calculated as , and we take the following initial guessed values:Table 2 evidently illustrates the supremacy behavior of MWM2 over MWM1 in terms of the estimated absolute error and in CPU time on the same number of iterations *n* = 7 for guesstimating all real roots of (34).

*Example 3. *(see [4]). Consider the functionThe problem describes the fractional alteration of nitrogen-hydrogen (NH) feed into ammonia at 250 atm pressure and C temperature. Since the (37) is of order four, it has four roots:The initial approximated value for (27) is taken asTable 3 evidently shows the supremacy behavior of MWM2 over MWM1 in terms of estimated absolute error and in CPU time on the same number of iterations for guesstimating all real and complex roots of (37). Minuscule alteration of nitrogen-hydrogen (NH) feed into ammonia lies between (0,1); therefore, our desire root is up to 1900 decimal places:Remaining other approximating roots are = −1.283 404 526*e* − 1457-8.93 219 631 521*e* − 1457i, = −8.21 745 235 223 462*e* − 1091 + 0i, and −1.2 834 045 268 801*e* − 1457 + 8.99 321 963152*e* − 1457*i*.

*Example 4. *(see [8]).Considerwith multiple exact roots:The initial estimations have been taken asFor distinct roots,Table 4 evidently shows the supremacy behavior of MWM2 over MWM1 in terms of estimated absolute error and in CPU time on the same number of iterations for guesstimating all real and complex roots of (41).

*Example 5. *(see [5]). The sourness of a soaked solution of magnesium-hydroxide (MgOH) in hydroelectric acid (HCl) is given byfor the cranium ion concentration . If we set , we obtain the following polynomial:with exact roots of (46), up to one decimal places. The initial estimates have been taken asTable 5 evidently illustrates the supremacy behavior of MWM2 over MWM1 in terms of estimated absolute error and in CPU time on the same number of iterations for guesstimating all real and complex roots of (46).

*Example 6. *(see [21]). In general, mechanical engineering, as well as the majority of other scientists, uses thermodynamics extensively in their research work. The following polynomial is used to relate the zero-pressure specific heat of dry air, , to temperature:The temperature that corresponds to specific heat of needs to be determined. Putting in (48), we havewith exact roots , , and . The initial estimations of (49) have been taken asTable 6 clearly illustrates the supremacy behavior of MWM1 over MWM2 in estimated absolute error and in CPU time on the same number of iterations for guesstimating all real and complex roots of (49).

#### 5. Conclusion

A new derivative-free family of inverse numerical methods of convergence order 2 for simultaneous estimations of all distinct and multiple roots of (1) was introduced and discussed in this paper. Tables 1–5 and Figure 1 clearly show that computational order of convergence of the proposed and existing methods are agreed with the theoretical results. Simulation time, from Figure 2, clearly indicates the supremacy of our newly proposed method MWM2 over existing Weierstrass method MWM1. The results of numerical test cases from Tables 1–5, CPU time, and residual error graph from Figure 3 demonstrated the effectiveness and rapid convergence of our proposed iterative method MWM2 as compared to MWM1.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

#### Data Availability

No data were used to support this study.

#### Disclosure

The statements made and views expressed are solely the responsibility of the authors.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Authors’ Contributions

All authors contributed equally and significantly in writing this article. All authors read and approved the final manuscript.