Next Article in Journal
The Black Hole Universe, Part I
Next Article in Special Issue
The Fractional Hilbert Transform of Generalized Functions
Previous Article in Journal
Distributed Multiple High-Speed Trains Consensus Control Based on Event-Triggered Mechanism
Previous Article in Special Issue
Numerical Investigation of Nonlinear Shock Wave Equations with Fractional Order in Propagating Disturbance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Composition Structures of Certain Fractional Integral Operators

by
Min-Jie Luo
1,*,† and
Ravinder Krishna Raina
2,†,‡
1
Department of Mathematics, College of Science, Donghua University, Shanghai 201620, China
2
Department of Mathematics, College of Technology & Engineering, Maharana Pratap University of Agriculture and Technology, Udaipur 313001, India
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Current address: 10/11, Ganpati Vihar, Opposite Sector 5, Udaipur 313002, India.
Symmetry 2022, 14(9), 1845; https://doi.org/10.3390/sym14091845
Submission received: 9 August 2022 / Revised: 28 August 2022 / Accepted: 2 September 2022 / Published: 5 September 2022
(This article belongs to the Special Issue Applied Mathematics and Fractional Calculus II)

Abstract

:
This paper investigates the composition structures of certain fractional integral operators whose kernels are certain types of generalized hypergeometric functions. It is shown how composition formulas of these operators can be closely related to the various Erdélyi-type hypergeometric integrals. We also derive a derivative formula for the fractional integral operator and some applications of the operator are considered for a certain Volterra-type integral equation, which provide two generalizations to Khudozhnikov’s integral equation (see below). Some specific relationships, examples, and some future research problems are also discussed.

1. Introduction

In 1978, Saigo [1] introduced his widely used fractional integral operators I α , β , η and J α , β , η (see Equations (16) and (17) below). Saigo’s operators involve the Gauss hypergeometric functions 2 F 1 as kernels and possess many properties (see, for example, Refs. [1,2,3,4,5]). Over the past few decades, Saigo’s operators have been applied in various branches of mathematics, especially in the Geometric Function Theory (see Refs. [6,7,8]). The symmetry of parameters of various hypergeometric functions injects more choice and flexibility into the theory of Generalized Fractional Calculus.
A natural question that arises is: Can an operator involving a generalized hypergeometric function p F q as kernel have such properties as Saigo’s operators? In this direction, some efforts have been made by some authors to find particular forms of operators. In 1987, Goyal and Jain [9] introduced two fractional integral operators I α h and K β λ , which involve the generalized hypergeometric functions p F q as kernels. Later, Goyal et al. [10,11] introduced two more general fractional integral operators involving the generalized hypergeometric function p F q and Srivastava’s polynomial S n m .
Although very general in form, the properties of the operators I α h and K β λ introduced by Goyal et al. are far less succinct than those of Saigo’s operators. For Saigo’s operators I α , β , η and J α , β , η , we have the following useful properties (see Refs. [12,13]):
I α , β , η x λ = Γ ( λ ) Γ ( λ β + η + 1 ) Γ ( λ β + 1 ) Γ ( λ + α + η + 1 ) x λ β
( α ) > 0 , ( λ ) > max 0 , ( β η ) 1
and
J α , β , η x λ = Γ ( β λ ) Γ ( η λ ) Γ ( λ ) Γ ( α + β + η λ ) x λ β
( α ) > 0 , ( λ ) < max ( β ) , ( η ) .
Under certain conditions, we also have the following composition properties (see Ref. [1], p. 140, Equations (2.22) and (2.23), see also Ref. [3]):
I α , β , η I γ , δ , α + η f = I α + γ , β + δ , η f ,
I α , β , η I γ , δ , η β γ δ f = I α + γ , β + δ , η γ δ f ,
J γ , δ , α + η J α , β , η f = J α + γ , β + δ , η f
and
J γ , δ , η β γ δ J α , β , η f = J α + γ , β + δ , η γ δ f .
However, it seems rather difficult to find properties for the operators I α h and K β λ similar to those given above by (1)–(6). Moreover, it is still unknown whether the corresponding generalized fractional derivatives of the forms (see Ref. [3], Equations (3.2) and (3.4))
I α , β , η f = d n d x n I α + n , β n , η n f and J α , β , η f = 1 n d n d x n J α + n , β n , η n f
can be defined for the operators I α h and K β λ .
Very recently, the authors [14] introduced two fractional integral operators I and J (see below Equations (12) and (13)) whose kernels involve a very special class of generalized hypergeometric function. The authors have to some extent overcome the limitations of the operators I α h and K β λ and obtained results similar to (1) and (2). Subsequently, some further results and applications related to I and J were discovered in the papers [15,16].
The aim of the present paper is to first establish for the operators I and J some results relating to the composition structures of the defined operators analogous to Formulas (3)–(7). We also consider defining the corresponding fractional derivative operators of these operators I and J . Finally, we shall consider some connections of our work with Khudozhnikov’s work [17] on Volterra-type integral equations.

2. Preliminaries

