Certain Results for the Twice-Iterated 2D q-Appell Polynomials

In this paper, the class of the twice-iterated 2D q-Appell polynomials is introduced. The generating function, series definition and some relations including the recurrence relations and partial q-difference equations of this polynomial class are established. The determinant expression for the twice-iterated 2D q-Appell polynomials is also derived. Further, certain twice-iterated 2D q-Appell and mixed type special q-polynomials are considered as members of this polynomial class. The determinant expressions and some other properties of these associated members are also obtained. The graphs and surface plots of some twice-iterated 2D q-Appell and mixed type 2D q-Appell polynomials are presented for different values of indices by using Matlab. Moreover, some areas of potential applications of the subject matter of, and the results derived in, this paper are indicated.

The q-analogues of a complex number α and of the factorial function are defined as follows: and where N is the set of positive integers. The q-binomial coefficients [ m s ] q are defined by The q-analogue of the classical derivative D f or d dt f of a function f at a point t ∈ C\{0} is defined by We also note that , where d dt denotes the classical ordinary derivative, (ii) D q a 1 f (t) + a 2 g(t) = a 1 D q f (t) + a 2 D q g(t), (iii) D q ( f g)(t) = f (qt)D q g(t) + g(t)D q f (t) = f (t)D q g(t) + D q f (t)g(qt), (vi) D q f (t) g(t) = g(t)D q f (t) − f (t)D q g(t) g(t)g(qt) = g(qt)D q f (t) − f (qt)D q g(t) g(t)g(qt) .
The two familiar q-analogues of the exponential function e t are given by and The above-defined q-exponential functions e q (t) and E q (t) satisfy the following properties: D q e q (t) = e q (t), D q E q (t) = E q (qt), (8) e q (t)E q (−t) = E q (t)e q (−t) = 1.
The class of Appell polynomials was introduced and characterized completely by Appell [7] in 1880. Further, Throne [8], Sheffer [9] and Varma [10] studied this class of polynomials from different points of view. For some subsequent and recent developments associated with the Appell polynomials, one may refer to the works [11][12][13][14].
In the year 1954, Sharma and Chak [15] introduced a q-analogue of the Appell polynomials and called this sequence of polynomials as q-Harmonics. Later, in the year 1967, Al-Salam [16] introduced the class of the q-Appell polynomials {A m,q (x)} ∞ m=0 and studied some of their properties. Some characterizations of the q-Appell polynomials were presented by Srivastava [17] in the year 1982. These polynomials arise in numerous problems of applied mathematics, theoretical physics, approximation theory and many other branches of the mathematical sciences [7,[18][19][20]. The polynomials A m,q (x) (of degree m) are called q-Appell polynomials, provided that they satisfy the following q-differential equation: Recently, Keleshteri and Mahmudov [21] introduced the 2D q-Appell polynomials (2DqAP) A m,q (x 1 , x 2 ) ∞ m=0 which are defined by means of the following generating function: where We write A m,q := A m,q (0, 0), where A m,q denotes the 2D q-Appell numbers. For x 2 = 0, the 2DqAP A m,q (x 1 , x 2 ) reduce to the q-Appell polynomials A m,q (x) (see, for example, [16,17,22]), that is, where A m,q (x) are defined by and A m,q given by A m,q := A m,q (0) denotes the q-Appell numbers. The explicit form of the 2DqAP A m,q (x 1 , x 2 ) in terms qAP A m,q (x) is given as follows (see [21]): The function A q (t) may be called the determining function for the set A m,q (x 1 , x 2 ). Based on suitable selections for the function A q (t), the following different members belonging to the family of the 2D q-Appell polynomials A m,q (x 1 , x 2 ) can be obtained: [23,24]), that is, and B m,q given by B m,q := B m,q (0, 0) denotes the 2D q-Bernoulli numbers. [23,24]), that is, and E m,q given by E m,q := E m,q (0, 0) denotes the 2D q-Euler numbers.
We recall here that, in a recent paper, Khan and Riyasat [26] introduced the twice-iterated q-Appell polynomials A [2] m,q (x) which are defined by means of the following generating function: In this paper, the class of the twice-iterated 2D q-Appell polynomials is introduced by means of generating functions, recurrence relations, partial q-difference equations, and series and determinant expressions. Further, several results are obtained for the members corresponding to this polynomial family. In Section 2, the twice-iterated 2D q-Appell polynomials are introduced by means of the generating functions and series definition. Also, the recurrence relation and q-difference equations involving the twice-iterated 2D q-Appell polynomials are derived. In Section 3, a determinant expression for the twice-iterated 2D q-Appell polynomials is established. In Section 4, the determinant expressions and some other properties of the members belonging to the family of the twice-iterated 2D q-Appell polynomials are obtained. Section 5 provides several graphical representations and surface plots associated with several members of families of q-polynomials which have investigated in this paper. Finally, in Section 6, we present some concluding remarks and observations.

