On the stability analysis of numerical schemes for solving non-linear polynomials arises in engineering problems

: This study shows the link between computer science and applied mathematics. It conducts a dynamics investigation of new root solvers using computer tools and develops a new family of single-step simple root-finding methods. The convergence order of the proposed family of iterative methods is two, according to the convergence analysis carried out using symbolic computation in the computer algebra system CAS-Maple 18. Without further evaluations of a given nonlinear function and its derivatives, a very rapid convergence rate is achieved, demonstrating the remarkable computing e ffi ciency of the novel technique. To determine the simple roots of nonlinear equations, this paper discusses the dynamic analysis of one-parameter families using symbolic computation, computer animation, and multi-precision arithmetic. To choose the best parametric value used in iterative schemes, it implements the parametric and dynamical plane technique using CAS-MATLAB @ R 2011 b . The dynamic evaluation of the methods is also presented utilizing basins of attraction to analyze their convergence behavior. Aside from visualizing iterative processes, this method illustrates not only iterative processes but also gives useful information regarding the convergence of the numerical scheme based on initial guessed values. Some nonlinear problems that arise in science and engineering are used to demonstrate the performance and e ffi ciency of the newly developed method compared to the existing method in the literature


Introduction
One of the most critical areas of applied mathematics are iterative techniques to solve nonlinear equations.Their importance is due to their wide applications in different fields, such as engineering, computer science, physics, biology, chemistry, banking, business, digital signal processing, control theory, insurance, social science, and many others.For instance, [1][2][3][4] evidence some of these applications.Computer visual representation, which converts numerical data into 2D or 3D images, unquestionably offers a modern and unconventional solution to some mathematical topics.It gives a fresh perspective and provides a deep knowledge of the characteristics and behavior of iterative processes for finding polynomial roots and allocating them (e.g., [5][6][7][8][9][10][11]).The development of digital computer hardware and software has substantially impacted many scientific disciplines and is frequently utilized in almost all scientific domains, including mathematics.The usage of advanced computational tools, such as symbolic computation and computer graphics, is necessary to tackle many complex problems in applied mathematics and other scientific disciplines.Thanks to the rapid development of computer graphics over the past 20 years, scientists from various fields can now use cutting-edge computer techniques to study and visualize phenomena that are not easily observable.Architecture and engineering, biomedical applications (X-ray computed tomography (CT), magnetic resonance imaging (MRI), mathematics, physics, business and management graphics, education and learning, user interface design, image syntheses, mapping and cartography, video and multimedia technology (animation, move, etc.), and business are just a few examples where high-quality computer visualization provides new and remarkable insights.Many problems in the areas above remain unsolvable by humans.Computer graphics allow a quick and practical link between machines and people.The technology and mathematics that computer graphics provide to a wide range of applications in theory and practice, and they indeed have great potential in the future (see e.g., [12][13][14][15][16]).The initial guess value affects the classical Newton-like method's performance, particularly for large-scale nonlinear problems.For this reason, many iterations of Newton's technique have been devised.For example, Knoll et al. [17] suggested the arc-length approaches, and Lemieux et al. [18] offered the Jacobian-Free Newton-Krylov method.Wu et al. [19] provided novel contributions to Newton's research.Qureshi et al. [20] presents the trapezoidal second-order Newton's, He et al. [21] presented the variational Newton-type scheme, Saheya et al. [22] employed the rational approximation function to suggest a Newton-like method, Abbasbandy [23] proposed an improved version of Newton's scheme by using a modified Adomian decomposition technique.Motivated by the above mentioned work, this study develops a novel family of Newton-type iterative algorithms using a single parameter for locating the simple root of nonlinear equations.The parameter value influences the newly developed family.The best parameter value for that family of methods should be chosen based on the method that has a larger region of convergence at starting points where the classical methods diverge.To achieve this aim, we deal with the corresponding complex dynamical system, determining fixed and critical points, analyzing asymptotic convergence, and identifying values that decompose the corresponding rational function.The main contribution of this work is reflected as follows: • Construction of a family of numerical schemes for finding roots of scalar nonlinear equations.
• Using CAS-Maple to verify the convergence order of the proposed method.
• Finding fixed and critical points utilizing complex dynamical analysis.
• Dynamical and parametric planes are generated using critical points.
• Computational tools are used to analyze the proposed technique's convergence rate, efficiency, stability, and applicability.• Analyzing how the suggested technique can be applied in various research disciplines.Some engineering applications are considered to show improved efficiency and stability compared to other classical techniques discussed in the literature.This research aims to construct a novel family of iterative methods using the computer algebra system CAS-Maple, Mathematica, MATLAB, Scilab, GNU Octave, Python, and many others (see e.g., [24,25] ).It uses parameter and dynamical planes to apply these techniques for picking the best parameter value for the iterative scheme.This paper is organized as follows: In section 2, the proposed method is introduced.In section 3, the stability of the proposed method using complex dynamical analysis is discussed.In section 4, the efficiency of this method is illustrated by applications.The last section is the conclusion of this article.

Construction of numerical methods and their convergence analysis
This section discusses some well-known methods in the literature.One of the most ancient and popular numerical methods for approximating single roots of a nonlinear equation is Newton's method [26] (NM), which has local quadratic convergence.

Special cases:
Now we construct some special cases of the family of methods described by (2.5).We construct the following numerical schemes by choosing β = 2 in (2.5).

Complex dynamical analysis of the family of numerical methods
This section discusses the stability of the numerical scheme Eq (2.16) against the backdrop of the complex dynamics.Iterative method Eq (2.16) leads to the following rational map: Recalling some basic concepts of this theory (detailed information can be found in [29][30][31][32][33][34][35][36][37][38]).Taking a rational function • A strange fixed point if it not associated to any root of non-linear equation (2.1).
An attracting point x * ∈ C defines a basin of attraction ℜ(x * ) as the set of starting points whose orbit tends to x * .
The iterative method BB1 holds the scaling theorem which allows a suitable change of co-ordinate to reduce of dynamics of iteration of general maps to the study of a specific family of iteration of similar maps.The one point iterative method (2.16) has a universal Julia set if a rational map exist which conjugates by the Mobius transformation.
, where a, b ∈ ℜ.The Möbius transformation given by: which we consider as a map from C ∪ ∞.Then, we have

□
Similarly, we can get the following conclusions.
Theorem 3. A map f g (x) (rational) arising from BB2, BB3, AC, NM, and YO applied to f ) The fixed points of (2.22) are x = −1, x = 0, x = 1, and x = ∞.We calculate f ′ g (x, α) for the stability of fixed points of BB1 as: (2.28) From (2.28) it is clear that x = 0 and ∞ are always super attractive fixed points but the stability of other fixed points are depends on α, which we present here.The operator f ′ g (x, α) for x = −1 gives Analyzing (2.29), as α −→ ±∞ we obtain horizontal asymptotes for f ′ g (−1, α) = 1 and vertical asymptotes for α = 0 (Newton method).
Theorem 4. The strange point x = −1 is known as: 4 , then x = −1 is an attractor and it can be super attractor for α = −1.

