# Development of a New Multi-step Iteration Scheme for Solving Non-Linear Models with Complex Polynomiography

Authors: Amanullah Soomro; Amir Naseem; Sania Qureshi; Nasr Al Din Ide
Journal: Complexity
 (2022)

---

## Abstract

The appearance of nonlinear equations in science, engineering, economics, and medicine cannot be denied. Solving such equations requires numerical methods having higher-order convergence with cost-effectiveness, for the equations do not have exact solutions. In the pursuit of efficient numerical methods, an attempt is made to devise a modified strategy for approximating the solution of nonlinear models in either scalar or vector versions. Two numerical methods of second-and sixth-order convergence are carefully merged to obtain a hybrid multi-step numerical method with twelfth-order convergence while using seven function evaluations per iteration, resulting in the efficiency index of about 1.4262. The convergence is also ascertained theoretically, and the asymptotic error constant is computed. Furthermore, various numerical methods of varying orders are used to compare the numerical results obtained with the proposed hybrid method in approximate solution, number of iterations, absolute error, absolute functional value, and the machine time measured in seconds. The obtained results outperformed the chosen methods when applied models arising from physical and natural fields were solved. Finally, to observe the convergence graphically, some complex polynomials are plotted as polynomiographs, wherein the rapid convergence of the proposed method is guaranteed.


---

## Body

## 1. Introduction

Computing the approximate zeros of the nonlinear scalar and vector functions is one of the most important and interesting research areas in the modern age. There are many applications of the root-finding methods in different disciplines of science as well as in arts and economics. With the help of several mathematical techniques, a variety of complex problems in different applied sciences can be modulated in the form of nonlinear equations and then can be readily solved via different root-finding techniques. A root-finding method in mathematics and computer technology is a method for finding zeros, commonly known as ”roots” of continuous functions. A zero of a real-valued or a complex-valued functionf, is a value r such that fr=0. Mostly, the root-finding techniques give approximations to zeros, expressed either as floating-point integers or as tiny isolating intervals, or discs for real or complex roots, because the zeros of a function cannot be calculated accurately with available analytical techniques.Most numerical approaches for root-finding rely on the iteration process, which generates a series of discrete points that should converge towards the root as a limit. These iteration schemes start with one or more estimations of the root as initial inputs, and every iteration of the process generates a more accurate estimation of the exact root [1–4].Since iterations must be ended at some point, these approaches yield estimation to the root rather than a precise solution. Many approaches compute successive approximations by considering an auxiliary function on the values that came before them. Therefore, the limit is a fixed point of the functionf, which is selected to have the solution of the original equation as fixed points and to converge to these fixed points quickly. In pursuing accurate and efficient roof-finding techniques, several techniques have been proposed in the past and current literature from the Newton method through the techniques proposed in the ongoing year.Perhaps, the most commonly used root-finding method is the Newton Rahpson methodN2 [5] with quadratic convergence. Its computational step is shown below that uses two function evaluations: 1 for the function fx itself and 1 for the first-order derivative f′x:(1)xi+1=xi−fxif′xi,i=0,1,2,where f′xi≠0.The researchers in [6] devised a modified version of an existing algorithm aiming at the removal of first-order derivatives. They came up with a two-step method with fourth-order convergence abbreviated by N4. One of the advantages of the algorithm was the use of only four function evaluations per iteration, as depicted in the following computational scheme:(2)yi=xi−fxiκxi,xi+1=yi−fyiκyi−f2yifxi.yi2f3yi,i=0,1,…,where(3)κxi=fxi+fxi−fxifxi,κyi=fyi+fyi−fyifyi,fxi,yi=κyi−κxifyi−fxi.The researchers in [7] devised a three-step method having eighth-order convergence as denoted by W8. Despite being eighth-order convergent, one of the advantages of the algorithm was the use of only four function evaluations per iteration, as depicted in the following computational scheme:(4)yi=xi−fxif′xi,zi=xi−fxif′xi4fxi2−5fxifyi−fyi24fxi2−9fxifyi,i=0,1,2,…,xi+1=zi−fzif′xi1+4fzifxi8fyi4fxi−11fyi+1+fzifyi.A ninth-order convergence for a method as denoted byN9 was proved in [8]. The algorithm merely needs 5 function evaluations per iteration, as depicted in the following computational scheme:(5)yi=xi−fxif′xi,zi=yi−1+fyifxi2fyif′yi,i=0,1,2,…,xi+1=zi−1+2fyifxi2+2fzifyifzif′yi.Another higher-order three-step iterative method as denoted byP10 is in [9]. The method is shown to be three-step that requires only 6 function evaluations per iteration, as depicted in the following computational scheme:(6)yi=xi−fxif′xi,zi=yi−fyif′yi,xx+1=zi−fzif′zi+3fzif′yi5f′yif′zi−f′2yi.forWe have also chosen an iterative process with eleventh-order convergence as denoted byN11 in [10]. The process consists of four steps and requires 7 function evaluations per iteration, as depicted in the following computational scheme:(7)yi=xi−fxif′xi,zi=yi−1+fyifxi2fyif′yi,wi=zi−1+2fyifxi2+2fzifyifzif′yi,i=0,1,2,…,xi+1=wi−fwif′wi.Finally, a three-step iterative method with twelfth-order convergence as denoted byN12 is in [11, 12]. The method consists of 5 function evaluations per iteration, as depicted in the following computational scheme:(8)yi=xi−fxif′xi,zi=yi−xi−yifyifxi−2fyi,i=0,1,2,….,xi+1=zi−fzif′zi1−0.5fzif′zi2.


