Next Article in Journal
Reduced Multiplicative (BURA-MR) and Additive (BURA-AR) Best Uniform Rational Approximation Methods and Algorithms for Fractional Elliptic Equations
Next Article in Special Issue
On the Nonlinear Integro-Differential Equations
Previous Article in Journal
Degenerate Derangement Polynomials and Numbers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Iterative Methods for Solving Nonlinear Equations in Quantum Calculus

by
Gul Sana
1,
Pshtiwan Othman Mohammed
2,*,
Dong Yun Shin
3,*,
Muhmmad Aslam Noor
1 and
Mohammad Salem Oudat
4
1
Department of Mathematics, COMSATS University Islamabad, Islamabad 44000, Pakistan
2
Department of Mathematics, College of Education, University of Sulaimani, Sulaimani 46001, Kurdistan Region, Iraq
3
Department of Mathematics, University of Seoul, Seoul 02504, Korea
4
Administrative Sciences, Accounting and Finance, Applied Science University, Al Eker 1048, Bahrain
*
Authors to whom correspondence should be addressed.
Submission received: 10 May 2021 / Revised: 19 June 2021 / Accepted: 21 June 2021 / Published: 25 June 2021

Abstract

:
Quantum calculus (also known as the q-calculus) is a technique that is similar to traditional calculus, but focuses on the concept of deriving q-analogous results without the use of the limits. In this paper, we suggest and analyze some new q-iterative methods by using the q-analogue of the Taylor’s series and the coupled system technique. In the domain of q-calculus, we determine the convergence of our proposed q-algorithms. Numerical examples demonstrate that the new q-iterative methods can generate solutions to the nonlinear equations with acceptable accuracy. These newly established methods also exhibit predictability. Furthermore, an analogy is settled between the well known classical methods and our proposed q-Iterative methods.

1. Introduction

In most scientific and engineering applications, the problem of finding the solution of nonlinear equations have become an active area of research. Many researchers have explored various order iterative methods to find solutions of the nonlinear equations using various techniques such as homotopy perturbation technique, variational iterative methods and decomposition technique, for details, see [1,2,3,4,5,6,7,8,9,10,11]. Firstly, Traub [12] initiated the study of the iterative methods for the solution of the nonlinear equations and introduced a basic quadratic convergent Newton iterative method for the solution of the nonlinear equations, which have much significance in the literature. Later on, in order to improve efficiency and local order of convergence of the Newtons method, Cordero and Torregrosa [13], Frontini and Sormani [14], Hasanov [15], Weerakoon and Fernando [16] and Ozban [17] have presented different modifications of the Newton’s method using quadrature rules. Daftardar-Gejji and Jafari [18] have used different modifications of Adomian decomposition method [19] and suggested a simple technique that does not need derivative evaluation of the Adomian polynomial, which is major advantage of using this technique over Adomian decomposition method. Saqib and Iqbal [20] and Ali et al. [21,22] have used this decomposition technique and developed a family of iterative methods with better efficiency and convergence order for solving the nonlinear equations. This study shifts the paradigm of determining higher order iterative methods for solving the nonlinear equations towards the q-analogue of the iterative methods in the q-calculus.
In the last quarter of the 20th century, q-calculus appeared as an amalgamation of mathematics and physics (see [23,24,25,26,27,28,29]), and much consideration has been given by many researchers because of its wide range of applications in many fields of mathematics such as combinatorics, theory of relativity, mechanics, number theory and orthogonal polynomials. Firstly, Jackson [30] introduced the q-Taylor’s formula. Then, Jing and Fan [31] derived q-Taylor’s formula with its q-remainder by using the q-differentiation approach and established results on the q-remainder in the q-Taylor’s formula. Ernest [32] presented the four different q-Taylor’s formulas along with q integral remainder. Prashant et al. [33] have used the q-Taylor’s formula and investigated the q-analogue of the iterative methods, particularly the q-analogue of generalized Newton Raphson method and the q-analogue of the Newton Raphson method for the solution of algebraic transcendental equations and compared the accuracy of the results obtained by the classical methods. Many linear and nonlinear models appearing in science and engineering problems can be modeled by using the q differential equations. Jafari et al. [34] have adopted Daftardar decomposition technique for solving the q difference equations and also determined the convergence of the method.
In this study, we determine the q-analogue of the iterative methods proposed and suggested by Noor and Noor [35] and Ullah et al. [36] with the help of the q-Taylor’s series and decomposition technique [14].
Now, we recall some of the basic results in the area of the q-calculus [37] for 0 < q < 1 that will support the development of our proposed q-iterative methods for the solution of the nonlinear equations.
Let q∈(0,1) the q-integer be defined as:
[ n ] q = 1 + q + q 2 + + q n 1 = 1 q n 1 q for n = 1 , 2 ,
[ n ] q = n for q = 1 .
For q factorial and for 0 k n , the q binomials are defined as:
[ n ] q ! = [ n ] q [ n 1 ] q [ 1 ] q [ 0 ] q ! = 1 , n k q = [ n ] q ! [ k ] q ! [ n k ] q ! .
Definition 1
(see [37]). The q-derivative for real valued continuous function f ( x ) is defined as follows
D q ( x ) = d d x q f ( x ) = f ( q x ) f ( x ) q x x , q 1 .
When q 1 , then the q-derivative is reduced to the standard derivative. Furthermore, the q-derivative can be represented as ( D q f ) ( x ) , and it is known as the Jackson Derivative. The higher order q-derivative for the function f ( x ) is given as
D q 0 f = f , D q n f = D q D q n 1 for n = 1 , 2 , 3 , .
Definition 2
(see [37]). The q-derivative of product and quotient of function f ( x ) and g ( x ) is defined as follows
D q ( f ( x ) g ( x ) ) = g ( x ) D q f ( x ) + f ( q x ) D q g ( x ) = g ( q x ) D q f ( x ) + f ( x ) D q g ( x ) , D q f ( x ) g ( x ) = g ( x ) D q f ( x ) f ( x ) D q g ( x ) g ( q x ) g ( x ) such that g ( x ) g ( q x ) 0 .
Definition 3
(see [30,31,32]). Let f ( x ) be a continuous function on some interval ( a , b ) and c [a,b] then Jackson q-Taylor’s formula is given as:
f ( x ) = k = 1 D q k ( x c ) k [ k ] ! x ( a , b ) ,
where,
( x c ) 0 = 1 , ( x c ) k = i = 0 k 1 ( x c q ) i , k N ,
where 0 < q < 1 and D q , D q 2 , are all q-derivatives.
The rest of this article is organized as follows: in Section 2, the structures of the q-iterative methods will be designed by proposing the Lists A–C. In Section 3, we deal with the convergence analysis of the proposed q-iterative methods, and it is established that these methods have the same order of convergence as the classical methods for q = 1 . In Section 4, we present some of the examples to check the efficacy and performance of these methods. Furthermore, the comparisons of the results obtained by q-iterative methods with the previously known iterative methods will be discussed in the same section. Section 5 explores a general form of the q-iterative method based on the proposed iterative methods. Finally, the findings of our article are given in Section 6.