.40)
If α 1, 1/4, then observe that cr1 = 1 cr2 and cr1 = cr2 = 1 only when α = −1, while cr1 = cr2 = −1 only when α = 0. Figure 2 presents the behavior of the fixed and critical points for real of α.Fixed points −1, 0, and 1 are represented by orange, black, and green dotted lines, while critical points (cr1, cr2) are represented by red and blue lines, respectively.Theorem 6.The only member of BB1 whose operator is always conjugated to the rational map x 2 is the element corresponding to α = 0.

Parametric planes:
We generate a parametric plane by taking α over a mesh of 250, 000 points in a complex plane.Critical points are chosen as initial guessed values and then iterate until it converges to any fixed points.Taking 10 −3 as a tolerance.We paint the complex parameter value in red if it converges , and black in other cases.The parametric planes for the critical points are shown in Figure 4.

Dynamical planes:
Dynamical planes are generated as parametric planes.To generate them, the real and imaginary parts of the starting estimations are represented as two axes over a mesh of 250 × 250, and the same stopping criteria as in the parametric plane.The study assigns different colors to indicate which root the method converges to, and black in other cases.The dynamical planes of the iterative scheme BB1, YO, NM, and AC for different parametric values are shown in Figures 5-7.In Table 1, Els-Time denotes the amount of time in seconds that has passed, Start-Points denotes the number of starting points (62500.00 in a square), Average-It denotes the average number of iterations, C-Points denotes the number of converging points, and D-Points denotes the number of diverging points utilized to form dynamical planes (the basins of attractions [39,40]).Table 1  In Table 3, Els-Time denotes the amount of time in seconds that has passed, Start-Points denotes the number of starting points (62500.00 in a square), Average-It denotes the average number of iterations, C-Points denotes the number of converging points, and D-Points denotes the number of diverging points utilized to form dynamical planes (the basins of attractions).

Engineering applications
This section discusses an engineering application to check the efficiency and performance of the newly developed methods compared to existing methods, i.e., the AC, YO, and NM methods.
Example 1. Current-Voltage Diode [41] A manufacturing company provided a tunnel diode with voltage-current output: The tunnel diode is connected with resistor R and voltage source E. By applying Kirchhoff's voltage law, we find a relation for the steady current as follows: For a given E=1.4 and R =19, we find x by setting f(x)=0, 14 − x 19 = x 3 − 1.5x 2 + 0.6x, The exact roots of (2.41) are The initial estimates of f 1 (x) have been taken as x 2 = 0.5, x 3 = 0.7.
Table 4 clearly shows the dominance behavior of BB1-BB3 over the AC, YO, and NM iterative methods in approximated absolute error on various number of iterations n for the non-linear function f 1 (x).Various initial guess values are used in order to approximate all roots of f 1 (x).A spring develops a non-linear restoring force given by where F represents force in kg force and x is displacement in cm.Using value of constants r 1 = 5.25kg/cm, r 3 = 0.60kg/cm, r 5 = 0.0118kg/cm in (2.42), we get a non-linear equation in terms of the displacement x f 2 (x) = 0.0118x 5 + 0.6x 3 + 5.25x. (2.43) The exact roots of (2.43) are The initial estimates for simultaneous determination of all roots of non-linear function f 2 (x) has been taken as: (0) x 2 = −3.3i, x 3 = 3.3i, x 4 = −6.2i, x 5 = 6.2i.
Table 6 clearly shows the dominance behavior of BB1-BB3 over the AC, YO, and NM iterative methods in approximated absolute error on various number of iterations n for the non-linear function f 3 (x).Various initial guess values are used in order to approximate all roots of f 3 (x).

