Next Article in Journal
Lepton Flavor Universality Tests in Semileptonic bc Decays
Next Article in Special Issue
Linear Arrangement of Euler Sums with Multiple Argument
Previous Article in Journal
Anthropometric Formulas Repurposed to Predict Body Fat Content from Ultrasound Measurements of Subcutaneous Fat Thickness
Previous Article in Special Issue
Convolution Theorem for (p,q)-Gamma Integral Transforms and Their Application to Some Special Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Fractional Operators Involving the Incomplete Mittag-Leffler Matrix Function and Its Applications

1
College of Mathematics and System Sciences, Xinjiang University, Urumqi 830046, China
2
Department of Mathematics, Faculty of Science, Al-Azhar University, Assiut 71524, Egypt
3
Mathematics Department, College of Science, King Khalid University, Abha 61413, Saudi Arabia
*
Authors to whom correspondence should be addressed.
Symmetry 2024, 16(8), 963; https://doi.org/10.3390/sym16080963
Submission received: 15 May 2024 / Revised: 20 June 2024 / Accepted: 26 June 2024 / Published: 29 July 2024
(This article belongs to the Special Issue Research in Special Functions)

Abstract

:
In this study, we derive multiple incomplete matrix Mittag-Leffler (ML) functions. We systematically investigate several properties of these incomplete matrix ML functions, which include some general properties and distinct representations of integral transforms. We further study the properties of the Riemann–Liouville fractional integrals and derivatives related to the incomplete matrix ML functions. Additionally, some interesting special cases of this work are highlighted. Finally, we establish the solution to the kinetic equations as an application.

1. Introduction

Gosta Mittag-Leffler (ML), a Swedish mathematician, introduced the given function in 1903 [1].
E α ( z ) = n = 0 z n Γ ( α n + 1 ) .
The well-known Gamma function, represented by Γ , is involved in the definition of the ML function, E α ( z ) , for a complex variable z. The variable α in this case is non-negative. When α = 1 , this function corresponds to the standard exponential function, which is a straightforward extension of the exponential function. It is estimated between the hypergeometric function 1 1 z and the pure exponential function for 0 < α < 1 . Its application in a number of disciplines, such as engineering, biology, chemistry, physics, and applied sciences, has increased its importance over the past 20 years. Whenever fractional order differential equations or fractional order integral equations are obtained, the ML function naturally appears. Understanding its characteristics and applications in a variety of domains has been made possible by its generalisation, as Wiman first investigated in 1905. Re ( z ) and Im ( z ) denote the real part and imaginary part, respectively, of a complex number z.
The two parametric Mittag-Leffler function is defined by [2]
E α , β ( z ) = n = 0 z n Γ ( α n + β ) ,
where α , β C with Re ( α ) > 0 and Re ( β ) > 0 .
In 1971, Prabhakar [3] introduced the function E α , β γ ( z ) in the form of
E α , β γ ( z ) = n = 0 ( γ ) n z n Γ ( α n + β ) n ! ,
where α , β , γ C and Re ( α ) > 0 , Re ( β ) > 0 , Re ( γ ) > 0 , where ( γ ) n is the Pochhammer symbol, as
( γ ) n = γ ( γ + 1 ) . . . ( γ + n 1 ) = Γ ( γ + n ) Γ ( γ ) , n 1 , 1 , n = 0 .
Thus, it is clear that the following special cases holds.
Shukla and Prajapati [4] introduced the ML function generalized in 2007, as follows:
E α , β γ , k ( z ) = n = 0 ( γ ) n k z n Γ ( α n + β ) n ! ,
where α , β , γ C and Re ( α ) > max { 0 , Re ( k ) 1 } , Re ( β ) > 0 , Re ( k ) > 0 .
After that, many other authors investigated and explored several properties and applications of the generalized ML function in the solution of fractional order integral and fractional order differential equations (see, e.g., [5,6,7,8,9]). Later on, a detailed and comprehensive review was provided by Haubold et al. (see [10,11,12]).
The solutions for special matrix differential equations have become accessible in recent years due to their special function with matrix parameters. These are systems of differential equations as these matrix differential equations are each carried out by a corresponding scalar special function. A corresponding scalar special function may fulfil the system of equations that represent the other results for special matrix functions, including generating functions, series definitions, recurrence relations, etc. (see [13,14,15,16,17,18]).
Additionally, the ML function with matrix arguments is equally as useful as its scalar version, and it can be utilized successfully in a wide range of applications, such as control theory and other related fields, the efficient and stable solution of systems of fractional differential equations (FDEs), the estimation of the solution for particular multiterm FDEs, and other related fields (see [19,20]).
At the present, the generalized version allows for matrices as arguments and provides a powerful tool for solving fractional differential and integral equations in matrix form. In 2018, Garrappa et al. [21] computed the Mittag-Leffler function with matrix arguments, with some applications in fractional calculus, and in 2023, Pal et al. [22] introduced a special matrix analog of the four-parameter Mittag-Leffler function.
The main objective of this work is to develop various aspects of incomplete matrix ML functions, which include certain fundamental properties and various representations of integral transformations, and to analyse incomplete ML functions with matrix parameters by means of the incomplete Pochhammer symbol. We investigate some of the features of the Riemann–Liouville fractional integrals and derivatives associated with the incomplete matrix ML functions. We additionally highlight some interesting particular examples of our major results. In the conclusion, incomplete matrix ML functions are utilized as an application to address the kinetic equations.
This work follows the following outline: in Section 2, we look at some of the most common definitions and basic applications of matrix arguments in special functions. A matrix variation of the incomplete ML functions is presented and its convergence is investigated in Section 3. In Section 4, representations based on the Euler-Beta and Laplace transforms for the incomplete matrix ML functions are presented, in addition to some important theorems. Utilizing incomplete matrix ML functions, we acquire several interesting characteristics of the fractional calculus operators in Section 5. The implementation of ML functions in the fractional kinetic equation solution will be discussed in Section 6. Finally, some concluding remarks are collected in Section 7.