2. Construction of the q -Iterative Methods

In this section, some new different order multi-step q-iterative methods are constructed by considering the Taylor’s series in the q-calculus. Here, we consider the nonlinear equation
f ( x ) = 0 ( x R ) .
Suppose that α is the root of Equation (8) and λ is an initial guess in the neighborhood of α . By the same technique used in [28], we rewrite the nonlinear Equation (8) as the coupled system of equations by using the Taylor’s series in the neighborhood of λ in the q-calculus:
f ( λ ) + ( x λ ) D q f ( λ ) + G ( x ) = 0 ,
G ( x ) = f ( x ) f ( λ ) ( x λ ) D q f ( λ ) .
Since f ( x ) = 0 , the relation (10) can be written as:
x = λ f ( λ ) D q f ( λ ) G ( x ) D q f ( λ ) ,
x = c + M q ( x ) ,
where
c λ f ( λ ) D q f ( λ ) ,
and
M q ( x ) G ( x ) D q f ( λ )
is a nonlinear operator and c is treated as constant.
It is noted that if we consider x 0 as an initial guess, then from (10), we have
f ( x 0 ) = G ( x 0 ) .
It is wort mentioning that Equation (15) plays a very significant role in the development of new multi-step q-iterative methods. Now, we establish a sequence of higher order iterative methods implementing the decomposition technique presented by Daftardar-Gejji and Jafari [18]. The main idea behind the implementation of this technique is to find out the solution of q-type functional Equation (12) in terms of infinite series:
x = i = 0 x i .
Now, we decompose the operator M q ( x ) defined in (14), such as:
M q ( x ) = M q ( x 0 ) + i = 1 M q j = 0 i x j M q j = 0 i 1 x j .
From the Equations (12), (16) and (17), we have
i = 1 x i = c + M q ( x 0 ) + i = 1 M q j = 0 i x j M q j = 0 i 1 x j ,
which generates the following iterative scheme
x 0 = c , x 1 = M q ( x 0 ) , x 2 = M q ( x 0 + x 1 ) M q ( x 0 ) , x n + 1 = M q j = 0 n x j M q j = 0 n 1 x j , n = 1 , 2 , .
It follows that
x 1 + x 2 + + x n + 1 = M q x 0 + x 1 + x 2 + + x n ,
and
x = c + i = 1 x i .
Theorem 1
(see [34]). If M is a contraction mapping, the defined series in (16) is absolutely convergent.
Proof. 
If M is a contraction mapping, we see that
M q ( x ) M q ( y ) α x y 0 < α < 1 ,
then in view of (19), we have
x n + 1 = M q ( x 1 + x 2 + + x n ) M q ( x 1 + x 2 + + x n 1 ) α x n α n x 0 n = 0 , 1 , 2 , .
So, the series x= i = 0 x i converges uniformly and absolutely to the solution of Equation (12) (see [38]). It is noted that x is approximated by
X n = x 0 + x 1 + x 2 + + x n ,
and thus lim x X n = x .
This completes the proof. □
Our iterative techniques proceed with the following algorithms:
List A: 
From (19), we have for n = 0 :
x X 0 = x 0 = c = λ f ( λ ) D q f ( λ ) .
This formulation suggests the following iterative scheme for solving the nonlinear Equation (8). Now, for the given initial guess x 0 , the approximate solution is computed by the iterative representation
x n + 1 = x n f ( x n ) D q f ( x n )
This represents the Newton’s method in the q-calculus which has quadratic convergence. It is proved by Singh et al. [33].
By replacing the value of the q-derivative in (23), we have
x n + 1 = x n f ( x n ) f ( q x n ) f ( x n ) ( q 1 ) x n
= x n ( q 1 ) x n f ( x n ) f ( q x n ) f ( x n ) , n = 0 , 1 , .
This method resembles the method of secants (chords).
Now, with the help of (10) and (19), we get
x 1 = M q ( x 0 ) = G ( x 0 ) D q f ( λ ) = f ( x 0 ) D q f ( λ ) .
List B: 
From (19), we have for n = 1 :
x X 1 = x 0 + x 1 = x 0 + M q ( x 0 ) .
By using (22) and (26), we have
x = λ f ( λ ) D q f ( λ ) ) f ( x 0 ) D q f ( λ ) .
This formulation suggests the following iterative scheme for solving the nonlinear Equation (8).
For the given initial guess x 0 , the approximate solution is computed by the following iterative method:
y n = x n f ( x n ) D q f ( x n ) ,
x n + 1 = y n f ( y n ) D q f ( x n ) , n = 0 , 1 , 2 , .
This is q-analogue of Chun method [3], which has cubic convergence for q = 1 . The error term for this algorithm is computed in Theorem (2).
By using (10), (19) and (26), we can obtain
x 1 + x 2 = M q ( x 0 + x 1 ) = G ( x 0 + x 1 ) D q f ( λ ) = f ( x 0 + x 1 ) f ( λ ) ( x 0 + x 1 λ ) D q f ( λ ) D q f ( λ ) = f ( x 0 + x 1 ) D q f ( λ ) f ( x 0 ) D q f ( λ ) .
List C: 
Again, from (19), we have for n = 2 :
x X 2 = x 0 + x 1 + x 2 = x 0 + M q ( x 0 ) + M q ( x 0 + x 1 ) M q ( x 0 ) = c + M q ( x 0 + x 1 ) .
By using (22) and (31), we get
x = λ f ( λ ) D q f ( λ ) f ( x 0 ) D q f ( λ ) f ( x 0 + x 1 ) D q f ( λ ) .
This formulation allows us to suggest the following iterative method for solving the nonlinear Equation (8).
For the given initial guess x 0 , the approximate solution is computed by the following iterative method:
y n = x n f ( x n ) D q f ( x n ) ,
z n = y n f ( y n ) D q f ( x n ) ,
x n + 1 = z n f ( z n ) D q f ( x n ) , n = 0 , 1 , 2 , .
This is q-analogue of convergent iterative method was investigated by Ullah et al. [36]. Furthermore, it has fourth order convergence for q = 1. The Error equation for this algorithm is computed in Theorem (2).

3. Convergence Analysis