Twice-Iterated 2D q-Appell Polynomials
In order to introduce the twice-iterated 2D q-Appell polynomials (2I2DqAP), we consider two different sets of the 2D q-Appell polynomialsȦ m,q (x 1 , whereȦ whereÄ The generating function for the 2I2DqAP is asserted by the following result. Theorem 1. The generating function for the twice-iterated 2D q-Appell polynomials A [2] m,q (x 1 , x 2 ) is given bẏ Proof. By expanding the first q-exponential function e q (x 1 t) in the left-hand side of the Equation (20) and then replacing the powers of x, that is, x 0 1 , x 1 , x 2 1 , · · · , x m 1 by the polynomialsÄ 0,q (x 1 ),Ä 1,q (x 1 ), A 2,q (x 1 ), · · · ,Ä m,q (x 1 ) in the left-hand side and x 1 byÄ 1,q (x 1 ) in the right-hand side of the resultant equation, we havė Moreover, by summing up the series in the left-hand side and then using the Equation (24) in the resulting equation, we geṫ Finally, denoting the resulting 2I2DqAP in the right-hand side of the above equation by A [2] m,q (x 1 , x 2 ), that is,Ȧ m,q (Ä 1,q (x 1 ), x 2 ) = A [2] m,q (x 1 , x 2 ), (28) the assertion (25) of Theorem 1 is proved.
m,q (x 1 , x 2 ) reduce to the twice-iterated q-Appell polynomials (see [26]) such that A It is also noted that We next give the series definition for the 2I2DqAP A [2] m,q (x 1 , x 2 ) by proving the following result.
m,q (x 1 , x 2 ) are given by the following series expression: Proof. In view of the Equations (21) and (22), the Equation (25) can be written as follows: which, on using the Cauchy product rule, gives Equating the coefficients of like powers of t in both sides of the above equation, we arrive at the assertion (31) of Theorem 2.
We now state and prove the following result.
Proof. Using the Equations (7) and (19) in the left-hand side of the generating function (25), we get which, on applying the Cauchy product rule in the left-hand side, yields Finally, equating the coefficients of like powers of t on both sides of this last equation, we obtain the assertion (35) of Theorem 3.

Remark 3. By taking x
Remark 4. The following statements are equivalent: In order to derive the q-recurrence relations and the q-difference equations for the twice-iterated 2D q-Appell polynomials by using the lowering operators that are, in fact, the q-derivative operator D q , we first prove the following lemma. Lemma 1. The twice-iterated 2D q-Appell polynomials A [2] m,q (x 1 , x 2 ) satisfy the following operational relations: and q s(s−1) 2 Proof. In view of the Equation (25), the proof of the above lemma requires a direct use of the identity (5). We, therefore, skip the details involved.
We now derive the q-recurrence relations for the 2I2DqAP A [2] m,q (x 1 , x 2 ).
m,q (x 1 , x 2 ) satisfy the following linear homogeneous recurrence relation: where Proof. Consider the following generating function: By taking the q-derivative of the Equation (47) partially with respect to t, we get Thus, upon factorizing G q (qx 1 , x 2 , t) occurring in the left-hand side and multiplying both sides of the identity (48) by t, we find that In view of the assumption (46) and the Equation (47), the Equation (49) becomes which, on using the Cauchy product rule, gives Finally, upon equating the coefficients of like powers of t on both sides of the above equation and dividing both sides of the resulting equation by [m] q , we get the assertion (45) of Theorem 4.
We now state and prove the following result.
Theorem 5. The following recurrence relation for the twice-iterated 2D q-Appell polynomials A [2] m,q (x 1 , x 2 ) holds true: Proof. We first use the Equation (47)  We next derive the q-difference equations which are satisfied by the twice-iterated 2D q-Appell polynomials.
Theorem 6. The twice-iterated 2D q-Appell polynomials A [2] m,q (x 1 , x 2 ) are the solutions of the following q-difference equations: Proof. The proof of the assertions (53) and (54) of Theorem 6 would follow directly upon using the Equations (43) and (44), respectively, in the recurrence relation (45).
In the next section (Section 3 below), the determinant forms for the 2I2DqAP are established.

The Twice-Iterated 2D q-Appell Polynomials from the Determinant Viewpoint
One of the important aspects of the study of any polynomial system is to find its potentially useful determinant representation. Recently, Keleshteri and Mahmudov [21] introduced the determinant definitions for the q-Appell polynomials and the 2D q-Appell polynomials. These polynomials are useful in finding the solutions of some general linear interpolation problems and can also be used for computational purposes. Khan and Riyasat [26], on the other hand, established the determinant expressions for the twice-iterated q-Appell polynomials. This fact provides motivation for us to establish the determinant definitions and the determinant expressions for the twice-iterated 2D q-Appell polynomials 2I2DqAP by proving the following result. Theorem 7. The 2I2DqAP A [2] m,q (x 1 , x 2 ) of degree m are defined by A [2] 0,q (x 1 , A [2] m,q (x 1 , Proof. Consider A [2] m,q (x 1 , x 2 ) as a sequence of the 2I2DqAP defined by the Equation (25). Also letȦ m,q and B m,q be two numerical sequences (the coefficients of the q-Taylor series expansions of functions) such thaṫ also satisfying the following condition:Ȧ q (t)ˆȦ q (t) = 1.
On using the Cauchy product rule for the two-series productȦ q (t)ˆȦ q (t), we geṫ that is, Next, upon multiplying both sides of the Equation (25) byˆȦ q (t), we geṫ Further, in view of the Equations (22), (58) and (59), the above Equation (62) becomes Now, on using the Cauchy product rule for the two series in the right-hand side of the Equation (63), we obtain the following infinite system for the unknowns A [2] m,q (x 1 , x 2 ): Obviously, the first equation of the system (64) leads to our first assertion (55). The coefficient matrix of the system (64) is lower triangular, so this helps us to obtain the unknowns A [2] m,q (x 1 , x 2 ) by applying the Cramer rule to the first m + 1 equations of the system (64). Accordingly, we can obtain where m ∈ N. Thus, upon expanding the determinant in the denominator and taking the transpose of the determinant in the numerator, we get Finally, after m circular row exchanges, that is, after moving the jth row to the (j + 1)st position for j = 1, 2, 3, · · · , m − 1, we arrive at our assertion (56) of Theorem 7. Theorem 8. The following identity for the 2I2DqAP A [2] m,q (x 1 , x 2 ) holds true: Proof. Expanding the determinant in the Equation (56) with respect to the (m + 1) st row, we get Next, by applying the same argument for the last determinant, we find that Again, we apply the same technique recursively until we arrive at the following consequence: Finally, upon summing up the series in the left-hand side of the Equation (68), we arrive at the assertion (67) of Theorem 8. Corollary 1. The following identity for the 2DqAPÄ n,q (x 1 , x 2 ) holds true:

Several Members of the Twice-Iterated 2D q-Appell Polynomials
During the last two decades, much research work has been done for different members of the family of the q-Appell polynomials and the 2D q-Appell polynomials. By making suitable selections for the functionsȦ q (t) andÄ q (t), the members belonging to the family of the twice-iterated 2D q-Appell polynomials A [2] k,q (x 1 , x 2 ) can be obtained. The 2D q-Bernoulli polynomials B m,q (x 1 , x 2 ), the 2D q-Euler polynomials E m,q (x 1 , x 2 ) and the 2D q-Genocchi polynomials G m,q (x 1 , x 2 ) are important members of the 2D q-Appell family. Therefore, in this section, we first introduce the 2D q-Euler based Bernoulli polynomials (2DqEBP) E B m,q (x 1 , x 2 ) and the 2D q-Genocchi based Bernoulli polynomials (2DqGBP) G B m,q (x 1 , x 2 ) by means of their respective generating functions and series definitions. We then explore other properties of these members.

The 2D q-Euler-Bernoulli Polynomials
Since, for A q (t) = 2 e q (t) + 1 and A q (t) = t e q (t) − 1 , the 2DqAP A m,q (x 1 , x 2 ) reduce to the 2DqEP E m,q (x 1 , x 2 ) and the 2DqBP B m,q (x 1 , x 2 ), respectively. Therefore, for the same choices of A q (t), that is, the 2I2DqAP reduce to 2DqEBP E B m,q (x 1 , x 2 ) and are defined by means of generating functions as follows: 2t The 2DqEBP E B m,q (x 1 , x 2 ) of degree m are defined by the following series: The following relation between the 2DqEBP E B m,q (x 1 , x 2 ) and the qEBP E B m,q (x 1 ) holds true: which, for x 2 = 1, yields The 2DqEBP E B m,q (x 1 , x 2 ) satisfy the following recurrence relation: Further, by taking B 0,q = 1, in the Equation (56), we obtain the determinant definition of the 2DqEBP E B m,q (x 1 , x 2 ) as given below.

Definition 1.
The 2D q-Euler-Bernoulli polynomials E B m,q (x 1 , x 2 ) of degree m are defined by where E m,q (x 1 , x 2 ) (m ∈ N 0 ) are the 2D q-Euler polynomials of degree m.

The 2D q-Genocchi-Bernoulli Polynomials
Since, for the 2DqAP A m,q (x 1 , x 2 ) reduce to the 2DqGP G m,q (x 1 , x 2 ) and the 2DqBP B m,q (x 1 , x 2 ), respectively. Therefore, for the same choices of A q (t), that is, the 2I2DqAP reduce to 2DqGBP G B m,q (x 1 , x 2 ) and are defined by means of generating functions as follows: The 2DqGBP G B m,q (x 1 , x 2 ) of degree m are defined by the following series: The following relation between the 2DqGBP G B m,q (x 1 , x 2 ) and the qGBP G B m,q (x 1 ) holds true: which, for x 2 = 1, gives      Tables 2 and 3.

Concluding Remarks and Observations
As long ago as 1910, Jackson [27] studied the q-definite integral of an arbitrary function f (t), which is defined as follows: We note also that Applying the double q-integral to both sides of the Equation (42), that is, m−1,q (t 1 , qt 2 ) d q t 1 d q t 2 = m,q (t 1 , t 2 ) d q t 1 d q t 2 , we have [m] q m−1,q (t 1 , qt 2 ) d q t 1 d q t 2 = m,q (t 1 , x 2 ) − A [2] m,q (t 1 , 0) d q t 1 .
Further investigations along the lines presented in this paper, which are associated with the various recent generalizations and extensions of the Apostol type Bernoulli, Euler and Genocchi polynomials introduced by, for example, Srivastava et al. (see [29,30]) may be worthy of consideration by the targeted readers.