Conclusions
We developed a family of iterative techniques with two parameters for locating the roots of nonlinear equations.Due to the lengthy and complicated mathematical formulas used in the convergence study, we used symbolic computation utilizing several programs created in the Mathematica computer system.We implemented a different method to evaluate the root solver efficiency utilizing parametric planes (Figures 1-4) and dynamical planes (Figures 5-7) to select the best parameter value.The stability region and dynamical behavior of the strange fixed points were discussed.Root-finding methods can be analyzed graphically using this approach concerning various starting points.Well-built software will offer a visual knowledge of the convergence behavior of iterative methods for establishing comparison standards to show the basins of attraction (Figures 8-10).However, it will also consider various qualitative factors, such as CPU time and the convergence region.The numerical results show that the newly developed technique outperforms the current method regarding the residual error, CPU-time, and iterations.In the future, we will construct higher-order efficient, optimal, and stable iterative methods for finding simple as well as all distinct and multiple roots of Eq (2.1).We may explore the applications of these results in other areas.For example, the dynamics nonlinear inventory management system, numerical solution of fractional SIR epidemiological model, stability and optimal control strategies for a novel epidemic model of COVID-19.

Figure 1 .
Figure 1.Zone of stability of strange fixed points of numerical scheme BB1.Left-zone of stability of strange fixed point x = −1; Right-zone of stability of strange fixed point x = 1.

Figure 3 :
Figure 3: The union of respective stability function of the all strange fixed points of BB1-BB3, AC, and YO, respectively, from left to right.

Figure 3 .
Figure 3.The union of the respective stability function of the all strange fixed points of BB1-BB3, AC, and, YO respectively.

( a )
Details of the parametric plane for the critical points cr1.(b) Details of the parametric plane for the critical points cr2.

Figure 4 .
Figure 4. (a-d): The parametric planes of the numerical scheme BB1 correspond to the critical points cr1 and cr2.

Figure 7 .
Figure 7. Iterative algorithms (BB1-BB3, AC, YO, NM) with dynamical planes are demonstrated from top left to bottom right end for various values of the parameter.
clearly shows that our recently developed methods BB1-BB3 outperform the existing iterative methods in terms of CPU-Time, Average-It, T-Points, C-Points, and D-Points.In Table 2, Els-Time denotes the amount of time in seconds that has passed, Start-Points denotes the number of starting points (62500.00 in a square), Average-It denotes the average number of iterations, C-Points denotes the number of converging points, and D-Points denotes the number of diverging points utilized to form dynamical planes (the basins of attractions).Table 2 clearly shows that our recently developed methods BB1-BB3 outperform existing iterative methods in terms of CPU-Time, Average-It, T-Points, C-Points, and D-Points.

Figure 8 .
Figure 8.The basin of attraction for the numerical schemes BB1-BB3, AC, YO, and NM technique, respectively, is shown from left to right (for polynomial equation used in example 1).

Figure 9 .
Figure 9.The basin of attraction for the numerical schemes BB1-BB3, AC, YO, and NM technique, respectively, is shown from left to right (for polynomial equation used in example 2).

Figure 10 .
Figure 10.The basin of attraction for the numerical schemes BB1-BB3, AC, YO, and NM technique, respectively, is shown from left to right (for polynomial equation used in example 3).

Table 3
clearly shows that our recently developed methods BB1-BB3 outperform the existing iterative methods in terms of CPU-Time, Average-It, T-Points, C-Points, and D-Points.

Table 1 .
Convergence divergence analysis of numerical schemes for f 1 (x).

Table 2 .
Convergence divergence analysis of numerical schemes for f 2 (x).

Table 3 .
Convergence divergence analysis of numerical schemes for f 3 (x).

Table 4 .
Simple roots for polynomial equations f 1 (x) used in engineering application 1.

Table 5
clearly shows the dominance behavior of BB1-BB3 over the AC, YO, and NM iterative methods in approximated absolute error on various number of iterations n for the non-linear function f 2 (x).Various initial guess values are used in order to approximate all roots of f 2 (x).

Table 5 .
Simple roots for polynomial equations f 2 (x) used in engineering application 2.
for the hydronium ion concentration [H 3 O + ] .If we set x = 10 4 [H 3 O + ] , we obtain the following nonlinear equation f 3

Table 6 .
Simple roots for polynomial equations f 3 (x) used in engineering application 3.