2. Basic Definitions and Preliminaries

Consider C h , which denotes a complex vector space of dimension h, and let C h × h be the set of all square complex matrices of order h. For any matrix T in C h × h , σ ( T ) represents the spectrum of T, which is the set containing all eigenvalues of T, referred to as T, then
a ( T ) = max { Re ( z ) : z σ ( T ) } , b ( T ) = min { Re ( z ) : z σ ( T ) } ,
where a ( T ) = b ( T ) and a ( T ) represent the spectral abscissa of T. The square matrix T becomes positive stable if and only if b ( T ) > 0 , and we refer to it as [16]
e t T e t a ( T ) s = 0 h 1 ( T h 1 2 t ) s s ! ; h 1 ; a n d t R
and
n T = e ( ln n ) T n a ( T ) s = 0 h 1 ( T h 1 2 ln n ) s s ! ; h 1 ; a n d t R .
The Beta matrix function is described similarly (see, e.g., [23,24,25]), where T and U are positive stable and commuting matrices in C h × h , so that the matrices T + n I , U + n I , and T + U + n I are invertible for any integer n 0 . Now
B ( T , U ) = 0 1 t T I ( 1 t ) U I d t = Γ ( T ) Γ ( U ) [ Γ ( T + U ) ] 1 ,
where
Γ ( T ) = 0 e t t T I d t , t T I = exp ( ( T I ) ln t ) ,
is the matrix function of Gamma (see, for example, [26]). In addition, if for all integers n 0 , we have
T + n I is an invertible matrix ,
then, Γ ( T ) is proven to be invertible and its inverse is denoted by Γ 1 ( T ) . This leads to the emergence of the Pochhammer symbol with a matrix argument, as detailed in references such as [27].
( T ) n = T ( T + I ) ( T + ( n 1 ) I ) = Γ ( T + n I ) Γ 1 ( T ) , n 1 , I , n = 0 .
As described in Abdalla’s work [24], the Gamma matrix functions can be decomposed into two incomplete Gamma matrix functions.
γ ( T , z ) = 0 z e t t T I d t ,
Γ ( T , z ) = z e t t T I d t
and
Γ ( T ) = γ ( T , z ) + Γ ( T , z ) .
where z C and T are a positive stable matrix in C h × h .
Consider that T represents any stable positive matrix within a complex system C h × h . In this context, for every non-negative x, the Pochhammer incomplete matrix ( T ; x ) m , along with its complement [ T ; x ] m , as discussed in [28], is described as follows:
( T ; x ) m = γ ( T + m I , x ) Γ 1 ( T )
and
[ T ; x ] m = Γ ( T + m I , x ) Γ 1 ( T ) ,
satisfy the decomposition formula, respectively.
( T ; x ) m + [ T ; x ] m = ( T ) m .
The Laplace transform of the original ϕ (t) is [29]
ϕ ¯ ( h ) = L ϕ ( t ) ( h ) = 0 e h t ϕ ( t ) d t , R e ( h ) > 0 ,
where ϕ ¯ ( h ) represents the Laplace transform of ϕ ( t ) .
Moreover, the Euler-Beta transform (see [29]) of the function f ( z ) is given by
β f ( z ) ; a , b = 0 1 z a 1 ( 1 z ) b 1 f ( z ) d z .
Definition 1 
([22]). Let T and U be two positive stable matrices in C h × h , such that U satisfys the condition (11). Then, the matrix ML function is as follows:
E α , U β , T ( z ) = Γ 1 ( T ) m = 0 Γ ( T + β m I ) Γ 1 ( U + α m I ) z m m ! = m = 0 ( T ) β m Γ 1 ( U + α m I ) z m m ! ,
where α , β R + .

3. Incomplete Matrix ML Functions

