On 2-Variables Konhauser Matrix Polynomials and Their Fractional Integrals

: In this paper, we ﬁrst introduce the 2-variables Konhauser matrix polynomials; then, we investigate some properties of these matrix polynomials such as generating matrix relations, integral representations, and ﬁnite sum formulae. Finally, we obtain the fractional integrals of the 2-variables Konhauser matrix polynomials.


Introduction
Special functions play a very important role in analysis, physics, and other applications, and solutions of some differential equations or integrals of some elementary functions can be expressed by special functions. In particular, the family of special polynomials is one of the most useful and applicable family of special functions. The Konhauser polynomials which were first introduced by J.D.E. Konhauser [1] include two classes of polynomials Y α n (x; k) and Z α n (x; k), where Y α n (x; k) are polynomials in x and Z α n (x; k) are polynomials in x k , α > −1 and k ∈ Z + . Explicit expressions for the polynomials Z α n (x; k) are given by Z α n (x; k) = Γ(α + kn + 1) n! n ∑ r=0 (−1) r n r x kr Γ(α + kr + 1) , (1) where Γ(·) is the classical Gamma function and for the polynomials Y α n (x; k), Carlitz [2] subsequently showed that where (a) n is Pochhammer's symbol of a as follows: (a) n = a(a + 1)(a + 2) . . . (a + n − 1), n ≥ 1, 1, It is easy to verify that the polynomials Y α n (x; k) and Z α n (x; k) are biorthogonal with respect to the weight function w(x) = x α e −x over the interval (0, ∞), which means where α > −1, k ∈ Z + and δ ij is the Kronecker delta. The Laguerre polynomials L α n (x) are defined as (see, e.g., [3]) For p, q ∈ N, we can define the general hypergeometric functions of p-numerator and q-denominator by such that β j = 0, −1, −2, . . . ; j = 1, 2, . . . , q. Then, according to [3], we can rewrite L α n (x) as For k = 1, we note that the Konhauser polynomials (1) and (2) reduce to the Laguerre Polynomials L α n (x) and their special cases; when k = 2, the case was encountered earlier by Spencer and Fano [4] in certain calculations involving the penetration of gamma rays through matter and was subsequently discussed in [5].
On the other hand, the matrix theory has become pervasive to almost every area of mathematics, especially in orthogonal polynomials and special functions. The special matrix functions appear in the literature related to statistics [6], Lie theory [7], and in connection with the matrix version of Laguerre, Hermite, and Legendre differential equations and the corresponding polynomial families (see, e.g., [8][9][10]). In the past few years, the extension of the classical Konhauser polynomials to the Konhauser matrix polynomials of one variable has been a subject of intensive studies [11][12][13][14]. Recently, many authors (see, e.g., [15][16][17][18]) have proposed the generating relations of Konhauser matrix polynomials of one variable from the Lie algebra method point of view and found some properties of Konhauser matrix polynomials of one variable via the Lie algebra technique; they also obtained operational identities for Laguerre-Konhauser-type matrix polynomials and their applications for the matrix framework.
The purpose of the present work is to introduce and study 2-variables Konhauser matrix polynomials and find the hypergeometric matrix function representations; we try to establish some basic properties of these polynomials which include generating matrix functions, finite sum formulae, and integral representations, and we will also discuss the fractional integrals of the 2-variables Konhauser matrix polynomials.
The rest of this paper is structured as follows. In the next section, we give basic definitions and previous results to be used in the following sections. In Section 3, we introduce the definition of 2-variables Konhauser matrix polynomials for parameter matrices A and B and some generating matrix relations involving 2-variables Konhauser matrix polynomials deriving the integral representations. Finally, we provide some results on the fractional integrals of 2-variables Konhauser matrix polynomials in Section 4.

