Next Article in Journal
Non-Equilibrium Thermodynamic Foundations of the Origin of Life
Next Article in Special Issue
Nonlocal ψ-Hilfer Generalized Proportional Boundary Value Problems for Fractional Differential Equations and Inclusions
Previous Article in Journal / Special Issue
A Note on a Coupled System of Hilfer Fractional Differential Inclusions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generalized Fractional Integrals Involving Product of a Generalized Mittag–Leffler Function and Two H-Functions

1
Department of Mathematics, Poornima University, Jaipur 303905, India
2
Department of Mathematics, Poornima College of Engineering, Jaipur 302022, India
3
Department of Mathematics, Anand International College of Engineering, Jaipur 303012, India
4
Nonlinear Dynamics Research Center (NDRC), Ajman University, Ajman 346, United Arab Emirates
*
Author to whom correspondence should be addressed.
Foundations 2022, 2(1), 298-307; https://doi.org/10.3390/foundations2010021
Submission received: 14 December 2021 / Revised: 19 February 2022 / Accepted: 1 March 2022 / Published: 11 March 2022
(This article belongs to the Special Issue Recent Advances in Fractional Differential Equations and Inclusions)

Abstract

:
The objective of this research is to obtain some fractional integral formulas concerning products of the generalized Mittag–Leffler function and two H-functions. The resulting integral formulas are described in terms of the H-function of several variables. Moreover, we give some illustrative examples for the efficiency of the general approach of our results.

1. Introduction and Preliminaries