## 2. Formulation of the Proposed Method

It has been observed in the current literature that new modified root-finding techniques are being proposed because of increasing the efficiency of the existing ones. In search of such algorithms, some researchers have merged two existing iterative methods of convergence orderm and n to obtain a hybrid method with convergence order mn. In this respect, the convergence order is improved. Nonetheless, the computational aspect was ignored, resulting in an increased number of function evaluations in most newly modified hybrid approaches. For example, authors in [13] proposed an iterative third-order method with five function evaluations required per iteration, including another algorithm presented in the same research work with a fourth-order three-step method that requires eight function evaluations. Likewise, authors in [14,15] have employed an excessive number of first-order derivatives, leading to high computational effort and machine time. The primary concern of the present work is to propose a hybrid method with possible higher-order convergence with the minimum number of function evaluations so that the computational cost in terms of arithmetic operations and CPU time be reduced. The proposed method comes from Newton’s method m=2 and a three-step method n=6 in [16,17], leading to produce the proposed approach with convergence order mn=12 while using just seven function evaluations per iteration. The four-step proposed method results in the following computational steps, whose flowchart is depicted in Figure 1:(9)wi=xi−fxif′xi,yi=wi−fwif′wi,zi=yi−fyif′yi,xi+1=yi−fyi+fzif′yi.Figure 1
Flow chart for the proposed four-step hybrid method given in (9).fori=0,1,2,….The methods as mentioned earlier, including the one proposed as the four-step hybrid method in the present article, are compared in Figure2 with each other based on efficiency index p1/n, order of convergence p, and the number of function evaluations n used by each method per iteration.Figure 2
Comparison of the methods under consideration on the basis of efficiency index, order of convergence, and number of function evaluations.


## 3. Convergence Analysis of the Proposed Method

This section has been divided into two sections wherein the convergence of the proposed four-step hybrid method in both scalar and vector form has been discussed in detail. It is noted that the twelfth-order convergence is theoretically verified in each case with the aid of Taylor’s expansion.


### 3.1. Convergence Theory with Scalar Form

In this subsection, we theoretically prove the local order of convergence for the proposed method given in (9).Theorem 1.
Suppose thatα∈Q be the required simple root for a sufficiently differentiable function f:Q⊆ℝ⟶ℝ within an open real interval Q. Then, the proposed four-step numerical method (9) possesses twelfth-order convergence with the error equation given by:(10)ζi+1=2f22f111ζi12+Oζi13,where fi=fiα and ζi=xi−α,i=1,2,…Proof.
Supposeα be a simple root of fxi=0, where xi be the i th approximation for the root by the proposed method (9), and ζi=xi−α be the error term in variable x at the i th iteration step. Employing the single real variable Taylor’s series in [9] for fxi around α, we obtain(11)fxi=f1ζi+f22!ζi2+f33!ζi3+f44!ζi4+Oζi5.
Similarly, using the Taylor’s series for1/f′xi around α, we obtain(12)1f′xi=f1−1−f2ζif12+ζi2f1−f32f1+f22f12+ζi3f1−f46f1+f2f32f12+f3f1−2f22f22f13+Oζi4.
Multiplying (11) and (12), we obtain(13)fxif′xi=−2f4f12−7f3f2f1+6f23ζi412f13−2f3f1−3f22ζi36f12−ζi2f22f1+ζi.
Substituting (13) in the first step of (9), we obtain(14)ηi=2f4f12−7f3f2f1+6f23ζi412f13+2f3f1−3f22ζi36f12+ζi2f22f1,where ηi=wi−α. Using the Taylor’s series for fwi around α, we obtain(15)fwi=f1ηi+f22!ηi2+f33!ηi3+f44!ηi4+Oηi5.
Similarly, using the Taylor’s series for1/f′wi around α, we obtain(16)1f′wi=f1−1−f2ηif12+ηi2f1−f32f1+f22f12+ηi3f1−f46f1+f2f32f12+f3f1−2f22f22f13+Oηi4.
Multiplying (15) and (16), we obtain(17)fwif′wi=−2f4f12−7f3f2f1+6f23ηi412f13−2f3f1−3f22ηi36f12−ηi2f22f1+ηi.
Substituting (17) in the second step of (9), we obtain(18)νi=2f4f12−7f3f2f1+6f23ηi412f13+2f3f1−3f22ηi36f12+ηi2f22f1.where νi=yi−α. Using the Taylor’s series for fyi around α, we obtain(19)fyi=f1νi+f22νi2+f33!νi3+f44!νi4+Oνi5.
Similarly, using the Taylor’s series for1/f′yi around α, we obtain(20)1f′yi=f1−1−f2νif12+νi2f1−f32f1+f22f12+νi3f1−f46f1+f2f32f12+f3f1−2f22f22f13+Oνi4.
Multiplying (19) and (20), we obtain(21)fyif′yi=−2f4f12−7f3f2f1+6f23νi412f13−2f3f1−3f22νi36f12−νi2f22f1+νi.
Substituting (21) in the first step of (9), we obtain(22)ρi=2f4f12−7f3f2f1+6f23νi412f13+2f3f1−3f22νi36f12+νi2f22f1,where ρi=zi−α. Using the Taylor’s series for fzi around α, we obtain(23)fzi=f1ρi+f22!ρi2+f33!ρi3+f44!ρi4+Oρi5.
Substituting (19), (20), and (23) in the fourth step of (9), one obtains(24)ζi+1=νi3f33f1+f24νi52f14+f32νi512f12−f32ρi2νi22f13+f42ρi2νi32f14−νi3f222f12−f2ρi22f1+νi2f22f1+νi4f46f1+νi4f232f13+ρiνi3f23f13+f22ρi2νi2f12−ρiνi2f22f12+ρiνi3f46f1+f3νi6f236f14+ρif2νif1+ρiνi2f32f1−2f22νi5f33f12+f3νi6f436f12−f32νi6f26f12−7νi4f2f312f12+f2νi5f412f12+f2ρi2νi3f412f12+f2ρi2νi2f34f12−ρiνi3f2f3f12−f22ρi2νi3f32f13−ρi.
Finally, using , and (14) and (18), (22) for the above equation, we obtain(25)ζi+1=2f22f111ζi12+Oζi13.
Hence, the twelfth-order convergence of the proposed method P12 given by (9) for the nonlinear functions in single variable fx=0 is proved.