This section is comprised of the convergence analysis of the q-iterative methods determined by Lists B and C in the previous section.
Theorem 2.
Let I R be an open interval and f : I R be a differentiable function. If α I is a simple root of f ( x ) = 0 and x 0 is sufficiently close to α, then the convergence order of multi-step methods determined byLists BandChave convergence of the order at least three and four, respectively, and we write it as [ 3 ; q ] and [ 4 ; q ] , where q represents q-calculus. Furthermore, it satisfies the error equations
e n + 1 = 2 b 2 2 e n 3 + ( 7 b 2 b 3 9 b 2 3 ) e n 4 + O ( e n 5 ) , e n + 1 = 4 b 2 3 e n 4 + ( 20 b 3 b 2 2 26 b 2 4 ) e n 5 + O ( e n 6 ) .
Proof. 
Expanding f ( x n ) and D q f ( x n ) in terms of q-Taylor’s series about α to get
f ( x n ) = D q f ( α ) e n + 1 [ 2 ] ! D q 2 f ( α ) e n 2 + 1 [ 3 ] ! D q 3 f ( α ) e n 3 + ,
f ( x n ) = D q f ( α ) e n + b 2 e n 2 + b 3 e n 3 + b 4 e n 4 + ,
D q f ( x n ) = D q f ( α ) 1 + 2 b 2 e n + 3 b 3 e n 2 + 4 b 4 e n 3 + ,
where
b m = D q m f ( α ) [ m ] ! D q f ( α ) for m = 2 , 3 , and e n = x n α .
In view of the Equations (38) and (39), we get
f ( x n ) D q f ( x n ) = e n b 2 e n 2 2 b 3 2 b 2 2 e n 3 + O e n 4 .
Now, by substituting (41) into (29), we get
y n = x n f q ( x n ) D q f ( x n ) = α + b 2 e n 2 + 2 b 3 2 b 2 2 e n 3 + O e n 4 .
Expanding f ( y n ) in terms of Taylor’s series about α to get
f ( y n ) = D q f ( α ) b 2 e n 2 + 2 b 3 2 b 2 2 e n 3 + 3 b 4 7 b 2 b 3 + 5 b 2 3 e n 4 + O e n 5 .
From (39) and (43), we have
f ( y n ) D q f ( x n ) = b 2 e n 2 + ( 2 b 3 4 b 2 2 ) e n 3 + ( 3 b 4 14 b 2 b 3 + 13 b 2 3 ) e n 4 + O ( e n 5 ) .
By applying (42) and (44) into (30), we get the error term of the List B:
x n + 1 = α + 2 b 2 2 e n 3 + 7 b 2 b 3 9 b 2 3 e n 4 + O e n 5 , e n + 1 = 2 b 2 2 e n 3 + 7 b 2 b 3 9 b 2 3 e n 4 + O e n 5 , z n = α + 2 b 2 2 e n 3 + 7 b 2 b 3 9 b 2 3 e n 4 + O e n 5 .
Expanding f ( z n ) in terms of Taylor’s series about α to get
f ( z n ) = D q f ( α ) { 2 b 2 2 e n 3 + 7 b 2 b 3 9 b 2 3 e n 4 + 10 b 2 b 4 + 6 c 3 2 44 b 3 b 2 2 + 30 b 2 4 e n 5 + O e n 6 } .
By the help of (39) and (46), we can deduce
f ( z n ) D q f ( x n ) = 2 b 2 2 e n 3 + 7 b 2 b 3 13 b 2 3 e n 4 + 10 b 2 b 4 + 6 b 3 2 64 b 3 b 2 2 + 56 b 2 4 e n 5 + O e n 6 .
By applying (45) and (47) into (34), we get the error term of the List C:
x n + 1 = α + 4 b 2 3 e n 4 + 20 b 3 b 2 2 26 b 2 4 e n 5 + O e n 6 ,
e n + 1 = 4 b 2 3 e n 4 + 20 b 3 b 2 2 26 b 2 4 e n 5 + O ( e n 6 ) .
This completes the proof. □

4. Numerical Examples and Comparison Results