The Fractional calculus operators having various special functions have been used for modeling systems in turbulence and fluid dynamics, stochastic dynamic systems, thermonuclear fusion, image processing, nonlinear biological system, and quantum mechanics. For example, Beleanu et al. [1] described a fractional sub-equation method to solve fractional differential equations. Purohit and Kalla [2] studied the solutions of generalized fractional partial differential equations in quantum mechanics. Mathai et al. [3,4,5] examine the various application of H- function by using fractional calculus for the various problem of physics.
Especially the diverse applications of fractional calculus motivated us to establish, here two image formula for the product of two H- functions and generalized Mittag–Leffler functions involving left and right sided fractional operator of Saigo–Meada [6], the resulting integral formulas are described in terms of the H-function of several variables. Similarly in next two image formula established with the help of the product of two H-functions and generalized multi-index Mittag–Leffler functions involving left and right sided fractional operator of Saigo–Meada [6], the resulting integral formulas are described in terms of the H-function of several variables. By virtue of the unified nature of generalized Mittag–Leffler functions involved in our results, a large number of new and known results are shown to follow as special cases of our main results but here we presented two of them.
First of all, we recall Mittag–Leffler function and its generalizations.
In 1903, Mittag–Leffler [7] defined a function in term of a power series:
E ρ ( y ) = l = 0 y n Γ ( ρ l + 1 ) , ( ρ > 0 , y C ) .
Then in 1905, Wiman [8] has given a two-index generalization of this function as:
E ρ , χ ( y ) = l = 0 y n Γ ( ρ l + χ ) , ( ρ > 0 , χ > 0 , y C ) .
Further, Prabhakar [9] presented the generalizing series representation of (2) as:
E ρ , χ ϖ ( y ) = l = 0 ( ϖ ) l Γ ( ρ l + χ ) y n l ! , ( ( ρ ) > 0 , ( χ ) > 0 , ρ , χ , ϖ , y C ) .
In addition, Kilbas et al. [10] defined an extension of (3) as:
E ϖ [ ( ρ , χ ) n ; y ] = E ϖ [ ( ρ 1 , χ 1 ) , , ( ρ n , χ n ) ; y ] = l = 0 ( ϖ ) l i = 1 n Γ ( ρ i l + χ i ) y l l ! .
If χ i R ( χ i 0 ) , ρ i C ( i = 1 , , n ) and y C , then
(i)
if i = 0 n χ i > 0 , then the above extended Wright function is an entire function of y.
(ii)
if i = 0 n χ i > 0 and either | y | < i = 1 n | χ i | χ i or y = i = 1 n | χ i | χ i , i = 1 n ( ρ ) > ( ϖ ) + n / 2 , then the series in l = 0 ( ϖ ) l i = 1 n Γ ( ρ i l + χ i ) y l l ! is absolutely convergent.
Then, after some time, generalization of multi-index Mittag–Leffler function was studied by Saxena and Nishimoto [11] in 2010, and defined in the following manner:
E ϖ , r [ ( ρ , χ ) n ; y ] = E ( ρ i , χ i ) n ϖ , r [ y ] = l = 0 ( ϖ ) r l i = 1 n Γ ( ρ i l + χ i ) y l l ! , ( ( ρ i ) > 0 , ρ j , χ i , ϖ , r , y C , i = 1 , 2 , , n ; i = 1 m ρ i > max { 0 , ( r ) 1 } ) .
In the present paper, we use the generalized Mittag–Leffler function, defined by Saxena and Kalla [12] as:
E ( ρ i ) , τ ϖ i ( t χ 1 ( d c t ) δ 1 , , t χ m ( d c t ) δ m ) = E ( ρ 1 , , ρ m ) , τ ( ϖ 1 , , ϖ m ) t χ 1 ( d c t ) δ 1 , , t χ m ( d c t ) δ m = l 1 , , l m = 0 ( ϖ 1 ) l 1 ( ϖ m ) l m ( t χ 1 ( d c t ) δ 1 ) l 1 ( t χ m ( d c t ) δ m ) l m Γ ( τ + ρ 1 l 1 + + ρ m l m ) ( l 1 ) ! ( l m ) ! , ( ϖ i , ρ i , τ . y , t , c , d , δ i , χ i C , ( ρ i ) > 0 , i = 1 , , m ) .
Further, the generalized multi-index Mittag–Leffler function E ( ρ j , χ j ) n ( ϖ , c ) ; r ( t χ ( d c t ) δ ; p ) , defined as [13]:
E ( ρ j , χ j ) n ( ϖ , e ) ; r ( t τ ( d c t ) δ ; p ) = l = 0 B p ( ϖ + r l , e ϖ ) B ( ϖ , e ϖ ) ( e ) r l i = 1 n Γ ( ρ i l + χ i ) ( t τ ( d c t ) δ ) l l ! , ( ρ i , χ i , ϖ , r , t , c , d , e , δ C , r 0 , ( e ) > ( ϖ ) > 0 , ( χ i ) > 0 , ( i = 1 , 2 , , n ) ; ( i = 1 m ρ i ) > max { 0 , ( r ) 1 } ) .
Here, B p is extended Beta function, define as B p ( x , y ) = 0 1 t x 1 ( 1 t ) y 1 e p t ( t 1 ) , ( ( x ) > 0 , ( y ) > 0 , ( p ) > 0 .)
B is classical Beta function, see [14].
Further, we are recalling fractional integral operator of arbitrary order involving Appell function F 3 in the kernel defined and studied by Saigo and Maeda [6], as given below:
Let κ > 0 , y > 0 , ( κ ) > 0 , ζ , ζ , ϑ , ϑ , κ C , then
( I 0 , y ζ , ζ , ϑ , ϑ , κ f ) ( y ) = y ζ Γ ( κ ) 0 y ( y t ) κ 1 t ζ F 3 ζ , ζ , ϑ , ϑ ; κ ; 1 t y , 1 y t f ( t ) d t ,
( I y , ζ , ζ , ϑ , ϑ , κ f ) ( y ) = y ζ Γ ( κ ) y ( t y ) κ 1 t ζ F 3 ζ , ζ , ϑ , ϑ ; κ ; 1 y t , 1 t y f ( t ) d t ,
where
F 3 ( ζ , ζ , ϑ , ϑ ; κ ; z , ω ) = m , n = 0 ( ζ ) m ( ζ ) n ( ϑ ) m ( ϑ ) n ( κ ) m + n z m m ! ω n n ! ( max { | z | , | ω | } < 1 ) .
The Fox’s H-function is a generalized hypergeometric function, defined by means of the Mellin-Barnes type contour integral (see [15]),
H r , s m , n z a i , ζ i 1 , r b j , ϑ i 1 , s = H r , s m , n z a 1 , ζ 1 , , a r , ζ r b 1 , ϑ 1 , , b s , ϑ s = 1 2 π i L θ s k z k d k ,
where, for convenience,
θ s k = i = 1 m Γ b i ϑ i k i = 1 n Γ 1 a i ζ i k i = n + 1 r Γ a i ζ i k i = m + 1 s Γ 1 b i ϑ i k .
The H-function of several variables is defined as (see [15]):
H y 1 , , y k = H r , s : r 1 , s 1 ; ; r k , s k 0 , n : m 1 , n 1 ; ; m k , n k y 1 . y k a i ; ζ i , ζ i ( k ) 1 , r : ( c i , κ i ) 1 , r 1 ; ; c i ( k ) , κ i ( k ) 1 , r k b i ; ϑ i , , ϑ i ( k ) 1 , s : ( d i , δ i ) 1 , s 1 ; ; d i ( k ) , δ i ( k ) 1 , s k = 1 ( 2 π i ) k L 1 L k ϕ ( ω 1 , , ω k ) θ 1 ( ω 1 ) θ k ( ω k ) y 1 ω 1 , , y k ω k d ω 1 d ω k ,
here
ϕ ( ω 1 , , ω k ) = i = 1 n Γ ( 1 a i + j = 1 k ζ i ( j ) ω j ) i = n + 1 r Γ ( a i j = 1 k ζ i ( j ) ω j ) i = 1 s Γ ( 1 b i + j = 1 k ϑ i ( j ) ω j ) ,
and
θ i ( ω i ) = i = 1 n j Γ ( 1 c i ( j ) + κ i ( j ) ω i ) i = 1 m j Γ ( d i ( j ) δ i ( j ) ω j ) i = n j + 1 r j Γ ( c i ( j ) κ i ( j ) ω j ) i = m j + 1 s j Γ ( 1 d i ( j ) + δ j ( i ) ω j ) ,
here for all j { 1 , , k } .
Here for convenience a i ; ζ i , , ζ i ( k ) 1 , r abbreviate a 1 ; ζ 1 , , ζ 1 ( k ) , , a r ; ζ r , , ζ r ( k ) . While c i ( j ) , κ i ( j ) 1 , r j abbreviate the array of r j pairs of parameters: c 1 ( j ) , κ 1 ( j ) , , c r j ( j ) , κ r j ( j ) , j { 1 , , k } , and so on.
Suppose, as usual, that the parameters: a i , ( i = 1 , , r ) ; c i ( j ) , ( i = 1 , , r j )   b i , ( i = 1 , , s ) ; d i ( j ) , ( i = 1 , , s j ) here for all j { 1 , , k } , are complex number and the associated coefficients ζ i ( j ) , ( i = 1 , , r ) ; κ i ( j ) , ( i = 1 , , r j ) , ϑ i ( j ) , ( i = 1 , , s ) ; δ i ( j ) , ( i = 1 , , s j ) , j { 1 , , k } , are positive real numbers such that:
Λ j : = i = 1 r ζ i ( j ) i = 1 s ϑ i ( j ) + i = 1 r j κ i ( j ) i = 1 s j δ i ( j ) 0 ,
j : i = n + 1 r ζ i ( j ) i = 1 s ϑ i ( j ) + i = 1 n j κ i ( j ) i = n j + 1 r k κ i ( j ) + i = 1 m j δ i ( j ) i = m j + 1 s j δ i ( j ) > 0 ,
where the integer n , r , s , m j , n j , r k , s k are constrained by the inequalities 0 n r , s 0 , 1 m j s k and 0 n j r j , for all j { 1 , , k } and the equality in (17) holds true for suitably restricted values of the complex variables y 1 , , y k .
The multiple Maline-Barnes contour integral representing the multivariate H function converge absolutely, under the condition (17) when: | a r g ( y j ) | < 1 2 j π for all j { 1 , , k } . The point y j = 0 for j { 1 , , k } and varius exceptional parameter values being tacitly excluded.
Remark 1.
The special case of (13) when r = 1 reduces to the Fox’s H-function.
To prove our main results we required the following Lemma [6]:
Lemma 1.
Let ζ , ζ , ϑ , ϑ , κ C ; if R e ( κ ) > 0 and R e ( ρ ) > m a x [ 0 , R e ( ζ + ζ + ϑ κ ) , R e ( ζ ϑ ) ] , then
I 0 , y ζ , ζ , ϑ , ϑ , κ y ρ 1 = y ρ ζ ζ + κ 1 Γ ( ρ ) Γ ( ρ + κ ζ ζ ϑ ) Γ ( ρ + ϑ ζ ) Γ ( ρ + κ ζ ζ ) Γ ( ρ + κ ζ ϑ ) Γ ( ρ + ϑ ) .
Lemma 2.
Let ζ , ζ , ϑ , ϑ , κ C ; if R e ( κ ) > 0 and R e ( ρ ) < 1 + m i n [ ( ϑ ) , R e ( ζ + ζ + ϑ κ ) , R e ( ζ + ϑ κ ) ] , then
I y , ζ , ζ , ϑ , ϑ , κ y ρ 1 = y ρ ζ ζ + κ 1 Γ ( 1 + ζ + ζ κ ρ ) Γ ( 1 + ζ + ϑ κ ρ ) Γ ( 1 ϑ ρ ) Γ ( 1 ρ ) Γ ( 1 + ζ + ζ + ϑ κ ρ ) Γ ( 1 + ζ ϑ ρ ) .