### 3.2. Convergence Theory with Vector Form

(26)wi=xi−f′xi−1f′xi,yi=wi−f′wi−1f′wi,zi=yi−f′yi−1f′yi,xi+1=yi−f′yi−1fyi+fziii=0,1,2,3….fi=Fiα,i=1,2,….Proof.
Supposeα be a simple root of fxi=0, where xi be the i th approximation for the root by the proposed method (9), and ζi=xi−α be the error term in variable x at the i th iteration step. Employing the multi variable Taylor’s series given in the theorem [9] for fxi around α, we obtain(27)fxi=f1ζi+f22!ζi2+f33!ζi3+f44!ζi4+Oζi5.
Again, using the Taylor’s expansion for the inverted Jacobian matrixf′xi−1 around α, we obtain(28)f′xi−1=f1−1−f2ζif12+ζi2f1−f32f1+f22f12+ζi3f1−f46f1+f2f32f12+f3f1−2fi2f22f13+Oζi4.
Multiplying (27) and (28), we obtain(29)f′xi−1fxi=−2f4f12−7f3f2f1+6f23ζi412f13−2f3f1−3f22ζi36f12−ζi2f22f1+ζi.
Substituting (29) in the first step of (9), we obtain(30)ηi=2f4f12−7f3f2f1+6f23ζi412f13+2f3f1−3f22ζi36f12+ζi2f22f1,where ηi=wi−α. Using the Taylor’s series for fwi around α, we obtain(31)fwi=f1ηi+f22!ηi2+f33!ηi3+f44!ηi4+O,Again, using the Taylor’s expansion for the inverted Jacobian matrix f′wi−1 around α, we obtain(32)f′wi−1=f1−1−f2ηif12+ηi2f1−f32f1+f22f12+ηi3f1−f46f1+f2f32f12+f3f1−2f22f22f13+Oηi4.
Multiplying (31) and (32), we obtain(33)f′wi−1fwi=−2f4f12−7f3f2f1+6fi3ηi412f13−2f3f1−3f22ηi36f12−ηi2f22f1+ηi.
Substituting (33) in the second step of (9), we obtain(34)νi=2f4f12−7f3f2f1+6f22ηi412f13+2f3f1−3f22ηi36f12+ηi2f22f1,where νi=yi−α. Using the Taylor’s expansion for fyi around α, we obtain(35)fyi=f1νi+f22νi2+f33!νi3+f44!νi4+Oνi5.
Again, using the Taylor’s expansion for the inverted Jacobian matrixf′yi−1 around α, we obtain(36)f′yi−1=f1−1−f2νif12+νi2f1−f32f1+f22f12+νi3f1−f46f1+f2f32f12+f3f1−2f22f22f13+Oνi4.
Multiplying (35) and (36), we obtain(37)f′yi−1fyi=−2f4f12−7f3f2f1+6f23νi412f13−2f3f1−3f22νi36f12−νi2f22f1+νi.Substituting (37) in the first step of (9), we obtain(38)ρi=2f4f12−7f3f2f1+6f23νi412f13+2f3f1−3f22νi36f12+νi2f22f1.where ρi=zi−α. Using the Taylor’s series for fzi around α, we obtain(39)fzi=f1ρi+f22!ρi2+f33!ρi3+f44!ρi4+Oρi5.
Substituting (35), (36), and (39) in the fourth step of (9), one obtains(40)ζi+1=νi3f33f1+f24νi52f14+f32νi512f12−f32ρi2νi22f13+f24ρi2νi32f14−νi3f222f12−f2ρi22f1+νi2f22f1+νi4f46f1+νi4f232f13+ρiνi3f23f13+f22ρi2νi2f12−ρiνi2f22f12+ρiνi3f46f1+f3νi6f236f14+ρif2νif1+ρiνi2f32f1−2f22νi5f33f13+f3νi6f436f12−f23νi6f26f13−7νi4f2f312f12+f2νi5f412f12+f2ρi2νi3f412f12+f2ρi2νi2f34f12−ρiνi3f2f3f12−f22ρi2νi3f32f13−ρi
Finally, using , and (30) and (34), (38) for the above equation, we obtain(41)ζi+1=2f22f111ζi12+Oζi13.
Hence, the twelfth-order convergence of the proposed multi-step (four-step) hybrid method P12 mentioned by (9) for the nonlinear functions in multi-variable fx=0 case is proved. □