This section elaborates on the efficacy of algorithms introduced in this paper with the support of examples. All the numerical experiments are performed with Intel (R) Core [TM] 2 × 2.1 GHz, 12 GB of RAM, and all the codes are written in Maple. We use ε = 10 15 and obtain an approximated simple root rather than the exact based on the exactness ε of the computer.
For the computational work, we use the following stopping criteria:
( i ) | x n + 1 x n | < ε a n d ( i i ) | f ( x n + 1 ) | < ε
Abbreviation CAG is used for classical iterative method and QAG for the q-analogue of classical iterative method, and term div is used for divergence of method.
Recall the classical List 2.2 in [35] ( CAG 1 ), defined by
y n = x n f ( x n ) f ( x n ) and x n + 1 = y n f ( y n ) f ( x n ) ( n = 0 , 1 , 2 , ) .
and the classical List 2.3 in [36] ( CAG 2 ), defined by
y n = x n f ( x n ) f ( x n ) , z n = y n f ( y n ) f ( x n ) , x n + 1 = z n f ( z n ) f ( x n ) , for n = 0 , 1 , 2 , .
For simplicity, we denote the iterative Lists B and C by QAG 1 and QAG 2 , respectively. The computational results are presented in Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6 to elaborate the performance and efficacy of our q iterative methods that is the main motivation of transformation of the classical methods towards the q-iterative methods.
For simplicity, initially in Examples 1–3, we check the performance of q-iterative methods with the classical methods for different values of q up to three iterations. Similarly, we can check the performance of the q-iterative methods for the different values of q for the rest of the iterations until we achieve the desired accuracy.
Example 1
(see [3]). We consider the nonlinear equation:
f ( x ) = x e x 2 sin 2 x + 3 cos x + 5 .
The exact solution for this example is x = 1.207647827130919 . We take x 0 = 2 as an initial guess. Table 1 shows the computation of x i and f ( x i ) for i = 1 , 2 , 3 and different values of q by using QAG 1 .
Proceeding in the way of Table 1, we get x = 1.2076478271 for different values of q, which is the required solution. One can observe from Table 1 that more accurate values of x i s can be obtained when q approaches towards one and for which f ( x i ) tend towards zero. The values of f ( x 1 ) = 2.746138e + 01, f ( x 2 ) = 5.942719e + 00, f ( x 3 ) = 5.488326e − 01 calculated by QAG 1 at q = 0.9999 are closer to zero as compared to the values f ( x 1 ) = 2.747477e + 01, f ( x 2 ) = 5.949309e + 00, f ( x 3 ) = 5.505169e − 01 calculated by CAG 1 . Furthermore, Equation (51) converges towards the root x 5 = 1.2076478271 for q = 0.9999 and f ( x ) = 4.435401e − 12. Table 2 shows the computation of x i and f ( x i ) for i = 1, 2, 3 for different values of q by using the List C ( QAG 2 ) .
Proceeding in the way of Table 2, we get the required solution x = 1.2076478271. From Table 2, we see that if q approaches one, we can obtain more accurate values of x i s for which f ( x i ) s tends towards zero, where i = 1 , 2 , 3 . It is also observe that the values of f ( x 1 ) = 2.117115e + 01, f ( x 2 ) = 2.896627e + 00, f ( x 3 ) = 3.622758e − 02 at q = 0.9999 calculated by QAG 2 are closer to zero as compared to the values of f ( x 1 ) = 2.118172e + 01, f ( x 2 ) = 2.900487e + 00, f ( x 3 ) = 3.644647e − 02 calculated by CAG 2 . Furthermore, Equation (51) converges to the root x 5 = 1.2076478271 for q = 0.9999 and f ( x ) = 1.424874e − 20.
Example 2
(see [39] (Population growth model)). Consider the nonlinear equation
1564000 = 1000000 e λ + 435000 λ e λ 1 .
This equation appears in the mathematical modeling of the growth of population over short periods of time, where λ denotes the constant birth rate of population whose value needs to determined.
For computational work, we take x 0 = 1.5 as an initial estimate. The solution of this example approximated to 16 decimal digits is 0.1009979296857498. In Table 3, we compute the values of x 1 , x 2 , x 3 and f ( x 1 ) , f ( x 2 ) , f ( x 3 ) for different values of q by using QAG 1 .
Proceeding as the way of Table 3, we get the required solution x = 0.1009979296857498 . From Table 3, It can easily observe that we obtain more accurate values of x 1 , x 2 , x 3 when q approaches one and for which f ( x 1 ) , f ( x 2 ) , f ( x 3 ) tend towards zero. The values of f ( x 1 ) = 6.591112e + 05, f ( x 2 ) = 2.736999e + 04, f ( x 3 ) = 4.839200e + 00 at q = 0.9999 computed by QAG 1 are closer to zero as compared to the other values of f ( x 1 ) = 6.592007e + 05, f ( x 2 ) = 2.738286e + 04, f ( x 3 ) = 4.850247e + 00 computed by CAG 1 . Furthermore, Equation (52) converges towards the root x 5 = 0.1009979297 in the fifth iteration for q = 0.9999 and f ( x 5 ) = 4.545278e − 22.
Meantime, in Table 4, we compute the values of x 1 , x 2 , x 3 and f ( x 1 ) , f ( x 2 ) , f ( x 3 ) for different values of q by using QAG 2 .
We can observe from Table 4 that we get more accurate values of x 1 , x 2 , x 3 when q approaches one and for which f ( x 1 ) , f ( x 2 ) , f ( x 3 ) tend towards zero. In addition, the values of f ( x 1 ) = 4.206350e + 05, f ( x 2 ) = 2.169407e + 03, f ( x 3 ) = 3.799508e − 06 calculated by QAG 2 at q = 0.9999 are closer to zero as compared to the values f ( x 1 ) = 4.207000e + 05, f ( x 2 ) = 2.171013e + 03, f ( x 3 ) = 3.859455e − 06 calculated by CAG 2 . Furthermore, Equation (52) converges towards the root x 4 = 0.1009979296857498 in the fourth iteration for q = 0.9999 and f ( x ) = 4.103556e − 22.
Example 3
(see [40]). Consider the van der Waal’s equation
f ( x ) = 0.986 x 3 5.181 x 2 + 9.067 x 5.289 .
This equation is used to interpret the real and ideal gas behavior that has been converted to the non-linear form after choosing the appropriate values of the parameters. Its exact solution is x = 1.92984624284786221849. Here, we take x 0 = 3.10.
By proceeding as the above table, we can obtain the required solution x = 1.9298462428. From Table 5, we can observe that we get more accurate values of x 1 , x 2 , x 3 for which f ( x 1 ) , f ( x 2 ) , f ( x 3 ) tend towards zero when q approaches one. The values of f ( x 1 ) = 4.378340e − 01, f ( x 2 ) = 7.835843e − 02, f ( x 3 ) = 1.283785e − 02 at q = 0.9999 computed by QAG 1 are closer to zero as compared to the values of f ( x 1 ) = 4.379863e − 01, f ( x 2 ) = 7.842636e − 0, f ( x 3 ) = 1.286153e − 02 computed by CAG 1 . Furthermore, Equation (53) converges towards the root x 7 = 1.92984624284786221850 in the seventh iteration for q = 0.9999 and f ( x ) = 3.081547e − 18.
The computational results obtained from Table 6 illustrate the accuracy of the values of x 1 , x 2 , x 3 when q approaches one and for which f ( x 1 ) , f ( x 2 ) , f ( x 3 ) tend towards zero. Moreover, the values of f ( x 1 ) = 3.105124e − 01, f ( x 2 ) = 3.856881e − 02, f ( x 3 ) = 3.623214e − 03 at q = 0.9999 computed by QAG 2 , are closer to zero as compared to f ( x 1 ) = 3.106207e − 01, f ( x 2 ) = 3.860481e − 02, f ( x 3 ) = 3.632441e − 03 computed by CAG 2 . Furthermore, Equation (53) converges towards the root x 6 = 1.92984624284786221849 in the sixth iteration for q = 0.9999 and f ( x ) = 2.935505e − 21.

4.1. Error Analysis and Application of the q-Iterative Methods

Error is considered as the difference between a true value and an estimate (see [39]), or an approximation, it can easily be observed from the numerical values that if we compute the errors of q-iterative methods than they fluctuate for different values of q. The error decreases when q approaches to the extreme values between 0 and 1. In view of this result, the q-iterative methods are calculated for a large value of q = 0.99999 , which will approximate the ordinary iterative methods.
List B ( QAG 1 )
Nonlinear equationTrue SolutionApproximate solution at q = 0.9999Error
Equation (51)1.207647827130918927011.207647827130918927007.187741e-21
Equation (52)0.100997929685749788950.100997929685749788953.394397e-28
Equation (53)1.929846242847862218491.929846242847862182833.565288e-17
List C ( QAG 2 )
Nonlinear equationTrue SolutionApproximate solution at q = 0.9999Error
Equation (51)1.207647827130918927011.207647827130918927017.016521e-22
Equation (52)0.100997929685749788950.100997929685749788953.064831e-28
Equation (53)1.929846242847862218491.929846242847862218453.396321e-20

4.2. Comparison of the Classical and q-Analogue of Iterative Methods

Here, we check the robustness and efficiency of our new iterative methods by considering some of the nonlinear equations. Furthermore, we compare the standard Newton’s method ( NM ), fourth order Chun method ( CM ) (see [3]), CAG 1 and CAG 2 with our new iterative methods QAG 1 and QAG 2 . In Table 7, we display the number of iterations ( IT ), the approximate root x n , the value f ( x n ) and δ be the distance between two successive estimations. It is important to mention that in order to get better computational results of the q-iterative methods, we take the value of q = 0.9999 .
f 1 ( x ) = x e x 2 sin 2 x + 3 cos ( x ) + 5 , f 2 ( x ) = 0.1056885297 x 2 + e x e x 2 sin ( x ) , f 3 ( x ) = 0.986 x 3 5.181 x 2 + 9.067 x 5.289 , f 4 ( x ) = 1,000,000 e x + 435,000 e x 1 x 1,564,000 , f 5 ( x ) = x 3 + 4 x 2 10 , f 6 ( x ) = x 2 e x 3 x + 2 , f 7 ( x ) = x 2 ( 1 x ) 5 , f 8 ( x ) = x 3 9 x 2 + 24 x 20 1 / 3 + e x / 2 , f 9 ( x ) = x 3 3 x 2 3 + log ( x + 1 ) .
Table 7 compares the solutions obtained by using the classical and our the q-iterative methods. The results show that our q-analogue iterative methods QAG 1 and QAG 2 give the same results as the classical methods CAG 1 and CAG 2 . The functions f 8 ( x ) , f 9 ( x ) are not differentiable at x = 2 and x = 3 , respectively, when we choose x = 2 and x = 3 as initial guesses for f 8 ( x ) and f 9 ( x ) , respectively. Then, the iterative methods: NM , CM , CAG 1 and CAG 2 fail, while the new iterative methods QAG 1 and QAG 2 are applicable and give rapid convergent results.
Remark 1
(see [12]). The efficiency index is considered as P 1 m , where P represents the order of the method and m is the total number of function evaluations per iteration necessary by the method.
  • Efficiency index of CAG 1 is 3 1 3 = 1.442250 .
  • Efficiency index of CAG 2 is 4 1 4 = 1.414214 .
  • Efficiency index of QAG 1 is 3 1 2 = 1.442250 .
  • Efficiency index of QAG 2 is 4 1 4 = 1.414214 .