In this paper, the symbols N , R + , and C denote the set of natural, positive real, and complex numbers, respectively. The Pochhammer symbol a k is defined by
a k : = Γ a + k Γ a = 1 ( k = 0 ; a C \ { 0 } ) a a + 1 a + k 1 ( k N ; a C ) .
In addition, we shall use the convention of writing the finite sequence of parameters a 1 , , a p by ( a p ) and the product of p Pochhammer symbols by ( ( a p ) ) k ( a 1 ) k ( a p ) k , where an empty product p = 0 is treated as unity.
We are particularly interested in the generalized hypergeometric function r + p F r + q of the form
r + p F r + q ( a p ) , ( f r + m r ) ( b q ) , ( f r ) ; z : = k = 0 ( ( a p ) ) k ( ( b q ) ) k ( ( f r + m r ) ) k ( ( f r ) ) k z k k ! ,
where m 1 , , m r N . The conditions of convergence of (8) follow easily from the usual definition of the generalized hypergeometric function; see Ref. [18], p. 62 and Ref. [19], p. 30. Several recent results concerning this particular type of generalized hypergeometric function have been obtained in Ref. [20] (see also Ref. [21]).
For convenience, we put
m : = m 1 + + m r ,
and let σ j ( 0 j m ) be determined by the generating relation
j = 1 r ( f j + x ) m j = j = 0 m σ m j x j .
Obviously, σ j ’s depend only on f j ( 1 j r ). Additionally, we define A k ( 0 k m ) by
A k = j = k m j k σ m j , A 0 = ( f 1 ) m 1 ( f r ) m r , A m = 1 ,
where the notation j k denotes the Stirling number of the second kind.
Definition 1
([14], p. 423, Definition 1.1). Let x , h , ν R + , δ , a , b , f 1 , , f r C and m 1 , , m r N . Also, let ( μ ) > 0 and φ be a suitable complex-valued function defined on R + . Then the fractional integral of the first kind of a function φ is defined by
I φ x I μ ; a , b : ( f r + m r ) h ; ν , δ : ( f r ) φ x : = ν x δ ν ( μ + h ) Γ ( μ ) 0 x x ν s ν μ 1 r + 2 F r + 1 a , b , μ , ( f r + m r ) ( f r ) ; 1 s ν x ν φ s s ν h + ν 1 d s ,
and the fractional integral of the second kind of a function φ x is defined by
J φ x J μ ; a , b : ( f r + m r ) h ; ν , δ : ( f r ) φ x : = ν x ν h + ν 1 Γ ( μ ) x s ν x ν μ 1 r + 2 F r + 1 a , b , μ , ( f r + m r ) ( f r ) ; 1 x ν s ν φ s s δ ν ( μ + h ) d s .
When r = 0 , we obtain
I μ ; a , b h ; ν , δ φ x = ν x δ ν ( μ + h ) Γ ( μ ) 0 x x ν s ν μ 1 2 F 1 a , b μ ; 1 s ν x ν φ s s ν h + ν 1 d s
and
J μ ; a , b h ; ν , δ φ x = ν x ν h + ν 1 Γ ( μ ) x s ν x ν μ 1 2 F 1 a , b μ ; 1 x ν s ν φ s s δ ν ( μ + h ) d s .
Some properties of the operators (12) and (13) have been presented in Refs. [14,16]. Further, the operators I μ ; a , b h ; ν , δ and J μ ; a , b h ; ν , δ have the following special cases:
(a)
For h = 0 , ν = 1 and δ = 0 in (14) and (15), we obtain
I μ ; a , b 0 ; 1 , 0 φ x = 2 I 0 + μ a , b φ x and I μ ; a , b 0 ; 1 , 0 φ x = 4 I μ a , b φ x ,
where 2 I 0 + μ a , b and 4 I μ a , b are two of the four operators introduced by Grinko and Kilbas [22].
(b)
When h = 0 , ν = 1 , δ = β , μ = α , a = α + β and b = η in (14) and (15), then we obtain Saigo’s fractional integral operators
I α , β , η φ x = I α ; α + β , η 0 ; 1 , β φ x = x β α Γ ( α ) 0 x x s α 1 2 F 1 α + β , η α ; 1 s x φ s d s ( α ) > 0
and
J α , β , η φ x = J α ; α + β , η 0 ; 1 , β φ x = 1 Γ ( α ) x s x α 1 2 F 1 α + β , η α ; 1 x s φ s s β α d s ( α ) > 0 .
(c)
When a = b = 0 , it is not difficult to observe that I μ ; a , b h ; ν , δ and J μ ; a , b h ; ν , δ contain the Erdélyi–Kober operators (see Ref. [19], p. 105 and Ref. [23], p. 322)
I + ; ν , h μ f x = I μ ; 0 , 0 h ; ν , 0 f x = ν x ν ( μ + h ) Γ ( μ ) 0 x x ν s ν μ 1 f s s ν h + ν 1 d s = 1 Γ ( μ ) 0 1 1 u μ 1 f x u 1 / ν u h d u ( μ ) > 0 , ν , h R +
and
I ; ν , h μ f x = J μ ; 0 , 0 h 1 + 1 / ν ; ν , 0 f x = ν x ν h Γ ( μ ) x s ν x ν μ 1 f s s ν ( 1 μ h ) 1 d s = 1 Γ ( μ ) 1 u 1 μ 1 f x u 1 / ν u μ h d u ( μ ) > 0 , ν , h R +
as special cases. The operators obtained by letting ν = 1 in (18) and (19) are usually denoted by I η , α + and K η , α , respectively (see Ref. [19], p. 106).
The operators defined above by (12) and (13) were previously studied in Refs. [14,16] in the space X c p ( c R , 1 p ) of those complex-valued Lebesgue measurable functions φ on R + for which φ X c p < , where
φ X c p : = 0 | u c φ u | p d u u 1 / p .
It follows at once that X 1 / p p = L p R + . For convenience, we define
c 1 t : = 1 + h + t ν and c 2 t : = c 1 δ 1 t ν .
The following lemma gives some useful properties of the operators I and J relating to the norm defined in (20).
Lemma 1.
Let φ X c p .
(i) 
If ( μ ) > 0 and c 1 c + min 0 , ( μ a b m ) > 0 , then the operator I is bounded from X c p into X c + δ , and
I φ X c + δ p C 1 φ X c p .
(ii) 
If ( μ ) > 0 and c 2 c + min 0 , ( μ a b m ) > 0 , then the operator J is bounded from X c p into X c + δ , and
J φ X c + δ p C 2 φ X c p .
(iii) 
If ( μ ) > 0 and c 1 c + min 0 , ( μ a b ) > 0 , then the operator I μ ; a , b h ; ν , δ is bounded from X c p into X c p , and
x δ I μ ; a , b h ; ν , δ φ X c p C 1 * φ X c p .
(iv) 
If ( μ ) > 0 and c 2 c + min 0 , ( μ a b ) > 0 , then the operator J μ ; a , b h ; ν , δ is bounded from X c p into X c p , and
x δ J μ ; a , b h ; ν , δ φ X c p C 2 * φ X c p .
(v) 
If ( μ ) > 0 and c 1 c > 0 , then the operator I + ; ν , h μ is bounded from X c p into X c p , and
I + ; ν , h μ φ X c p C 1 * * φ X c p .
(vi) 
If ( μ ) > 0 and ν h + c > 0 , then the operator I ; ν , h μ is bounded from X c p into X c p , and
I ; ν , h μ φ X c p C 2 * * φ X c p .
Proof. 
The results (i) and (ii) are established in Ref. [14], p. 437, Theorem 3.1.
On the other hand, the results (iii) and (iv) are the corollaries of (i) and (ii) (see also Ref. [16], p. 614).
Finally, the results (v) and (vi) follow immediately from (iii) and (iv). These results are consistent with the classical ones. It may be noted that if we set c = 1 / p in (v) and (vi), then the operator I + ; ν , h μ is bounded in L p R + provided that ( μ ) > 0 and h > 1 + 1 / p ν , and the operator I ; ν , h μ is bounded in L p R + provided that ( μ ) > 0 and h > 1 / p ν (see Ref. [19], p. 107, Lemma 2.28 and Ref. [23], p. 323).    □
It should be particularly emphasized here that the operators I and J are quite different from the multiple Erdélyi–Kober fractional integral operators (see Ref. [4], p. 11, see also Refs. [24,25]), though some special cases of I and J when r = 0 (e.g., Saigo’s operators) can be expressed as multiple Erdélyi–Kober fractional integral operators. The cases that r = 0 are very special because Meijer’s G-function G 2 , 2 2 , 0 [ σ ] and 2 F 1 [ 1 σ ] have the following relationship (see [4], p. 18, Equation (1.1.18))
G 2 , 2 2 , 0 σ γ 1 + δ 1 , γ 2 + δ 2 γ 1 , γ 2 = σ γ 2 1 σ δ 1 + δ 2 1 Γ ( δ 1 + δ 2 ) 2 F 1 γ 2 + δ 2 γ 1 , δ 1 δ 1 + δ 2 ; 1 σ
for σ < 1 . However, there is no such relationship between G m , m m , 0 [ σ ] and r + 2 F r + 1 [ 1 σ ] . A slightly more general case than (21) will lead us to the Marichev–Saigo–Maeda fractional integral operators (see Refs. [26,27]), which are also very different from our operators I and J . In addition, the operators I and J cannot be regarded as special cases of G-transform studied in Ref. [28]. Since the kernels of I and J are not of Sonine’s type, they cannot be included in the theory developed very recently by Luchko (see Ref. [29]).

3. The Main Results

3.1. Composition Formulas