## 3.1. Convergence Theory with Scalar Form

In this subsection, we theoretically prove the local order of convergence for the proposed method given in (9).Theorem 1.
Suppose thatα∈Q be the required simple root for a sufficiently differentiable function f:Q⊆ℝ⟶ℝ within an open real interval Q. Then, the proposed four-step numerical method (9) possesses twelfth-order convergence with the error equation given by:(10)ζi+1=2f22f111ζi12+Oζi13,where fi=fiα and ζi=xi−α,i=1,2,…Proof.
Supposeα be a simple root of fxi=0, where xi be the i th approximation for the root by the proposed method (9), and ζi=xi−α be the error term in variable x at the i th iteration step. Employing the single real variable Taylor’s series in [9] for fxi around α, we obtain(11)fxi=f1ζi+f22!ζi2+f33!ζi3+f44!ζi4+Oζi5.
Similarly, using the Taylor’s series for1/f′xi around α, we obtain(12)1f′xi=f1−1−f2ζif12+ζi2f1−f32f1+f22f12+ζi3f1−f46f1+f2f32f12+f3f1−2f22f22f13+Oζi4.
Multiplying (11) and (12), we obtain(13)fxif′xi=−2f4f12−7f3f2f1+6f23ζi412f13−2f3f1−3f22ζi36f12−ζi2f22f1+ζi.
Substituting (13) in the first step of (9), we obtain(14)ηi=2f4f12−7f3f2f1+6f23ζi412f13+2f3f1−3f22ζi36f12+ζi2f22f1,where ηi=wi−α. Using the Taylor’s series for fwi around α, we obtain(15)fwi=f1ηi+f22!ηi2+f33!ηi3+f44!ηi4+Oηi5.
Similarly, using the Taylor’s series for1/f′wi around α, we obtain(16)1f′wi=f1−1−f2ηif12+ηi2f1−f32f1+f22f12+ηi3f1−f46f1+f2f32f12+f3f1−2f22f22f13+Oηi4.
Multiplying (15) and (16), we obtain(17)fwif′wi=−2f4f12−7f3f2f1+6f23ηi412f13−2f3f1−3f22ηi36f12−ηi2f22f1+ηi.
Substituting (17) in the second step of (9), we obtain(18)νi=2f4f12−7f3f2f1+6f23ηi412f13+2f3f1−3f22ηi36f12+ηi2f22f1.where νi=yi−α. Using the Taylor’s series for fyi around α, we obtain(19)fyi=f1νi+f22νi2+f33!νi3+f44!νi4+Oνi5.
Similarly, using the Taylor’s series for1/f′yi around α, we obtain(20)1f′yi=f1−1−f2νif12+νi2f1−f32f1+f22f12+νi3f1−f46f1+f2f32f12+f3f1−2f22f22f13+Oνi4.
Multiplying (19) and (20), we obtain(21)fyif′yi=−2f4f12−7f3f2f1+6f23νi412f13−2f3f1−3f22νi36f12−νi2f22f1+νi.
Substituting (21) in the first step of (9), we obtain(22)ρi=2f4f12−7f3f2f1+6f23νi412f13+2f3f1−3f22νi36f12+νi2f22f1,where ρi=zi−α. Using the Taylor’s series for fzi around α, we obtain(23)fzi=f1ρi+f22!ρi2+f33!ρi3+f44!ρi4+Oρi5.
Substituting (19), (20), and (23) in the fourth step of (9), one obtains(24)ζi+1=νi3f33f1+f24νi52f14+f32νi512f12−f32ρi2νi22f13+f42ρi2νi32f14−νi3f222f12−f2ρi22f1+νi2f22f1+νi4f46f1+νi4f232f13+ρiνi3f23f13+f22ρi2νi2f12−ρiνi2f22f12+ρiνi3f46f1+f3νi6f236f14+ρif2νif1+ρiνi2f32f1−2f22νi5f33f12+f3νi6f436f12−f32νi6f26f12−7νi4f2f312f12+f2νi5f412f12+f2ρi2νi3f412f12+f2ρi2νi2f34f12−ρiνi3f2f3f12−f22ρi2νi3f32f13−ρi.
Finally, using , and (14) and (18), (22) for the above equation, we obtain(25)ζi+1=2f22f111ζi12+Oζi13.
Hence, the twelfth-order convergence of the proposed method P12 given by (9) for the nonlinear functions in single variable fx=0 is proved.


