Non-Lipschitz heterogeneous reaction with a p-Laplacian operator

1.   Introduction

It is widely acknowledged that the study of the nonlinear equations characterized as $f(\tilde u) = 0$ can be used to investigate a broader range of problems that occur in physical sciences. Because of its significance, many scientists have investigated numerous different order multistep methods to explore the solutions of the nonlinear equations using diversified approaches; such as variational iterative methods, homotopy perturbation method, homotopy analysis method, and the decomposition techniques, for details see ^{[1,2,3,4,5,6,7,8,9]}. These developed approaches are of varying order of convergence. Firstly, Traub ^[16] has initiated the study of repetitious schemes for solving nonlinear equations who developed a central quadratic convergent Newton iterative method which has much importance in literature.

Later on, to increase the practical usefulness and efficiency index of Newton's method, its various rectifications have been presented by many researchers (see ^{[17,18,19,20]}). Daftardar-Gejji and Jafari ^[21] have proposed a straightforward approach that does not necessitate the derivative evaluation of the Adomian polynomial by making different modifications in the Adomian decomposition method ^[1]. Moreover, this technique helps us to write the nonlinear equation as a combination of both linear and nonlinear components, and plays a remarkable role in developing different iterative schemes to estimate the solution of the nonlinear equations. Saqib and Iqbal ^[22] have determined the fourth and fifth-order convergent iterative methods for computing roots of the nonlinear equations by using a modified decomposition approach and presented some test examples to check the efficacy and performance of the newly established iterative methods. Ali et al. ^[23] have established a new class of the iterative methods by implementing the technique ^[21] and testified the validity of these schemes by considering some mathematical models. Alharbi et al. ^[2] have introduced some new and efficient iterative methods and implemented a decomposition technique along with an auxiliary function. Variational iteration method (VIM) is another effective tool that is employed to develop effective iterative methods for getting approximate, converging solutions of the nonlinear equations. Based on VIM, Naseem et al. ^[17] have investigated a new class of iterative methods that are superior in convergence and efficient as compared to other methods. They also elaborated the behavior and dynamical aspects of the suggested iterative schemes by using polynomiographs. The $q-$calculus is materialized as the composition of Physics and Mathematics in the last twenty-five years of the XX century (see ^{[24,25,26,27,43]}). Many researchers have been designated considerable thought due to its diversified choice of the utilization in mathematical spheres; such as mechanics, theory of relativity, basic hypergeometric function, quantum, and number theory. Firstly, the $q-$analogue of derivative and the $q-$Taylor's formula were introduced by Jackson ^[29]. Then, by using the differentiation technique in the $q-$calculus, Jing and Fan ^[30] have presented some modifications in Taylor's formula together with its remainder in the sphere of the $q-$calculus. They compared the $q$-Taylor's formula and the ordinary Taylor's formula and found signified results on the $q-$remainder.

Ernst ^[31] has investigated some novel modifications in the $q-$Taylor's formula with the help of $q-$ integration by parts and developed its varient formulations. Firstly, some novel recursive schemes under the $q-$analysis were suggested and analyzed by Singh et al. ^[32] and introduced some varied forms of the $q-$iterative schemes by opting several values of the $q-$parameter. The stability and reliability of the $q-$iterative methods are checked by presenting comparative analysis of several nonlinear algebraic equations with some classical methods. The $q-$difference equation plays a vital role in the realm of the $q-$calculus. For solving the partial differential equations in the $q-$calculus, Jafari et al. ^[33] have applied an iterative method called the Daftardar-Jafari decomposition method. It is demonstrated that the proposed procedure's computational outcome converges to the true solution of the $q-$difference equations subject to specific constraints. The study of $q-$integro differential equation with three criteria was investigated by Abdeljawad and Samei ^[12] and checked its solution existence by applying the $q-$calculus. Sadik and Orie ^[11] have introduced a convenient and efficient method based on $q-$calculus known as $q$-differential transform method for solving partial $q$-differential equations. The solution obtained by this method is expressed in terms of convergent power series and the validity of this method is checked by computing several examples. Liang and Samei ^[13] have determined the existence of solutions for non-linear problems regular and singular fractional $q-$differential equation subject to certain constraints. They have presented some results with the support of numerical examples and by applying definitions of the fractional $q$-derivative of Riemann–Liouville & Caputo type. Many real-life problems can be modeled in the form of $q-$fuzzy differential equations. Noeiaghdam et al. ^[15] have introduced two fuzzy numerical methods based on the generalized Hukuhara $q-$differentiability named as the fuzzy $q-$Euler's and the local $q-$Taylor's expansion method for solving $q-$fuzzy initial value differential equations. In an attempt of transformation of the classical results towards the $q-$calculus considered by Srivastava et al. ^[45], the two subclasses of normalized analytic functions are investigated by using various operators of $q-$calculus and fractional $q-$calculus in the complex $z-$plane. Sana et al. ^[10] have transformed the classical iterative methods over the $q-$iterative methods and presented a comparative analysis of these methods with the classical methods. They also presented the generalized formulation of new methods and test their reliability, effectiveness & convergence speed via various numerical examples.