We conclude that the efficiency indexes calculated by QAG 1 and QAG 2 are the same as those calculated by CAG 1 and CAG 2 .

5. Generalization of the Iterative Scheme in q -Calculus

In this section, based on our previous results, we determine the generalized q-iterative scheme. By adding the values of x i s ( i = 1 , 2 , , n ) in (19), we obtain
x 1 + x 2 + + x n = x 0 + M q x 0 + x 1 + + x n 1 .
From (22), we have
x 0 = λ f ( λ ) D q f ( λ ) , x 0 + x 1 = x 0 f ( x 0 ) D q f ( λ ) = x 0 + M q ( x 0 ) , x = x 0 + x 1 + x 2 = x 0 + x 1 f ( x 0 + x 1 ) D q f ( λ ) = x 0 + M q ( x 0 + x 1 ) .
Now, if x is approximated by
x = x 0 + x 1 + + x n = x 0 + M q ( x 0 + x 1 + + x n 1 ) x = x 0 + x 1 + + x n 1 f ( x 0 + x 1 + + x n 1 ) D q f ( λ ) x = x 0 + x 1 + + x n = x 0 + x 1 + + x n 1 f ( x 0 + x 1 + + x n 1 ) D q f ( λ ) .
Therefore, (54) gives the following iterative scheme:
x n = f ( x 0 + x 1 + + x n 1 ) D q f ( λ ) ( n = 1 , 2 , 3 , ) ;
x 0 = λ f ( λ ) D q f ( λ ) , x 0 + x 1 = x 0 f ( x 0 ) D q f ( λ ) , x 0 + x 1 + x 2 = x 0 + x 1 f ( x 0 + x 1 ) D q f ( λ ) , x 0 + x 1 + x 2 + + x n 1 + x n = x 0 + x 1 + x 2 + + x n 1 f ( x 0 + x 1 + x 2 + + x n 1 ) D q f ( λ ) .
This formulation allows us to suggest the following generalized q iterative method:
y 0 = x m , y 1 = y 0 f ( y 0 ) D q f ( y 0 ) , y 2 = y 1 f ( y 1 ) D q f ( y 0 ) , y 3 = y 2 f ( y 2 ) D q f ( y 0 ) , y n + 1 = y n f ( y n ) D q f ( y 0 ) , x m + 1 = y n + 1 .
The order of convergence of the iterative scheme is n + 2 for n = 0 , 1 , 2 , and the number of functional evaluations is also n + 2 as well.

6. Conclusions

The main target of this article is to introduce novel algorithms for solution of the nonlinear equations in the context of the q-calculus. The new algorithms are introduced by using the Daftardar–Jafari decomposition technique. The comparison of these newly established algorithms with the classical methods reflects that the proposed q-iterative methods are reliable and best alternatives to the already known algorithms. The computational results conclude that the q-analogue of the iterative methods for solving the algebraic nonlinear equations generate the same results as the classical methods, but convergence rate towards approaching the root is higher than convergence rate suggested by the classical methods. Furthermore, the errors associated with the proposed methods are comparatively lesser by the appropriately chosen value of q being close to one. The difficulty in this method, which needs future investigation, is that we need to estimate the value of q ( 0 , 1 ) .

Author Contributions