Theorem 1.
Assume that φ X c p . Let
λ 1 λ a m , λ 2 λ b m and p m λ a b m ,
where m is given by (9). Let ( ϑ m ) be the nonvanishing zeros of the parametric polynomial Q m t defined by
Q m t = k = 0 m 1 k A k λ 1 k λ 2 k t k a + k m k b + k m k · 3 F 2 k m , k + t , p m a + k , b + k ; 1 ,
where A k ( 0 k m ) is defined in (11). Then for ( γ ) > 0 , ( μ ) > 1 / p > 0 ,
h + min 0 , ( γ + μ a b m ) > γ + p m + ( ρ c ) / ν
and h + 1 + min 0 , ( μ λ p m ) > ( c + ρ ) / ν , we have
I μ ; λ 1 , λ 2 : ( f r + m r ) h ; ν , δ : ( f r ) I γ ; p m , λ μ h γ p m ρ / ν ; ν , ρ φ x = I γ + μ ; a , b : ( ϑ m + 1 ) h γ p m ρ / ν ; ν , δ + ρ : ( ϑ m ) φ x
where I μ ; a , b : ( f r + m r ) h ; ν , δ : ( f r ) and I μ ; a , b h ; ν , δ are defined by (12) and (14), respectively.
Proof. 
Denote the left-hand side of (24) by Φ x . Then by interchanging the order of integration, we obtain
Φ x = ν x δ ν ( μ + h ) Γ ( μ ) 0 x x ν s ν μ 1 r + 2 F r + 1 λ 1 , λ 2 , μ , ( f r + m r ) ( f r ) ; 1 s ν x ν s ν h + ν 1 · ν s ν ( h p m ) Γ ( γ ) 0 s s ν t ν γ 1 2 F 1 p m , λ μ γ ; 1 t ν s ν φ t t ν ( h γ p m ) ρ + ν 1 d t d s = ν 2 x δ ν ( μ + h ) Γ ( μ ) Γ ( γ ) 0 x φ t t ν ( h γ p m ) ρ + ν 1 Δ 1 t d t ,
where
Δ 1 t : = t x s ν h + ν 1 ν h + ν p m x ν s ν μ 1 s ν t ν γ 1 · r + 2 F r + 1 λ 1 , λ 2 , μ , ( f r + m r ) ( f r ) ; 1 s ν x ν 2 F 1 p m , λ μ γ ; 1 t ν s ν d s .
We shall tackle Equation (24) and leave the verification of the validity of interchanging the order of integration in (25) at the end of the proof.
Letting s ν = x ν u x ν t ν in (26), we have
Δ 1 t = 1 ν x ν p m x ν t ν μ + γ 1 0 1 u μ 1 1 u γ 1 1 1 t ν x ν u p m · r + 2 F r + 1 λ 1 , λ 2 , μ , ( f r + m r ) ( f r ) ; 1 t ν x ν u 2 F 1 p m , λ μ γ ; ( 1 u ) ( 1 t ν / x ν ) 1 u ( 1 t ν / x ν ) d u .
The right-hand side of (27) can be evaluated by using an Erdélyi-type integral established by Luo and Raina [21]. For ( γ ) > ( μ ) > 0 and z C [ 1 , ) , Luo and Raina proved that (Ref. [21], p. 482, Theorem 3.2)
m + 2 F m + 1 a , b , ϑ m + 1 γ , ϑ m ; z = Γ ( γ ) Γ ( μ ) Γ ( γ μ ) 0 1 t μ 1 1 t γ μ 1 1 t z p m · r + 2 F r + 1 λ 1 , λ 2 , ( f r + m r ) μ , ( f r ) ; z t 2 F 1 p m , λ μ γ μ ; 1 t z 1 t z d t ,
where λ 1 , λ 2 and p m are given by (22) and ϑ m are the nonvanishing zeros of the parametric polynomial defined in (23). We note that the parametric polynomial is independent of parameter γ , and thus we may replace γ by γ + μ (without changing the values of λ 1 , λ 2 , p m and Q m t ) in (28) to get
m + 2 F m + 1 a , b , ϑ m + 1 γ + μ , ϑ m ; z = Γ ( γ + μ ) Γ ( μ ) Γ ( γ ) 0 1 t μ 1 1 t γ 1 1 t z p m · r + 2 F r + 1 λ 1 , λ 2 , ( f r + m r ) μ , ( f r ) ; z t 2 F 1 p m , λ μ γ ; 1 t z 1 t z d t ,
where min { ( γ ) , ( μ ) } > 0 .
Using the Erdélyi-type integral (29) in (27), we obtain
Δ 1 u = Γ ( μ ) Γ ( γ ) Γ ( γ + μ ) x ν p m x ν t ν μ + γ 1 m + 2 F m + 1 a , b , γ + μ , ϑ m + 1 ϑ m ; 1 t ν x ν .
Finally, substituting (30) into (25), we get
Φ x = ν x δ ρ Γ ( γ + μ ) x ν μ + γ + ( h γ p m ρ / ν ) 0 x x ν t ν μ + γ 1 · m + 2 F m + 1 a , b , γ + μ , ϑ m + 1 ϑ m ; 1 t ν x ν φ t t ν ( h γ p m ρ / ν ) + ν 1 d t = I γ + μ ; a , b : ( ϑ m + 1 ) h γ p m ρ / ν ; ν , δ + ρ : ( ϑ m ) φ x ,
which is the desired right-hand side of (24).
Now, we validate the interchanging of the integration. It is sufficient to show that
I = 0 x x ν s ν ( μ ) 1 r + 2 F r + 1 λ 1 , λ 2 , μ , ( f r + m r ) ( f r ) ; 1 s ν x ν s ν + ( ν p m ) 1 Δ 2 s d s < ,
where
Δ 2 s = 0 s s ν t ν ( γ ) 1 2 F 1 p m , λ μ γ ; 1 t ν s ν | φ t | t ν h ( ν γ + ν p m + ρ ) + ν 1 d t = 1 ν s ν h ( ν p m + ρ ) 0 1 1 u ( γ ) 1 2 F 1 p m , λ μ γ ; 1 u · φ s u 1 / ν u h ( γ + p m + ρ / ν ) 1 d u .
Note that (see Ref. [18], p. 63, Theorem 2.1.3 and [30], p. 387)
2 F 1 a , b c ; 1 z = O 1 , c a b > 0 ; O z ( c a b ) , c a b < 0 ; O log z , a + b = c ; O z ( c a b ) + O 1 , c a b = 0 , c a + b
as z 0 + , so for each s, we have
Δ 2 s D 1 · s ν h ( ν p m + ρ ) 0 1 1 u ( γ ) 1 · u h ( γ + p m + ρ / ν ) + min { 0 , ( γ + p m λ + μ ) } 1 φ s u 1 / ν d u ,
where D 1 is a positive number. In view of the definition of the Erdélyi–Kober operator (18), we have
Δ 2 s D 2 · s ν h ( ν p m + ρ ) F s ,
where D 2 : = D 1 Γ ( γ ) ( ( γ ) > 0 ) and
F s : = I + ; ν , h ( γ + p m + ρ / ν ) + min { 0 , ( γ + p m λ + μ ) } 1 ( γ ) | φ | s .
From Lemma 1, we have F X c p , since φ X c p and
h + min 0 , ( γ + p m λ + μ ) > γ + p m + ( ρ c ) / ν .
For the generalized hypergeometric function p + 1 F p [ z ] , we have (see, for example Ref. [31], p. 149)
p + 1 F p a 1 , , a p + 1 b 1 , , b p ; 1 z = O 1 , ( ψ p ) > 0 ; O z ( ψ p ) , ( ψ p ) < 0 ; O log z , ψ p = 0 ,
as z 0 + , where ψ p : = = 1 p b = 1 p + 1 a . Therefore, for each x R + , we find that
I D 2 D 3 x ν min { 0 , ( μ λ 1 λ 2 ) m } · 0 x x ν s ν ( μ ) 1 s ν min { 0 , ( μ λ 1 λ 2 ) m } + ν h + ν ( ρ ) F s d s s D 2 D 3 x ν min { 0 , ( μ λ 1 λ 2 ) m } F X c p · 0 x x ν s ν p ( μ ) p s p ν min { 0 , ( μ λ 1 λ 2 ) m } + p ν ( h + 1 ) p ( ρ ) p c 1 d s 1 / p D 2 D 3 v 1 / p x ν ( μ ρ / ν ) + ν h c F X c p · 0 1 1 u p ( μ ) p u p min { 0 , ( μ λ 1 λ 2 ) m } + p ( h + 1 ) p ( ρ ) / ν p c / ν 1 d s 1 / p < ,
where D 3 is a positive number, 1 / p + 1 / p = 1 , p ( μ ) p + 1 > 0 and
min 0 , ( μ λ 1 λ 2 ) m + h + 1 > ( c + ρ ) / ν .
Thus, Fubini’s theorem is applicable and the proof is complete.    □
Remark 1.
When r = 0 , we can set h = 0 , ν = 1 , μ = α , δ = λ b α and ρ = a + b λ γ in (24) to get
I α ; λ a , λ b 0 ; 1 , λ b α I γ ; a + b λ , λ α 0 ; 1 , a + b λ γ φ x = I α ; λ b , λ a 0 ; 1 , λ b α I γ ; a + b λ , λ α 0 ; 1 , a + b λ γ φ x = I γ + α ; a , b 0 ; 1 , a α γ φ x .
By comparing it with (16), we find that (33) is equivalent to the identity
I α , λ b α , a λ I γ , a + b λ γ , α λ φ x = I γ + α , a γ α , b φ x .
If we let further a = β + γ + δ + α , b = γ + δ η and λ = β η + γ + δ + α , then (34) reduces to (4).
Theorem 2.
Assume that φ X c p . Let λ 1 , λ 2 , and p m be defined in (22). Let ( ϑ m ) be the nonvanishing zeros of the parametric polynomial Q m t defined in (23). Then for ( γ ) > 0 , ( μ ) > 1 / p > 0 ,
h + 1 + ( ρ + δ ) / ν p m γ + min 0 , ( γ + p m λ + μ ) + ( c 1 ) / ν > 0
and 1 + h + ( 1 + c ) / ν + min 0 , ( μ λ p m ) + ( ρ + δ ) / ν > 0 , we have
J μ ; λ 1 , λ 2 : ( f r + m r ) h ; ν , δ : ( f r ) J γ ; p m , λ μ h γ p m + δ / ν ; ν , ρ φ x = J γ + μ ; a , b : ( ϑ m + 1 ) h p m γ ; ν , δ + ρ : ( ϑ m ) φ x ,
where J μ ; a , b : ( f r + m r ) h ; ν , δ : ( f r ) and J μ ; a , b h ; ν , δ are defined by (13) and (15), respectively.
Proof. 
Denote the left-hand side of (35) by Ψ s . Then, following a similar procedure as described in the proof of Theorem 1, we have
Ψ s = ν x ν h + ν 1 Γ ( μ ) x s ν x ν μ 1 r + 2 F r + 1 λ 1 , λ 2 , μ , ( f r + m r ) ( f r ) ; 1 x ν s ν s δ ν ( μ + h ) · { ν s ν h + ν 1 Γ ( γ ) s ν ( γ + p m ) + δ s t ν s ν γ 1 2 F 1 p m , λ μ γ ; 1 s ν t ν · φ t t ρ ν ( h p m ) δ d t } d s = ν 2 x ν h + ν 1 Γ ( μ ) Γ ( γ ) x φ t t ρ ν ( h p m ) δ Δ 3 t d t ,
where
Δ 3 t = x t s ν 1 ν ( μ + γ + p m ) s ν x ν μ 1 t ν s ν γ 1 · r + 2 F r + 1 λ 1 , λ 2 , μ , ( f r + m r ) ( f r ) ; 1 x ν s ν 2 F 1 p m , λ μ γ ; 1 s ν t ν d s .
Letting
s = t x t ν + ( x ν t ν ) u 1 / ν ,
so that
d s = 1 ν t x t ν x ν t ν + ( x ν t ν ) u 1 1 / ν d u and u = t ν ( x ν s ν ) s ν ( x ν t ν ) 0 , 1 ,
we have
Δ 3 t = 1 ν t ν x ν μ + γ 1 t ν μ x ν ( γ + p m ) 0 1 u μ 1 1 u γ 1 1 1 x ν t ν u p m · r + 2 F r + 1 λ 1 , λ 2 , μ , ( f r + m r ) ( f r ) ; 1 x ν t ν u 2 F 1 p m , λ μ γ ; ( 1 u ) ( 1 x ν / t ν ) 1 u ( 1 x ν / t ν ) d u .
The use of Erdélyi-type integral (29) gives
Δ 3 t = 1 ν Γ ( μ ) Γ ( γ ) Γ ( γ + μ ) t ν x ν μ + γ 1 t ν μ x ν ( γ + p m ) m + 2 F m + 1 a , b , ϑ m + 1 γ + μ , ϑ m ; 1 x ν t ν ,
and thus (36) becomes
Ψ s = ν x ν ( h γ p m ) + ν 1 Γ ( μ + γ ) x t ν x ν μ + γ 1 m + 2 F m + 1 a , b , ϑ m + 1 γ + μ , ϑ m ; 1 x ν t ν · φ t t ( ρ + δ ) ν ( μ + γ + h γ p m ) d t = J γ + μ ; a , b : ( ϑ m + 1 ) h p m γ ; ν , δ + ρ : ( ϑ m ) φ x ,
where ( ϑ m ) are the nonvanishing zeros of the parametric polynomial (23).
As in the proof of Theorem 1, we verify the validity of interchanging the order of integration by checking the finiteness of the integral
I = x s ν x ν ( μ ) 1 r + 2 F r + 1 λ 1 , λ 2 , μ , ( f r + m r ) ( f r ) ; 1 x ν s ν s ν 1 ( ν μ + ν γ + ν p m ) Δ 4 s d s ,
where
Δ 4 s = s t ν s ν ( γ ) 1 2 F 1 p m , λ μ γ ; 1 s ν t ν φ t t ( ρ + δ ν p m ) ν h d t = 1 ν s 1 ν ( ρ + δ ν γ ν p m ) ν h 1 u 1 ( γ ) 1 2 F 1 p m , λ μ γ ; 1 1 u · | φ s u 1 / ν | u 1 / ν 1 ( ρ / ν + δ / ν p m ) h d u .
Using (31) gives
Δ 4 s D 4 · s 1 ν ( ρ + δ ν γ ν p m ) ν h 1 u 1 ( γ ) 1 · φ s u 1 / ν u ( γ ) h + ( γ ) + 1 / ν 1 ( ρ / ν + δ / ν p m ) min { 0 , ( γ + p m λ + μ ) } d u ,
where D 4 is a positive number. Thus we have
Δ 4 s D 4 · s 1 ν ( ρ + δ ν γ ν p m ) ν h G s ,
where D 5 : = D 4 Γ ( γ ) ( ( γ ) > 0 ) and
G s : = I ; ν , h 1 / ν + 1 + ( ρ / ν + δ / ν p m γ ) + min { 0 , ( γ + p m λ + μ ) } ( γ ) | φ | s X c p .
Then from (32) we have
I D 5 D 6 x ν min { 0 , ( μ λ 1 λ 2 ) m } · x s ν x ν ( μ ) 1 s ν min { 0 , ( μ λ 1 λ 2 ) m } ν h ( ρ + δ + ν μ ) 1 G s d s s D 5 D 6 x ν min { 0 , ( μ λ 1 λ 2 ) m } G X c p · x s ν x ν p ( μ ) p s p ν min { 0 , ( μ λ 1 λ 2 ) m } p ν h p ( ρ + δ + ν μ ) p p c 1 d s 1 / p D 5 D 6 x ν ν h ( ρ + δ ) 1 c 1 / p G X c p · 1 u 1 p ( μ ) p u p min { 0 , ( μ λ 1 λ 2 ) m } p h p ( ( ρ + δ ) / ν + μ ) p ( 1 + c ) / ν 1 d s 1 / p < .
This completes the proof.   □
Remark 2.
When r = 0 , we can set h = 0 , ν = 1 , δ = γ + λ a b and ρ = a + b λ γ in (35) to get
J μ ; λ b , λ a 0 ; 1 , γ + λ a b I γ , a + b λ γ , μ λ φ x = J γ + μ ; a , b a + b λ γ ; 1 , 0 φ x .
Letting further a = μ + γ in (37), we have
I μ , λ μ b , γ + μ λ I γ , μ + b λ , μ λ φ x = J γ + μ ; γ + μ , b μ + b λ ; 1 , 0 φ x = I ; 1 , μ λ γ + μ φ x = K μ λ , γ + μ φ x .
Additionally, by putting b = β + λ μ in (38) and then letting λ = μ η in the resulting equation we get the following clearer form
I μ , β , γ + η I γ , β , η φ x = K η , γ + μ φ x ,
which is a special case of (5) when δ = β . It does not seem possible to deduce (5) by merely specializing the parameters in (35). Therefore, it should be interesting to find a composition formula from (35) which may include (5) or (6) as particular cases.
As depicted in Theorems 1 and 2, the study of the composition structure of the operators I and J rests heavily on the existence of a suitable Erdélyi-type integral, because we derive (24) and (35) from the Erdélyi-type integral (29). However, there may possibly be an alternative approach by which the Erdelyi-type integral may be obtained from a known composition structure [1] (see also Refs. [22,32]). Such an approach may be of special interest since our operators involve the generalized hypergeometric function r + 2 F r + 1 and the methodology may lead to some new results.