## 3.2. Convergence Theory with Vector Form

(26)wi=xi−f′xi−1f′xi,yi=wi−f′wi−1f′wi,zi=yi−f′yi−1f′yi,xi+1=yi−f′yi−1fyi+fziii=0,1,2,3….fi=Fiα,i=1,2,….Proof.
Supposeα be a simple root of fxi=0, where xi be the i th approximation for the root by the proposed method (9), and ζi=xi−α be the error term in variable x at the i th iteration step. Employing the multi variable Taylor’s series given in the theorem [9] for fxi around α, we obtain(27)fxi=f1ζi+f22!ζi2+f33!ζi3+f44!ζi4+Oζi5.
Again, using the Taylor’s expansion for the inverted Jacobian matrixf′xi−1 around α, we obtain(28)f′xi−1=f1−1−f2ζif12+ζi2f1−f32f1+f22f12+ζi3f1−f46f1+f2f32f12+f3f1−2fi2f22f13+Oζi4.
Multiplying (27) and (28), we obtain(29)f′xi−1fxi=−2f4f12−7f3f2f1+6f23ζi412f13−2f3f1−3f22ζi36f12−ζi2f22f1+ζi.
Substituting (29) in the first step of (9), we obtain(30)ηi=2f4f12−7f3f2f1+6f23ζi412f13+2f3f1−3f22ζi36f12+ζi2f22f1,where ηi=wi−α. Using the Taylor’s series for fwi around α, we obtain(31)fwi=f1ηi+f22!ηi2+f33!ηi3+f44!ηi4+O,Again, using the Taylor’s expansion for the inverted Jacobian matrix f′wi−1 around α, we obtain(32)f′wi−1=f1−1−f2ηif12+ηi2f1−f32f1+f22f12+ηi3f1−f46f1+f2f32f12+f3f1−2f22f22f13+Oηi4.
Multiplying (31) and (32), we obtain(33)f′wi−1fwi=−2f4f12−7f3f2f1+6fi3ηi412f13−2f3f1−3f22ηi36f12−ηi2f22f1+ηi.
Substituting (33) in the second step of (9), we obtain(34)νi=2f4f12−7f3f2f1+6f22ηi412f13+2f3f1−3f22ηi36f12+ηi2f22f1,where νi=yi−α. Using the Taylor’s expansion for fyi around α, we obtain(35)fyi=f1νi+f22νi2+f33!νi3+f44!νi4+Oνi5.
Again, using the Taylor’s expansion for the inverted Jacobian matrixf′yi−1 around α, we obtain(36)f′yi−1=f1−1−f2νif12+νi2f1−f32f1+f22f12+νi3f1−f46f1+f2f32f12+f3f1−2f22f22f13+Oνi4.
Multiplying (35) and (36), we obtain(37)f′yi−1fyi=−2f4f12−7f3f2f1+6f23νi412f13−2f3f1−3f22νi36f12−νi2f22f1+νi.Substituting (37) in the first step of (9), we obtain(38)ρi=2f4f12−7f3f2f1+6f23νi412f13+2f3f1−3f22νi36f12+νi2f22f1.where ρi=zi−α. Using the Taylor’s series for fzi around α, we obtain(39)fzi=f1ρi+f22!ρi2+f33!ρi3+f44!ρi4+Oρi5.
Substituting (35), (36), and (39) in the fourth step of (9), one obtains(40)ζi+1=νi3f33f1+f24νi52f14+f32νi512f12−f32ρi2νi22f13+f24ρi2νi32f14−νi3f222f12−f2ρi22f1+νi2f22f1+νi4f46f1+νi4f232f13+ρiνi3f23f13+f22ρi2νi2f12−ρiνi2f22f12+ρiνi3f46f1+f3νi6f236f14+ρif2νif1+ρiνi2f32f1−2f22νi5f33f13+f3νi6f436f12−f23νi6f26f13−7νi4f2f312f12+f2νi5f412f12+f2ρi2νi3f412f12+f2ρi2νi2f34f12−ρiνi3f2f3f12−f22ρi2νi3f32f13−ρi
Finally, using , and (30) and (34), (38) for the above equation, we obtain(41)ζi+1=2f22f111ζi12+Oζi13.
Hence, the twelfth-order convergence of the proposed multi-step (four-step) hybrid method P12 mentioned by (9) for the nonlinear functions in multi-variable fx=0 case is proved. □


