On the Two-Variable Analogue Matrix of Bessel Polynomials and Their Properties

: In this paper, we explore a study focused on a two-variable extension of matrix Bessel polynomials. We initiate the discussion by introducing the matrix Bessel polynomials involving two variables and derive specific differential formulas and recurrence relations associated with them. Additionally, we present a segment detailing integral formulas for the extended matrix Bessel polynomials. Lastly, we introduce the Laplace–Carson transform for the two-variable matrix Bessel polynomial analogue.


Introduction
In modern mathematics, special functions play a crucial role in various disciplines.Special versions of these functions have proven invaluable in various fields, including probability theory, computer science, mathematical physics, engineering, and many other areas (see [1,2]).
Bessel polynomials are important due to their natural occurrence in seemingly unrelated situations.For example, they appear in the solution of the wave equation in spherical polar coordinates (see [3]), in network synthesis and design (see [4]), in the analysis of the Student t-distribution (as shown in [5]), and in the development of a matrix approach suitable for solving differential equations with multiple orders and fractions of orders.In addition, Bessel polynomials play a role in the representation of the energy spectrum functions for a family of isotropic turbulence fields (see [6]).
In 1949, Krall and Frink [3] presented an article on what they termed Bessel polynomials.Within this work, they introduced the elementary Bessel polynomials in the following manner: where 2 F 0 is the Gauss hypergeometric function of a two-numerator as The generalized Bessel polynomials, denoted as Y ϵ (µ, ν; x) and extending the Bessel polynomials Y ϵ by introducing two parameters, are defined based on the research by Krall and Frink [3] as follows: ). ( where (µ) k is the Pochhammer symbol, defined by They also obtained certain recurrence relations between these polynomials and one generating function for the Bessel polynomials proper and also provided some of their qualities including orthogonality and the linkages between them and Bessel functions.
The Laplace-Carson transform of the function G(s) for all s ≥ 0 is defined as [12] L where L is Laplace-Carson transform operator.The field of generalized special matrix functions has undergone significant development in recent years.This interest is due to various reasons.Focusing on practical applications, one realizes that the use of new classes of special matrix functions in certain physical problems has led to solutions that are difficult to obtain using traditional analytical and numerical methods.In this context, Hermite, Chebyshev, Jacobi, Laguerre, and Gegenbauer matrix polynomials have been introduced and studied (see, e.g., [13][14][15][16][17]).
The objective of this paper is to present a novel two-variable analogue, denoted as (z, w), and to derive specific outcomes related to the two-variable matrix Bessel polynomials Y (θ,ϑ) ϵ (z, w).Additionally, we explore applications involving the Laplace-Carson transform of functions.
The structure of this article is outlined as follows: In Section 2, we give a brief introduction to certain matrix functions that are important for the further development of the article.In Section 3, we introduce a novel extension of Bessel matrix polynomials called Y (θ,ϑ) ϵ (z, w) and state several theorems on recurrence relations and the derivation formula for this extension.Section 4 deals with various integral formulas for the Bessel matrix polynomials Y (θ,ϑ) ϵ (z, w).In Section 5, we give theorems on the Laplace-Carson transform of functions containing the new extension of the matrix Bessel polynomials Y (θ,ϑ) ϵ (z, w).The concluding remarks and future work are presented in Section 6.

Some Definitions and Notations
Let C r be the r-dimensional complex vector space and C r×r denote all square complex matrices of order r.A positive stable matrix θ in C r×r is denoted as Re(λ) > 0 for all λ ∈ σ(θ), where σ(θ) is the set of all Eigenvalues of θ and θ n ̸ = 0, then we call a matrix polynomial in z of degree n.
In the matrix complex space C r×r , we use I to represent the identity matrix and O for the zero matrix.
The spectrum of a matrix θ in C r×r is defined as the set of all its Eigenvalues and is denoted by σ(θ).If g(z) and h(z) are holomorphic functions defined on an open set D ⊆ C and θ is a matrix in C r×r , such that σ(θ) ⊂ D, then the commutative property g(θ)h(θ) = h(θ)g(θ) holds ( [18,19]).
Laguerre matrix polynomials may be defined in [16] by The Whittaker matrix function W θ,−θ+I/ϵ+2 (z) is defined [26] by the matrix function W θ,−θ+I/ϵ+2 (z) defined [26] by We also see that the Whittaker matrix functions and Laguerre's matrix polynomials are the fundamental components of the generalized Bessel matrix polynomials.In fact, we have and The integral representation is an immediate result of ( 12) as The orthogonality of the generalized Bessel matrix polynomials on the unit circle with regard to the weight matrix function ((cf.[26]): This is satisfied by the associated matrix nonhomogeneous equation For n ̸ = ϵ, we have Now, a few features of the Pochhammer symbol are needed in our present study.

The Matrix Bessel Polynomial
In this section, we introduce the matrix Bessel polynomial Y (θ,ϑ) ϵ (z, w) of two variables and discuss some important basic properties of it as follows.
Definition 1.Let θ and ϑ be positive stable matrices in C r×r satisfying the condition (9), then the matrix Bessel polynomials of two variables Y Remark 1.The following lists some particular results on matrix Bessel polynomials: 2. Also, from (27), we can deduce that 3.
If we put w = 0 in (27), we find that 4.
If putting w = 0 and ϑ = O in (27), we obtain the relation: Theorem 1. Suppose that θ and ϑ are matrices in C r×r satisfying the condition (9) such that ϵ, s, and r are non-negative integer numbers and where all matrices are commutative, then we have Proof.Using the R.H.S. of (32), we determine that By using (24), we have Putting ϵ − r − s = v and rearranging the terms, we obtain This complete the proof.
Theorem 2. If we assume that θ and ϑ are matrices in C r×r satisfying the condition (9) such that ϵ, p, q, and r are non-negative integer numbers and where all matrices are commutative, then we have that Proof.By using the definition in (27), we obtain Now, we find that By using (20), we have and this finish the proof.
Theorem 3. Let parameters θ and ϑ be positive stable matrices in C r×r satisfying the condition (9), where ϵ and r are non-negative integer numbers and where all matrices are commutative, then we obtain Proof.The Bessel matrix polynomials are expanded in double series using (27) in this case, and it is sufficient to show that the coefficient of z v w r is the same on both sides of Equation (34).Form the definition in (27), we have the L.H.S. as then, we find the coefficient of z v w r given by Also, we see the R.H.S. of Equation (34) as follows: The coefficient of z v w r of the R.H.S is We find this from (36) as Hence, from (35) and (37), the coefficient of z v w r is the same on both sides of Equation (34), and this complete the proof.

Some Integrals Involving the Matrix Bessel Polynomial
In this section, our attention is directed towards presenting integral representations for the matrix Bessel polynomial Y (θ,ϑ) m (z, w) in the form of theorems.
Theorem 4. Let θ, ϑ, and R be positive stable matrices in C r×r satisfying the condition (9) and where all matrices are commutative, then we obtain Where a, b > 0 Proof.By using the definition (27) in the L.H.S., we obtain By substituting p = v + r, rearranging the terms, and using (1), we obtain This concludes the proof.
Theorem 5.For θ andϑ, R and K are matrices in C r×r such that RK = KR and satisfying the condition (9), then we obtain where a, b > 0 Proof.By using the L.H.S. of ( 39) and (27), we find that By putting v + r = p and rearranging the terms, we have This complete the proof.
Theorem 6.If θ and ϑ are positive stable matrices in C r×r such that θϑ = ϑθ and satisfy the condition (9), then the matrix Bessel polynomials Y (θ,ϑ) ϵ (z, w) satisfy the following integral representation as where a, b > 0.
Proof.By using the L.H.S. of (40) and (27), we obtain This leads to the assertion (40).
Theorem 7. Let θ and ϑ be positive stable matrices in C r×r such that θϑ = ϑθ satisfies the condition (9), then we have where a, b > 0 Proof.By using the L.H.S. of ( 41) and (27), we obtain Substituting z = tu, we obtain and this complete the proof.

The Laplace-Carson Matrix Transform
In this section, we introduce the Laplace-Carson transform of the extended matrix Bessel polynomial.Initially, we provide the definition of the Laplace-Carson transform for matrix functions.Definition 2. Let F(u, s) be a function defined on the collection of all the positive stable matrices observed in C r×r , then the two-dimensional Laplace-Carson transform is provided by: such that the integral on the right side of the equation exists.
Theorem 8. Suppose that θ and ϑ are in C r×r satisfying the condition (9) and where all matrices are commutative, then we obtain the Laplace-Carson transform as where Re(u) > 0, Re(s) > 0.
Proof.By using Definition 2, we find that where K = (θ + I) ϵ (ϑ + I) ϵ .Putting uz = t and sw = n, we have which is the required proof.
Theorem 10.Let θ and ϑ be matrices in C r×r , satisfying the condition (9) and where all matrices are commutative, then we obtain where A = (θ + I) ϵ (ϑ + I) ϵ , and the Humbert matrix function ψ 1 (θ, ϑ; R, R ′ ; z, w) of two complex variables is given in [29] as This concludes the proof of Theorem 10.

Conclusions
The article examines a two-variable counterpart of matrix Bessel polynomials and investigates specific differential formulas and recurrence relations associated with them.Part of the integral formula for this new extension of matrix Bessel polynomials is also presented.In addition, we presented the Laplace-Carson transform for the analogous matrix Bessel polynomial with two variables.Future research efforts could be devoted to unveiling further properties and features of these polynomials.This could include the exploration of extended and generalized forms, as well as integral representations.The analysis of these facets may contribute to a deeper understanding of the polynomials and their behavior.