3.2. Derivative Formula

In this section we derive a derivative formula involving the fractional integral operator (12).
We introduce here some notations describing necessary concepts that would be used in the sequel. Let ( ξ m ) be the nonvanishing zeros of the parametric polynomial Q m t of degree m defined by
Q m t = j = 0 m σ m j k = 0 j j k ( b ) k ( t ) k ( μ b t ) m k ,
where the σ j ( 0 j m ) are determined by the generating relation (10). We define the parametric polynomial Q ˜ m t by
Q ˜ m t = j = 0 m σ ˜ m j k = 0 j j k ( μ b m ) k ( t ) k ( b + m n t ) m k ,
where σ ˜ j ( 0 j m ) are determined by the generating relation
j = 1 m ( ξ j + x ) = j = 0 m σ ˜ m j x j .
Theorem 3.
For ( μ ) > n ( n N ) , we have
n x n x δ + ν ( μ a + h ) I μ ; a , b : ( f r + m r ) h ; ν , δ : ( f r ) φ x = ν n x δ + ν ( μ n a + h ) I μ n ; a , b n : ( η m + 1 ) h ; ν , δ : ( η m ) φ x ,
where ( η m ) are the nonvanishing zeros of the parametric polynomial Q ˜ m t given by (40).
Proof. 
Using the Euler-type transformation due to Miller and Paris [20], p. 305, Theorem 3
r + 2 F r + 1 a , b , μ , ( f r + m r ) ( f r ) ; x = 1 x a m + 2 F m + 1 a , μ b m , μ , ( ξ m + 1 ) ( ξ m ) ; x x 1 ,
we have
x δ + ν ( μ a + h ) I μ ; a , b : ( f r + m r ) h ; ν , δ : ( f r ) φ x = ν Γ ( μ ) 0 x x ν s ν μ 1 m + 2 F m + 1 a , μ b m , μ , ( ξ m + 1 ) ( ξ m ) ; 1 x ν s ν φ s s ν ( h a ) + ν 1 d s ,
where ( ξ m ) are the nonvanishing zeros of the parametric polynomial Q m t defined by (39). By making use of the Leibniz integral rule, we obtain
x x δ + ν ( μ a + h ) I μ ; a , b : ( f r + m r ) h ; ν , δ : ( f r ) φ x = ν Γ ( μ ) x 0 x x ν s ν μ 1 m + 2 F m + 1 a , μ b m , μ , ( ξ m + 1 ) ( ξ m ) ; 1 x ν s ν φ s s ν ( h a ) + ν 1 d s = ν Γ ( μ ) 0 x x x ν s ν μ 1 m + 2 F m + 1 a , μ b m , μ , ( ξ m + 1 ) ( ξ m ) ; 1 s ν x ν · φ s s ν ( h a ) + ν 1 d s .
Taking into account the formula [33], p. 442, Equation (51)
n z n z σ 1 p F q ( a p ) ( b q 1 ) , σ ; z = ( σ n ) n z σ n 1 p F q ( a p ) ( b q 1 ) , σ n ; z ,
we have
x x δ + ν ( μ a + h ) I μ ; a , b : ( f r + m r ) h ; ν , δ : ( f r ) φ x = ν 2 Γ ( μ 1 ) 0 x x ν s ν ( μ 1 ) 1 · m + 2 F m + 1 a , μ b m , μ 1 , ( ξ m + 1 ) ( ξ m ) ; 1 x ν s ν φ s s ν ( h a ) + ν 1 d s .
Next, differentiating n times, we obtain
n x n x δ + ν ( μ a + h ) I μ ; a , b : ( f r + m r ) h ; ν , δ : ( f r ) φ x = ν 1 + n Γ ( μ n ) 0 x x ν s ν μ n 1 · m + 2 F m + 1 a , μ b m , μ n , ( ξ m + 1 ) ( ξ m ) ; 1 x ν s ν φ s s ν ( h a ) + ν 1 d s .
By applying the Euler-type transformation (43) again, we get
n x n x δ + ν ( μ a + h ) I μ ; a , b : ( f r + m r ) h ; ν , δ : ( f r ) φ x = ν 1 + n x ν a Γ ( μ n ) 0 x x ν s ν μ n 1 m + 2 F m + 1 a , b n , μ n , ( η m + 1 ) ( η m ) ; 1 s ν x ν φ s s ν h + ν 1 d s = ν n x δ + ν ( μ n a + h ) I μ n ; a , b n : ( η m + 1 ) h ; ν , δ : ( η m ) φ x ,
where the sequence of parameters ( η m ) are the nonvanishing zeros of the parametric polynomial Q ˜ m t of degree m given by (40). This completes the proof of (42).    □
Before proceeding further, we consider here a simple example.
Example 1.
When r = 1 and m = m 1 = 1 , f 1 = f and η 1 = η in (42), we get
n x n x δ + ν ( μ a + h ) I μ ; a , b : f + 1 h ; ν , δ : f φ x = ν n x δ + ν ( μ n a + h ) I μ n ; a , b n : η + 1 h ; ν , δ : η φ x ,
where η is the nonvanishing zero of the parametric polynomial
Q ˜ 1 t = j = 0 1 σ ˜ 1 j k = 0 j j k ( μ b 1 ) k ( t ) k ( b + 1 n t ) 1 k = σ ˜ 1 0 0 ( b + 1 n t ) + σ ˜ 0 1 0 ( b + 1 n t ) + σ ˜ 0 1 1 ( μ b 1 ) t = σ ˜ 1 ( b + 1 n ) + [ σ ˜ 0 ( μ b 1 ) σ ˜ 1 ] t .
Therefore, η can be expressed as
η = σ ˜ 1 ( b + 1 n ) σ ˜ 1 σ ˜ 0 ( μ b 1 ) .
It follows from (41) that σ ˜ 0 = 1 and σ ˜ 1 = ξ , where ξ is the nonvanishing zero of the parametric polynomial
Q 1 t = σ 1 ( μ b ) + [ σ 0 b σ 1 ] t .
From (10), we have σ 0 = 1 and σ 1 = f and thus ξ can be written as ξ = f ( μ b ) / ( f b ) . Hence,
η = f ( μ b ) ( b + 1 n ) f + b ( μ b 1 ) ,
wherein we note that η depends on n.
It may be observed that the Euler-type transformation (43) is used twice, so we need to be careful while finding special cases of Theorem 3.
(i)
By letting b = n ( n N ) in (42) and noting that m + 2 F m + 1 -function in (42) reduces to 1, we get
n x n x δ + ν ( μ a + h ) I μ ; a , n : ( f r + m r ) h ; ν , δ : ( f r ) φ x = ν 1 + n x ν a Γ ( μ n ) 0 x x ν s ν μ n 1 φ s s ν h + ν 1 d s = ν n x ν ( μ n a + h ) I + ; ν , h μ n φ x ,
where I + ; ν , h μ n denotes the Erdélyi–Kober type fractional integral defined by (18).
In fact, letting b = n changes the parametric polynomials Q m t and Q ˜ m t defined by (39) and (40), respectively. However, if the new polynomials, say Q m * t and Q ˜ m * t , also have nonvanishing zeros, denoted by ( ξ m * ) and ( η m * ) respectively, then (47) holds true. To illustrate here, let us set b = n in Example 1, then Q 1 t becomes Q 1 * t = f ( μ n ) + ( n f ) t with ξ * = f ( μ n ) / ( f n ) its nonvanishing zero and Q ˜ 1 t becomes Q ˜ 1 * t = ξ * + ( μ n 1 ξ * ) t . The nonvanishing zero of Q ˜ 1 * t is
η * = f ( μ n ) f + n ( μ n 1 ) ( f 0 , μ n ) .
Therefore, we obtain from (46) that
n x n x δ + ν ( μ a + h ) I μ ; a , n : f + 1 h ; ν , δ : f φ x = ν n x δ + ν ( μ n a + h ) I μ n ; a , 0 : η * + 1 h ; ν , δ : η * φ x = ν n x ν ( μ n a + h ) I + ; ν , h μ n φ x .
We also observe that the subsitution b = n may always reduce the right-hand side of (42) to a Erdélyi–Kober type integral.
(ii)
When r = 0 , then in view of (14) and (42), we simply obtain
n x n x δ + ν ( μ a + h ) I μ ; a , b h ; ν , δ φ x = ν n x δ + ν ( μ n a + h ) I μ n ; a , b n h ; ν , δ φ x .
Further, if h = 0 , ν = 1 , δ = β , a = α + β , b = η + n and μ = α + n in (48), we then have
n x n x n I α + n ; α + β , η + n 0 ; 1 , β φ x = I α ; α + β , η 0 ; 1 , β φ x .
In addition, in view of (16) and the relation
x n I α + n ; α + β , η + n 0 ; 1 , β φ x = x β α Γ α + n 0 x x s α + n 1 2 F 1 α + β , η + n α + n ; 1 s x φ s d s = I 0 , x α + n , β n , η n φ x ,
we note the following interesting and remarkable relation:
n x n I 0 , x α + n , β n , η n φ x = I 0 , x α , β , η φ x ,
which serves as the definition of Saigo’s generalized fractional derivative (see Ref. [3], p. 8, Equation (3.2)).