Conceptualization, G.S., P.O.M. and M.A.N.; methodology, G.S., P.O.M. and D.Y.S.; software, P.O.M., D.Y.S. and M.A.N.; validation, P.O.M., D.Y.S. and M.S.O.; formal analysis, P.O.M. and D.Y.S.; investigation, P.O.M.; resources, P.O.M. and M.A.N.; data curation, G.S., P.O.M. and D.Y.S.; writing—original draft preparation, G.S. and D.Y.S.; writing—review and editing, G.S., P.O.M. and M.S.O.; visualization, D.Y.S.; supervision, P.O.M., M.A.N. and D.Y.S. All authors have read and agreed to the final version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to express our sincere gratitude to the anonymous referees for their helpful comments that will help to improve the quality of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mohammed, P.O.; Machado, J.A.T.; Guirao, J.L.G.; Agarwal, R.P. Adomian decomposition and fractional power series solution of a class of nonlinear fractional differential equations. Mathematics 2021, 9, 1070. [Google Scholar] [CrossRef]
  2. Mohammed, P.O.; Alqudah, M.A.; Hamed, Y.S.; Kashuri, A.; Abualnaja, K.M. Solving the modified regularized long wave equations via higher degree B-spline algorithm. J. Funct. Space 2021, 2021, 5580687. [Google Scholar]
  3. Chun, C. Iterative methods improving Newton’s method by the decomposition method. Comput. Math. Appl. 2005, 50, 1559–1568. [Google Scholar] [CrossRef] [Green Version]
  4. Noor, M.A.; Waseem, M.; Noor, K.I.; Ali, M.A. New iterative technique for solving nonlinear equations. Appl. Math. Comput. 2015, 265, 1115–1129. [Google Scholar]
  5. Hamasalh, F.K.; Mohammed, P.O. Generalized quartic fractional spline interpolation with applications. Int. J. Open Probl. Compt. Math. 2015, 8, 67–80. [Google Scholar] [CrossRef]
  6. Solaiman, O.S. Two new efficient sixth order iterative methods for solving nonlinear equations. J. King. Saud Uni. 2019, 31, 701–705. [Google Scholar] [CrossRef]
  7. Sana, G.; Noor, M.A.; Noor, K.I. Some multistep iterative methods for nonlinear equation using quadrature rule. Int. J. Anal. Appl. 2020, 18, 920–938. [Google Scholar]
  8. Alqudah, M.A.; Mohammed, P.O.; Abdeljawad, T. Solution of singular integral equations via Riemann–Liouville fractional integrals. Math. Prob. Eng. 2020, 2020, 1250970. [Google Scholar] [CrossRef]
  9. Wu, J.; Yuan, J.; Gao, W. Analysis of fractional factor system for data transmission in SDN. Appl. Math. Nonlinear Sci. 2019, 4, 191–196. [Google Scholar] [CrossRef] [Green Version]
  10. Kurt, A.; Şenol, M.; Tasbozan, O.; Chand, M. Two reliable methods for the solution of fractional coupled Burgers’ equation arising as a model of Polydispersive sedimentation. Appl. Math. Nonlinear Sci. 2019, 4, 523–534. [Google Scholar] [CrossRef] [Green Version]
  11. Touchent, K.A.; Hammouch, Z.; Mekkaoui, T. A modified invariant subspace method for solving partial differential equations with non-singular kernel fractional derivatives. Appl. Math. Nonlinear Sci. 2020, 5, 35–48. [Google Scholar] [CrossRef]
  12. Traub, J.F. Iterative Methods for Solution of Equations; Prentice-Hall: Englewood Cliffs, NJ, USA, 1964. [Google Scholar]
  13. Cordero, A.; Torregrosa, J.R. Variants of Newton’s method using fifth-order quadrature formulas. Appl. Math. Comput. 2007, 190, 686–698. [Google Scholar]
  14. Frontini, M.; Sormani, E. Some variants of Newtons method with third order convergence. Appl. Math. Comput. 2003, 140, 419–426. [Google Scholar]
  15. Hasanov, V.I.; Ivanov, I.G.; Nedzhibov, G. A new modification of Newton method. Appl. Math. Eng. 2002, 27, 278–286. [Google Scholar]
  16. Weerakoon, S.; Fernando, T.G.I. A variant of Newton’s method with accelerated third-order convergence. Appl. Math. Lett. 2000, 13, 87–93. [Google Scholar] [CrossRef]
  17. Ozban, A.Y. Some New variants of Newton’s method. Appl. Math. Lett. 2004, 17, 677–682. [Google Scholar] [CrossRef] [Green Version]
  18. Daftardar-Gejji, V.; Jafari, H. An iterative method for solving nonlinear functional equations. J. Math. Anal. Appl. 2006, 316, 753–763. [Google Scholar] [CrossRef] [Green Version]
  19. Adomian, G. Nonlinear Stochastic Systems and Applications to Physics; Springer Science & Business Media: Dordrecht, The Netherlands, 1989; Volume 46. [Google Scholar]
  20. Saqib, M.; Iqbal, M. Some multi-step iterative methods for solving nonlinear equations. Open J. Math. Sci. 2017, 1, 25–33. [Google Scholar] [CrossRef]
  21. Ali, F.; Aslam, W.; Ali, K.; Anwar, M.A.; Nadeem, A. New family of iterative methods for solving nonlinear models. Discret. Dyn. Nat. Soc. 2018, 2018, 9619680. [Google Scholar] [CrossRef] [Green Version]
  22. Ali, F.; Aslam, W.; Khalid, I.; Nadeem, A. Iteration methods with an auxiliary function for nonlinear equations. J. Math. 2020, 2020, 7356408. [Google Scholar] [CrossRef]
  23. Ernst, T. A New Notation for q-Calculus a New q-Taylor’s Formula; UUDM Report; Department of Mathematics, Uppsala University: Uppsala, Sweden, 1999; pp. 1–28. [Google Scholar]
  24. Koelink, E. Eight lectures on quantum groups and q-special functions. Rev. Colomb. Mat. 1996, 30, 93–180. [Google Scholar]
  25. Alqudah, M.A.; Kashuri, A.; Mohammed, P.O.; Abdeljawad, T.; Raees, M.; Anwar, M.; Hamed, Y.S. Hermite-Hadamard integral inequalities on coordinated convex functions in quantum calculus. Adv. Differ. Equ. 2021, 2021, 264. [Google Scholar] [CrossRef]
  26. Eryılmaz, A. Spectral analysis of q-sturm-liouville problem with the spectral parameter in the boundary condition. J. Funct. Spaces 2012, 2012, 736437. [Google Scholar] [CrossRef] [Green Version]
  27. Erzan, A. Finite q-differences and the discrete renormalization group. Phys. Lett. A 1997, 4–6, 235–238. [Google Scholar] [CrossRef]
  28. He, J.H. A new iteration method for solving algebraic equations. Appl. Math. Comput. 2003, 135, 81–84. [Google Scholar] [CrossRef]
  29. Koornwinder, T.H.; Swarttouw, R.F. On q-analogues of the Fourier and Hankel transforms. Trans. Am. Math. Soc. 1992, 333, 445–461. [Google Scholar]
  30. Jackson, F.H. A q-form of Taylors formula. Mess. Math. 1909, 38, 62–64. [Google Scholar]
  31. Jing, S.C.; Fan, H.Y. q-Taylor’s Formula with its q remainder. Commun. Theor. Phys. 1995, 23, 117–120. [Google Scholar] [CrossRef]
  32. Ernst, T. A method for q-calculus. J. Nonlinear Math. Phys. 2003, 10, 487–525. [Google Scholar] [CrossRef] [Green Version]
  33. Singh, P.; Mishra, P.K.; Pathak, R.S. q-iterative methods. IOSR J. Math. 2013, 9, 6–10. [Google Scholar] [CrossRef]
  34. Jafari, H.; Johnston, S.J.; Sani, S.M.; Baleanu, D. A decomposition method for solving q-difference equations. Appl. Math. Inf. Sci. 2015, 9, 2917–2920. [Google Scholar]
  35. Noor, M.A.; Noor, K.I. Three step iterative methods for nonlinear equations. Appl. Math. Comput. 2006, 183, 322–327. [Google Scholar]
  36. Ullah, M.Z.; Ahmad, F.; Jbbar, M.A.A.A. A correction note on three-step iterative methods for nonlinear equations and generalization of method. J. Mod. Methods Num. Anal. 2014, 5, 10–16. [Google Scholar] [CrossRef] [Green Version]
  37. Kac, V.; Cheung, P. Quantum Calculus; Springer: New York, NY, USA, 2002. [Google Scholar]
  38. Cherruault, Y. Convergence of Adomians method. Kybernetes 1989, 18, 31–38. [Google Scholar] [CrossRef]
  39. Burden, R.L.; Faires, J.D. Numerical Analysis; PWS Publishing Company: Boston, MA, USA, 2001. [Google Scholar]
  40. Waals, V.D.; Diderik, J. Over de Continuiteit van den Gasen Vloeistoftoestand (on the Continuity of the Gas and Liquid State). Ph.D. Thesis, University of Leiden, Leiden, The Netherlands, 1873. [Google Scholar]
Table 1. Calculation of x i and f ( x i ) for i = 1 , 2 , 3 and different values of q by using QAG 1 .
Table 1. Calculation of x i and f ( x i ) for i = 1 , 2 , 3 and different values of q by using QAG 1 .
q x 1 f ( x 1 ) x 2 f ( x 2 ) x 3 f ( x 3 )
1.021.72502069313.025329e+011.44748352537.379764e+001.25186467819.612675e-01
1.011.71461186922.883931e+011.43058591396.635881e+001.24206816597.368529e-01
0.99991.70399727272.746138e+011.41384444325.942719e+001.23362429685.488326e-01
0.981.68279737612.489016e+011.38197104964.732768e+001.22100482072.767866e-01
0.971.67200531792.366777e+011.36658772974.195734e+001.21649759031.821333e-01
0.991.69349715182.615875e+011.39779061745.316156e+001.22670619333.983936e-01
0.961.66112199702.249026e+011.35167792883.701888e+001.21306828671.109783e-01
0.951.65014773702.135624e+011.33727822673.248313e+001.21058738725.995948e-02
0.91.59388586021.628822e+011.27415719561.498487e+001.20749406973.121695e-03
0.81.47373141298.631573e+001.20378186547.805411e-021.20707828991.155594e-02
0.71.34265986413.415231e+001.20987453414.537035e-021.20834004931.407187e-02
0.61.20213585291.110139e-011.20432724206.709762e-021.20568583063.972595e-02
0.51.05516399182.509831e+000.80002700365.058134e+005.39210252872.284390e+13
0.41.26182152341.196531e+001.24499904778.032748e-011.23951895696.795706e-01
Table 2. Calculation of x i and f ( x i ) for i = 1 , 2 , 3 and different values of q by using QAG 2 .
Table 2. Calculation of x i and f ( x i ) for i = 1 , 2 , 3 and different values of q by using QAG 2 .
q x 1 f ( x 1 ) x 2 f ( x 2 ) x 3 f ( x 3 )
1.021.66949913162.339181e+011.35319606263.751021e+001.21274902001.043912e-01
1.011.65896040012.226275e+011.33914001003.305714e+001.21077715676.384855e-02
0.99991.64830762552.117115e+011.32560852212.896627e+001.20942701953.622758e-02
0.991.63786051592.014678e+011.31301927252.532342e+001.20858079011.897265e-02
0.981.62730475491.915594e+011.30100273942.198535e+001.20807374848.654896e-03
0.971.61674709841.820720e+011.28970821271.896564e+001.20780486523.189788e-03
0.961.60618791891.729859e+011.27915065091.624142e+001.20768480527.509703e-04
0.951.59562678201.642815e+011.26934118911.379186e+001.20764652972.634743e-05
0.91.54273008991.258434e+011.23162176515.049568e-011.20760686618.317613e-04
0.81.43514916766.832250e+001.20612948703.076338e-021.20772955211.659826e-03
0.71.32234547502.800732e+001.20663487702.053914e-021.20783277673.756890e-03
0.61.20297034719.432425e-021.20979431094.373048e-021.20673901111.843054e-02
0.51.07809191742.195871e+001.41722303456.079244e+001.20811623629.518872e-03
0.40.94953375333.745734e+000.56419974506.473432e+00divdiv
Table 3. The computed values of x i and f ( x i ) for i = 1 , 2 , 3 and different values of q by using QAG 1 .
Table 3. The computed values of x i and f ( x i ) for i = 1 , 2 , 3 and different values of q by using QAG 1 .
q x 1 f ( x 1 ) x 2 f ( x 2 ) x 3 f ( x 3 )
1.020.51427916966.772056e+050.12320182983.004429e+040.10100350737.468324e+00
1.010.50971378126.681808e+050.12221214572.869170e+040.10100245056.053356e+00
0.99990.50510562376.591112e+050.12124417012.736999e+040.10100154374.839200e+00
0.990.50059165346.502654e+050.12032566452.611695e+040.10100079323.834194e+00
0.980.49603515996.413748e+050.11942780752.489313e+040.10100015732.982748e+00
0.970.49148187336.325288e+050.11855949522.371057e+040.10099962922.275633e+00
0.960.48693191066.237276e+050.11772019152.256843e+040.10099919521.694573e+00
0.950.48238538616.149710e+050.11690935832.146590e+040.10099884321.223128e+00
0.90.45970816775.718559e+050.11326332281.651858e+040.10099796334.506777e-02
0.80.41467012304.889359e+050.10779505219.130505e+030.10099810602.360760e-01
0.70.37012145184.103351e+050.10431322614.446092e+030.10099844056.839589e-01
0.60.32611909903.358939e+050.10231013771.758120e+030.10099837685.987359e-01
0.50.28269230502.654100e+050.10132835474.425044e+020.10099811732.511882e-01
0.40.23984332661.986490e+050.10097777532.698620e+010.10099791292.244417e-02
Table 4. The computed values of x i and f ( x i ) for i = 1 , 2 , 3 and different values of q by using QAG 2 .
Table 4. The computed values of x i and f ( x i ) for i = 1 , 2 , 3 and different values of q by using QAG 2 .
q x 1 f ( x 1 ) x 2 f ( x 2 ) x 3 f ( x 3 )
1.020.38362458134.338095e+050.10287226792.511938e+030.10099792973.473304e-05
1.010.37985721734.272299e+050.10274143352.336453e+030.10099792971.377622e-05
0.99990.37606755364.206350e+050.10261687642.169407e+030.10099792973.799508e-06
0.990.37236804284.142198e+050.10250185972.015174e+030.10099792973.969817e-07
0.980.36864646344.077891e+050.10239247371.868507e+030.10099792977.033782e-08
0.970.36494037244.014077e+050.10228960621.730594e+030.10099792977.550798e-08
0.960.36124987873.950754e+050.10219295731.601031e+030.10099792971.396345e-06
0.950.35757508883.887920e+050.10210223631.479425e+030.10099792975.200839e-06
0.90.33944023763.581029e+050.10172801129.779077e+020.10099792967.367781e-05
0.80.30440203403.002834e+050.10128219503.806792e+020.10099792953.115139e-04
0.70.27106675002.470282e+050.10108646591.185540e+020.10099792943.461354e-04
0.60.23948527511.981025e+050.10101602632.423142e+010.10099792961.702203e-04
0.50.20968436221.532571e+050.10099845487.031754e-010.10099792979.707395e-06
Table 5. The calculated values of x i and f ( x i ) for i = 1 , 2 , 3 and different values of q by using QAG 1 .
Table 5. The calculated values of x i and f ( x i ) for i = 1 , 2 , 3 and different values of q by using QAG 1 .
q x 1 f ( x 1 ) x 2 f ( x 2 ) x 3 f ( x 3 )
0.99992.52025891474.378340e-012.19495469687.835843e-022.02017482571.283785e-02
0.992.51151053204.227815e-012.18313731077.175670e-022.00875783771.058966e-02
0.982.50250587154.076474e-012.17098043836.533731e-021.99716096478.509070e-03
0.972.49332597753.925895e-012.15858827875.917107e-021.98552081326.617240e-03
0.952.47441097803.627241e-012.13302698164.760752e-021.96224255753.387207e-03
0.92.42344830302.897601e-012.06326269472.326221e-021.90801087941.646041e-03
0.82.29935952031.537076e-011.86345636913.701777e-031.90495397391.839852e-03
0.72.12334178174.361666e-021.73020415264.077219e-031.75845458624.282942e-03
0.61.86388675173.689049e-031.88510912082.899285e-031.90244866391.992358e-03
0.51.56333904259.364898e-031.57917220077.962836e-031.59317689791.567623e-03
0.41.38933198184.831340e-021.43420009103.331323e-021.46792079972.455640e-02
Table 6. The calculated values of x i and f ( x i ) for i = 1 , 2 , 3 and different values of q by using QAG 2 .
Table 6. The calculated values of x i and f ( x i ) for i = 1 , 2 , 3 and different values of q by using QAG 2 .
q x 1 f ( x 1 ) x 2 f ( x 2 ) x 3 f ( x 3 )
1.022.45419064393.324888e-012.12949716234.612880e-021.97931325585.685628e-03
1.012.44654602543.215032e-012.11999490264.228580e-021.97165116174.607177e-03
0.99992.43872473273.105124e-012.11030295363.856881e-021.96413220143.623214e-03
0.992.43095732662.998435e-012.10070674363.508499e-021.95704286102.760648e-03
0.982.42300625792.891735e-012.09091187723.172520e-021.95025018771.991531e-03
0.972.41494523142.786122e-012.08100814262.852389e-021.94393965801.325652e-03
0.962.40676942112.681611e-012.07098682092.547916e-021.93825767827.649241e-04
0.952.39847333262.578211e-012.06083646252.258868e-021.93339498043.134112e-04
0.92.35493749812.077991e-012.00741639671.033870e-021.93086937108.898415e-05
0.82.25245657791.158895e-011.85144309784.009256e-031.93366728463.380165e-04
0.72.10599494693.698096e-021.80461335984.487988e-031.86791477773.563896e-03
0.61.86578720513.631368e-031.89536880052.393342e-031.91501644121.168978e-03
0.51.56935218508.800824e-031.59119714237.070040e-031.60987117735.951690e-03
0.41.42477373463.613315e-021.47484901032.300092e-021.51053458771.619046e-02
Table 7. The methods and their IT , x n , f ( x n ) and δ .
Table 7. The methods and their IT , x n , f ( x n ) and δ .
f 1 ( x ) = x e x 2 sin 2 x + 3 cos x + 5 , x 0 = 2
Methods IT x n f ( x n ) δ = | x n x n 1 |
NM 81.20764782713091895.538944e-204.261204e-11
CM61.20764782713091891.063857e-174.175854e-10
CAG 1 61.20764782713091896.780624e-341.948602e-12
QAG 1 51.20764782713091894.435401e-127.441078e-05
CAG 2 51.20764782713091891.009799e-371.383883e-10
QAG 2 51.20764782713091891.424874e-201.175311e-10
f 2 ( x ) = 0.1056885297 x 2 + e x e x 2 sin ( x ) , x 0 = 1.2
NM 90.31706177457295715.220965e-211.571531e-10
CM70.31706177457295712.444123e-464.508430e-16
CAG 1 60.31706177457295751.510838e-171.782637e-06
QAG 1 60.31706177457301571.964445e-151.741348e-06
CAG 2 50.31706177457295713.494619e-191.009576e-05
QAG 2 50.31706177457295712.502232e-191.000084e-05
f 3 ( x ) = 0.986 x 3 5.181 x 2 + 9.067 x 5.289 , x 0 = 3.10
NM 101.92984624284786221.815641e-255.866916e-13
CM71.92984624284786216.121465e-182.408838e-09
CAG 1 71.92984624284786221.803102e-283.036789e-10
QAG 1 71.92984624284786223.081547e-182.566058e-11
CAG 2 61.92984624284786227.761545e-405.605978e-11
QAG 2 61.92984624284786222.935505e-212.073804e-11
f 4 ( x ) = 1,000,000 e x + 435,000 e x 1 x 1,564,000 , x 0 = 1.5
NM 70.10099792968574987.697779e-311.104193e-18
CM60.10099792968574982.569791e-334.511244e-20
CAG 1 50.10099792968574985.620211e-452.113407e-17
QAG 1 50.10099792968574984.545278e-221.496790e-17
CAG 2 40.10099792968574983.875494e-412.882365e-12
QAG 2 40.10099792968574984.103556e-222.837594e-12
f 5 ( x ) = x 3 + 4 x 2 10 , x 0 = 1
NM 51.20764782713091893.662513e-212.126976e-11
CM411.20764782713091895.538944e-204.261204e-11
CAG 1 41.36523001341409683.000665e-381.557777e-13
QAG 1 41.36523001341409686.937438e-209.377190e-13
CAG 2 31.36523001341409682.672414e-271.361253e-07
QAG 2 31.36523001341409686.648073e-191.349188e-07
f 6 ( x ) = x 2 e x 3 x + 2 , x 0 = 2
NM 50.25753028543986083.439576e-279.869210e-14
CM50.25753028543986089.348485e-203.638191e-10
CAG 1 40.25753028543986081.218871e-415.694586e-14
QAG 1 40.25753028543986089.378610e-254.284785e-14
CAG 2 30.25753028543986089.720474e-292.979348e-07
QAG 2 30.25753028543986081.658292e-232.973073e-07
f 7 ( x ) = x 2 ( 1 x ) 5 , x 0 = 0.9
NM 50.34595481584824206.819645e-176.158927e-09
CM60.34595481584824204.839061e-221.160085e-11
CAG 1 40.34595481584824209.410608e-421.327460e-14
QAG 1 30.34595481584824207.943736e-171.736512e-05
CAG 2 30.34595481584824201.606415e-382.055302e-10
QAG 2 30.34595481584824201.855349e-231.990689e-10
f 8 ( x ) = x 3 9 x 2 + 24 x 20 1 / 3 + e x / 2 , x 0 = 2
NM Fail
CM Fail
CAG 1 Fail
QAG 1 40.96942644858323262.465462e-157.651474e-06
CAG 2 Fail
QAG 2 40.96942644858323147.378707e-251.231459e-10
f 9 ( x ) = x 3 3 x 2 3 + log ( x + 1 ) , x 0 = 3
NM Fail
CM Fail
CAG 1 Fail
QAG 1 62.69251767626217181.218786e-166.562094e-10
CAG 2 Fail
QAG 2 52.69251767626217173.027817e-174.453735e-07
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sana, G.; Mohammed, P.O.; Shin, D.Y.; Noor, M.A.; Oudat, M.S. On Iterative Methods for Solving Nonlinear Equations in Quantum Calculus. Fractal Fract. 2021, 5, 60. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5030060

AMA Style

Sana G, Mohammed PO, Shin DY, Noor MA, Oudat MS. On Iterative Methods for Solving Nonlinear Equations in Quantum Calculus. Fractal and Fractional. 2021; 5(3):60. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5030060

Chicago/Turabian Style

Sana, Gul, Pshtiwan Othman Mohammed, Dong Yun Shin, Muhmmad Aslam Noor, and Mohammad Salem Oudat. 2021. "On Iterative Methods for Solving Nonlinear Equations in Quantum Calculus" Fractal and Fractional 5, no. 3: 60. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5030060

Article Metrics

Back to TopTop