Motivated and inspired by the research going on in this direction, we have restructured some new multistep iterative methods for computing zeros of the nonlinear equations in the context of the $q-$calculus. First, we find some new $q-$analogues of the iterative methods initiated and advanced by Shah and Noor ^[9]. Then, to obtain the needed results, we rephrase the supposed nonlinear equation accompanied by an auxiliary function and apply the $q-$Taylor's formula. For the best implementation of the results and the derivations of recursive schemes, we utilize the decomposition methodology ^[21] under the $q-$paradigm. It is essential to mention that the new suggested algorithms can reduce the number of computing costs compared to conventional iterative methods while good numerical accuracy is maintained by appropriately choosing the parameter $q\in [0, 1]$. Now, we recollect some of the basic ideas in the $q-$analysis ^[34] that are prerequisites and reinforce the construction of our novel $q-$iterative schemes for computing solutions of the nonlinear equations. Let the $q-$integer, for $q$ $\in$(0, 1) is described such as:

For $0\leq p \leq m$, the $q-$factorial and the $q-$binomials are defined as:

Definition 1 (see ^[34]). A $q-$analogue of classical exponential function $e_q^{\tilde u}$ is defined as

The derivative of the classical exponential function remains unchanged under differentiation. The $q-$analogue of exponential function also remains the same in the $q-$calculus such as:

Definition 2 (see ^[34]). Let $f(\tilde u)$ is a real valued continuous function and its $q-$derivative is prescribed as follows:

where $(D_qf)(\tilde u)$ represents $q-$derivative is known as Jackson derivative. It reduces to the standard derivative when $q$ approaches to one. The $q-$derivative with higher-order for the function $f(\tilde u)$ is prescribed as:

Definition 3 (see ^[34]). The $q-$derivative of product and quotient of function $f(\tilde u)$ and $g(\tilde u)$ is defined as follows:

Definition 4 (see ^[29,30,31]). Let $f(\tilde u)$ is a continuous function defined on an interval $(k, l)$ and $c\in$ $[k, l]$. Then, the $q-$Taylor's formula for the function $f(\tilde u)$ instigated by Jackson is explained as:

where

$D_q, D^2_q, \ldots$ are all $q-$derivatives, where $0 < q < 1$.

2.   Evolvement of the iterative schemes

This section deals with the construction of some novel iterative schemes by employing Taylor's formula and Daftardar-Jafari decomposition technique ^[21] in the paradigm of the quantum calculus.

2.1.   Main results

We consider the nonlinear algebraic equation of the general form:

Let $g(\tilde u)$ be an auxiliary function. Suppose $\kappa$ is an initial guess in the neighbourhood of $\beta$ which is the simple root of nonlinear equation (2.1).

Using $q-$Taylor's series about $\kappa$ and the technique of He ^[28], we rewrite the nonlinear equation (2.2) as a parallel coupled system of the equation :

Eq (2.4) can be rewritten as:

where

and

The term $N_q(\tilde u)$ is treated as nonlinear and $c$ as a constant.

Let $\tilde u_0$ be an initial guess then from relation (2.4), we can easily compute a key equation that is helpful in the development of new $q-$iterative methods :

We now carry out a decomposition technique primarily due to Daftardar-Gejji and Jafari ^[21], known as the Daftardar-Jafari decomposition technique, to set up arrangements of higher-order iterative methods. The central idea behind using this methodology is to seek the solution of the $q-$basic equation (2.6) in the series form.

Now, we deteriorate the nonlinear operator $N_q(\tilde u)$ which is defined in (2.8) such as:

From the equations (2.6), (2.10) and (2.11), we have

finally, we obtain the following iterative procedure:

It follows that

and

Note that $\tilde u$ is approximated by

and thus $\lim\limits_{x\to\infty}\tilde u_n = \tilde u$.

Theorem 2.1 (see ^[33]). If $N_q$ is a contraction, then the series specified in $(2.10)$ is absolutely convergent.

Proof. Let $N_q$ is a contraction mapping, then by definition we can write:

then in view of (2.13), we have

then the series $\tilde u$ = $\sum_{k = 0}^{\infty}\tilde u_k$ is uniformly and absolutely convergent to an answer of the equation (2.6) (see ^[36]).

This completes the proof.

Now, we construct the following iterative schemes to find the solution of the nonlinear algebraic equation (2.1)

$\bf Algorithm \ A:$ From (2.13), we have for $n = 0:$

This composition permits us to put forward the subsequent recursive approach for solving the nonlinear equation (2.1), and the iterative schema computes the approximate solution $\tilde u_{ n+1}$ for a given starting guess $\tilde u_0$:

This represents main $q-$analogue of iterative scheme which is prospected by He ^[28] and Shah ^[9]. This main iterative scheme helps us to generate different $q-$algorithms for solving the nonlinear equation (2.1). Now, with the help of (2.4) and (2.13), we get:

$\bf Algorithm \ B:$ From (2.13), we have for $n = 1:$

By using (2.16) and (2.18), we have

This composition permits us to put forward the subsequent iterative approach for solving the nonlinear equation (2.1), and the iterative schema computes the approximate solution $\tilde u_{ n+1}$ for a given starting guess $\tilde u_0$:

This is $q-$analogue of Algorithm 2.2 which is investigated by ^[9] and the error equation of Algorithm B is determined in Theorem 3.1.

By using (2.4), (2.13) and (2.18), we can obtain

$\bf Algorithm \ C:$ Now, considering from (2.13), we have for $n = 2:$

By using (2.16) and (2.23), we get

This composition permits us to put forward the subsequent iterative approach for solving the nonlinear equation (2.1), and the iterative schema computes the approximate solution $\tilde u_{ n+1}$ for a given starting guess $\tilde u_0$:

This is $q-$analogue of Algorithm 2.3 which is investigated by ^[9] and the error equation of Algorithm C is determined in Theorem 3.1.

Algorithm A, Algorithm B, Algorithm C are the main and general iterative schemes that are used to generate some new algorithms by considering different choices of auxiliary functions that are the main attractiveness of modification of this technique. To convey the idea, we consider the following auxiliary functions.

$\bf Case :$ Let $g(\tilde u_n) = e^{-\beta \tilde u_n}$ and $D_qg(\tilde u_n) = -\beta e^{-\beta \tilde u_n}$. Using these values we obtain the following iterative methods for the solving nonlinear equations.

$\bf Algorithm \ D:$ For a given $\tilde u _0$ (initial guess), approximate solution $\tilde u _{n+1}$ is computed by the following iterative scheme:

$\bf Algorithm \ E:$ For a given $\tilde u _0$ (initial guess), approximate solution $\tilde u _{n+1}$ is computed by the following iterative scheme:

$\bf Algorithm \ F:$ For a given $\tilde u _0$ (initial guess), approximate solution $\tilde u _{n+1}$ is computed by the following iterative scheme:

To the best of our knowledge, the new schemes Algorithm D, Algorithm E, Algorithm F appear to be new ones.

3.   Convergence analysis

In this part, the order of convergence of the primary $q-$iterative methods made out by $\textbf{Algorithm A}$, $\textbf{Algorithm B}$, and $\textbf{ Algorithm C}$ is investigated. In the same approach, the rest of the iterative procedures can be established.

Theorem 3.1. Let $f:E\subset R\rightarrow$ R be a differentiable function, where $E$ is an open interval in R. If $\tilde u _0$ is sufficiently close to $\beta\in E$ which is the root of $f(\tilde u) = 0$ then the iterative methods Algorithm A, Algorithm B and Algorithm C are convergent algorithms of order at least $2, 3, 4$ respectively and we format it as follows: $[2;q]$, $[3;q]$ and $[4;q]$, where parameter $q$ corresponds to the quantum calculus. Error equations for these newly established algorithms are given as:

Proof. Let $f$ is adequately differentiable function and $\beta$ is root of $f(\tilde u)$. Now, expanding $f(\tilde u_n)$ and $D_qf(\tilde u_n)$ in the $q-$Taylor's series about $\beta$ we obtain

By expanding $g(\tilde u_n)$, $D_qg(\tilde u_n)$ in the $q-$Taylor's series, we obtain

where

By expanding $f(\tilde u_n)g(\tilde u_n)$, $f(\tilde u_n)D_qg(\tilde u_n)$, $D_qf(\tilde u_n)g(\tilde u_n)$, in the $q-$Taylors series about $\beta$, we obtain

From (3.7), (3.9), (3.8), we get

Now, using (3.10) into (2.17), we get the error term of the Algorithm A:

Choosing (3.12), we have

By expanding $f(\tilde v_n)$ in the $q-$Taylor's series about $\beta$ and using (3.14), we have

From (3.4), (3.9), (3.8) and (3.15) we have

Using (3.14), (3.16) into (2.22), we obtain the error term for the Algorithm B

By expanding $\tilde w_n$, $f(\tilde w_n)$ in terms of the $q-$Taylor's series about $\beta$

From (3.4), (3.9), (3.8) and (3.19) we have

Using (3.18) and (3.20) into (3.12), we obtain the error equation for the Algorithm C:

Equation (3.20) shows the error equation for the Algorithm C and has at least fourth-order convergence. It is noted that Algorithm C is the main iterative scheme and all other schemes investigated from this scheme are at least fourth-order convergent.