4. Relationship with Khudozhnikov’s Work

In a very short paper, Khudozhnikov [17] considered in a certain class of integrable functions the following Volterra-type integral equation
a x x s γ 1 Γ ( γ ) 3 F 2 α , β , γ , ε + m ε ; 1 x s φ s d s = g x ,
where 0 < ( γ ) < 1 , m N and 0 < a x b < + . By using some known formulas from Ref. [33], Khudozhnikov obtained the following result [17], p. 79, Equation (2).
Theorem 4
(Khudozhnikov). The Volterra-type integral Equation (49) can be reduced to the following system of differential and integral equations:
k = 0 m m k ( α ) k ( β ) k ( ε ) k ( x ) k y ( m k ) ( x ) = g ( x ) x α + β γ , a x ( x s ) γ + m 1 Γ ( γ + m ) 2 F 1 γ α , γ β γ + m ; 1 x s s α + β γ φ s d s = y x ,
with initial conditions y a = y a = = y ( m 1 ) a = 0 .
In Ref. [17], Khudozhnikov briefly mentioned that the result can be generalized to those equations involving the generalized hypergeometric functions p + 1 F p , p F p and p 1 F p . However, he did not give possible forms of the generalizations or the formulas to be used. In fact, the most likely generalization requires use of a generalized Euler-type transformation, which is not included in Ref. [33]. Therefore, we think that the question of finding a generalization of Theorem 4 is still open.
In this section, we first propose a generalization of Theorem 4. We then consider a Volterra-type integral equation generated by the operator I defined by (12) and obtain an analogue of Khudozhnikov’s theorem.