## 4. Polynomiography

Polynomiography is a process that integrates mathematics and art to create a new type of visual art. The produced graphics result from algorithmic visualization of iterative approaches for solving a polynomial equation. This term was first introduced by Dr. Bahman Kalantari at the start of the 21st century [18]. Dr. Bahman Kalantari’s study on polynomial root-finding, which is an old and traditional discipline that continues to find new implications with each generation of mathematicians and scientists, inspired the concepts of Polynomiography. Dr. Kalantari invented the term ”polynomiography,” which is a mixture of the word ”polynomial” with the suffix ”graphy.” A ”polynomiograph” is a separately produced picture resulting from Polynomiography. It is defined as ”An iterative procedure for producing two-dimensional colored pictures (polynomiographs) that represent polynomials.”In recent years, researchers worked in the field of Polynomiography along with its implementations in other fields. In [19], the authors introduced a new mathematical art with the help of Newton–Ellipsoid method. Gdawiec et al. in [20], presented the visual analysis of Newton’s method with fractional-order derivatives. The authors employed the techniques of coloring by roots and coloring via iterations to study the convergence and dynamical aspects of the processes visualized by polynomiographs.Naseem et al. in [21] presented some new graphical art with the help of newly suggested ninth-order iteration schemes. Scot et al. in [22], presented the basin of attraction for several methods and examined its dependence on their convergence orders. In [23], the authors introduced a new family of eighth-order methods and then drew their basins of attractions by assuming different polynomials. Finally, in [24], the authors generated some new fractal patterns by combing two root-finding methods. The obtained fractal patterns were diverse and had many applications in the textile and ceramic industries.We use a rectangleR∈ℂ along with the dimension −2,2×−2,2, accuracy ε=1×10−3 and the max. no. of iterations T=20 to create the polynomiographs over the complex plane C through the computer software by taking multiple complex polynomials. The color black is allocated to the spots where the method failed to converge. The partitioning of R determines the pixel density of the created visual representations; for example, if we partition the rectangle R into a grid of 2000×2000, the plotting polynomiographs will then have better resolution.For drawing graphical objects in the complex plane, we use the four complex polynomials listed below:(42)q1t=t3−1,q2t=3t3+2t2−t+1,q3t=t4−1,q4t=t2+1t2−2.For coloring the iterations, we employ the colormap given in Figure3.Figure 3
The colormap for drawing polynomiographs.Problem 1.
Polynomiographs for the Polynomialq1t Through Various Methods
In this example, we investigate and compare the dynamical results obtained through different iteration schemes with our presented method by considering the cubic polynomialt3−1 which possesses three distinct simple zeros:1, −1/2−3/2i, −1/2+3/2i. We executed all the methods to achieve the simple zeros of the considered polynomials and the results can be visualized in Figure 4.Figure 4
Polynomiographs forq1t using different methods.
(a)(b)(c)(d)(e)(f)(g)(h)Problem 2.
Polynomiographs for the Polynomialq2t Through Various Methods
In this example, we investigate and compare the dynamical results obtained through different iteration schemes with our presented method by considering the 3rd-degree polynomial3t3+2t2−t+1 which possesses three distinct simple zeros:(43)−20031690,8673344±15043251i.
We executed all methods to achieve the simple zeros of the considered polynomials and the results can be seen in Figure5.Figure 5
Polynomiographs forq2t using different methods.
(a)(b)(c)(d)(e)(f)(g)(h)Problem 3.
Polynomiographs for the Polynomialq3t Through Various Methods
In this example, we investigate and compare the dynamical results obtained through different methods with our presented method by considering the 4th-degree polynomialt4−1 which possesses the following simple zeros:1, −1, i, and −i. We executed all the methods to achieve the simple zeros of the considered polynomials and the results are given in Figure 6.Figure 6
Polynomiograph forq3t using different methods.
(a)(b)(c)(d)(e)(f)(g)(h)Problem 4.
Polynomiographs for the Polynomialq4t Through Various Methods
In this example, we investigate and compare the dynamical results obtained through different iteration methods with our presented method by considering 4th-degree polynomialt2−1t2−2 which possesses the simple zeros:1, −1, 2, and −2. We executed all methods to achieve the simple zeros of the considered polynomials and the results can be visualized in Figure 7.
In all given examples, a detailed graphical analysis of the designed algorithm has been provided via polynomigraphs. For plotting polynomiographs on the complex plane, we take two cubic-degree polynomials namely:q1t, and q2t and two quatric-degree polynomials represented by q3t, and q4t respectively. The plotted graphs tell us about the convergence speed and the iterations performed by the method for drawing these objects. The second characteristic is the dynamics of the iteration scheme. In each polynomiographs, the individual root has been denoted by the blue colored dot. The black colored zones denote the divergence area or deficiency of the method through which the polynomiographs has been plotted. The darker or brighter zones in the provided polynomiographs showing less iterations performed to approximate the solution. One can easily observe the superiority of the proposed method over the others by examining the more darker or brighter zones of the polynomiographs drawn by the suggested method.Figure 7
Polynomiograph forq4t using different methods.
(a)(b)(c)(d)(e)(f)(g)(h)