Preliminaries
In this section, we give the brief introduction related to Konhauser matrix polynomials and recall some previously known results.
Let C N×N be the vector space of N-square matrices with complex entries; for any matrix A ∈ C N×N , its spectrum σ(A) is the set of all eigenvalues of A, A square matrix A ∈ C N×N is said to be positive stable if and only if β(A) > 0. Furthermore, the identity matrix and the null matrix or zero matrix in C N×N will be symbolized by I and 0, respectively. If Φ(z) and Ψ(z) are holomorphic functions of the complex variable z, which are defined as an open set Ω of the complex plane and A is a matrix in C N×N with σ(A) ⊂ Ω, then, from the properties of the matrix functional calculus [27,28], we have Let A be a positive stable matrix in C N×N . Then, Γ(A) is well defined as where t A−I = exp((A − I) ln t). Then, the matrix Pochhammer symbol (A) n of A is denoted as follows (see, e.g., [29][30][31]): The Laguerre matrix polynomials are defined by Jódar et al. [8] L where A ∈ C N×N is a matrix such that −k ∈ σ(A), ∀k ∈ Z + , (A + I) k are given by Equation (12) and λ is a complex number with Re(λ) > 0.
are matrices such that B j + kI are invertible for all integers k ≥ 0, the generalized hypergeometric matrix functions are defined as [32] It follows that, for λ = 1 in (13), we have For commuting matrices A i , B i , C i , D i , E i and F i in C N×N , we define the Kampé de Fériet matrix series as [32] F m 1 ,n 1 ,l 1 m 2 ,n 2 ,l 2 where A abbreviates the sequence of matrices A 1 , ...., A m 1 , etc. and D i + kI, E i + kI and F i + kI are invertible for all integers k ≥ 0.
If A ∈ C N×N is a matrix satisfying the condition and λ is a complex numbers with Re(λ) > 0, we recall the following explicit expression for the Konhauser matrix polynomials (see, e.g., [11]) and which are biorthogonal with respect to matrix weight function w(x) = x A e −λx over the interval (0, ∞).

2-Variables Konhauser Matrix Polynomials
In this section, we first introduce the 2-variables Konhauser matrix polynomials with parameter matrices A and B; then, we get the hypergeometric matrix function representations, generating matrix functions, finite summation formulas, and related results for the 2-variables Konhauser matrix polynomials.

Hypergeometric Representation
Now, by using (16) and (20), we obtain the hypergeometric matrix function representations where ∆(k; A) abbreviates the array of k parameters such that and F 1 k,l is defined in (16).

Remark 2.
If A ∈ C N×N is a matrix satisfying the condition (17), letting B = 0 and y = 0 in (22), we obtain where Z (A,λ) n (x; k) are Konhauser matrix polynomials in [11] and 1 F k is hypergeometric matrix function of 1-numerator and k-denominator defined in (14).

−nI,
where L A n (x) are the Laguerre's matrix polynomials defined in (15).

Generating Matrix Relations for the 2-Variables of Konhauser Matrix Polynomials
Generating matrix relations always play an important role in the study of polynomials, first, we give some generating matrix relations for the 2-variables of Konhauser matrix polynomials as follows: Theorem 1. Letting A, B ∈ C N×N be matrices satisfying the condition (17), we obtain the explicit formulae of matrix generating relations for the 2-variables Konhauser matrix polynomials as follows: where 0 F k and 0 F l are hypergeometric matrix functions of 0-numerator and k, l-denominator as (14), ∆(k; A + I) and ∆(l; B + I) are defined as (23), and the short line "−" means that the number of parameters is zero. (20)

Proof. From Equation
by using This completes the proof.
For a matrix E in C N×N , we can easily obtain the following generating relations for the 2-variables Konhauser matrix polynomial similar to Theorem 1 where F 1 k,p are defined in Equation (16), ∆(k; A + I) and ∆(l; B + I) are defined as Equation (23).

Corollary 1.
Letting A, B ∈ C N×N be matrices satisfying the condition (17), the following generating matrix relations of the 2-variables Konhauser matrix polynomials hold: where 0 F k and 0 F l are hypergeometric matrix functions of 0-numerator and k,l-denominator as (14).

Corollary 2.
Letting A, B, and E be matrices in C N×N satisfying the condition (17), we give explicit formulae of matrix generating relations for the 2-variables Konhauser matrix polynomials as follows: Considering the double series, Now, by making use of Theorem 1, we find Here, Equation (33) may be regarded as a double generating matrix relations for (20).
we find generating matrix relations of the Laguerre's matrix polynomials.

Some Properties of the 2-Variables Konhauser Matrix Polynomials
For the finite sum property of the 2-variables Konhauser matrix polynomials Z (A,B,λ,ρ) n (x, y, k, l), we get the generating relations together as follows: The integral representations for the 2-variables Konhauser matrix polynomials are derived in the following theorem.
where c 1 , c 2 are the paths around the origin in the positive direction, beginning at and returning to positive infinity with respect for the branch cut along the positive real axis.
Proof. The right side of the above formulae are deformed into and using the integral representation of the reciprocal Gamma function, which are given in [34] 1 Γ(z) = 1 2πi c e t t −z dt, where c is the path around the origin in the positive direction, beginning at and returning to positive infinity with respect for the branch cut along the positive real axis. Thus, from Equation (40), we obtain the following integral matrix functional