4.1. A Generalization of Khudozhnikov’s Theorem

Let us consider the Volterra-type integral equation
a x x s γ 1 Γ ( γ ) r + 2 F r + 1 α , β , γ , ( f r + m r ) ( f r ) ; 1 x s φ s d s = g x ,
where 0 < ( γ ) < 1 , m N and 0 < a x b < + . Obviously, (50) reduces to (49) when r = 1 , f 1 = ε and m 1 = m .
By using a lemma due to Miller and Paris [20], p. 298, Lemma 4, and the classical Euler transformation [18], p. 68, Equation (2.2.7), we can express the r + 2 F r + 1 -function as a finite sum of 2 F 1 -functions given by
r + 2 F r + 1 α , β , γ , ( f r + m r ) ( f r ) ; x = k = 0 m A k A 0 ( α ) k ( β ) k ( γ ) k 2 F 1 α + k , β + k γ + k , ; x x k = 1 x γ α β k = 0 m A k A 0 ( α ) k ( β ) k ( γ ) k 2 F 1 γ α , γ β γ + k , ; x x 1 x k .
Then (50) can be written as
g x = a x x s γ 1 Γ ( γ ) r + 2 F r + 1 α , β , γ , ( f r + m r ) ( f r ) ; 1 x s φ s d s = k = 0 m A k A 0 ( α ) k ( β ) k x k x γ α β a x x s γ + k 1 Γ ( γ + k ) 2 F 1 γ α , γ β γ + k ; 1 x s s α + β γ φ s d s .
Let
y x : = a x x s γ + m 1 Γ ( γ + m ) 2 F 1 γ α , γ β γ + m ; 1 x s s α + β γ φ s d s .
In view of the derivative Formula (44), we have
m k x m k x s γ + m 1 Γ ( γ + m ) 2 F 1 γ α , γ β γ + m ; 1 x s = x s γ + k 1 Γ ( γ + k ) 2 F 1 γ α , γ β γ + k ; 1 x s ,
and therefore
y ( m k ) x = a x x s γ + k 1 Γ ( γ + k ) 2 F 1 γ α , γ β γ + k ; 1 x s s α + β γ φ s d s
and y ( m k ) a = 0 for k = 1 , , m 1 . Now (52) can be expressed as
k = 0 m A k A 0 ( α ) k ( β ) k x k y ( m k ) x = x α + β γ g x .
The above steps concerning the integral Equation (50) therefore yield the following theorem.
Theorem 5.
The Volterra-type integral Equation (50) can be reduced to the following system of differential and integral equations:
k = 0 m A k A 0 ( α ) k ( β ) k x k y ( m k ) ( x ) = g ( x ) x α + β γ , a x ( x s ) γ + m 1 Γ ( γ + m ) 2 F 1 γ α , γ β γ + m ; 1 x s s α + β γ φ s d s = y x ,
with initial conditions y a = y a = = y ( m 1 ) a = 0 , where A k ( 0 k m ) is defined in (11).
To show that Theorem 5 contains Khudozhnikov’s result as a special case, we only need to prove that
A k A 0 = m k 1 ε k .
Our calculations require some basics on the theory of combinatorics.
When r = 1 , f 1 = ε and m 1 = m , we get
A 0 = ( ε ) m and A k = j = k m j k σ ^ m j ,
where σ ^ m j is generated by
( ε + x ) m = j = 0 m σ ^ m j x j .
We need in fact to find an explicit expression for σ ^ m j . By using the Chu–Vandermonde identity [18], p. 70, we have
( ε + x ) m = k = 0 m m k ( ε ) m k ( x ) k .
Recall that
( x ) k = j = 0 k ( 1 ) k j s ( k , j ) x j = j = 0 k k j x j ,
where s ( k , j ) is the Stirling number of the first kind and the symbol k j is usually used to denote the unsigned Stirling number of the first kind (see Ref. [34], p. 239). Substituting (57) into (56) and then interchanging the order of summation, we obtain
( ε + x ) m = k = 0 m m k ( ε ) m k j = 0 k k j x j = j = 0 m k = j m m k ( ε ) m k k j x j .
Comparing (58) with (55), it follows that
σ ^ m j = k = j m m k ( ε ) m k k j ,
and combining (54) with (59) and taking into account the index factorization
[ k j m ] [ j m ] = [ k j m ] = [ k m ] [ k j ] ,
we obtain
A k = j = k m j k = j m m ( ε ) m j = = k m m ( ε ) m j = k j k j = 1 k 1 ! = k m m k ( ε ) m 1 1 = 1 k 1 ! = 0 m k m + k + k k ( ε ) m k 1 + k 1 = m k = 0 m k m k ( ε ) m k k = m k ε m ε k ,
where we have used the familiar convoluation identity (see, for example Ref. [34], p. 240)
j = k j k j = k 1 ! k 1 ! ( k 1 ) .
Evidently, (60) is equivalent to (53).