## 5. Numerical Simulations: Real-world Scenarios

This section of the paper discusses the real-life applications by applying the newly proposed hybrid method. We also present a numerical comparison with other existing most frequently used methods: N2, N4, W8, N9, P10, N11, and N12, whose computational steps are shown in the introductory section above. In each applied model, we set the tolerance to beε=10−100 as the stopping criterion of the iterative process of every method under consideration: xN−xe<ε. Two additional methods with fifth- and sixth-order of convergence are also included for the simulations of a six-dimensional model chosen from the field of neurophysiology based on the reason that some of the methods under discussion in the above introductory section did not prove to be valuable candidates when it comes to the simulations of the nonlinear models presented in the system or vector form.Problem 5.
The Plank’s radiation law in physics explains the spectral density of radiation emitted by a black body in thermal equilibrium at the temperatureT and the condition that there must not be a flow of energy between the body and its surroundings. In other words, the law is introduced to determine the amount of energy density in a black body based on isothermal properties. Moreover, it is sometimes used to estimate the maximal radiation’s wavelength. As described in [25], the maximal wavelength of the radiation may be written in the form of the following nonlinear equation in scalar version:(44)f1x=exp−x+x5−1,where x stands for the maximal wavelength. The exact solution of the above equation is as follows: 0.0.
The Plank’s radiation model described in (44) is simulated with several numerical algorithms while assuming two different initial guesses. The maximum number of iterations in each case is set to be N=4. It can be observed in Table 1 that the accuracy is maximum for P12 at the cost of a slightly higher amount of CPU time, regardless of the initial estimate’s location. The eleventh-order method N11 did not converge with the initial estimate taken to be x0=0.49 while it converges second to P12 when the initial guess x0=0.75 is taken. With this second initial guess, the CPU time consumption also slightly increases for N4 and N11.Table 1
Comparison of several methods with the same number of iterationsN=4 for the Plank’s radiation model given in Problem 5.
MethodxN−xef1xNTimexN−xef1xNTimef1x,x0=0.49f1x,x0=0.75N23.80e-047.24e-084.70e-021.91e-021.85e-044.70e-02N46.34e-521.25e-011.25e-019.98e-356.19e-1399.40e-02W84.32e-2541.24e-20303.10e-022.98e-1126.44e-8966.30e-02N91.67e-2087.90e-18721.60e-022.82e-359.01e-3134.60e-02P107.38e-4581.00e-30003.20e-022.18e-2472.78e-24694.70e-02N11divergence––1.18e-3051.00e-30004.70e-02N125.05e-3179.90e-25347.80e-021.12e-2105.92e-16839.40e-02P123.71e-7370.00e+009.30e-021.34e-3752.00e-30009.30e-02Problem 6.
Fraction Conversion of Nitrogen-Hydrogen to Ammonia [26]. This nonlinear scalar problem depicts the fractional conversion of nitrogen-hydrogen to ammonia and has appeared in several research works conducted in the past and recent literature. In this experiment, we set the pressure value to be 250 atmospheric pressure while the temperature is set to 500 degrees Celsius. In terms of a nonlinear function, the model mentioned above has the following polynomial form:(45)f2x=x4−7.79075x3+14.744x2+2.511x−1.674.
It has been identified in the recent work [27] that one of its positive real roots lies in the open unit interval 0,1 which is estimated to be 0.2777595428.
For this nonlinear model, the numerical simulations are shown in Table2 while the number of iterations for each method under consideration is set to be 7. Two different initial guesses, that is, x0=0.5 (near to the exact solution) and x0=0.95 (away from the exact solution), are chosen. It is seen that the fourth-order method N4 converges towards some other solution for the first initial guess while the method abbreviated as N11 failed after three iterations while the same is the case for the W8 method, but the method W8 produced the correct approximate solution till four iterations and failed after that. The Newton method N2 has the comparatively most significant absolute error at the fourth iteration compared to other methods. Nonetheless, the proposed hybrid method, in addition to a few more methods, always converged towards the required solution. More so, the proposed method has achieved the minor error tolerance with a reasonable amount of time. Hence, it can be concluded that the initial location of the estimate does not matter much when it comes to the proposed hybrid method devised in this article.Table 2
Comparison of several methods with the same number of iterationsN=4 for the Problem 6.
MethodxN−xef2xNTimexN−xef2xNTimef2x,x0=0.5f2x,x0=0.95N22.75e-076.60e-131.60e-021.51e-061.99e-111.60e-02N4other sol.––1.29e-023.53e-049.30e-02W85.37e-3335.04e-26576.20e-029.41e-2094.51e-16639.40e-02N91.86e-4880.00e+001.60e-027.91e-4070.00e+003.10e-02P103.21e-7830.00e+003.10e-024.19e-6910.00e+003.10e-02N11failed––failed––N128.40e-3342.10e-26631.60e-022.04e-2242.51e-17884.60e-02P122.18e-13730.00e+001.50e-021.16e-12130.00e+007.80e-02Problem 7.
Application of mechanical engineering [28]: There are several fields wherein the use of thermodynamics is extensively required. The particular areas include mechanical, civil, mechatronics, electronic, chemical engineering, and many others. The fourth-order polynomial is used to show a relation between the zero-pressure specific heat of dry air cp (kJ/kgK) to temperature x:(46)f3x=1.9520×10−14x4−9.5838×10−11x3+9.721510−8x2+1.671×10−4x−0.99403−1.2,where cp=1.2 is used.
As described in [28], the above nonlinear model given in terms of fourth-order polynomial has two real distinct roots given as: r1=1126.009751 and r2=−1289.950382. It can be observed in Table 3 that each method converges for the initial guesses chosen to determine the approximate solution of the above model. The eleventh-order method N11 performed better from an accuracy viewpoint, followed by the proposed hybrid method. Looking at the CPU values, it is clear that the methods N11 and P12 take the same amount of time, thereby being equally time-efficient for this particular fourth-order polynomial. Moreover, the absolute functional values are identical for both methods, including others such as N9 and P10.Table 3
Comparison of several methods with the same number of iterationsN=3 for the Problem 7.
MethodxN−xef3xNTimexN−xef3xNTimef3x,x0=1126.0f3x,x0=−1286N21.83e-182.61e-433.20e-021.84e-072.25e-203.10e-02N42.44e-476.07e-1947.80e-022.10e-078.04e-341.57e-01W83.41e-3319.35e-26716.30e-021.53e-1727.94e-14027.80e-02N92.22e-4121.00e-30003.10e-023.22e-1887.03e-17154.60e-02P102.80e-5191.00e-30003.10e-022.57e-2414.58e-24373.20e-02N115.53e-16591.00e-30001.60e-022.56e-7613.00e-30004.70e-02N122.42e-3315.66e-26723.10e-023.37e-1591.69e-12933.20e-02P121.72e-7481.00e-30001.60e-021.06e-3453.00e-30004.70e-02Problem 8.
Neurophysiology Application: As a final experiment, we consider a six-dimensional nonlinear system first proposed in [29] and later was used by several researchers for the simulation purpose of their newly developed algorithms. See, for example, [30], and some cited references therein. The nonlinear model consists of the following six equations:(47)x12+x32=1,x22+x42=1,x5x33+x6x43=c1,x5x13+x6x23=c2,x5x1x32+x6x42x2=c3,x5x12x3+x6x22x4=c4.
The constantsci in the above model can be randomly chosen. In our experiment, we consideredci=0,i=1,…,4.
The simulations for the above neurophysiology application system are shown in Table4 wherein the two columns represent the absolute error at the last iteration and the CPU time consumption. Two more methods, including fifth-order Halley’s (HM5) in [6] and a sixth-order Hameer-Mujtaba method (HM6) in [17], is used for the simulations of the above system. It is evident from Table 4 that the accuracy is much higher for the proposed approach in comparison to other competitive methods, while the CPU time is also reasonable.Table 4
Comparison of several methods with the same number of iterationsN=5 while the initial guess is set to be Fx10,x20,x30,x40,x50,x60=1.8,2.6,3.5,1.3,1.0,1.1 for the Problem 8.
Methodε=xN−xeTimeN28.77e-024.70e-02HM58.22e-942.34e-01HM61.56e-2421.57e-01P102.37e-16672.03e-01P125.02e-28231.57e-01


## 6. Concluding Remarks and Future Directions

A new four-step nonlinear method for solvingfx=0 type models is introduced in this research work with twelfth-order convergence, and seven function evaluations are required per iteration. The theoretical order of convergence for the proposed hybrid method is proved under both scalar and vector cases, along with an asymptotic error constant. Comparison with various existing numerical methods discloses the better performance of the proposed approach when the absolute errors, the absolute functional value computed at the last iteration, and the time of machine in seconds are taken into consideration. The proposed method brings out the slightest absolute error regardless of initial conditions chosen for the simulations of the nonlinear models that belong to real-world scenarios from science and engineering. The rapid convergence of the proposed hybrid method is confirmed with the aid of polynomiography when the method is applied to some complex-valued polynomials.Declarations



---

*Source: 2596924-2022-06-27.xml*