In this research article, we construct a family of derivative free simultaneous numerical schemes to approximate all real zero of non-linear polynomial equation. We make a comparative analysis of the newly constructed numerical schemes with a well-known existing simultaneous method for determining all the distinct real zeros of polynomial equations using computer algebra system Mat Lab. Lower bound of convergence of simultaneous schemes is calculated using Mathematica. Global convergence property of the numerical schemes is presented by taking random starting initial approximation and their convergence history are graphically presented. Some real life engineering applications along with some higher degree polynomials are considered as numerical test problems to show performance and efficiency of the derivative free family of numerical methods with comparison of an existing method of same order in literature. Local computational order of convergence, CPU time, graph of computational order of convergence and residual error graphs elaborate efficiency, robustness and authentication of the suggested family of numerical methods in its domain.
Polynomialssimultaneous iterative methodsrandom initial guesseslower boundlocal computational orderCAS-mathematica and mat labIntroduction
One of the most primal problem of science and engineering is locating the zeros of polynomial of degree k with arbitrary real coefficient.
f(x)=∑i=1kaixi,ak≠0.
Let ξ1,…,ξk denote all the simple zeros of Eq. (1). According to Abel’s impossibility theorem [1] “There is no solution in radical to general polynomial with arbitrary co-efficient of degree five or higher” we therefore look toward numerical schemes to approximate zeros of polynomial Eq. (1).
There exist a lot of numerical schemes in literature which approximate single zeros at a time (see, e.g., [2–4]). Here, we consider the following family of numerical schemes [5]:
The numerical scheme Eq. (2) approximate single zero of Eq. (1) at a time.
Beside these single roots finding methods, mathematicians and engineers are interested in simultaneous numerical schemes which approximate all roots simultaneously. More detail on simultaneous methods, their global convergence and parallel implementation on computer algebra system (CAS) and stability are found in [6,7] and reference cite there in [8–14].
Therefore, the main aim of this research article is to construct a derivative free numerical scheme which approximates all real zero of Eq. (1). Using CAS-Mathematica, we find the lower bound of convergence to verify convergence order theoretically. Computational order of convergence [15] and convergence history for random initial approximations are graphed to show the efficiency and performance of numerical schemes as compared to other existing methods of same order. Log of residual graph, graphs of computational order of convergence and local computational order of convergence [16] support the global convergence behavior of our newly constructed numerical scheme for estimating all real zeros of Eq. (1).
Construction of Numerical Scheme
Corresponding to numerical schemes
y(n)=x(n)-f(x(n))f′(x(n)),
for approximating all zeros of Eq. (1), the numerical method is
This numerical method is well known Weierstrass method [17] for approximating all zeros of polynomial Eq. (1) having local quadratic convergence. Using in analogical way as above, we convert Eq. (2) into simultaneous method as:
where yi(n)=xi(n)-f(xi(n))∏j≠ij=1k(xi(n)-xj(n)) and β∈R. Thus, we have constructed here, a new derivative free family of numerical schemes (abbreviated as MD).
Nedzhibov et al. [18] in 2005, present the following cubic convergence derivative free family of simultaneous numerical schemes (abbreviated as ND) as:
Here, we prove the convergence order of the suggested derivative free family of numerical schemes.
Theorem: Let algebraic polynomial Eq. (1) has k number of simple zerosξ1,…,ξkand for sufficiently close initial guessesx1(0),…,xk(0)of the zeros, then for arbitrary real parameterβ,the numerical scheme MD has third order convergence.
Proof: let ∈i=xi(n)-ξi, ∈i′=yi(n)-ξi, ∈i′′=xi(n+1)-ξi be the errors in xi(n), yi(n) ,xi(n+1) respectively and
Ði=∏j≠ij=1k(xi(n)-ζjxi(n)-xj(n)).
Iterative schemes Eq. (5) can be written as:
xi(n+1)=yi(n)-wi(xi(n))(Ģi+2Ģi2+βĢi3),
where Ģi=f(yi(n))f(xi(n)) and wi(xi(n))=f(xi(n))∏j≠ij=1k(xi(n)-xj(n)).
If we express Eq. (1) as f(xi(n))=(x1(n)-ξ1)…(xk(n)-ξk)=∈i∏j≠ij=1k(xi(n)-ξj), then we have:
f(yi(n))=(y1(n)-ξ1)…(yi(n)-ξi).
Substitution in Eq. (5), we have:
f(yi(n))=∏j=1k(xi(n)-wi(xi(n)-ξj)=(ξi-wi(xi(n)))∏j≠ij=1k(xi(n)-wi(xi(n))-ξj)=∈i(1-Ði)∏j≠ij=1k(xi(n)-wi(xi(n))-ξj).
holds. Using ÐiRi=∏j≠ij=1k(xi(n)-wi(xi(n))ξjxi(n)-xj(n)), we have:
ÐiRi-1=∑t≠ikxt(n)-wi(xi(n))-ξtxi(n)-xt(n)∏j≠it-1xi(n)-wi(xi(n))-ξt(xi(n)-xj(n)),=∑t≠ik∈t-∈iÐixi(n)-xt(n)∏j≠it-1xi(n)-wi(xi(n))-ξt(xi(n)-xj(n))=O(|∈|).
Using CAS-Mathematica for Finding Lower Bound of Convergence
Consider
f(x)=(x-φ1)(x-φ2)(x-φ3)
and the first component H(x^{(n)}) of iterative scheme Eq. (5) to find zeros of Eq. (15), simultaneously. In order to find lower bound of convergence, we have to express the differential of an operator H(x^{(n)}) in terms of their partial derivate of its component as Hi(x):
∂H1(x)∂x1∂H1(x)∂x2∂H1(x)∂x3∂2H1(x)∂x12∂2H1(x)∂x1∂x2∂2H1(x)∂x22∂2H1(x)∂x2∂x3∂3H1(x)∂x13∂3H1(x)∂x12∂x2∂3H1(x)∂x1∂x22∂3H1(x)∂x23∂3H1(x)∂x22∂x3⋮⋮⋮⋮⋮⋯
and to so on.
The lower bound of the convergence is obtained until the first non-zero element of row is found zero (see [19]). The Mathematica program is given for each of the considered method as:
Here, some numerical examples from [20,21] are considered with and without random initial approximation to estimate all real zeros of polynomial equation of higher degree. All the computations are performed using Mat Lab R@2011 with 64 digits floating point arithmetic. We take ∈=10-30 as tolerance and use the following stopping criteria
ei=|xi(n+1)-xi(n)|<∈,
where e_{i} represents the absolute error. We compare our iterative schemes MD with ND of the same convergence order. In all numerical calculations, we take β=2.
Example 1: Consider a stirred tank reactor (SCRT) in which two items A and R are feed in reactor at Q and q-Q rate respectively. A complex reaction
a+r→b;b+r→c;c+r→r;c+r→e, develops the following function
hc2.98(x+2.25)(x+1.45)(x+2.85)2(x+4.35)=1.
We obtained an algebraic polynomial equation of degree 4 by taking ℏc=0 in Eq. (16)f1(x)=x4+11.5x3+47.49x2+83.06325x+51.32366875=0
with exact roots:
ζ1=-1.45,ζ2=-2.85,ζ3=-2.85,ζ4=-4.45
Convergence history and computational order of convergence graph (Figs. 1, 4, 7 and 10) of numerical schemes ND, MD are obtained by taking the following random initial guessed valued in our computer program i.e.,
X1=[0.032601;0.5612;0.88187;0.66918],
where X1=[xi(0),i=1,…,4]. Using random initial guessed value X1, the iterative scheme MD converges to exact zeros after 100 iteration by consuming 16.848368 s CPU time while ND converges after 155 iterations and consumes 37.249147 s.
Shows convergence history of numerical scheme MD for polynomial equation <italic>f</italic><sub>1</sub>(<italic>x</italic>)Shows convergence history of numerical scheme MD for polynomial equation <italic>f</italic><sub>2</sub>(<italic>x</italic>)
<!--<pag>\vskip13pt</pag>-->
Shows convergence history of numerical scheme MD for polynomial equation <italic>f</italic><sub>3</sub>(<italic>x</italic>)
<!--<pag>\vskip13pt</pag>-->
Shows convergence history of numerical scheme ND for polynomial equation <italic>f</italic><sub>1</sub>(<italic>x</italic>)
<!--<pag>\vskip13pt</pag>-->
Shows convergence history of numerical scheme ND for polynomial equation <italic>f</italic><sub>2</sub>(<italic>x</italic>)
<!--<pag>\vskip13pt</pag>-->
Shows convergence history of numerical scheme ND for polynomial equation <italic>f</italic><sub>3</sub>(<italic>x</italic>)
<!--<pag>\vskip10pt</pag>-->
Shows computational order of convergence of numerical scheme MD for polynomial equation <italic>f</italic><sub>1</sub>(<italic>x</italic>)
<!--<pag>\vskip10pt</pag>-->
Shows computational order of convergence of numerical scheme MD for polynomial equation <italic>f</italic><sub>2</sub>(<italic>x</italic>)
<!--<pag>\vskip13pt</pag>-->
Shows computational order of convergence of numerical scheme MD for polynomial equation <italic>f</italic><sub>3</sub>(<italic>x</italic>)
<!--<pag>\vskip13pt</pag>-->
Shows computational order of convergence of numerical scheme ND for polynomial equation <italic>f</italic><sub>1</sub>(<italic>x</italic>)
<!--<pag>\vskip13pt</pag>-->
Convergence rates increase by taking the following initial guessed value:
x1(0)=-1.0,x2(0)=-1.1,x3(0)=-2.2,x4(0)=-3.9.
The numerical results are presented in Tabs. 1–6. In all Tabs. 1–6, CO presents convergence order, CPU, presents CPU time and local computational order of convergence (LCOC) by σ.
Example 2: Consider
f2(x)=x5-3x4-23x3+51x2+94x-120
with exact roots:
ζ1=1,ζ2=-2,ζ3=3,ζ4=-4,ζ5=5.
For convergence history and computational order of convergence graph (Figs. 2, 5, 8 and 11) of numerical schemes ND, MD, we used the following initial guessed valued in our computer program i.e.,
X2=[0.81472;0.90579;0.12699;0.91338;0.63236]
where X2=[xi(0),i=1,…,5]. Using random initial guessed values X2, the iterative scheme MD converges to exact zeros after 92 iteration by consuming 27.056098 s CPU time while ND converges after 98 iterations and consumes 31.444833 s.
Convergence rates increase by taking the following initial guessed value:
x1(0)=0.9,x2(0)=-1.9,x3(0)=2.9,x4(0)=-3.9,x5(0)=4.9
with exact roots:
ζ1=-1,ζ2=1,ζ3=3,ζ4=4,ζ5=7,ζ6=5,ζ7=8.
For convergence history and computational order of convergence graph (Figs. 3, 6, 9 and 12) of numerical schemes ND, MD are obtained by taking the following initial guessed valued in our computer program i.e.,
X3=[0.77029;0.35022;0.66201;0.41616;0.84193;0.83193;0.83292;0.25644],
where X3=[xi(0),i=1,…,7]. Using random initial guessed value X3, the iterative scheme MD converges to exact zero after 19 iterations by consuming 13.7714 s CPU time while ND converges after 23 iterations and consumes 18.120062 s.
Simultaneous determination of all zeros of polynomial <italic>f</italic><sub>1</sub>(<italic>x</italic>)
Shows computational order of convergence of numerical scheme ND for polynomial equation <italic>f</italic><sub>2</sub>(<italic>x</italic>)
<!--<pag>\vskip13pt</pag>-->
Shows computational order of convergence of numerical scheme ND for polynomial equation <italic>f</italic><sub>3</sub>(<italic>x</italic>)
<!--<pag>\vskip13pt</pag>-->
Shows residual fall for iterative method MD and ND for polynomial <italic>f</italic><sub>1</sub>(<italic>x</italic>) respectively
<!--<pag>\vskip13pt</pag>-->
Shows residual fall for iterative method MD and ND for polynomial <italic>f</italic><sub>2</sub>(<italic>x</italic>) respectively
<!--<pag>\vskip13pt</pag>-->
Shows residual fall for iterative method MD and ND for polynomial <italic>f</italic><sub>3</sub>(<italic>x</italic>) respectively
<!--<pag>\vskip13pt</pag>-->
Convergence rates increase by taking the following initial guessed value:
x1(0)=-0.9,x2(0)=0.9,x3(0)=2.5,x4(0)=3.9,x5(0)=6.9,x6(0)=4.5,x7(0)=7.9.
Conclusions
Here, we have developed a family of derivative free method for approximating all real zeros of polynomial. Lower bound of convergence of iterative methods MD and ND are calculated using CAS-Mathematica. Using Mat Lab, we graph convergence history and computational order of convergence. From Tabs. 1–6 and Figs. 1–15, we observe that our method MD is much better in terms of convergence history, computational order of convergence, numerical results, log of residual and local computational order of convergence as compared to ND method.
Funding Statement: The authors received no specific funding for this study.
Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.
ReferencesM. I.Rosen, “Niels hendrik abel and equations of the fifth degree,” O. S.Solaiman and I.Hashim, “An iterative scheme of arbitrary odd order and its basins of attraction for nonlinear systems,” S. A.Sariman and I.Hashim, “New optimal newton-householder methods for solving nonlinear equations and their dynamics,” C.Chun, “Some variants of king’s fourth-order family of methods for nonlinear equations,” N.Rafiq, S.Akram, N. A.Mir and M.Shams, “Study of dynamical behavior and stability of iterative methods for nonlinear equation with application in engineering,” Y. M.Chu, N.Rafiq, M.Shams, S.Akram, N. A.Miret al., “Computer methodologies for the comparison of some efficient derivative free simultaneous iterative methods for finding roots of non-linear equations,” P. D.Proinov and S. I.Ivanov, “On the convergence of halley’s method for multiple polynomial zeros,” O. S.Solaiman, S. A. A.Karim and I.Hashim, “Dynamical comparison of several third-order iterative methods for nonlinear equations,” N. A.Mir, M.Shams, N.Rafiq, S.Akram and M.Rizwan, “Derivative free iterative simultaneous method for finding distinct roots of polynomial equation,” P. D.Proinov and M. T.Vasileva, “On a family of weierstrass-type root-finding methods with accelerated convergence,” M.Shams, N. A.Mir, N.Rafiq, A. O.Almatroud and S. A.Akram, “On dynamics of iterative technique for nonlinear equation with application in engineering,” S. C.Chapra, V. K.Kyncheva, V. V.Yotov and S. I.Ivanov, “Convergence of newton, halley and chebyshev iterative methods as methods for simultaneous determination of multiple polynomial zeros,” S. I.Cholakov, “Local and semi local convergence of Wang-Zheng’s method for simultaneous finding polynomial zeros,” F.Soleymani and D. K. R.Babajee, “Computing multiple zeros using a class of quartically convergent methods,” S.Akram, M.Shams, N.Rafiq and N. A.Mir, “On the stability of weierstrass type method with king’s correction for finding all roots of non-linear function with engineering application,” N. A.Mir, M.Shams, S.Akram and R.Ahmed, “On the family of simultaneous method for finding distinct as well as multiple roots of nonlinear equation,” G. H.Nedzhibov, “Convergence of the modified inverse weierstrass method for simultaneous approximation of polynomial zeros,” M.Shams, N.Rafiq, B.Ahmad and N. A.Mir, “Inverse numerical iterative technique for finding all roots of non-linear equations with engineering applications,” M.Shams, N.Rafiq, N. A.Mir, B.Ahmad, S.Abbasiet al., “On computer implementation for comparison of inverse numerical schemes,” M. R.Farmer, “