4.2. A Variant of Khudozhnikov’s Theorem

A comparison of the fractional integral operator I with Equations (49) and (50) inspire us to consider the following integral equation
ρ x x ν s ν μ 1 Γ ( μ ) r + 2 F r + 1 a , b , μ , ( f r + m r ) ( f r ) ; 1 s ν x ν φ s s ν h + ν 1 d s = g x ,
where 0 < ( μ ) < 1 and 0 < ρ s x < .
Using the Euler-type transformation (43), then Equation (61) can be converted into
ρ x x ν s ν μ 1 Γ ( μ ) m + 2 F m + 1 a , μ b m , μ , ( ξ m + 1 ) ( ξ m ) ; 1 x ν s ν φ s s ν ( h a ) + ν 1 d s = x ν a g x ,
where ( ξ m ) are nonvanishing zeros of the parametric polynomial Q m t of degree m given by (39).
By using the same lemma of Miller and Paris [20], p. 298, Lemma 4 and the Euler transformation [18], p. 68, Equation (2.2.7) or else using Equation (51), we can express (as in the proof of Theorem 5) the m + 2 F m + 1 -function as a finite sum of 2 F 1 -functions given by
m + 2 F m + 1 a , μ b m , μ , ( ξ m + 1 ) ( ξ m ) ; x = k = 0 m A k ( μ ) k 2 F 1 a + k , μ b m + k , μ + k , ; x x k = 1 x b a + m k = 0 m A k ( μ ) k 2 F 1 μ a , b + m , μ + k , ; x x 1 x k ,
where
A k : = ( a ) k ( μ b m ) k ξ 1 ξ m j = k m j k σ ˜ m j
and σ ˜ j ( 0 j m ) are generated by (41). With the help of (63), the integral Equation (62) can then be written as
k = 0 m A k ( x ν ) k ρ x x ν s ν μ + k 1 Γ ( μ + k ) 2 F 1 μ a , b + m , μ + k , ; 1 x ν s ν φ s s ν ( h b m ) + ν 1 d s = x ν ( b + m ) g x .
By making use of (44), we obtain
m k z m k z μ + m 1 2 F 1 μ a , b + m μ + m ; z = ( μ ) m ( μ ) k z μ + k 1 2 F 1 μ a , b + m μ + k ; z ,
and thus
m k x m k ( x ν s ν ) μ + m 1 Γ ( μ + m ) 2 F 1 μ a , b + m μ + m ; 1 x ν s ν = ν m k ( x ν s ν ) μ + k 1 Γ ( μ + k ) 2 F 1 μ a , b + m μ + k ; 1 x ν s ν .
Substituting (66) into (65), we get
k = 0 m A k ν k ( x ν ) k ρ x m k x m k x ν s ν μ + m 1 Γ ( μ + m ) 2 F 1 μ a , b + m , μ + m , ; 1 x ν s ν φ s s ν ( h b m ) + ν 1 d s = ν m x ν ( b + m ) g x .
Finally, using the Leibniz integral rule and simplifying the resulting formula by the Pfaff transformation [18], p. 68, Equation (2.2.6), we obtain
k = 0 m A k ν k ( x ν ) k m k x m k { x ν ( a μ ) ρ x x ν s ν μ + m 1 Γ ( μ + m ) 2 F 1 μ a , μ b , μ + m , ; 1 s ν x ν φ s s ν ( h + μ b a m ) + ν 1 d s } = ν m x ν ( b + m ) g x .
If
y x = x ν ( a μ ) ρ x x ν s ν μ + m 1 Γ ( μ + m ) 2 F 1 μ a , μ b , μ + m , ; 1 s ν x ν φ s s ν ( h + μ b a m ) + ν 1 d s ,
then the above details concerning the integral equation (61) may be put in the following theorem.
Theorem 6.
The Volerra-type integral Equation (61) can be reduced to the following system of differential and integral equations:
k = 0 m A k ν k ( x ν ) k y ( m k ) x = ν m x ν ( b + m ) g x , x ν ( a μ ) ρ x x ν s ν μ + m 1 Γ ( μ + m ) 2 F 1 μ a , μ b , μ + m , ; 1 s ν x ν φ s s ν ( h + μ b a m ) + ν 1 d s = y x ,
with initial conditions y ( ρ ) = y ( ρ ) = = y ( m 1 ) ( ρ ) = 0 , where A k ( 0 k m ) is given by (64).

5. Conclusions

In this paper, some composition formulas of I and J defined by (12) and (13) are obtained by making use of a Erdélyi-type integral. We find a derivative formula, which in the future may enable us to define a new fractional derivative operator. Finally, we generalize Khudozhnikov’s work on Volterra-type integral equation and find its relationship with our operator I .
Considering the obtained properties of the operators I and J , we briefly mention here some problems that deserve further study.
(i)
Since only two composition formulas for I and J are found in the present work, which is still a very small number compared to the number of the composition formulas of Saigo’s operators I α , β , η and J α , β , η , it may be worthwhile if additional composition structures can be discovered for the operators I and J . The exploration in this direction may also lead us to new discoveries related to the Erdélyi-type integrals;
(ii)
The present work together with our previous papers [14,16] have established many fundamental properties of I and J . For further possible work, some new properties and problems may be worthy of attention in view of the classical books [4,23] on the subject and some recent review articles contained, for example, in Ref. [35]. In particular, it may be worthwhile to first focus on the problem of finding a reasonable analogue of the well known limit case formula, viz. lim α 0 I a + α φ x = φ x concerning the Riemann–Liouville fractional integral operator (see Ref. [23], p. 51, Theorem 2.7).

Author Contributions

Conceptualization, M.-J.L. and R.K.R.; methodology, M.-J.L. and R.K.R.; writing—original draft preparation, M.-J.L. and R.K.R.; writing—review and editing, M.-J.L. and R.K.R.; funding acquisition, M.-J.L. All authors have read and agreed to the published version of the manuscript.

Funding