2. Main Results

In the present section, we introduced fractional integrals involving the product of generalized Mittag–Leffler function and two Fox’s H-functions. These Integrals are defined in form of generalized multivariate H-function.
Theorem 1.
Let ζ , ζ , ϑ , ϑ , κ , μ , η , δ i , ϖ i , ρ i , χ i , v 1 , v 2 , z 1 , z 2 , c , d , y C and τ , α 1 , α 2 > 0 , ( ρ i ) > 0 , i = 1 , , m . Then the following relation holds:
( I 0 , y ζ , ζ , ϑ , ϑ , κ ( t μ 1 ( d c t ) η E ( ρ i ) , τ ϖ i t χ 1 ( d c t ) δ 1 , , t χ m ( d c t ) δ m H r 1 , s 1 m 1 , n 1 z 1 t α 1 ( d c t ) v 1 a i , A i 1 , r 1 b i , B i 1 , s 1 H r 2 , s 2 m 2 , n 2 z 2 t α 2 ( d c t ) v 2 c i , C i 1 , r 2 d i , D i 1 , s 2 ) ) ( y ) = d η y μ ζ ζ + κ 1 l 1 , , l m = 0 ( ϖ ) l 1 ( ϖ ) l m Γ ( τ + ρ 1 l 1 + + ρ m l m ) c d δ l l 1 y χ 1 l 1 1 l 1 ! c d δ m l m y χ m l m 1 l m ! H 4 , 4 : r 1 , s 1 ; r 2 , s 2 ; 0 , 1 0 , 4 : m 1 , n 1 ; m 2 , n 2 ; 1 , 0 z 1 y α 1 z 2 y α 2 c d y E 1 : a i , A i 1 , r 1 ; c i , C i 1 , r 2 ; E 2 : b i , B i 1 , s 1 ; d i , D i 1 , s 2 ; 0 , 1 ,
where E 1 = ( 1 η δ 1 l 1 δ m l m ; v 1 , v 2 , 1 ) , ( 1 μ χ 1 l 1 χ m l m ; α 1 , α 2 , 1 ) , ( 1 μ κ + ζ + ζ + ϑ χ 1 l 1 χ m l m ; α 1 , α 2 , 1 ) , ( 1 μ + ζ ϑ χ 1 l 1 χ m l m ; α 1 , α 2 , 1 ) and E 2 = ( 1 η δ 1 l 1 δ m l m ; v 1 , v 2 , 0 ) , ( 1 μ κ + ζ + ζ χ 1 l 1 χ m l m ; α 1 , α 2 , 1 ) , ( 1 μ κ + ζ + ϑ χ 1 l 1 χ m l m ; α 1 , α 2 , 1 ) , ( 1 μ ϑ χ 1 l 1 χ m l m ; α 1 , α 2 , 1 ) . and, satisfying the following condition
(i) 
| a r g z 1 | < 1 2 1 π , 1 > 0 , where
1 = i = 1 m 1 B i + i = 1 n 1 A i i = m 1 + 1 s 1 B i i = n 1 + 1 r 1 A i .
(ii) 
| a r g z 2 | < 1 2 2 π , 2 > 0 , where
2 = i = 1 m 2 D i + i = 1 n 2 C i i = m 2 + 1 s 2 D i i = n 2 + 1 r 2 C i .
(iii) 
| a b y | < 1 , also we have
R e ( μ ) + α 1 min 1 i m 1 R e ( b i B i ) + α 2 min 1 i m 2 R e ( d i D i ) > max [ 0 , R e ( ζ + ζ + ϑ κ ) , R e ( ζ ϑ ) ] ,
R e ( η ) + v 1 min 1 i m 1 R e ( b i B i ) + v 2 min 1 i m 2 R e ( d i D i ) > max [ 0 , R e ( ζ + ζ + ϑ κ ) , R e ( ζ ϑ ) ] .
Proof. 
Let 𝚥 1 be the left-hand side of (20), then by apply Equation (6), we get
𝚥 1 = ( I 0 , y ζ , ζ , ϑ , ϑ , κ ( t μ 1 ( d c t ) η l 1 , , l m = 0 ( ϖ 1 ) l 1 ( ϖ m ) l m ( t χ 1 ( d c t ) δ 1 ) l 1 ( t χ m ( d c t ) δ m ) l m Γ ( τ + ρ 1 l 1 + + ρ m l m ) ( l 1 ) ! ( l m ) ! H r 1 , s 1 m 1 , n 1 z 1 t α 1 ( d c t ) v 1 a i , A i 1 , r 1 b i , B i 1 , s 1 H r 2 , s 2 m 2 , n 2 z 2 t α 2 ( d c t ) v 2 c i , C i 1 , r 2 d i , D i 1 , s 2 ) ) ( y ) .
Using (11) to replace H functions in its Mellin-Barnes contour integral, we get
𝚥 1 = ( I 0 , y ζ , ζ , ϑ , ϑ , κ ( t μ 1 ( d c t ) η l 1 , , l m = 0 ( ϖ 1 ) l 1 ( ϖ m ) l m ( t χ 1 ( d c t ) δ 1 ) l 1 ( t χ m ( d c t ) δ m ) l m Γ ( τ + ρ 1 l 1 + + ρ m l m ) ( l 1 ) ! ( l m ) ! L 1 1 2 π i ϕ 1 ( ζ 1 ) ( z 1 ( d c t ) v 1 ) ζ 1 d ζ 1 L 2 1 2 π i ϕ 2 ( ζ 2 ) ( z 2 ( d c t ) v 2 ) ζ 2 d ζ 2 ) ) ( y ) .
By using the generalized binomial theorem, expanding the term ( d c t ) v 1 , we get:
( d c t ) v 1 = d v 1 k = 0 ( v 1 ) k k ! ( c t d ) k , | c t d | < 1 ,
similarly, we define the terms ( d c t ) v 2 and ( d c t ) η , and arranging the order of integral, we get
𝚥 1 = d η l 1 , , l m = 0 ( ϖ 1 ) l 1 ( ϖ m ) l m Γ ( τ + ρ 1 l 1 + + ρ m l m ) c d δ l l 1 1 l 1 ! c d δ m l m 1 l m ! 1 2 π i 3 L 1 ϕ 1 ( ζ 1 ) ( z 1 ) ζ 1 d v 1 ζ 1 d ζ 1 L 2 ϕ 2 ( ζ 2 ) ( z 2 ) ζ 2 d v 2 ζ 2 d ζ 2 L 3 Γ ( η + δ 1 l 1 δ m l m + v 1 ζ 1 + v 2 ζ 2 + ζ 3 ) Γ ( η + δ 1 l 1 δ m l m + v 1 ζ 1 + v 2 ζ 2 ) Γ ( 1 + ζ 3 ) ( a b ) ζ 3 d ζ 3 I 0 , y ζ , ζ , ϑ , ϑ , κ t χ 1 l 1 + + χ m l m + μ + α 1 ζ 1 + α 2 ζ 2 + ζ 3 1 ( y ) .
Now using the Equation (18), we get
𝚥 1 = d η l 1 , , l m = 0 ( ϖ 1 ) l 1 ( ϖ m ) l m Γ ( χ + ρ 1 l 1 + + ρ m l m ) c d δ l l 1 1 l 1 ! c d δ m l m 1 l m ! 1 2 π i 3 L 1 ϕ 1 ( ζ 1 ) ( z 1 ) ζ 1 d v 1 ζ 1 d ζ 1 L 2 ϕ 2 ( ζ 2 ) ( z 2 ) ζ 2 d v 2 ζ 2 d ζ 2 L 3 Γ ( η + δ 1 l 1 + + δ m l m + v 1 ζ 1 + v 2 ζ 2 + ζ 3 ) Γ ( η + δ 1 l 1 + + δ m l m + v 1 ζ 1 + v 2 ζ 2 ) Γ ( 1 + ζ 3 ) ( c d ) ζ 3 d ζ 3 Γ ( μ + χ 1 l 1 + + χ m l m + α 1 ζ 1 + α 2 ζ 2 ) Γ ( μ + χ 1 l 1 + + χ m l m + α 1 ζ 1 + α 2 ζ 2 + ζ 3 + κ ζ ζ ) Γ ( μ + χ 1 l 1 + + χ m l m + δ 1 ζ 1 + α 2 ζ 2 + ζ 3 + κ ζ ζ ϑ ) Γ ( μ + χ 1 l 1 + + χ m l m + α 1 ζ 1 + α 2 ζ 2 + ζ 3 κ ζ ϑ ) Γ ( μ + χ 1 l 1 + + χ m l m + δ 1 ζ 1 + α 2 ζ 2 + ζ 3 + ϑ ζ ) Γ ( μ + χ 1 l 1 + + χ m l m + α 1 ζ 1 + α 2 ζ 2 + ζ 3 + ϑ ) y ( μ + χ 1 l 1 + + χ m l m α 1 ζ 1 + α 2 ζ 2 + ζ 3 ζ ζ + κ 1 ) .
Hence the above equation can be written in form of R.H.S. of the Equation (20) by using the Equation (13) and we get our desired result. □
Theorem 2.
Let ζ , ζ , ϑ , ϑ , κ , μ , η , δ , v 1 , v 2 , z 1 , z 2 , c , d , ϖ i , χ i , ρ i , δ i C and τ , α 1 , α 2 > 0 , ( ρ i ) > 0 , i = 1 , , m . Then the relation holds:
( I y , ζ , ζ , ϑ , ϑ , κ ( t μ 1 ( d c t ) η E ( ρ i ) , τ ϖ i t χ 1 ( d c t ) δ 1 , , t χ m ( d c t ) δ m H r 1 , s 1 m 1 , n 1 z 1 t α 1 ( d c t ) v 1 a i , A i 1 , r 1 b i , B i 1 , s 1 H r 2 , s 2 m 2 , n 2 z 2 t α 2 ( d c t ) v 2 c i , C i 1 , r 2 d i , D i 1 , s 2 ) ) ( y ) = d η y μ ζ ζ + κ 1 l 1 , , l m = 0 ( ϖ ) l 1 ( ϖ ) l m Γ ( τ + ρ 1 l 1 + + ρ m l m ) c d δ l l 1 y χ i l 1 l 1 ! c d δ m l m y χ m l m 1 l m ! H 4 , 4 : r 1 , s 1 ; r 2 , s 2 ; 0 , 1 0 , 4 : m 1 , n 1 ; m 2 , n 2 ; 1 , 0 z 1 y α 1 z 2 y α 2 c d y F 1 : a i , A i 1 , r 1 ; c i , C i 1 , r 2 ; F 2 : b i , B i 1 , s 1 ; d i , D i 1 , s 2 ; 0 , 1 ,
where F 1 = ( 1 η δ 1 l 1 δ m l m ; v 1 , v 2 , 1 ) , ( 1 μ χ 1 l 1 χ m l m ; α 1 , α 2 , 1 ) , ( 1 + ζ + ζ + ϑ κ μ χ 1 l 1 χ m l m ; α 1 , α 2 , 1 ) , ( 1 + ζ ϑ μ χ 1 l 1 χ m l m ; α 1 , α 2 , 1 ) , and F 2 = ( 1 η δ 1 l 1 δ m l m ; v 1 , v 2 , 0 ) , ( 1 + ζ + ζ κ μ χ 1 l 1 χ m l m ; α 1 , α 2 , 1 ) , ( 1 + ζ + ϑ κ χ 1 l 1 χ m l m ; α 1 , α 2 , 1 ) , ( 1 ϑ μ χ 1 l 1 χ m l m ; α 1 , α 2 , 1 ) and the condition (i)–(iii) in (20) are also satisfied.
Proof. 
Let 𝚥 2 be the left-hand side of(25), using Equation (6) and (11), we get:
𝚥 2 = ( I y , ζ , ζ , ϑ , ϑ , κ ( t μ 1 ( d c t ) η l 1 , , l m = 0 ( ϖ 1 ) l 1 ( ϖ m ) l m ( t χ 1 ( d c t ) δ 1 ) l 1 ( t χ m ( d c t ) δ m ) l m Γ ( τ + ρ 1 l 1 + + ρ m l m ) ( l 1 ) ! ( l m ) ! L 1 1 2 π i ϕ 1 ( ζ 1 ) ( z 1 ( d c t ) v 1 ) ζ 1 d ζ 1 L 2 1 2 π i ϕ 2 ( ζ 2 ) ( z 2 ( d c t ) v 2 ) ζ 2 d ζ 2 ) ) ( y ) .
By using the generalized binomial theorem, expanding the term ( d c t ) v 1 , we get:
( d c t ) v 1 = d v 1 k = 0 ( v 1 ) k k ! ( c t d ) k , | c t d | < 1 ,
similarly, we define the terms ( d c t ) v 2 and ( d c t ) η , and arranging the order of integral, we get:
𝚥 2 = d η l 1 , , l m = 0 ( ϖ 1 ) l 1 ( ϖ m ) l m Γ ( τ + ρ 1 l 1 + + ρ m l m ) c d δ l l 1 1 l 1 ! c d δ m l m 1 l m ! 1 2 π i 3 L 1 ϕ 1 ( ζ 1 ) ( z 1 ) ζ 1 d v 1 ζ 1 d ζ 1 L 2 ϕ 2 ( ζ 2 ) ( z 2 ) ζ 2 d v 2 ζ 2 d ζ 2 L 3 Γ ( η + δ 1 l 1 + + δ m l m + v 1 ζ 1 + v 2 ζ 2 + ζ 3 ) Γ ( η + δ 1 l 1 + + δ m l m + v 1 ζ 1 + v 2 ζ 2 ) Γ ( 1 + ζ 3 ) ( c d ) ζ 3 d ζ 3 I y , ζ , ζ , ϑ , ϑ , κ t χ 1 l 1 + + χ m l m + μ + α 1 ζ 1 + α 2 ζ 2 + ζ 3 1 ( y ) .
Next by using the Equation (19), we get:
𝚥 2 = d η l 1 , , l m = 0 ( ϖ 1 ) l 1 ( ϖ m ) l m Γ ( τ + ρ 1 l 1 + + ρ m l m ) c d δ l l 1 1 l 1 ! c d δ m l m 1 l m ! 1 2 π i 3 L 1 ϕ 1 ( ζ 1 ) ( z 1 ) ζ 1 d v 1 ζ 1 d ζ 1 L 2 ϕ 2 ( ζ 2 ) ( z 2 ) ζ 2 d v 2 ζ 2 d ζ 2 L 3 Γ ( η + δ 1 l 1 + + δ m l m + v 1 ζ 1 + v 2 ζ 2 + ζ 3 ) Γ ( η + δ 1 l 1 + + δ m l m + v 1 ζ 1 + v 2 ζ 2 ) Γ ( 1 + ζ 3 ) ( c d ) ζ 3 d ζ 3 Γ ( 1 + ζ + ζ κ μ χ 1 l 1 χ m l m α 1 ζ 1 α 2 ζ 2 ζ 3 ) Γ ( 1 μ χ 1 l 1 χ m l m α 1 ζ 1 α 2 ζ 2 ζ 3 ) Γ ( 1 + ζ + ϑ κ μ χ 1 l 1 χ m l m α 1 ζ 1 α 2 ζ 2 ζ 3 ) Γ ( 1 + ζ + ζ + ϑ κ μ χ 1 l 1 χ m l m α 1 ζ 1 α 2 ζ 2 ζ 3 ) Γ ( 1 ϑ μ χ 1 l 1 χ m l m α 1 ζ 1 α 2 ζ 2 ζ 3 ) Γ ( 1 + ζ ϑ μ χ 1 l 1 χ m l m α 1 ζ 1 α 2 ζ 2 ζ 3 ) y ( μ + χ 1 l 1 + + χ m l m + α 1 ζ 1 + α 2 ζ 2 + ζ 3 ζ z e t a + κ 1 ) .
Hence, the above equation can be written in form of R.H.S of the Equation (25) by using Equation (13). Hence, we get our desired result. □
Theorem 3.
let ζ , ζ , ϑ , ϑ , κ , μ , η , δ , v 1 , v 2 , z 1 , z 2 , c , d , e , ϖ , χ i , ρ i C and χ , α 1 , α 2 > 0 , ( ρ i ) > 0 . Then the following result holds:
( I 0 , y ζ , ζ , ϑ , ϑ , κ ( t μ 1 ( d c t ) η E ( ρ i , χ i ) n ( ϖ , e ) : r ( t χ ( d c t ) δ ; p ) H r 1 , s 1 m 1 , n 1 z 1 t α 1 ( d c t ) v 1 a i , A i 1 , r 1 b i , B i 1 , s 1 H r 2 , s 2 m 2 , n 2 z 2 t α 2 ( d c t ) v 2 c i , C i 1 , r 2 d i , D i 1 , s 2 ) ) ( y ) = d η y μ ζ ζ + κ 1 e 2 p a , b = 0 L a L b β ( ϖ + r l + a + 1 , e ϖ + b + 1 ) β ( ϖ , e ϖ ) l = 0 ( e ) r l i = 1 n Γ ( χ i + ρ i l ) c d δ l y χ l 1 l ! H 4 , 4 : r 1 , s 1 ; r 2 , s 2 ; 0 , 1 0 , 4 : m 1 , n 1 ; m 2 , n 2 ; 1 , 0 z 1 y α 1 z 2 y α 2 c d y g 1 : a i , A i 1 , r 1 ; c i , C i 1 , r 2 ; g 2 : b i , B i 1 , s 1 ; d i , D i 1 , s 2 ; 0 , 1 ,
where g 1 = ( 1 η δ l ; v 1 , v 2 , 1 ) , ( 1 μ χ l ; α 1 , α 2 , 1 ) , ( 1 μ κ + ζ + ζ + ϑ χ l ; α 1 , α 2 , 1 ) , ( 1 μ + ζ ϑ χ l ; α 1 , α 2 , 1 ) and g 2 = ( 1 η δ l ; v 1 , v 2 , 0 ) , ( 1 μ κ + ζ + ζ χ l ; α 1 , α 2 , 1 ) , ( 1 μ κ + ζ + ϑ χ l ; α 1 , α 2 , 1 ) , ( 1 μ ϑ χ l ; α 1 , α 2 , 1 ) . and the condition (i)–(iii) in Theorem 1, are also satisfied.
Proof. 
The same argument as in the proof of Theorem 1 will establish the result in Theorem 3 by using the (7). So its proof details are omitted. □
Theorem 4.
let ζ , ζ , ϑ , ϑ , κ , μ , η , δ , v 1 , v 2 , z 1 , z 2 , c , d , e , ϖ , χ i , ρ i , δ , y C and χ , α 1 , α 2 > 0 , ( ρ i ) > 0 . Then the following result holds:
( I y , ζ , ζ , ϑ , ϑ , κ ( t μ 1 ( d c t ) η E ( ρ i , χ i ) n ( ϖ , e ) : r ( t χ ( d c t ) δ ; p ) H r 1 , s 1 m 1 , n 1 z 1 t α 1 ( d c t ) v 1 a i , A i 1 , r 1 b i , B i 1 , s 1 H r 2 , s 2 m 2 , n 2 z 2 t α 2 ( d c t ) v 2 c i , C i 1 , r 2 d i , D i 1 , s 2 ) ) ( y ) = d η y μ ζ ζ + κ 1 e 2 p a , b , l = 0 L a L b β ( ϖ + r l + a + 1 , e ϖ + b + 1 ) β ( ϖ , e ϖ ) l = 0 ( e ) r l i = 1 n Γ ( χ i + ρ i l ) 1 l ! c d δ l y χ l H 4 , 4 : r 1 , s 1 ; r 2 , s 2 ; 0 , 1 0 , 4 : m 1 , n 1 ; m 2 , n 2 ; 1 , 0 z 1 y α 1 z 2 y α 2 c d x h 1 : a i , A i 1 , r 1 ; c i , C i 1 , r 2 ; h 2 : b i , B i 1 , s 1 ; d i , D i 1 , s 2 ; 0 , 1 ,
where h 1 = ( 1 η δ l ; v 1 , v 2 , 1 ) , ( 1 μ χ l ; α 1 , α 2 , 1 ) , ( 1 + ζ + ζ + ϑ κ μ χ l ; α 1 , α 2 , 1 ) , ( 1 + ζ ϑ μ χ l ; α 1 , α 2 , 1 ) h 2 = ( 1 η δ l ; v 1 , v 2 , 0 ) , ( 1 + ζ + ζ κ μ χ l ; α 1 , α 2 , 1 ) , ( 1 + ζ + ϑ κ χ l ; α 1 , α 2 , 1 ) , ( 1 ϑ μ χ l ; α 1 , α 2 , 1 ) and the condition (i)–(iii) in Theorem 1 are also satisfied.
Proof. 
The same argument as in the proof of Theorem 2 will establish the result in Theorem 4 by using the (7). So, its proof details are omitted. □

3. Special Cases

Firstly, here, it is important to remark the fact that the Saigo-Maeda fractional integral and fractional derivative operators involved in Theorems 1–4 are unified ones in nature. Secondly, the product of Fox’s H-functions occurring in Theorems 1–4 can be suitably specialized to give a large number of useful functions, for example, the Bessel functions, Wright hypergeometric functions, and so on. Here, among a remarkably large number of possible special examples of the results in Theorems 1–4, we consider only the following two examples.
If we consider only single variable generalized Mittag–Leffler function E ( ρ ) , τ ϖ by putting i = 1 in (6) and Multivariate H-functions and then apply the Theorems 1 and 2, then we get our following special cases as described below:
Corollary 1.
Let ζ , ζ , ϑ , ϑ , κ , μ , η , δ , v 1 , v 2 , z 1 , z 2 , c , d , ϖ , χ , ρ , δ , y C and τ , α 1 , α 2 > 0 , ( ρ ) > 0 . Then the following result holds:
( I 0 , y ζ , ζ , ϑ , ϑ , κ ( t μ 1 ( d c t ) η E ( ρ ) , τ ϖ t χ ( d c t ) δ H r 1 , s 1 m 1 , n 1 z 1 t α 1 ( d c t ) v 1 a i , A i 1 , r 1 b i , B i 1 , s 1 H r 1 , s 2 m 2 , n 2 z 2 t α 2 ( d c t ) v 2 c i , C i 1 , r 2 d i , D i 1 , s 2 ) ) ( y ) = d η y μ ζ ζ + κ 1 l = 0 ( ϖ ) l Γ ( τ + ρ l ) c d δ l y χ l 1 l ! H 4 , 4 : r 1 , s 1 ; r 2 , s 2 ; 0 , 1 0 , 4 : m 1 , n 1 ; m 2 , n 2 ; 1 , 0 z 1 y α 1 z 2 y α 2 c d y E 1 : a i , A i 1 , r 1 ; c i , C i 1 , r 2 ; E 2 : b i , B i 1 , s 1 ; d i , D i 1 , s 2 ; 0 , 1 ,
where E 1 = ( 1 η δ l ; v 1 , v 2 , 1 ) , ( 1 μ χ l ; α 1 , α 2 , 1 ) , ( 1 μ κ + ζ + ζ + ϑ χ l ; α 1 , α 2 , 1 ) , ( 1 μ + ζ ϑ χ l ; α 1 , α 2 , 1 ) and E 2 = ( 1 η δ l ; v 1 , v 2 , 0 ) , ( 1 μ κ + ζ + ζ χ l 1 ; α 1 , α 2 , 1 ) , ( 1 μ κ + ζ + ϑ χ l ; α 1 , α 2 , 1 ) , ( 1 μ ϑ χ l ; α 1 , α 2 , 1 ) .
Proof. 
In the similar manner as the proof of Theorem 1, we can easily get our desired result. □
Corollary 2.
let ζ , ζ , ϑ , ϑ , κ , μ , η , δ , v 1 , v 2 , z 1 , z 2 , c , d , ϖ , χ , ρ , δ , y C and τ , α 1 , α 2 > 0 , R e ( ρ ) > 0 . Then the following result holds:
( I y , ζ , ζ , ϑ , ϑ , κ ( t μ 1 ( d c t ) η E ( ρ ) , τ ϖ t χ ( d c t ) δ H r 1 , s 1 m 1 , n 1 z 1 t α 1 ( d c t ) v 1 a i , A i 1 , r 1 b i , B i 1 , s 1 H r 2 , s 2 m 2 , n 2 z 2 t α 2 ( d c t ) v 2 c i , C i 1 , r 2 d i , D i 1 , s 2 ) ) ( y ) = d η y μ ζ ζ + κ 1 l = 0 ( ϖ ) l Γ ( χ + ρ l ) c d δ l y χ l 1 l ! H 4 , 4 : r 1 , s 1 ; r 2 , s 2 ; 0 , 1 0 , 4 : m 1 , n 1 ; m 2 , n 2 ; 1 , 0 z 1 y α 1 z 2 y α 2 c d x F 1 : a i , A i 1 , r 1 ; c i , C i 1 , r 2 ; F 2 : b i , B i 1 , s 1 ; d i , D i 1 , s 2 ; 0 , 1 ,
F 1 = ( 1 η δ l ; v 1 , v 2 , 1 ) , ( 1 μ χ l ; α 1 , α 2 , 1 ) , ( 1 + ζ + ζ + ϑ κ μ χ l ; α 1 , α 2 , 1 ) , ( 1 + ζ ϑ μ χ l ; α 1 , α 2 , 1 ) , F 2 = ( 1 η δ l ; v 1 , v 2 , 0 ) , ( 1 + ζ + ζ κ μ χ l ; α 1 , α 2 , 1 ) , ( 1 + ζ + ϑ κ χ l ; α 1 , α 2 , 1 ) , ( 1 ϑ μ χ l ; α 1 , α 2 , 1 ) and the condition (i)–(iii) in Theorem 1, are also satisfied.
Proof. 
In the similar manner as the proof of Theorem (2), we can easily get our desired result. □

4. Conclusions

The H- functions associated with fractional calculus have been recognized to play a fundamental role in the probability theory as well as in their applications, including non-Gaussian stochastic processes and phenomena of nonstandard (i.e., anomalous) relaxation and diffusion. Another hand the H function and the generalized Mittag–Leffler function reduce to hypergeometric function and polynomials, so it becomes more important from the application viewpoint. Therefore, fractional calculus formulae involving hypergeometric functions and polynomials play an important role in the theory of special function and mathematical physics [16,17,18]. Finally, we conclude that our results presented in this paper are new and important from an application point of view. In the future, we are also trying to find some basic applications and examples of those results presented here to different research regions.

Author Contributions

All authors contributed equally. All authors have read and agreed to the published version of the manuscript.

Funding

Shilpi Jain is very thankful to SERB (project number: MTR/2017/000194) for providing the necessary facility.

Conflicts of Interest

All authors are confirm that they have no Conflicts of Interest.

References

  1. Baleanu, D.; Diethelm, K.; Scalas, E.; Trujillo, J.J. Fractional Calculus Models and Numerical Methods; Series on Complexity, Nonlinearity and Chaos; World Scientific: Boston, MA, USA, 2012. [Google Scholar]
  2. Purohit, S.D.; Kalla, S.L. On fractional partial differential equations related to quantum mechanics. J. Phys. AMath. Theor. 2011, 44, 45–202. [Google Scholar] [CrossRef]
  3. Mathai, A.M.; Saxena, R.K.; Haubold, H.J. The H-Function: Theorey and Applications; Springer: New York, NY, USA, 2010. [Google Scholar]
  4. Mathai, A.M.; Saxena, R.K. Distribution of a product and the structural setup of densities. Ann. Math. Statist. 1969, 4, 1439–1448. [Google Scholar] [CrossRef]
  5. Agarwal, P. Fractional integration of the product of two H-functions and a general class of polynomials. Asian J. Appl. Sci. 2012, 5, 144–153. [Google Scholar] [CrossRef] [Green Version]
  6. Saigo, M.; Maeda, N. More Generalization of Fractional Calculus. In Transform Methods and Special Functions, Varna’96; Bulgarian Academy of Sciences: Varna, Bulgaria, 1996; pp. 386–400. [Google Scholar]
  7. Mittag–Leffler, G.M. Sur la nouvelle fonction Eξ(x). C. R. Acad. Sci. 1903, 137, 554–558. [Google Scholar]
  8. Wiman, A. U ber der Fundamentalsatz in der Theorie der Funktionen Eα(x). Acta Math. 1905, 29, 191–201. [Google Scholar] [CrossRef]
  9. Prabhakar, T.R. A singular integral equation with a generalized Mittag–Leffler function inthe kernel. Yokohama Math. J. 1971, 19, 7–15. [Google Scholar]
  10. Kilbas, A.A.; Saigo, M.; Trujullo, J.J. On the generalized Wright function. Frac. Calc. Appl. Anal. 2002, 5, 437–466. [Google Scholar]
  11. Saxena, R.K.; Nishimoto, K. N-fractional calculus of generalized Mittag–Leffler functions. J. Fract. Calc. 2010, 37, 43–52. [Google Scholar]
  12. Saxena, R.K.; Kalla, S.L.; Saxena, R. Multivariate analogue of generalized Mittag–Leffler function. Integral Transform. Spec. Funct. 2011, 22, 533–548. [Google Scholar] [CrossRef]
  13. Agarwal, P.; Suthar, D.; Jain, S.; Momani, S. The extended multi-index Mittag–Leffler functions and their properties connected with fractional calculus and integral transforms. Thai J. Math. 2021, 19, 4625. [Google Scholar]
  14. Choudhary, M.A.; Qadir, A.; Srivastva, H.M.; Paris, R.B. Extended hypergeometric and confluent hypergeometric functions. Appl. Math. Comput. 2004, 159, 589–602. [Google Scholar]
  15. Srivasta, H.M.; Gupta, K.C.; Goyal, S.P. The H Function of One and Two Variables with Applicatons; South Asian Publication: New Delhi, India; Madras, India, 1982. [Google Scholar]
  16. Baleanu, D.; Agarwal, P.; Purohit, S.D. Certain fractional integral formulas involving the product of generlized Bessel function. Sci. World J. 2014, 2014, 567132. [Google Scholar]
  17. Baleanu, D.; Mustafa, O.G. On the global existence of solutions to a class of fractional differential equations. Comput. Math. Appl. 2010, 59, 1835–1841. [Google Scholar] [CrossRef] [Green Version]
  18. Baleanu, D.; Mustafa, O.G.; Agarwal, R.P. On the solution set for a class of sequential fractional differential equations. J. Phys. A Math. Theor. 2010, 43, 385209. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Singh, P.; Jain, S.; Agarwal, P. Generalized Fractional Integrals Involving Product of a Generalized Mittag–Leffler Function and Two H-Functions. Foundations 2022, 2, 298-307. https://doi.org/10.3390/foundations2010021

AMA Style

Singh P, Jain S, Agarwal P. Generalized Fractional Integrals Involving Product of a Generalized Mittag–Leffler Function and Two H-Functions. Foundations. 2022; 2(1):298-307. https://doi.org/10.3390/foundations2010021

Chicago/Turabian Style

Singh, Prakash, Shilpi Jain, and Praveen Agarwal. 2022. "Generalized Fractional Integrals Involving Product of a Generalized Mittag–Leffler Function and Two H-Functions" Foundations 2, no. 1: 298-307. https://doi.org/10.3390/foundations2010021

Article Metrics

Back to TopTop