Remark 3.1. Based on the study of convergence analysis of proposed iterative methods, it can be easily observed that various order iterative methods can be developed by choosing appropriately multiple choices of the auxiliary function in Algorithm A, Algorithm B and, Algorithm C respectively.

then Algorithm A generates the following iterative method with the initial guess $\tilde u_0$.

$\bf Algorithm \ G:$ For a given initial guess $\tilde u_o$, approximate solution $\tilde u_{n+1}$ is computed by the following iterative scheme:

This is $q-$analogue of well known Halley method ^[32] which has cubic convergence i.e. [3, $q$], where $q$ represents the $q-$calculus. Now, again using the above stated specified value of an auxiliary function then Algorithm B and Algorithm C reduces to the following iterative procedures.

$\bf Algorithm \ H:$ For a given initial guess $\tilde u_o$, approximate solution $\tilde u_{n+1}$ is computed by the following iterative scheme:

This method is fourth-order convergent for solving nonlinear equations and appears to be a novel one.

$\bf Algorithm \ I:$ For a given initial guess $\tilde u_o$, approximate solution $\tilde u_{n+1}$ is computed by the following iterative scheme:

This method emerges as a new method that has fifth-order of convergence.

This completes the proof.

4.   Numerical examples and applications

This section discusses some nonlinear equations. With the support of these examples, we elaborate on the efficacy and performance of newly established methods initiated in this paper. The general algorithm for finding the estimated solution of the given nonlinear function is given as: in Algorithm A, Algorithm B, Algorithm C, we consider $\varepsilon = 10^{-100}$ as tolerance. We obtain an approximate solution relatively than the exact lean on the computational accuracy $\varepsilon$. We adopt the following stopping criterium for computational performance:

For convergence criteria, it is prerequisite that the space of two successive estimations for the zero must be less than $10^{-100}$. We make use of abbreviations $\mathsf {QG}$ & $\mathsf {CG}$ for the $q-$iterative methods and traditional iterative methods respectively. We symbolize Algorithm D, Algorithm E and Algorithm F by $\mathsf{QG}1$, $\mathsf{QG}2$ and $\mathsf{QG}3$ respectively and phrase $\it{div}$ served as the divergence of methods. We develop a comparative analysis between the standard Newton's method ($\mathsf{NM}$) ^[35], Chun method ($\mathsf{CM}$) ^[4], Noor method (NR)^[8], $\mathsf{CG}1$, $\mathsf{CG}2$ and $\mathsf{CG}3$ ^[9] and our newly proposed $q-$iterative methods $\mathsf{QG}1$, $\mathsf{QG}2$ and $\mathsf{QG}3$. The computational results of comparative analysis are presented in Tables (4, 8, 12, 14). We exhibit the number of iterations, the final estimated solution and the corresponding functional value by the symbols $\mathsf{IT}$, $\tilde u _n$ and $f(\tilde u_n)$, whereas, the distance in the middle of two successive estimates is shown by $\Delta$. It is necessary to mention that for the best implementation of results, we choose the value of $q = 0.9999$. We use Maple software to perform all the numerical computations.

Now, we recollect the classical Algorithm 2.1 ($\mathsf{CG}1$) in ^[9], elucidated as:

and the classical Algorithm 2.2 ($\mathsf{CG}2$) in ^[9], described as:

and the classical Algorithm 2.3 ($\mathsf{CG}3$) in ^[9], described as:

We present some examples of nonlinear equations (4.1–4.4) to illustrate the efficiency of the newly developed one-step, two-step and three-step iterative methods in this article. Firstly, for the sake of simplicity, we investigate the efficacy and credibility of the $q-$recursive schemes for multiple values of $q$ up to three iterations that can be extended to any number of iterations until we achieve the desired accuracy. The results in the Tables (1–3, 5–7, 9–11) demonstrate the calculations of $\tilde u_{i}$ and $f(\tilde u_{i})$, $i = 1, 2, 3$ by employing ${\mathsf{QG}}1$, ${\mathsf{QG}}2$, ${\mathsf{QG}}3$ for multiple values of $q$ and $\beta = 0.5$. We choose $\beta = 0.5$ for both $q$ and ordinary iterative methods.

Example 4.1 (see ^[4]). We assume the folllowing nonlinear equation$:$

We take $\tilde u_0 = -2$ as an initial guess for computational evaluations. The quantifiable outcomes for the equation (4.2) are calculated in Tables (1–3) by using $\mathsf{QG}1$, $\mathsf{QG}2$, $\mathsf{QG}3$ for multiple values of $q$ and $\beta = 0.5$. Following the steps of the Tables (1–3), we get the required solution of equation (4.2) i.e. $\tilde u = 1.2076478271$.

The results from the Table 1, elaborate that precise values of $\tilde u_i$s are achieved subject to the constraint $q$ approaces to one and the parallel functional values $f(\tilde u_i)$ tend to zero, where $i$ = 1, 2, 3. It is also noted that the values of $f(\tilde u_1) = 4.365516e+01, f(\tilde u_2) = 1.736630e+01, f(\tilde u_3) = 6.213934e+00$ computed by $\mathsf{QG}1$ at $q = 0.9999$ are nearer to zero in comparison with the values $f(\tilde u_1) = 4.367558e+01, f(\tilde u_2) = 1.738244e+01, f(\tilde u_3) = 6.223602e+00$ calculated by $\mathsf{CG}1$. And by choosing $q = 0.9999$ and $\beta$ = 0.5 the equation (4.2) converges to the root $\tilde u_8 = 1.2076478271$ and corresponding functional value is attained as $f(\tilde u_8) = 1.268198e-10$.

Columns in Table 2 display the more precise values of $\tilde u_1, \tilde u_2, \tilde u_3$ with the constraint that $q$ approaches one and for which $f(\tilde u_1), f(\tilde u_2), f(\tilde u_3)$ tend to zero. Furthermore, it is also noted that the values of $f(\tilde u_1) = 3.032134e+01, f(\tilde u_2) = 7.638732e+00, f(\tilde u_3) = 1.084323e+00$ at $q = 0.9999$ calculated by $\mathsf{QG}2$ exist closely in the neighbourhood of zero in contrast to the values of $f(\tilde u_1) = 3.033637e+01, f(\tilde u_2) = 7.646987e+00, f(\tilde u_3) = 1.087173e+00$ calculated by $\mathsf{CG}2$. Following the steps of Table 2 and by taking values $q = 0.9999$, $\beta = 0.5$, the equation (4.2) converges to the root $\tilde u_6 = 1.2076478271$ and corresponding function value is attained as $f(\tilde u_6) = 1.389449e-15$.

Table 3 illustrates the accuracy and precision of results for $\tilde u_i$, $i = 1, 2, 3$ whenever $q$ $\rightarrow$ one. One can also figure out that for $q = 0.9999$ computed values of $f(\tilde u_1) = 2.344843e+01$, $f(\tilde u_2) = 3.928572e+00$, $f(\tilde u_3) = 1.242253e-01$ by $\mathsf{QG}3$ give results nearer to zero in comparison to $f(\tilde u_1) = 2.345843e+01, f(\tilde u_2) = 3.933608e+00, f(\tilde u_3) = 1.248415e-01$ calculated by $\mathsf{CG}3$. Following the steps of Table 3 the equation (4.2) converges to the root $\tilde u_5 = 1.2076478271$ and $f(\tilde u_5) = 3.785784e-18$, for $q = 0.9999$ & $\beta$ = 0.5.

In Table 4, we compare our new $q-$iterative methods ($\mathsf{QG}1$, $\mathsf{QG}2$, $\mathsf{QG}3$) with some other methods to examine the reliability and effectiveness of the methods. The second column ($\mathsf{IT}$) in Table 4 exhibits the comparison of different iterative methods with newly established methods concerning to the number of iterations. It is clear from the computational results that proposed methods require less number of iterations compared to other methods to meet the stopping criteria (4.1) or number of iterations are the same in some cases when comparing with ($\mathsf{CG}1$, $\mathsf{CG}3$).

Example 4.2 (Van der Waal's Equation see ^[37]).

We consider the Van der Waal's equation representing the real and ideal behaviour of gass is prescribed as:

Eq (4.3) can be transformed to following nonlinear form:

After appropriately choosing the needed parameters and unknown constants we can find the following nonlinear function:

where the variable $\tilde u$ shows the volume of the gas. We take $\tilde u_0 = 3.10$ as an initial guess for computational evaluations. The mathematical computations for the equation (4.5) are calculated in Tables 5–7, for multiple values of $q$ and $\beta$ = 0.5. Following the steps of the Tables 5–7, we get the required solution for equation (4.5) i.e. $\tilde u = 1.9298462428.$

Continuing step by step and evaluating values, Table 5 illustrates the accuracy and precision of results for $\tilde u_i$, $i = 1, 2, 3$ whenever $q$ $\rightarrow$ one. One can also figure out that for $q = 0.9999$ computed values of $f(\tilde u_1)$ = 4.436290e-01, $f(\tilde u_2) = 1.030037e-01$, $f(\tilde u_3) = 2.540698e-02$ by adopting $\mathsf{QG}1$ give results more near to zero in comparison to $f(\tilde u_1) = 4.437270e-01$, $f(\tilde u_2) = 1.030704e-01$, $f(\tilde u_3) = 2.544189e-02$ calculated by $\mathsf{CG}1$. And for $q = 0.9999$ & $\beta$ = 0.5, equation (4.5) converges to the root $\tilde u_9 = 1.929846242847$ and $f(\tilde u_9) = 9.996931e-14$.

Table 6 investigate the precision and accuracy of the values of $\tilde u_1, \tilde u_2, \tilde u_3$ when the parameter $q$ $\rightarrow$ one parallel functional values $f(\tilde u_1), f(\tilde u_2), f(\tilde u_3)$ tend to zero. Moreover, one can also figure out that for $q = 0.9999$ enumerated values of $f(\tilde u_1) = 2.814977e-01, f(\tilde u_2) = 4.108624e-02, f(\tilde u_3) = 5.485842e-03$ at $q = 0.9999$ by adopting $\mathsf{QG}2$, are more near to zero in comparison with $f(\tilde u_1) = 2.815532e-01, f(\tilde u_2) = 4.111438e-02, f(\tilde u_3) = 5.496367e-03$ computed by $\mathsf{CG}2$. Also, following the steps of Table 6 and for $q = 0.9999$ the equation (4.5) converges towards the root $\tilde u_7 = 1.929846242847$ and $f(\tilde u_7) = 3.089402e-19$.

We calculate the values $\tilde u_i, f(\tilde u_i)$, where i = 1, 2, 3, in the Table 7 that illustrates the accuracy of values $\tilde u_i$ subject to constraint $q$ approaches one. Moreover, it is apparently observed that the values of $f(\tilde u_1) = 2.021924e-01, f(\tilde u_2) = 2.062325e-02, f(\tilde u_3) = 1.336893e-03$ at $q = 0.9999$ computed by $\mathsf{QG}3$, are near to zero in comparison with $f(\tilde u_1) = 2.022301e-01, f(\tilde u_2) = 2.063865e-02, f(\tilde u_3) = 1.340965e-03$ determined by $\mathsf{CG}3$. Also, the equation (4.5) converges to root $\tilde u_6 = 1.929846242847$ for $q = 0.9999$ and corresponding functional value is attained as $f(\tilde u_6) = 2.935505e-21$.

Table 8 presents the comparative analysis of our new $q-$iterative methods ($\mathsf{QG}1$, $\mathsf{QG}2$, $\mathsf{QG}3$) with some other classical methods. In Table 8 the efficiency of methods is displayed with respect to number of iterations using the stopping criterium (4.1). It is observed that new methods are comparable with other methods and need less number of iterations required by the other methods of the same order.

Example 4.3 (Motion of particle on an inclined plane see ^[35]). We consider the nonlinear model created due to the motion of a particle on an inclined plane whose inclination angle $\theta$ remodeled at a steady rate $\frac{d(\theta)}{dt} = \gamma < 0$.

We take $\tilde u_o = -1.2$ as an initial guess for computational evaluation. The numerical findings for equation (4.6) are calculated in the Tables (8–10) by using $\mathsf{QG}1$, $\mathsf{QG}2$, $\mathsf{QG}3$ for multiple values of $q$ and $\beta = 0.5$.

The numeric values in Table 9, illustrate that one can obtain more precise values of $\tilde u_i$ with the constraint $q$ $\rightarrow$ one and $f(\tilde u_i)$ attain zero value, where $1\leq i \leq 3$. It is also noted that the values of $f(\tilde u_1) = 1.562110e-01, f(\tilde u_2) = 5.384093e-02, f(\tilde u_3) = 1.734706e-02$ computed by $\mathsf{QG}1$ at $q = 0.9999$ exist in the neighbourhood of zero more nearly in comparison to the values $f(\tilde u_1) = 1.562383e-01, f(\tilde u_2) = 5.386087e-02, f(\tilde u_3) = 1.735750e-02$ calculated by $\mathsf{CG}1$. Following the steps of the Table 9, the equation (4.6) converges towards the root $\tilde u_9 = 0.3170617745$ and $f(\tilde u_9) = 8.400168e-13$, for $q = 0.9999$.

The results in Table 10, elaborate that one can obtain more precise values of $\tilde u_i$ with the constraint $q$ $\rightarrow$ one and $f(\tilde u_i)$ attain zero value, where $1\leq i \leq 3$. By choosing $q$ = 0.9999, it is also noted that the values of $f(\tilde u_1) = 9.698942e-02, f(\tilde u_2) = 1.973729e-02, f(\tilde u_3) = 3.342707e-03$ computed by $\mathsf{QG}2$ are nearer to zero in comparison to the values $f(\tilde u_1) = 9.700799e-02, f(\tilde u_2) = 1.974547e-02, f(\tilde u_3) = 3.345242e-03$ calculated by $\mathsf{CG}2$. Continuing the iterative procedure as presented in Table 10, the equation (4.6) converges to the root $\tilde u_6 = 0.3170617746$ for $q = 0.9999$ and $f(\tilde u_6) = 4.953192e-14$.

Columns in the Table 11 demonstrate that one can obtain more precise values of $\tilde u_i$ with the constraint that $q$ $\rightarrow$ 1 and $f(\tilde u_i)$ attain value zero, where $1\leq i \leq 3$. It is also noted that the values of $f(\tilde u_1) = 6.887648e-{02}, f(\tilde u_2) = 9.422356e-03, f(\tilde u_3) = 8.107219e-04$ computed by $\mathsf{QG}3$ at $q = 0.9999$ are more adjacent to zero in comparison to the values $f(\tilde u_1) = 6.889025e-02, f(\tilde u_2) = 9.426571e-03, f(\tilde u_3) = 8.115710e-04$ calculated by $\mathsf{CG}3$. Following the steps of Table 11 and for $q = 0.9999$, $\beta$ = 0.5, the equation (4.6) converges to the root $\tilde u_5 = 0.3170617746$ and parallel functional values are obtained as$f(\tilde u_5) = 4.850310e-{15}$.

The second column ($\mathsf{IT}$) in Table 12 illustrates the comparsion of different iterative methods with proposed methods in terms of number of iterations. It is clear from the computational results that new methods need less number of iterations as compared to other methods ($\mathsf{NM}$, $\mathsf{CM}$, $\mathsf{NR}$, $\mathsf{CG}1$) to meet the stopping criteria (4.1) or same in some cases when comparing with ($\mathsf{CG}2$, $\mathsf{CG}3$).

Remark 4.1. It is worthy to mention that when we evaluate the errors for the $q-$iterative schemes then it oscillate for various values of $q$. The error reduces when $q$ tends to the highest values uniting $0$ and $1$. Therefore, in Table 13 error for equations [(4.2), (4.5), (4.6)] are computed by using $q$ = 0.9999 and $\beta$ = 0.5 which will estimate the classical methods.

Example 4.4. (Algebraic and Transedental equations). This example comprises of a few nonlinear equations which help us to examine the reliability and effectiveness of our new $q-$iterative methods.

Some of these nonlinear equations are used by Chun^[4] & Singh et al.^[32] to validate the theoretical results. The last three numerical equations namely; $f_7(\tilde u)$, $f_{8}(\tilde u)$, $f_{9}(\tilde u)$ represent some real-world applications of nonlinear equations. These nonlinear equations are the transformations of some mathematical models that appeared in science and engineering. The first one nonlinear equation $f_7$ is generated as a solution of mathematical modeling of the growth of population over short periods of time that can be written as in the form of differential equation:

where $\lambda$ denotes the constant birth rate of population and $M(\tilde t)$ denotes the number in the population at time $\tilde t$, for details (see ^[35]). The second nonlinear equation $f_{8}$ represents the physical constraint problem of fractional modification in a chemical reactor. The variable $\tilde u$ illustrates a fractional conversion of certain kind in a chemical reactor problem (see ^[39]). The value of $\tilde u$ is chosen between [0, 1] because for negative values of $\tilde u$ the equation $f_{8}$ has no physical meaning. Therefore, within the limited reigion we have to select the initial guess carefully to find the real root of $f_{8}$. The third nonlinear equation $f_9$ is originated from the problem of the fraction conversion of nitrogen-hydrogen to ammonia which was investigated by ^[38]. The problem has the following form:

The values of temperature and pressure have been considered as 500℃ and 250 atm respectively which can be easily reduced to the equation $f_{9}$.

The Table 14 reflects the comparable outcomes of classical and $q-$iterative methods by implementing the stoping criterium (4.1). From Table 14, it is also noted that we obtain identical results for both $q-$analogue of iterative methods ($\mathsf{QG}1$, $\mathsf{QG}2$ and $\mathsf{QG}3$) and conventional classical iterative methods ($\mathsf{CG}1$, $\mathsf{CG}2$ and $\mathsf{CG}3$). We consider two functions that are non-differentiable at points $\tilde u$ = 2, 3. When, we choose $\tilde u = 2, 3$ as an initial guesses for $f_7(\tilde u)$ & $f_8(\tilde u)$ consectively then novel schemes $\mathsf{QG}1$, $\mathsf{QG}2$ and $\mathsf{QG}3$ are appropriately implemented and offer quick concurrent solutions whereas, the classical methods crash for these functions which is a major advantage of using $q-$iterative methods over the classical methods.

Remark 4.2 (see ^[16]). The efficiency index is calculated as $\rho^\frac{1}{k}$ where $k$ determines the required number of estimations per iteration imperatively applied to a step of a recursive method and $\rho$ symbolizes the convergence order of the method.

● Efficiency index of $\mathsf{CG}1$ is $2^\frac{1}{4} = 1.89207$.

● Efficiency index of $\mathsf{CG}2$ is $3^\frac{1}{5} = 1.245731$.

● Efficiency index of $\mathsf{CG}3$ is $4^\frac{1}{6} = 1.259921$.

● Efficiency index of $\mathsf{QG}1$ is $2^\frac{1}{4} = 1.89207$.

● Efficiency index of $\mathsf{QG}2$ is $3^\frac{1}{5} = 1.245731$.

● Efficiency index of $\mathsf{QG}3$ is $4^\frac{1}{6} = 1.259921$.

Finally, we come to an end that the efficiency indexes evaluated by $\mathsf{QG}1$, $\mathsf{QG}2$, $\mathsf{QG}3$ and $\mathsf{CG}1$, $\mathsf{CG}2$, $\mathsf{CG}3$ give identical outcomes.

General formulation of the $q$-iterative schemes

This section consists of some previous results which are used for the derivation of the generality of the $q-$iterative methods. Now, combining the entries $u_k s'$ $(\forall\; k = 1, 2, 3, \ldots, n)$ in (2.13), we acquire

From (2.16), we have

and

Now, if $\tilde u$ is approximated by

Therefore, (4.9) gives the following iterative scheme

This relation enables and allows us to propose the subsequent iterative method.

This is convergent generalized $q-$iterative scheme which has order of convergence $n+2$ for $n = 0, 1, 2, \ldots$ proof is simple and starightforward.

Remark 4.3. System of nonlinear equations emerges when several scientific and technological challenges are involved. Many integral equations, boundary value problems, minimization problems, and variational problems may also be reduced to the system of nonlinear equations (see ^[40,41,44]). We consider the system of nonlinear equations of the form:

where each function $f_j$, $j = 1, 2, \dots, n$ maps a vector $\tilde U = (\tilde u_1, \tilde u_2, ..., \tilde u_n)$ of the $n$-dimensional space $R^{n}$ to the real line $R$. The above system of n nonlinear equations in $n$-unknowns can also be represented as:

Where $F:R^n\rightarrow R^n$ be nonlinear mapping from $n$-dimensional real linear space $R^n$ into itself. The components $f_j, j = 1, 2, ..., n$, are the coordinate functions of $F$. The solution of the nonlinear system of equations in (4.12) may be defined as the process of finding a vector $U^* = (u_1^*, u_2^* \dots, u_n^*)^t$ such that $F(U^*) = 0$. We feel it worthwhile to mention that considering the methodology and idea of this article, one can present and analyze higher-order multistep iterative methods for solving a system of nonlinear equations (4.12). It is an open problem to broaden the concept and ideas of this study for solving the boundary value problems and associated issues. This is another recommendation for prospective research.

5.   Conclusions and observations

Study and formulation of numerical results in quantum calculus induce interest due to the high demand in mathematics and easy implementation. This manuscript introduces some novel iterative schemes to find the estimated solution of nonlinear equations with success in quantum calculus. The key motivation of proposing $q-$iterative schemes is to overcome differentiability and convergence issues while getting solutions of algebraic equations. These new iterative schemes are applicable for different choices of an auxiliary function and derived by considering the valuable Daftardar-Jafari decomposition technique. We develop the comparative analysis of newly proposed methods with the traditional iterative methods to demonstrate the performance and efficiency of $q-$iterative schemes. Moreover, it is shown that the numerical results obtained for both conventional and $q-$iterative methods remain identical. Also, the errors connected with the suggested schemes are relatively marginal by selecting the value of $q$ approaches to one. Hence, it is evident that the transformation of iterative methods in the $q-$calculus framework which we referred to as $q-$analogue of iterative schemes, is better than classical methods, and in limited cases when the parameter $q\rightarrow$ 1, these methods reduces to the classical iterative methods. The significant challenge of dealing with these schemes which necessitate more exploration is that to get results with high accuracy, we must estimate the value of $q$ in $(0, 1).$

Our utilization here of the $q-$calculus in the development of the iterative methods are supposed to promote and motivate major future breakthroughs in Mathematical analysis. It is noticed that in $(p, q)$ analysis the extra parameter $p$ is clearly redundant, Srivastava (see [42,p. 340] and [43,pp. 1511-1512]; see also the related recent works ^[45,46]) revealed that the so-called $(p, q)$ variations of the suggested $q-$results which are obtained by inconsequentially and trivially adding a redundant parameter $p$ as quite simple and insignificant modification of the standard $q-$calculus. Along these lines, while we reinforce and revitalize the $q-$results introduced in this paper, together with potential $q-$extensions of other similar developments in physical and engineering sciences, we do not encourage and support the so-called $(p, q)$-variations of the suggested $q-$results which are obtained by inconsequentially and trivially adding a redundant or superfluous parameter $p$.

Acknowledgments

The authors would like to thank the Rector, COMSATS University Islamabad, Islamabad Pakistan, for providing excellent research and academic environments.

Conflict of interest

The authors agree with the contents of the manuscript and there is no conflict of interest among the authors.