The research of the first author is supported by National Natural Science Foundation of China (No. 12001095).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank the referees for their comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Saigo, M. A remark on integral operators involving the Gauss hypergeometric functions. Math. Rep. Coll. Gen. Educ. Kyushu Univ. 1978, 11, 135–143. [Google Scholar]
  2. Saigo, M. On the Hölder continuity of the generalized fractional integrals and derivatives. Math. Rep. Coll. Gen. Educ. Kyushu Univ. 1980, 12, 55–62. [Google Scholar]
  3. Saigo, M. A generalization of fractional calculus and its applications to Euler-Darboux equation. RIMS Kokyuroku 1981, 412, 33–56. [Google Scholar]
  4. Kiryakova, V. Generalized Fractional Calculus and Applications; Pitman Research Notes in Mathematics Series No. 301; Longman Scientific and Technical: Harlow, UK, 1994. [Google Scholar]
  5. Naheed, S.; Mubeen, S.; Rahman, G.; Khan, A.Z.; Nisar, K.S. Certain integral and differential formulas involving the product of Srivastava’s polynomials and extended Wright function. Fractal Fract. 2022, 6, 93. [Google Scholar] [CrossRef]
  6. Dziok, J.; Raina, R.K.; Sokół, J. Applications of differential subordinations for norm estimates of an integral operator. Proc. R. Soc. Edinb. Sect. A Math. 2018, 148, 281–291. [Google Scholar] [CrossRef]
  7. Kiryakova, V. On two Saigo’s fractional integral operators in the class of univalent functions. Fract. Calc. Appl. Anal. 2006, 9, 159–176. [Google Scholar]
  8. Saigo, M.; Raina, R.K. On the fractional calculus operator involving Gauss’s series and its application to certain statistical distributions. Rev. Téc. Ing. Univ. Zulia 1991, 14, 53–62. [Google Scholar]
  9. Goyal, S.P.; Jain, R.M. Fractional integral operators and the generalized hypergeometric functions. Indian J. Pure Appl. Math. 1987, 18, 251–259. [Google Scholar]
  10. Goyal, S.P.; Jain, R.M.; Gaur, N. Fractional integral operators involving a product of generalized hypergeometric functions and a general class of polynomials. Indian J. Pure Appl. Math. 1991, 22, 403–411. [Google Scholar]
  11. Goyal, S.P.; Jain, R.M.; Gaur, N. Fractional integral operators involving a product of generalized hypergeometric functions and a general class of polynomials. II. Indian J. Pure Appl. Math. 1992, 23, 121–128. [Google Scholar]
  12. Saigo, M.; Saxena, R.K.; Ram, J. Certain properties of operators of fractional integration associated with Mellin and Laplace transformations. In Current Topics in Analytic Function Theory; Srivastava, H.M., Owa, S., Eds.; World Scientific: Singapore, 1992; pp. 291–304. [Google Scholar]
  13. Araci, S.; Rahman, G.; Ghaffar, A.; Azeema; Nisar, K.S. Fractional calculus of extended Mittag-Leffler function and its applications to statistical distribution. Mathematics 2019, 7, 248. [Google Scholar] [CrossRef]
  14. Luo, M.-J.; Raina, R.K. Fractional integral operators characterized by some new hypergeometric summation formulas. Fract. Calc. Appl. Anal. 2017, 20, 422–446. [Google Scholar] [CrossRef]
  15. Luo, M.-J.; Raina, R.K. On a multiple Čebyšev type functional defined by a generalized fractional integral operator. Tbil. Math. J. 2017, 10, 161–169. [Google Scholar] [CrossRef]
  16. Luo, M.-J.; Raina, R.K. The decompositional structure of certain fractional integral operators. Hokkaido Math. J. 2019, 48, 611–650. [Google Scholar] [CrossRef]
  17. Khudozhnikov, V.I. Integration of Volterra-type integral equations of the first kind with kernels containing some generalized hypergeometric functions. Matem. Mod. 1995, 7, 79. [Google Scholar]
  18. Andrews, G.E.; Askey, R.; Roy, R. Special Functions; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar]
  19. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; North-Holland Mathematics Studies; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  20. Miller, A.M.; Paris, R.B. Transformation formulas for the generalized hypergeometric function with integral parameter differences. Rocky Mt. J. Math. 2013, 43, 291–327. [Google Scholar] [CrossRef]
  21. Luo, M.-J.; Raina, R.K. Erdélyi-type integrals for generalized hypergeometric functions with integral parameter differences. Integral Transform. Spec. Funct. 2017, 28, 476–487. [Google Scholar] [CrossRef]
  22. Grinko, A.P.; Kilbas, A.A. On compositions of generalized fractional integrals. J. Math. Res. Expo. 1991, 11, 165–171. [Google Scholar]
  23. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach: New York, NY, USA, 1993. [Google Scholar]
  24. Kiryakova, V. On the origins of generalized fractional calculus. AIP Conf. Proc. 2015, 1690, 050007. [Google Scholar]
  25. Kiryakova, V. Fractional calculus operators of special functions? The result is well predictable! Chaos Solitons Fractals 2017, 102, 2–15. [Google Scholar] [CrossRef]
  26. Baleanu, D.; Agarwal, P. On generalized fractional integral operators and the generalized Gauss hypergeometric functions. Abstr. Appl. Anal. 2014, 2014, 630840. [Google Scholar] [CrossRef] [Green Version]
  27. Bansal, M.K.; Kumar, D.; Jain, R. A study of Marichev-Saigo-Maeda fractional integral operators associated with the S-generalized Gauss hypergeometric function. Kyungpook Math. J. 2019, 59, 433–443. [Google Scholar]
  28. Brychkov, Y.A.; Glaeske, H.-J.; Marichev, O.I. Factorization of integral transformations of convolution type. J. Math. Sci. 1985, 30, 2071–2094. [Google Scholar] [CrossRef]
  29. Luchko, Y. General fractional integrals and derivatives with the Sonine kernels. Mathematics 2021, 9, 594. [Google Scholar] [CrossRef]
  30. Olver, F.W.J.; Lozier, D.W.; Boisvert, R.F.; Clark, C.W. (Eds.) NIST Handbook of Mathematical Functions; Cambridge University Press: New York, NY, USA, 2010. [Google Scholar]
  31. Bühring, W. Generalized hypergeometric functions at unit argument. Proc. Am. Math. Soc. 1992, 114, 145–153. [Google Scholar] [CrossRef]
  32. Grinko, A.P.; Kilbas, A.A. On compositions of generalized fractional integrals and evaluation of definite integrals with Gauss hypergeometric functions. J. Math. Res. Expo. 1991, 11, 443–446. [Google Scholar]
  33. Prudnikov, A.P.; Brychkov, J.A.; Marichev, O.I. Integrals and Series. Volume 3: More Special Functions; Gould, G.G., Ed.; Gordon and Breach Science Publishers: New York, NY, USA, 1990. [Google Scholar]
  34. Chu, W. Disjoint convolution sums of Stirling numbers. Math. Commun. 2021, 26, 239–251. [Google Scholar]
  35. Hilfer, R.; Luchko, Y. Desiderata for Fractional Derivatives and Integrals. Mathematics 2019, 7, 149. [Google Scholar] [CrossRef] [Green Version]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Luo, M.-J.; Raina, R.K. On the Composition Structures of Certain Fractional Integral Operators. Symmetry 2022, 14, 1845. https://doi.org/10.3390/sym14091845

AMA Style

Luo M-J, Raina RK. On the Composition Structures of Certain Fractional Integral Operators. Symmetry. 2022; 14(9):1845. https://doi.org/10.3390/sym14091845

Chicago/Turabian Style

Luo, Min-Jie, and Ravinder Krishna Raina. 2022. "On the Composition Structures of Certain Fractional Integral Operators" Symmetry 14, no. 9: 1845. https://doi.org/10.3390/sym14091845

APA Style

Luo, M. -J., & Raina, R. K. (2022). On the Composition Structures of Certain Fractional Integral Operators. Symmetry, 14(9), 1845. https://doi.org/10.3390/sym14091845

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