In this section, we discuss a matrix variant of the incomplete ML functions and develop its convergence for | z | = 1 , and we demonstrate some properties and differential relations.
Definition 2. 
Suppose T and U are positive stable matrices in C h × h , such that U satisfies condition (11); we can now define a matrix incomplete ML functions as
E α , U β , T ( z ) = n = 0 Γ 1 ( U + n α I ) ( T ; x ) β n z n n ! = Γ 1 ( T ) n = 0 γ ( T + n β I , x ) Γ 1 ( U + n α I ) z n n !
and
E α , U β , T ( Z ) = n = 0 Γ 1 ( U + n α I ) [ T ; x ] β n z n n ! = Γ 1 ( T ) n = 0 Γ T + n β I , x Γ 1 ( U + n α I ) z n n ! ,
where α , β , x R + , and we can find that
E α , U β , T ( z ) + E α , U β , T ( Z ) = E α , U β , T ( z ) .
Remark 1. 
(i) 
For x = 0 in (22) and (23), we find the definition in (21)
(ii) 
By using (22) and (23), we find that
E α , U β , T ( z ) + E α , U β , T ( Z ) = E α , U β , T ( z ) .
Theorem 1. 
Suppose T and U are two positive constants if the following relation holds a ( U ) > b ( T ) in C h × h , then, the matrix of incomplete ML functions (22) converges completely for the given | z | = 1 .
Proof. 
Consider that relation a ( U ) > b ( T ) holds. Under this assumption, there exists a positive number λ , such as
a ( U ) b ( T ) = 2 λ .
Now, we can write
n 1 + λ Γ 1 ( T ) γ ( T + β n I , x ) Γ 1 ( U + α n I ) 1 n ! = n λ ( n 1 ) ! n T Γ 1 ( T ) γ ( T + β n I , x ) ( n 1 ) ! n T × n U ( n 1 ) ! Γ 1 ( U + α n I ) n U .
From (8) and a ( U ) = b ( U ) , we have
n T n U n b ( T ) a ( U ) u = 0 r 1 U r 1 / 2 ln n u u ! × u = 0 r 1 T r 1 / 2 ln n u u ! n 2 λ u = 0 r 1 max { U , T } r 1 / 2 ln n u u ! 2 .
Letting n in (27) and (28), for | z | = 1 , we have
lim n n 1 + λ Γ 1 ( T ) γ ( T + β n I , x ) Γ 1 ( U + α n I ) 1 n ! lim n n λ ( n 1 ) ! n T Γ 1 ( T ) γ ( T + β n I , x ) ( n 1 ) ! n T n U ( n 1 ) ! Γ 1 ( U + α n I ) n U lim n n λ ( n 1 ) ! n T Γ 1 ( T ) γ ( T + β n I , x ) ( n 1 ) ! n U ( n 1 ) ! Γ 1 ( U + α n I ) × u = 0 r 1 max { U , T } r 1 / 2 ln n u u ! 2 = 0 ,
on | z | = 1 the series (22) is convergent absolutely. □
Theorem 2. 
Assume T and U are positive stable matrices in C h × h and U + n I is invertible for every integer n 0 . The incomplete matrix ML function defined by (23), with α , β > 0 , satisfies
E α , U β , T ( z ) = U E α , U + I β , T ( z ) + α z d d z E α , U + I β , T ( z ) .
Proof. 
By using (23), we find that
U E α , U + I β , T ( z ) + α z d d z E α , U + I β , T ( z ) = U E α , U + I β , T ( z ) + α z d d z n = 0 Γ 1 ( U + ( n α + 1 ) I ) [ T ; x ] β n z n n ! = U E α , U + I β , T ( z ) + n = 0 Γ 1 ( U + ( n α + 1 ) I ) [ T ; x ] β n ( α n I + U U ) z n n ! = E α , U β , T ( z ) .
The Theorem 2 proof has become complete. □
Remark 2. 
For x = 0 in (30), we have
E α , U β , T ( z ) = U E α , U + I β , T ( z ) + α z d d z E α , U + I β , T ( z ) .
Theorem 3. 
Suppose U, T in C h × h is a positive stable matrix, such that U satisfies a condition (11). Then, the derivative of matrix of the incomplete ML function can be represented by
(i) 
d d z m E α , U β , T ( z ) = ( T ) β m E α , U + m α I β , T + β m I ( z ) ,
(ii) 
d d z m z U I E α , U β , T ( w z α ) = z U ( m + I ) I E α , U m I β , T ( w z α ) .
Proof. 
(i) By applying term-wise differentiation m times, (23) gives
d d z m E α , U β , T ( z ) = n = m Γ 1 ( U + n α I ) [ T ; x ] β n z n m ( n m ) ! = n = 0 Γ 1 ( U + ( α ( n + m ) I ) T ; x β ( n + m ) z n n ! = ( T ) β m n = 0 Γ 1 ( U + α ( n + m ) I ) [ T + β m I ; x ] β n z n n ! = ( T ) β m E α , U + m α I β , T + β m I ( z ) .
(ii) By using (23), we obtain
d d z m z U I E α , U β , T ( w z α ) = n = 0 Γ 1 ( U + α n I ) [ T ; x ] β n w n n ! d d z m z U + ( α n 1 ) I .
Differentiating term by term under the sign of summation, we find that
d d z m z U I E α , U β , T ( w z α ) = n = 0 Γ 1 ( U + ( α n m ) I ) [ T ; x ] β n w n n ! z U + ( α n m 1 ) I = z U ( m + I ) I E α , U m I β , T . ( w z α ) .
This complete the proof. □
Remark 3. 
(i) 
Taking α and β = 1 in (ii), we obtain
d d z m z U I 1 Γ 1 [ T ; x ] ; U ; w z = Γ 1 ( U m I ) Γ ( U ) z U ( m + I ) I 1 Γ 1 [ T ; x ] ; U m I ; w z .
(ii) 
Using x = 0 in Theorem 3, we obtain
(a) 
d d z m E α , U β , T ( z ) = ( T ) β m E α , U + m α I β , T + k m I ( z ) .
(b) 
d d z m z U I E α , U β , T ( w z α ) = z U ( m + I ) I E α , U m I β , T ( w z α ) .

4. Some Integral Transforms of Matrix Incomplete ML Function

We establish the Euler-Beta and Laplace transform representations for the incomplete matrix ML functions in this section. Utilizing incomplete Fox–Wright matrix functions, we first established incomplete matrix ML functions as follows:
Definition 3. 
Let T P and U q be positive stable matrices in C h × h , such that U q + n I is invertible for every integer n 0 ; thus, we can defne the incomplete Fox–Wright matrix function as follows
ψ q ( γ ) p ( T 1 , β 1 ; x ) , ( T 2 , β 2 ) . . . . ( T P , β P ) ( U 1 , α 1 ) . . . . ( U q , α q ) ; z
= n 0 γ ( T 1 + n β 1 I ; x ) l = 2 p Γ ( T l + n β l I ) j = 1 q Γ 1 ( U j + n α l I ) z n n ! ,
where p and q are finite positive integers.
Definition 4. 
Suppose T, U in C h × h is positive stable matrices, as U satisfies Equation (11). Using the incomplete Fox–Wright matrix, we can generate the incomplete matrix ML functions as follows:
E α , U β , T ( z ) = Γ 1 ( T ) ψ 1 ( γ ) 1 ( T , β ; x ) ( U , α ) ; z
and
E α , U β , T ( z ) = Γ 1 ( T ) ψ 1 ( γ ) 1 ( T , β ; x ) ( U , α ) ; z .
Theorem 4. 
Suppose A, T, and U in C h × h are positive stable matrices, such that U holds in Equation (11). Then, the Laplace transform representation of the incomplete matrix ML function satisfies
L z A I E α , U β , T w z m = s A Γ 1 ( T ) ψ 1 ( Γ ) 2 ( T , β ; x ) , ( A , m ) ( U , α ) ; w s m .
Proof. 
From (19) of the Laplace transform and utilizing Equation (23), we obtain
L z A I E α , U β , T w z m = 0 z A I e s z E α , U β , T w z m d z = 0 z A I e s z n = 0 Γ 1 ( U + α n I ) [ T ; x ] β m ω n z m n n ! d z = n = 0 Γ 1 ( U + α n I ) [ T ; x ] β m ω n n ! 0 z A + ( m n 1 ) I e s z d z = n = 0 Γ 1 ( U + α n I ) [ T ; x ] β m s A + m n I 1 Γ ( A + m n I ) ω n n ! = Γ 1 ( T ) s A ψ 1 ( Γ ) 2 ( T , β ; x ) , ( A , m ) ( U , α ) ; w s m .
Thus, the proof is completed. □
Remark 4. 
(i) 
If we put A = U , β = 1 , and m = α in Equation (36), we have
0 z U I e s z E α , U 1 , T ( w z α ) = s U 1 Γ 0 ( T ; x ) ; ; w s α .
(ii) 
If putting x = 0 in Equation (36), we find that
L z A I E α , U β , T w z m = Γ 1 ( T ) s A ψ 1 ( Γ ) 2 ( T , β ) , ( A , m ) ( U , α ) ; w s m .
Theorem 5. 
The Euler-Beta transform repression of incomplete matrix ML function is given as follows:
β E α , U β , T w z m ; A , B = Γ 1 ( T ) Γ ( B ) ψ 2 ( Γ ) 2 ( T , β ; x ) , ( A , m ) ; ( U , α ) , ( A + B , m ) ; w ,
where A, B, T, and U are positive stable matrices in C h × h , such that U + k I is invertible for all k 0 .
Proof. 
From Equation (20) of the Euler-Beta transform, one can easily obtain
β E α , U β , T w z m ; A , B = 0 z A I ( 1 z ) B I E α , U β , T ω z m d z = 0 1 z A I ( 1 z ) B I n = 0 Γ 1 ( U + n α I ) [ T ; x ] β n w n z m n n ! = n = 0 Γ 1 ( U + n α I ) [ T ; x ] β n w n n ! 0 1 z A + ( n m 1 ) I ( 1 z ) B I d z = n = 0 Γ 1 U + n α I [ T ; x ] β n Γ 1 ( A + B + n m I ) Γ ( B ) Γ ( A + n m ) w n n ! ,
which, in view of the definitions of ψ 2 ( Γ ) 2 in Definition 3, provide the proper representation (38). □
Remark 5. 
(i) 
If putting A = U and α = m in (36), we obtain
0 z U I ( 1 z ) B I E m , U β , T w z m d z = Γ ( B ) E m , U + B β , T ( w ) .

5. Fractional Calculus Operators with Incomplete Matrix ML Functions

We derive several interesting features of the fractional calculus operators with respect to incomplete matrix ML functions in this section.
For an operator of fractional order μ and x > 0 , the integral and derivatives of Riemann–Liouville with respect to Re ( μ ) > 0 are presented in the subsequent form (see [13,30])
( I a μ f ) ( x ) = D t μ 0 [ f ( t ) ] = 1 Γ ( μ ) a x ( x t ) μ 1 f ( t ) d t .
Moreover,
D a μ f ( x ) = I a n μ D n f ( x ) , D = d d x .
Bakhet [13] studied the fractional order integrals and derivatives using the operators (39) and (40) as follows
Definition 5. 
Let T be a stable positive matrix in C h × h with the properties Re ( μ ) > 0 and μ C . The fractional integrals of order μ in the Riemann–Liouville sense is defined as
I μ ( x T ) = 1 Γ ( μ ) 0 x ( x t ) μ 1 t T d t .
Lemma 1. 
Let T to be a positive stable matrix in C h × h , such that Re ( μ ) > 0 . Then, the Riemann–Liouville integrals fractional of order μ can be written as
I μ ( x T I ) = Γ ( T ) Γ 1 ( T + μ I ) x T + ( μ 1 ) I .
Theorem 6. 
Suppose T and U are a positive stable matrices in C h × h and Re ( μ ) > 0 ; then, the fractional integration of incomplete matrix ML function can be denoted as
I μ x U I E α , U β , T w x α = x U + ( μ 1 ) I E α , U + μ I β , T w x α .
Proof. 
From Definition (5) and (23),we obtain
I μ x U I E α , U T , k w x α = 1 Γ ( μ ) 0 x ( x t ) μ 1 t U I E α , U T , k w x α d t = m = 0 [ T ; x ] β m Γ 1 U + α m I β m m ! × I μ x U + ( α m 1 ) I .
By using Lemma 1, we obtain
I μ x U I E α , U β , T w x α = m = 0 [ T ; x ] β m Γ 1 U + α m I β m m ! × x U + ( α m + μ 1 ) I Γ U + α m I Γ 1 U + ( α m + μ I ) = m = 0 [ T ; x ] β m Γ 1 U + ( α m + μ I ) β m m ! × x U + ( α m + μ 1 ) I .
Hence the Theorem 6 is proved. □
Corollary 1. 
Suppose x = 0 is substituted in Theorem 6; then, the fractional integration of the Mittag-Leffler matrix function will be
I μ x U I E α , U β , T κ x α = x U + ( μ 1 ) I E α , U + μ I β , T κ x α ,
where Re ( μ ) > 0 .
Theorem 7. 
Suppose U and T are a positive stable matrices in C h × h and 0 < Re ( μ ) . Then, the fractional derivative of the incomplete matrix ML function denotes
D μ x U I E α , U β , T w x α = x U ( μ + 1 ) I E α , U μ I β , T w x α .
Proof. 
By using (40), we find
D μ x U I E α , U β , T w x α = d d x n I n μ [ x U I E α , U β , T w x α .
Using Theorem 6 gives
D μ x U I E α , U β , T w x α = d d x n x U ( n μ 1 ) I E α , U ( n μ ) I β , T w x α .
By using (ii) in Theorem 7, this immediately yields the desired proof. □
Corollary 2. 
When putting x = 0 in Theorem (7), we find the fractional derivative of matrix ML function as
D μ x U I E α , U β , T κ x α = x U ( μ + 1 ) I E α , U μ I β , T κ x α .
where, Re ( μ ) > 0 .

6. Application to the Solution of Fractional Kinetic Equation

The resolution of fractional kinetic equations (FKEs) has attracted significant attention from researchers in a number of applied scientific domains, which include engineering, dynamical systems, physics, and control systems. Its ability to support the development of mathematical models for an extensive variety of physical procedures and mathematical physics applications is the reason that it has attracted this attention. Kinetic equations (KEs) are crucial for mathematical physics and natural science as they explain the continuity of matter’s motion in a variety of astrophysical situations. As stated in references like [12,31], recent research has discovered various fractional calculus operators that enable extending and generalizing FKEs. Haubold and Mathai established a functional differential equation in [12], which relates the rate of reaction change, decomposition rate, and production rate. It appears as follows:
d N d t = d ( N t ) + p ( N t ) ,
where N = N ( t ) represents the reaction rate, d = d ( N ) represents the destruction rate, p = p ( N ) denotes the production rate, and N t satisfies the function defined by N t ( t * ) = N ( t t * ) , t * > 0 .
The differential equation given below is a particular case of (45) when spatial fluctuations or inhomogeneities in the value N ( t ) are ignored.
d N d t = m i N i ( t ) ,
The initial conditions specified for this differential equations are N i ( t = 0 ) = N 0 , where i represents the number of density of species at time t = 0 , m i > 0 . If term i is ignored and the specific kinetic Equation (46) is integrated, the following relation will be obtained
N ( t ) N 0 = m D t 1 0 N ( t ) .
Within this context, m 0 is a constant value, while D t 1 stands for the Riemann–Liouville integral operator with a degree of μ = 1 . The fractional kinetic equation (FKE) is redefined by Haubold and Mathai as follows, as detailed in [12].
N ( t ) N 0 = m μ D t μ 0 N ( t ) .
where D t μ 0 is given in (39).
In light of this, the solution for N ( t ) can be expressed as follows:
N ( t ) = N 0 r = 0 ( 1 ) r Γ ( μ r + 1 ) ( m t ) μ r = N 0 E μ ( m μ t μ ) ,
where E μ ( m μ t μ ) denotes the ML function (see [32,33]).
Furthermore, Saxena Kalla created an alternative FKE, as given in [31,34,35,36].
N ( t ) N 0 f ( t ) = m μ D t μ 0 N ( t ) , m > 0 , R e ( μ ) > 0 ,
here, m is a constant, N ( t ) represents the species amount, which is initially N 0 = N ( 0 ) at time level t = 0 ; however, f is an integrable function at interval ( 0 , ) .
As discussed in [34,37,38,39,40,41], a variety of research articles have appeared recently in this field utilising FKEs to solve different integral transforms, such as Fourier, Laplace, Sumudu, and Mellin transforms, including special functions and matrix functions.
In this theorem, we examined the solutions for FKEs that require the extension of incomplete matrix ML functions.
Theorem 8. 
Consider T, U, and R holds (11), where T, U, and R are positive stable matrices in C h × h . Then, for Re ( v ) > 0 and t , w C , the following generalized FK matrix equation of the incomplete matrix ML functions satisfies the following equation
N ( t ) I N 0 E α , U β , T ( w t ) = R ν D t v 0 N ( t )
has the solution
N ( t ) = N 0 n = 0 [ T ; x ] β n Γ 1 ( U + α n I ) × ( w t ) n E ν , β + 1 ( R ν t ν ) ,
where E v , β + 1 ( R v t v ) is called the generalized the Mittag–Leffler function (see [33]).
Proof. 
By using the Laplace transform and using incomplete matrix ML function, we have
N * ( s ) I = N 0 0 e s t Γ 1 ( U + α n I ) [ T ; x ] β n ( w t ) n d t R v s v N * ( s ) ,
where N * ( s ) = w { N ( t ) ; s } and
1 + R s v N * ( s ) = N 0 n = 0 [ T ; x ] β n Γ 1 ( U + n α I ) Γ ( ( n + 1 ) I ) w n s n + 1 ,
we find that
N * ( s ) I = N 0 n = 0 [ T ; x ] β n w n s n + 1 ; Γ 1 ( U + n α I ) 1 + R s v 1 = N 0 n = 0 T ; x 1 n Γ 1 ( U + n α I ) w n s n + 1 × h = 0 ( 1 ) h R s v h .
Employing the relation through the use of the inverse Laplace Transform
L 1 s v , t = t v 1 Γ ( v ) , ( R e ( v ) > 0 ) ,
we have,
N ( t ) I = L 1 { N * ( s ) , t } = N 0 n = 0 [ T ; x ] β n Γ 1 U + n α I × ( w t ) n × h = 0 ( 1 ) h ( R t ) ν h Γ ( ν h + β + 1 ) = N 0 n = 0 [ T ; x ] β n Γ 1 ( U + n α I ) × ( w t ) n × E ν , n + 1 R v t v .
The proof is now completed. □
Corollary 3. 
Consider T, U and R holds (11), where T, U, and R are positive stable matrices in C h × h . Then, for Re ( v ) > 0 and t C , the generalized FK matrix equation of matrix ML functions satisfies the following equation
N ( t ) I N 0 E α , U β , T ( w t ) = R ν D t v 0 N ( t )
has the solution
N ( t ) I = N 0 n = 0 ( T ) β n Γ 1 ( U + α n I ) × ( w t ) n E ν , β + 1 ( R ν t ν ) .

7. Conclusions

In conclusion, this paper has explored an incomplete Mittag-Leffler (ML) function where matrix arguments were introduced. Some properties of these functions, such as functional relations, convergent, integral formulas, and integral representations were investigated and the properties of the Riemann–Liouville fractional integrals and derivatives related to the incomplete matrix ML function were studied. Additionally, some interesting special cases of this work were highlighted. Also, we established a solution to the kinetic equations involving the incomplete matrix ML function. Ultimately, these theoretical advancements find practical applications, particularly in theorem of the incomplete matrix ML functions across diverse types, such as the k-incomplete matrix ML function. This research opens avenues for further exploration and development within this intricate field of study. Certain applications to other research subjects and investigation regarding other properties of these newly introduced functions are left to the authors and the interested researchers for future study.

Author Contributions

Conceptualization, A.B., S.H. and M.Z.; Formal analysis, A.B., S.H. and M.Z.; Funding acquisition, M.Z.; Investigation, A.B. and S.H.; Methodology, A.B., S.H. and M.Z.; Project administration, M.Z.; Validation, A.B. and S.H.; Writing original draft, A.B.; Writing the review and editing, A.B. and M.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by Deanship of Research and Graduate Studies at King Khalid University through Small Group Research Project under grant number RGP1/116/45.

Informed Consent Statement

Not applicable.

Data Availability Statement

No data are associated with this study.

Acknowledgments

The authors extend their appreciation to the Deanship of Research and Graduate Studies at King Khalid University for funding this work through Small Group Research Project under grant number RGP1/116/45.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Mittag-Leffler, G.M. Sur la nouvelle function Eα(x). C. R. Acad. Sci. Paris 1903, 137, 554–558. [Google Scholar]
  2. Abramowitz, M.; Stegun, I.A. (Eds.) Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables; Applied Mathematics Series 55, Ninth Printing; Reprint of the 1972 Edition; National Bureau of Standards: Washington, DC, USA, 1972; Dover Publications, Inc.: New York, NY, USA, 1992. [Google Scholar]
  3. Prabhakar, T.R. A singular integral equation with a generalized Mittag-Leffler function in the kernel. Yokohama Math. J. 1971, 19, 7–15. [Google Scholar]
  4. Shukla, A.K.; Prajapati, J.C. On a generalized Mittag-Leffler function and its properties. J. Math. Anal. Appl. 2007, 336, 797–811. [Google Scholar] [CrossRef]
  5. Kilbas, A.A.; Saigo, M.; Saxena, R.K. Generalized Mittag-Leffler function and generalized fractional calculus operators. Integral Transform. Spec. Funct. 2004, 15, 31–49. [Google Scholar] [CrossRef]
  6. Kilbas, A.A.; Saigo, M.; Saxena, R.K. Solution of Volterra integro-differential equations with generalized Mittag-Leffler function in the Kernels. J. Integral Equ. Appl. 2002, 14, 377–396. [Google Scholar] [CrossRef]
  7. Paneva-Konovska, J. The convergence of series in multi-index Mittag-Leffler functions. Integral Transform. Spec. Funct. 2012, 23, 207–221. [Google Scholar] [CrossRef]
  8. Parmar, R.K. A class of extended Mittag Leffler functions and their properties related to integral transforms and fractional calculus. Mathematics 2015, 3, 1069–1082. [Google Scholar] [CrossRef]
  9. Saxena, R.K.; Saigo, M. Certain properties of fractional calculus operators associated with generalized Mittag-Leffler function. Fract. Calc. Appl. Anal. 2005, 8, 141–154. [Google Scholar]
  10. Gorenflo, R.; Kilbas, A.; Mainardi, F.; Rogosin, S. Mittag- Leffler Function, Related Topics and Applications; Springer: Berlin/Heidelberg, Germany, 2014. [Google Scholar]
  11. Haubold, H.J.; Mathai, A.M.; Saxena, R.K. Mittag-Leffler function and their applications. J. Appl. Math. 2011, 2011, 1–51. [Google Scholar] [CrossRef]
  12. Haubold, H.J.; Mathai, A.M. The fractional kinetic equation and thermonuclear functions. Astrophys. Space Sci. 2000, 327, 53–63. [Google Scholar] [CrossRef]
  13. Bakhet, A.; Jiao, Y.; He, F.L. On the Wright hypergeometric matrix functions and their fractional calculus. Integral Transform. Spec. Funct. 2019, 30, 138–156. [Google Scholar] [CrossRef]
  14. Bakhet, A.; Zayed, M. Incomplete exponential type of R-matrix functions and their properties. AIMS Math. 2023, 8, 26081–26095. [Google Scholar] [CrossRef]
  15. Jódar, L.; Cortés, J.C. On the hypergeometric matrix Function. J. Comput. Appl. Math. 1998, 99, 205–217. [Google Scholar] [CrossRef]
  16. Jódar, L.; Cortés, J.C. Some properties oF Gamma and Beta matrix Functions. Appl. Math. Lett. 1998, 11, 89–93. [Google Scholar] [CrossRef]
  17. Goyal, R.; Agarwal, P.; Oros, I.G.; Jain, S. Extended Beta and Gamma Matrix Functions via 2-Parameter Mittag-LeFFler Matrix Function. Mathematics 2022, 10, 892. [Google Scholar] [CrossRef]
  18. Cuchta, T.; Grow, D.; Wintz, N. Discrete matrix hypergeometric functions. J. Math. Anal. Appl. 2023, 518, 126716. [Google Scholar] [CrossRef]
  19. Garrappa, R.; Popolizio, M. On the use of matrix functions for fractional partial differential equations. Math. Comput. Simul. 2011, 81, 1045–1056. [Google Scholar] [CrossRef]
  20. Popolizio, M. Numerical solution of multiterm fractional differential equations using the matrix Mittag-Leffler functions. Mathematics 2018, 1, 7. [Google Scholar] [CrossRef]
  21. Garrappa, R.; Popolizio, M. Computing the matrix Mittag-Leffler function with applications to fractional calculus. J. Sci. Comput. 2018, 77, 129–153. [Google Scholar] [CrossRef]
  22. Pal, A.; Jatav, V.K.; Shukla, A.K. Matrix analog of the four-parameter Mittag-Leffler function. Math. Meth. Appl. Sci. 2023, 46, 15094–15106. [Google Scholar] [CrossRef]
  23. Abdalla, M. Special matrix functions: Characteristics, achievements and future directions. Linear Multilinear Algebra 2020, 68, 1–28. [Google Scholar] [CrossRef]
  24. Abdalla, M. On the incomplete hypergeometric matrix functions. Ramanujan J. 2017, 43, 663–678. [Google Scholar] [CrossRef]
  25. Akel, M.; Bakhet, A.; Abdalla, M.; He, F. On degenerate gamma matrix functions and related functions. Linear Multilinear Algebra 2023, 71, 673–691. [Google Scholar] [CrossRef]
  26. Bakhet, A.; He, F.; Yu, M. On the matrix version of extended Bessel functions and its application to matrix differential equations. Linear Multilinear Algebra 2022, 70, 5661–5680. [Google Scholar] [CrossRef]
  27. Hidan, H.; Bakhet, A.; Abd-Elmageed, H.; Abdalla, M. Matrix-valued hypergeometric Appell-type polynomials. Electron. Res. Arch. 2022, 30, 2964–2980. [Google Scholar] [CrossRef]
  28. Zou, C.Y.; Bakhet, A.; He, F. On the Matrix Versions of Incomplete Extended Gamma and Beta Functions and Their Applications for the Incomplete Bessel Matrix Functions. Complexity 2021, 2021, 5586021. [Google Scholar] [CrossRef]
  29. Atangana, A. Fractional Operators with Constant and Variable Order with Application to Geo-Hydrology; Academic Press: Cambridge, MA, USA, 2018. [Google Scholar]
  30. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon & Breach: Yverdon, Switzerland, 1993. [Google Scholar]
  31. Saxena, R.K.; Kalla, S.L. On the solutions of certain fractional kinetic equations. Appl. Math. Comput. 2008, 199, 504–511. [Google Scholar] [CrossRef]
  32. Wiman, A. Aber den Fundamentalsatz in der Teorie der Funktionen Ea(x). Acta Math. 1905, 29, 191–201. [Google Scholar] [CrossRef]
  33. Diethelm, K. The Analysis of Fractional Differential Equations; Lecture Notes in Mathematics; Springer: Berlin/Heidelberg, Germany, 2010; Volume 2004. [Google Scholar]
  34. Habenom, H.; Suthar, D.L.; Gebeyehu, M. Application of Laplace Transform on Fractional Kinetic Equation Pertaining to the Generalized Galué Type Struve Function. Adv. Math. Phys. 2019, 2019, 5074039. [Google Scholar] [CrossRef]
  35. Suthar, L.D.; Purohit, S.D.; Araci, S. Solution of Fractional Kinetic Equations Associated with the p,q-Mathieu-Type Series. Discret. Dyn. Nat. Soc. 2020, 8645161. [Google Scholar] [CrossRef]
  36. Saichev, A.I.; Zaslavsky, G.M. Fractional kinetic equations: Solutions and applications. Chaos 1997, 7, 753–764. [Google Scholar] [CrossRef] [PubMed]
  37. Singh, G.; Agarwal, P.; Chand, M.; Jain, S. Certain fractional kinetic equations involving generalized k-Bessel function. Trans. A Razmadze Math. Inst. 2018, 172, 559–570. [Google Scholar] [CrossRef]
  38. Abdalla, M.; Akel, M. Contribution of Using Hadamard Fractional Integral Operator via Mellin Integral Transform for Solving Certain Fractional Kinetic Matrix Equation. Fractal Fract. 2022, 6, 305. [Google Scholar] [CrossRef]
  39. Hidan, M.; Akel, M.; Abd-Elmageed, H.; Abdalla, M. Solution of fractional Kinetic equations involving of extended (k, τ)- Gauss hypergeometric matrix functions. AIMS Math. 2022, 7, 14474–14491. [Google Scholar] [CrossRef]
  40. Nisar, K.S.; Purohit, S.D.; Mondal, S.R. Generalized fractional kinetic equations involving generalized Struve function of the first kind. J. King Saud. Univ. Sci. 2016, 28, 167–171. [Google Scholar] [CrossRef]
  41. Nisar, K.S.; Belgacem, F.B.M. Dynamic k-Struve Sumudu solutions for fractional kinetic equations. Adv. Differ. Equ. 2017, 2017, 340. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bakhet, A.; Hussain, S.; Zayed, M. On Fractional Operators Involving the Incomplete Mittag-Leffler Matrix Function and Its Applications. Symmetry 2024, 16, 963. https://doi.org/10.3390/sym16080963

AMA Style

Bakhet A, Hussain S, Zayed M. On Fractional Operators Involving the Incomplete Mittag-Leffler Matrix Function and Its Applications. Symmetry. 2024; 16(8):963. https://doi.org/10.3390/sym16080963

Chicago/Turabian Style

Bakhet, Ahmed, Shahid Hussain, and Mohra Zayed. 2024. "On Fractional Operators Involving the Incomplete Mittag-Leffler Matrix Function and Its Applications" Symmetry 16, no. 8: 963. https://doi.org/10.3390/sym16080963

APA Style

Bakhet, A., Hussain, S., & Zayed, M. (2024). On Fractional Operators Involving the Incomplete Mittag-Leffler Matrix Function and Its Applications. Symmetry, 16(8), 963. https://doi.org/10.3390/sym16080963

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop