Next Article in Journal
Hermite Interpolation Based Interval Shannon-Cosine Wavelet and Its Application in Sparse Representation of Curve
Previous Article in Journal
Picture Fuzzy MCDM Approach for Risk Assessment of Railway Infrastructure
Previous Article in Special Issue
The General Fractional Derivative and Related Fractional Differential Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Unified Approach to Fractional Calculus Images of Special Functions—A Survey

by
Virginia Kiryakova
Institute of Mathematics and Informatics, Bulgarian Academy of Sciences, 1113 Sofia, Bulgaria
Submission received: 17 November 2020 / Revised: 13 December 2020 / Accepted: 14 December 2020 / Published: 21 December 2020
(This article belongs to the Special Issue Fractional Integrals and Derivatives: “True” versus “False”)

Abstract

:
Evaluation of images of special functions under operators of fractional calculus has become a hot topic with hundreds of recently published papers. These are growing daily and we are able to comment here only on a few of them, including also some of the latest of 2019–2020, just for the purpose of illustrating our unified approach. Many authors are producing a flood of results for various operators of fractional order integration and differentiation and their generalizations of different special (and elementary) functions. This effect is natural because there are great varieties of special functions, respectively, of operators of (classical and generalized) fractional calculus, and thus, their combinations amount to a large number. As examples, we mentioned only two such operators from thousands of results found by a Google search. Most of the mentioned works use the same formal and standard procedures. Furthermore, in such results, often the originals and the images are special functions of different kinds, or the images are not recognized as known special functions, and thus are not easy to use. In this survey we present a unified approach to fulfill the mentioned task at once in a general setting and in a well visible form: for the operators of generalized fractional calculus (including also the classical operators of fractional calculus); and for all generalized hypergeometric functions such as p Ψ q and p F q , Fox H- and Meijer G-functions, thus incorporating wide classes of special functions. In this way, a great part of the results in the mentioned publications are well predicted and appear as very special cases of ours. The proposed general scheme is based on a few basic classical results (from the Bateman Project and works by Askey, Lavoie–Osler–Tremblay, etc.) combined with ideas and developments from more than 30 years of author’s research, and reflected in the cited recent works. The main idea is as follows: From one side, the operators considered by other authors are cases of generalized fractional calculus and so, are shown to be (m-times) compositions of weighted Riemann–Lioville, i.e., Erdélyi–Kober operators. On the other side, from each generalized hypergeometric function p Ψ q or p F q ( p q or p = q + 1 ) we can reach, from the final number of applications of such operators, one of the simplest cases where the classical results are known, for example: to 0 F q p (hyper-Bessel functions, in particular trigonometric functions of order ( q p ) ), 0 F 0 (exponential function), or 1 F 0 (beta-distribution of form ( 1 z ) α z β ) . The final result, written explicitly, is that any GFC operator (of multiplicity m 1 ) transforms a generalized hypergeometric function into the same kind of special function with indices p and q increased by m.

1. Introduction

Special functions (SF) have always been unavoidable tools for mathematicians, physicists, astronomers, applied scientists and engineers while looking to express and study (theoretically, in tables or by numerical algorithms) the solutions of treated mathematical models. On the other side, recently there has been an increased interest in fractional calculus (FC) and its applications, as evidence for which we refer the readers to the data in the survey by Machado–Kiryakova [1]. Fractional calculus is nowadays a favorite, and even a sort of fashionable research area, although the boom of publications and attempts to “fractalize” any kinds of integer order models can bring some threats to the prestige of this discipline, especially in cases of weak or wrong results and not adequate innovations. Let us mention also the phenomenon of hundreds of papers of the last few years (only a few of them can be cited here) dealing with “evaluation of FC images of SF”, most of which use the same standard techniques with changing only the particular special function (SF) and the particular case of the FC operator. Furthermore, it often happens that in such results the originals and the images are special functions of different kinds, or the images are not recognized as known special functions, and thus are not easy to use. In recent papers, such as [2,3,4,5], we share our criticism on this practice and show that all such results can be derived at once by following a general approach, based on ideas from older author’s works on generalized fractional calculus (GFC), since [6].
Here we try to collect the ideas, results and examples from our recent works on the subject. The survey starts with Preliminaries (Section 2) providing a short background on the considered SF and FC operators; followed by Section 3 with results for images of the generalized hypergeometric functions p F q and p Ψ q and their simpler cases under the operators the classical FC operators (Riemann–Liouville and Erdélyi–Kober integrals and derivatives of fractional order). Then, in Section 4 we present our unified approach for evaluation of GFC operators of arbitrary generalized hypergeometric functions ( p F q and p Ψ q ), resulting in the main Theorems 3 and 4. This allows to handle very wide classes of operators of generalized (m-tuple, m 1 ) fractional integration and differentiation and of considered special functions. In Section 5, Section 6 and Section 7 we consider specifications of these results for the Erdélyi–Kober, Saigo and Marichev–Saigo–Maeda (M-S-M) operators, that appear as cases of our GFC, resp. for m = 1 , m = 2 , m = 3 , give their images for the Wright generalized hypergeometric functions, and many illustrative examples for particular results by other authors. Section 8 considers more general cases of GFC operators with arbitrary multiplicity m 1 , as the multiple Gel’fond–Leontiev operators related to the multi-index Mittag–Leffler functions, and the hyper-Bessel operators related to the hyper-Bessel functions of Delerue. In Section 9 we comment on works of other authors on introducing some “new” special functions and show that these are again Wright generalized hypergeometric functions p Ψ q . Therefore, the various FC images they propose come as simple corollaries of our general results. To show the effectiveness of the proposed unified approach, in this survey we collected some 21 examples for FC images of SF, and referred to a long list of other authors’ works on the subject. Section 10 summarizes some conclusions.

2. Preliminaries

Here we provide a short and only necessary background on the considered classes of special functions (SF) and of operators of classical FC and of generalized fractional calculus (GFC), so as to explain the general ideas. All details on defining the single-valued branches of the considered functions, functional spaces, and necessary conditions on appearing parameters, can be found in our previous works, as cited, and for example in ([6], Section 5.5.i). Basically, we consider functions in the complex plane of the form f ( z ) = z μ f ˜ ( z ) ,   μ 0 ,   f ˜ ( z )   analytic   and   single   valued   in   Ω , where Ω is a starlike domain with respect to z = 0 , usually a disk Δ R : | z | < R . Most of the considered special functions are entire functions, or analytic ones in disks in C .
The results we consider are for the classes of so-called generalized hypergeometric functions (g.h.f) with Mellin–Barnes type integral representations, namely the Fox H-function, Meijer G-function and their most widely used cases of Wright g.h.f. p Ψ q and g.h.f. p F q . Even if our aim is to incorporate as large as possible classes of special functions, let us mention that other transcendental functions as the elliptic integrals, Lambert W-, Mathiew-, Zeta-, etc. functions are outside of our studies. Also, we emphasize on results for LHS integrals, although for the RHS ones similar techniques and results are applied; and consider Riemann–Liouville type fractional derivatives. For the Caputo-type differentiation operators, similar but different results will be exposed in a separate work.

2.1. Special Functions of Fractional Calculus

Under “classical” Special Functions (SF) we mean these “mathematical functions” and orthogonal polynomials of which the origin goes back to 18th and 19th centuries and are named after great mathematicians like Euler, Gauss, Riemann, Bessel, Kummer, Legendre, Laguerre. These “Special Functions of Mathematical Physics” appeared with the needs of applied sciences and serve as solutions of integer order (most commonly 2nd order) differential equations from models in mathematical physics. In the last two centuries it was observed that modeling of many phenomena of the physical and social world can be reflected much more adequately by means of differential equations of arbitrary fractional or higher integer orders, and the so-called special functions of fractional calculus (SF of FC) as providing tools for their explicit solutions became unavoidable tools in the hands of theoretical and applied scientists recognizing the power of fractional calculus (FC).
Recently, many handbooks and surveys appeared as dedicated not only to classical SF but also to the SF of FC, to mention some of them: Prudnikov–Brychkov–Marichev [7], Marichev [8], Srivastava–Gupta–Goyal [9], Kilbas–Srivastava–Trujillo [10], Podlubny [11], Kiryakova [6], Yakubovich–Luchko [12], Mathai–Haubold [13], Gorenflo–Kilbas–Mainardi-Rogosin [14]. Such a list cannot be full here, and for more sources see also the survey paper Machado–Kiryakova [1]. In the papers on the topic and in this survey, we limit ourselves to the Fox H-functions of one complex variable, as enough of a general level to expose the proposed approach.
Definition 1
(Ch. Fox 1960). see books such as [6,7,9,10], and earlier and latest ones) The Fox H-function is a generalized hypergeometric function, defined by means of the Mellin–Barnes type contour integral
H p , q m , n z ( a i , A i ) 1 p ( b j , B j ) 1 q = 1 2 π i L H p , q m , n ( s ) z s d s ,   w i t h   H p , q m , n ( s ) = j = 1 m Γ ( b j + B j s ) i = 1 n Γ ( 1 a i A i s ) j = m + 1 q Γ ( 1 b j B j s ) i = n + 1 p Γ ( a i + A i s ) ,
z 0 , where L is a suitable contour (of three possible types in C : L , L , ( γ i , γ + i ) ), the orders ( m , n , p , q ) are non negative integers so that 0 m q , 0 n p , the parameters A i > 0 , B j > 0 are positive, and a i , b j , i = 1 , , p ; j = 1 , , q can be arbitrary complex such that A i ( b j + l ) B j ( a i l 1 ) ,   l , l = 0 , 1 , 2 , ; i = 1 , , n ; j = 1 , , m . Note that the integrand H p , q m , n ( s ) with s s is the Mellin transform of the H-function (1).
The details on the properties of the Fox H-function can be found in many contemporary handbooks on SF such as [7,9,10], where its behavior is described in term of the denotations:
ρ = i = 1 p A i A i j = 1 q B j B j ;   Δ = k = 1 j B j i = 1 p A i ; μ = j = 1 q b j i = 1 p a i + p q 2 ;   a * = i = 1 n A i i = n + 1 p A i + j = 1 m B j j = m + 1 q B j .
Note that the H-function is an analytic function of z in circle domains | z | < ρ or outside them (or in sectors of them, or in the whole C ), depending on the above parameters and the contours.
If all A i = B j = 1 ,   i = 1 ,   , p ;   j = 1 ,   , q , the H-function H p , q m , n z ( a i , 1 ) 1 p ( b j , 1 ) 1 q reduces to the Meijer’s G-function (C.S. Meijer (1936), see details in ([15], Vol.1) and all above-mentioned books)
G p , q m , n z ( a i ) 1 p ( b j ) 1 q = 1 2 π i L G p , q m , n ( s ) z s d s = 1 2 π i L j = 1 m Γ ( b j + s ) i = 1 n Γ ( 1 a i s ) j = m + 1 q Γ ( 1 b j s ) i = n + 1 p Γ ( a i + s ) z s d s ,     z 0 .
Although simpler than (1), the G-function is yet enough general as it incorporates the Classical SF (known also as Named SF) and many elementary functions. See lists of examples, for example, in ([15], Vol.1), ([6], Appendix C).
Now, we attract the readers’ attention to the most typical examples of SF of FC, which are Fox H-functions but not reducible to Meijer G-functions in the general case (of irrational A j , B k ). These originate from works of Sir Edward Maitland (E.-M.) Wright in a series of his works (1935–1940).
Definition 2
(see, e.g., ([6,7,14], App.E)). The Wright generalized hypergeometric function p Ψ q ( z ) , called also Fox–Wright function (F-W g.h.f.) is defined as:
p Ψ q ( a 1 , A 1 ) , , ( a p , A p ) ( b 1 , B 1 ) , , ( b q , B q ) z = k = 0 Γ ( a 1 + k A 1 ) Γ ( a p + k A p ) Γ ( b 1 + k B 1 ) Γ ( b q + k B q )   z k k !
= H p , q + 1 1 , p z ( 1 a 1 , A 1 ) , , ( 1 a p , A p ) ( 0 , 1 ) , ( 1 b 1 , B 1 ) , , ( 1 b q , B q ) .
In terms of parameters (2), the p Ψ q -function is an entire function of z if Δ > 1 , while for Δ = 1 , it is an absolutely convergent series in the disk { | z | < ρ } , and also for | z | = ρ if R e ( μ ) > 1 / 2 , see, for example, [16].
If all A 1 = = A p = 1 ,   B 1 = = B q = 1 , the Wright g.h.f. reduces to the generalized hypergeometric p F q -function, which is a case of the G-function (3), see details in ([15], Vol.1):
p Ψ q ( a 1 , 1 ) , , ( a p , 1 ) ( b 1 , 1 ) , , ( b q , 1 ) z = c p F q ( a 1 , , a p ; b 1 , , b q ;   z ) = k = 0 ( a 1 ) k ( a p ) k ( b 1 ) k ( b q ) k z k k ! = G p , q + 1 1 , p z 1 a 1 , , 1 a p 0 , 1 b 1 , , 1 b q ;   where   c = i = 1 p Γ ( a i ) / j = 1 q Γ ( b j ) , ( a ) k : = Γ ( a + k ) / Γ ( a ) .
The Mittag-Leffler (M-L) function, introduced by G. Mittag-Leffler (1902–1905), with extended 2-parameters’ definition by R.P. Agarwal (1953), was presented yet in Bateman Project’s [15], Vol.3 (1954), in a chapter for “Miscellaneous Functions”. However, it was ignored for a long time in the books on special functions because the applied scientists suffered from a lack of tables for its Laplace transforms. Although appearing from studies not related to fractional calculus, nowadays the M-L function has become the most popular and most exploited SF of FC, honored to be the “Queen”-function of FC. See details, for example, in [14], also in [6,17,18].
Definition 3.
The Mittag-Leffler (M-L) functions E α and E α , β , are entire functions of order ρ = 1 / α and type 1, defined by the power series
E α ( z ) = k = 0 z k Γ ( α k + 1 ) ,   E α , β ( z ) = k = 0 z k Γ ( α k + β ) ,   α > 0 ,   β > 0 .
As “fractional index ( α > 0 ) analogs of the exponential and trigonometric functions that satisfy ODEs of 1st and 2nd order ( α = 1 , 2 ), the M-L functions serve as solutions of fractional order differential equations. A M-L type function with three indices, known as the Prabhakar function (T.R. Prabhakar, 1971) is also often studied and used, for details see [14,17,18,19], and other contemporary books and surveys on M-L type functions:
E α , β γ ( z ) = k = 0 ( γ ) k Γ ( α k + β ) z k k ! ,   α , β , γ C , Re α > 0 ,
where ( γ ) 0 = 1 , ( γ ) k = Γ ( γ + k ) / Γ ( γ ) denotes the Pochhammer symbol. For γ = 1 we get the M-L function E α , β , and if additionally β = 1 , it is E α .
These M-L type functions are simple cases of the Wright g.h.f. and of the H-function, namely:
E α , β ( z ) = 1 Ψ 1 ( 1 , 1 ) ( β , α ) z = H 1 , 2 1 , 1 z ( 0 , 1 ) ( 0 , 1 ) , ( 1 β , α ) ,
E α , β γ ( z ) = 1 Γ ( γ ) 1 Ψ 1 ( γ , 1 ) ( β , α ) z = H 1 , 2 1 , 1 z ( 1 γ , 1 ) ( 0 , 1 ) , ( 1 β , α ) .
A vector index extension of (7) appeared in the works by Luchko et al. (for example, Yakubovich-Luchko [12]) on operational calculus’ methods for some fractional order PDE. Under the name multi-index (multiple) M-L function, it was introduced by Kiryakova [20] using a different approach, via the Gelfond–Leontiev generalized integration and differentiation operators (see Section 8). Further, these functions are studied in detail by Kiryakova [21,22], by Kilbas–Koroleva–Rogosin [23], Paneva–Konovska [19], and many other followers.
Definition 4
(Kiryakova [21,22]). Let m > 1 be an integer, ( α 1 > 0 , α 2 > 0 , , α m > 0 ) and ( β 1 , β 2 , , β m ) be arbitrary real parameters. By means of these “multi-indices”, the multi-index Mittag-Leffler function (multi-M-L f.) is the entire function defined as:
E ( α i ) , ( β i ) ( z ) : = E ( α i ) , ( β i ) ( m ) ( z ) = k = 0 z k Γ ( α 1 k + β 1 ) Γ ( α m k + β m )
= 1 Ψ m ( 1 , 1 ) ( β i , α i ) 1 m z = H 1 , m + 1 1 , 1 z ( 0 , 1 ) ( 0 , 1 ) , ( 1 β i , α i ) 1 m .
Under weakened restrictions on the α’s not obligatory to be all nonnegative, the study was extended by Kilbas et al; see Kilbas–Koroleva–Rogosin [23].
The basic properties and results for the functions (9) and long lists of their examples, all of them having wide applications in solutions of integer- and fractional-order models, are provided in our previous papers like Kiryakova [21,22,24]. Let us shortly mention particular cases like: for m = 1 , we have the classical M-L function E α , β with all its particulars (error-, incomplete gamma-, Rabotnov, etc., functions); and for m > 1 (many of them treated in the examples in next sections) (see ([15], Vol.2, Section 7.5.4, Section 7.5.5), [25]):
the   W r i g h t   f u n c t i o n / B e s s e l - M a i t l a n d   f u n c t i o n   ϕ ( κ , ν + 1 ; z ) : = J ν κ ( z ) = E ( κ , 1 ) , ( ν + 1 , 1 ) ( 2 ) ( z ) = 1 Ψ 1 ( 1 , 1 ) ( ν + 1 , κ ) , ( 1 , 1 ) z = 0 Ψ 1 ( ν + 1 , κ ) z ; the   M a i n a r d i   f u n c t i o n   M ( z ; β ) = ϕ ( β , 1 β ; z ) = E ( β , 1 ) , ( 1 β , 1 ) ( 2 ) ( z ) ; and   its   examples ,   as   M ( z ; 1 / 2 ) = 1 / π exp ( z 2 / 4 ) ,   and   A i r y   f .   M ( z ; 1 / 3 ) = 3 2 / 3 A i ( z / 3 1 / 3 ) ; P a t h a k s   g e n .   W r i g h t - B e s s e l - L o m m e l f .   J ν , λ μ ( z ) = = ( z / 2 ) ν + 2 λ E ( 1 / μ , 1 ) , ( ν + λ + 1 , λ + 1 ) ( 2 ) ( z 2 / 4 ) ; of   which   ,   for   μ = 1 ,   the   L o m m e l   f . appears   as   s 2 λ + ν 1 , ν ( z ) = const · J ν , λ 1 ( z ) , and ,   thus   its   particular   case ,   the   S t r u v e f .   H ν ( z ) ; D z r b a s h j a n s   f u n c t i o n   with   2 × 2   parameters   Φ 1 / α 1 , 1 / α 1 ( z ; β 1 , β 2 ) : = E ( α 1 , α 2 ) , ( β 1 , β 2 ) ( 2 ) ( z ) ; if   all   α i = 1 , i = 1 , m : h y p e r - B e s s e l   f .   o f   D e l e r u e   J ν 1 , , ν m ( m ) ( z ) ,   as   multi - index   ext .   of   the   B e s s e l   f . ;   e t c .
Recently, in Kilbas–Koroleva–Rogosin [23] the definition (9) has been extended for arbitrary values of the α ’s parameters. Paneva-Konovska introduced and studied generalizations of the Prabhakar function (8) by means of three sets of parameters ( α 1 > 0 , α 2 > 0 , , α m > 0 ) , ( β 1 , β 2 , , β m ) , ( γ 1 , γ 2 , , γ m ) , called 3 m -parametric M-L functions, see, for example, [19,26] and references therein. Multivariate and matrix extensions of the M-L and multi-index M-L functions are also explored.
In another survey paper, Kiryakova [27], we are exposing many other details on the theory of the SF of FC, in the sense of Wright generalized hypergeometric functions p Ψ q and multi-index Mittag-Leffler functions, and provide an extensive list of their particular cases, studied in theoretical and applicable aspects by various authors.
Remark 1.
The techniques of the Mellin transform
M f ( z ) ; s : = F * ( s ) = 0 f ( z ) z s 1 d t
is one of the main tools to evaluate integrals and various integral transforms of special functions, including their images under operators of FC. After some classical publications of previous centuries, the main contribution to this approach is due to Marichev [8]. He proposed a natural but wide ranged scheme, based on the contour integral representations of Mellin–Barnes type for the H- and G-functions, like (1) and (3). Note that the integrands H p , q m , n and G p , q m , n are their Mellin transforms (of variable s s ) are fractions of products of 2 × 2 groups of Gamma-functions, and each special function being a special case of the generalized hypergeometric functions, has a particular representation of that kind. For example ([10], (1.11.24)):
M p Ψ q ( a i , A i ) 1 p ( b j , B j ) 1 a z ; s = Γ ( s ) i = 1 p Γ ( a i s A i ) j = 1 q Γ ( b j s B j ) .
For variations of results, one can use in addition the relations (see, e.g., in ([28], (2.6)–(2.8))):
M f ( λ z ) ; s = λ s F * ( s + γ ) ,   λ > 0 ;   M z γ   f ( z ) ; s = F * ( s + γ ) ;   M f ( z μ ) ; s = 1 μ F * ( s μ ) ,   μ > 0 .
Examples for the use of the Mellin transform in this respect are given (among many others works) in: Luchko and Kiryakova ([28], Section 4) (general scheme and examples with the M-L and Wright functions), Agarwal, Rogosin, and Trujillo [29] and Paneva-Konovska and Kiryakova [30] (images for multi-index M-L functions and their particular cases).

2.2. Operators of Generalized Fractional Calculus

In fractional calculus (FC), meant as a theory of the integration and differentiation of arbitrary (including fractional, not obligatorily integer) order, there are several almost equivalent definitions for “fractional” integrals and derivatives, applied in various functional spaces. Here we are interested in evaluating FC operator images of special functions, defined by power series, most of which are entire functions, or at least analytic ones inside/outside disks in a complex plane. Therefore, we restrict our statements to such functions, although they hold also for spaces of weighted continuous or Lebesgue integrable functions on the real half-line.
For the basic background on FC theory and related topics as SF, integral transforms, generalizations and applications, we can refer to the books (among many others, for a longer list see, e.g., Machado and Kiryakova [1]), such as those by: Samko, Kilbas, and Marichev [31], Podlubny [11], Kilbas, Srivastava, and Trujillo [10], Yakubovich and Luchko [12], including the author’s one, Kiryakova [6] and a recent one, Sandev and Tomovski [32].
We state results for the Riemann-Liouville (R-L) operator for integration R δ of order δ > 0 , the corresponding R-L fractional derivative D δ , and its counterpart * D δ in the Caputo sense, that is only the left-hand sided operators of FC (and skip similar details for the Weyl-type, right-hand sided operators).
The main operator of fractional integration we consider is the Erdélyi–Kober operator (E-K) of integration of order δ > 0 , depending on two additional parameters γ R and β > 0 ,
I β γ , δ f ( z ) = 1 Γ ( δ ) 0 1 σ γ ( 1 σ ) δ 1   f ( z σ 1 β ) d σ = z β ( γ + δ ) Γ ( δ ) 0 z ( z β ξ β ) δ 1 ξ β γ f ( ξ ) d ( ξ β ) ,
note it is the identity for δ = 0 . Especially for functional spaces of weighted analytic functions of the form f ( z ) = z μ f ˜ ( z ) ,   μ 0 (see beginning of Section 2), to be preserved by this operator, we require γ > 1 μ β , in addition to δ 0 , β > 0 . This operator, more general than the R-L integral, and having many more applications, was introduced in Sneddon’s works, such as [33], and is considered in books ([6,10,31], Ch.2), and recently in many other works on fractional order models. The Erdélyi–Kober-type fractional integrals, or briefly Erdélyi–Kober integrals, of the form
I f ( z ) = z δ 0 I β γ , δ f ( z ) ,   with   δ 0 0 .
are basic in our studies, and are called classical fractional integrals, and we consider their commutable compositions that are presented as our generalized fractional integrals, Kiryakova [6,34].
The Erdélyi–Kober operator (13) reduces to the R-L operator of integration for γ = 0 , β = 1 , δ 0 = δ ,
R 0 + , z δ f ( z ) : = R δ f ( z ) = z δ I 1 0 , δ f ( z ) ;   and   conversely ,   I 1 γ , δ f ( z ) = z γ δ R δ z γ f ( z ) .
Note that some authors often refer to the Erdélyi–Kober integral (12) as Euler integral transformation, when they are to handle various integral transforms of special functions.
The fractional order derivative of R-L type corresponding to the E-K integral (12), called E-K fractional derivative D β γ , δ , is an extension of the R-L fractional derivative D δ f ( z ) : = d d z n   R n δ   f ( z ) . Instead of ( d / d z ) n , a suitably chosen auxiliary differential operator D n of integer order is used, a polynomial of the Euler differential operator ( z   d / d z ) . It has been introduced and studied in the works of Kiryakova and Luchko et al., ([6], Ch.2) and ([12], Ch.3) and in the next ones, as [35],
D β γ , δ f ( z ) = D n I β γ + δ , n δ f ( z ) = j = 1 n 1 β z d d z + γ + j I β γ + δ , n δ f ( z ) ,   n 1 < δ n ,   n N .
The more formal representation ([6], Ch.1, Equation (1.6.7))
D β γ , δ f ( z ) = z γ D δ z γ + δ f ( z 1 / β ) z z β
serves to provide a better understanding on the structure and nature of (15).
The Caputo-type analogs of the R-L and E-K fractional derivatives are defined in the same way but with exchanged order of the nonnegative order integration and the integer order differentiation, see, for example, [36], also [35], namely,
* D β γ , δ f ( z ) = I β γ + δ , n δ D n f ( z ) .
The notion of generalized fractional integration operators was introduced by S. Kalla (1969–1979), who suggested the common form of such operators (see details and references in [37]),
I f ( z ) = 0 1 Φ ( σ ) σ γ f ( z σ ) d σ = z γ 1 0 z Φ ( ξ z ) ξ γ f ( ξ ) d ξ ,
where Φ ( σ ) is an arbitrary continuous (analytical) function for which the integral makes sense, most commonly a special function as the Bessel, Gauss, G- or H-function. The operators of such generalized fractional calculus (GFC) are expected to include, in particular, these of the classical FC and should satisfy the main axioms for the FC theory.
Note that for a rather general or rather narrow choice of the special function Φ , only some formal operational rules for the generalized fractional integrals (17) can be provided. Therefore, in our generalized fractional calculus (GFC), Kiryakova [6], the suitable choice of the kernel-functions Φ as G m , m m , 0 - and H m , m m , 0 -functions was crucial. In that case, the generalized fractional integrals can be decomposed into commutative products of operators of classical FC (Erdélyi–Kober operators). Thus, the tools of the special functions and the wide usage of the classical FC are combined into a GFC with developed detailed theory and many established applications.
Definition 5
(Kiryakova [6]). The multiple E-K integral (of multiplicity m > 1 ), is defined by means of the real parameters’ sets ( δ 1 0 , δ m 0 ) —multi-order of integration, ( γ 1 , , γ m ) — multi-weight; and ( β 1 > 0 , , β m > 0 ) —additional multi-parameter, as:
I ( β k ) , m ( γ k ) , ( δ k ) f ( z ) : = 0 1 H m , m m , 0 σ ( γ k + δ k + 1 1 β k , 1 β k ) 1 m ( γ k + 1 1 β k , 1 β k ) 1 m f ( z σ ) d σ ,
if k = 1 m δ k > 0 ; and as the identity operator: I ( β k ) , m ( γ k ) , ( 0 , , 0 ) f ( z ) = f ( z ) , if δ 1 = δ 2 = = δ m = 0 .
Note that the above kernel H m , m m , 0 -function is an analytic function in the unit disk and H m , m m , 0 ( σ ) 0 for | σ | > 1 (Kiryakova, [6]). Specially for functional spaces of weighted analytic functions of the form f ( z ) = z μ f ˜ ( z ) , μ 0 (see beginning of Section 2), to be preserved by this operator, we require γ k > 1 μ β k , in addition to δ k 0 , β k > 0 .
If all the β ’s are equal: β 1 = β 2 = = β m = β > 0 , then (18) has a simpler representation where the kernel is a G m , m m , 0 -function of Meijer, which is also analytic in unit disk and G m , m m , 0 ( σ ) 0 for | σ | > 1 ,
I ( β , , β ) , m ( γ k ) , ( δ k ) f ( z ) : = I β , m ( γ k ) , ( δ k ) f ( z ) = 0 1 G m , m m , 0 σ ( γ k + δ k ) 1 m ( γ k ) 1 m f ( z σ 1 / β ) d σ = k = 1 m I β γ k , δ k f ( z ) .
The operators of the form
I ˜ f ( z ) = z δ 0 I ( β k ) , m ( γ k ) , ( δ k ) f ( z ) ,   I ˜ f ( z ) = z δ 0 I β , m ( γ k ) , ( β k ) f ( z ) ,   with   δ 0 0 ,
are both referred to shortly as generalized fractional integrals of multi-order ( δ 1 , , δ m ) .
The following decomposition property is proved in [6], etc. (see, e.g., decomposition Th.5.2.1 in [6]). It is important because the GFC integrals (18) and (19) can be represented not only by using the kernel Fox H- and G-functions, but also by means of the repeated integral representations for the commutative product of classical E-K operators (12):
I ( β k ) , m ( γ k ) , ( δ k ) f ( z ) : = k = 1 m I β k γ k , δ k f ( z ) = 0 1 0 1 k = 1 m ( 1 σ k ) δ k 1 σ k γ k Γ ( δ k ) f z σ 1 1 / β 1 σ m 1 / β m d σ 1 d σ m .
In the book Kiryakova [6] and subsequent papers, we have provided the operational properties of the operators (18) and (19) as semigroup property, formal inversion formula, reduction to identity or to the conventional integration operators for special parameters’ choice. This is to justify their names as operators of GFC.
Following the idea of how the R-L and E-K fractional derivatives are defined, we have proposed the definition of the corresponding generalized fractional derivatives. To this end, the auxiliary differential operator D η , a polynomial of z ( d d z ) of degree η 1 + + η m , is used:
D η = r = 1 m j = 1 η r 1 β r z d d z + γ r + j ,   with   η k : = [ δ k ] + 1 , for   noninteger   δ k , δ k , for   integer   δ k ,   k = 1 , , m .
Definition 6
(Kiryakova ([6], Ch.1,Ch.5), [34,35]). The multiple (m-tuple) Erdélyi-Kober fractional derivative of R-L type of multi-order δ = ( δ 1 0 , , δ m 0 ) is defined by means of the differ-integral operator:
D ( β k ) , m ( γ k ) , ( δ k ) f ( z ) : = D η I ( β k ) , m ( γ k + δ k ) , ( η k δ k ) f ( z ) = D η 0 1 H m , m m , 0 σ ( γ k + η k + 1 1 β k , 1 β k ) 1 m ( γ k + 1 1 β k , 1 β k ) 1 m f ( z σ ) d σ .
Analogously, the Caputo-type generalized fractional derivative has been introduced in Kiryakova and Luchko [35], as
* D ( β k ) , m ( γ k ) , ( δ k ) f ( z ) = I ( β k ) , m ( γ k + δ k ) , ( η k δ k ) D η f ( z ) .
For all equal β ’s: β 1 = = β m = β > 0 , the R-L and Caputo-type “derivatives” corresponding to the generalized fractional integral (19) has a simpler form with Meijer G-function in the kernel:
D β , m ( γ k ) , ( δ k ) f ( z ) = D η I β , m ( γ k + δ k ) , ( η k δ k ) f ( z ) = r = 1 m j = 1 η r 1 β z d d z + γ r + j I β , m ( γ k + δ k ) , ( η k δ k ) f ( z ) , and   * D β , m ( γ k ) , ( δ k ) f ( z ) = I β , m ( γ k + δ k ) , ( η k δ k ) D η f ( z ) .
Under generalized (multiple, multi-order) fractional derivatives of the R-L type, resp. of the Caputo type, we have in mind all the differ-integral/integro-differential operators of the form
D ˜ f ( z ) = D ( β k ) , m ( γ k ) , ( δ k ) z δ 0 f ( z ) = z δ 0 D ( β k ) , m ( γ k δ 0 β ) , ( δ k ) f ( z ) ,   * D ˜ f ( z ) = * D ( β k ) , m ( γ k ) , ( δ k ) z δ 0 f ( z )   with   δ 0 0 .
A basic formula for the image of a power function in the general case of (18) and (19) (say from Kiryakova [6]) reads as
I ( β k ) , m ( γ k ) , ( δ k ) { z p } = c p z p ,   with   c p = i = 1 m Γ ( γ i + 1 + p / β i ) Γ ( γ i + δ i + 1 + p / β i ) ,   δ k 0 ,   p > β k ( γ k + 1 ) , k = 1 , , m ,
and a similar one holds for the generalized fractional derivatives, both analogous to the same formulas for the classical Erdélyi–Kober operators. These results are in the base of the standard techniques applied by other authors for the evaluation of FC operators of special functions, in various particular cases. Using (27), in particular for p = k = 0 , 1 , 2 , , n , , then interchanging the integration and summation of the power series for a particular special function, the authors of mentioned papers obtain a new power series to be recognized as another special function (in the successful cases, or the result is useless, left just as such series or as some p Ψ q -function). Our general result states as follows.
Theorem 1
(Kiryakova, since 1988, see, e.g., ([6,38], Ch.5)). Let the conditions β k ( γ k + 1 ) > μ ,   δ k 0 ,   β k > 0 ,   k = 1 , , m , be satisfied for the parameters of the multiple E-K integral (18). Then, it preserves the class of weighted analytic functions f ( z ) in a disk Δ R , denoted by H μ ( Δ R = { | z | < R } ) :
f ( z ) = z μ k = 0 a k z k = z μ ( a 0 + a 1 z + ) H μ ( Δ R )   with   R = lim sup k | a k | k 1 .
Namely, the images of such functions have the same form:
I ( β k ) , m ( γ k ) , ( δ k ) f ( z ) = z μ k = 0 a k b k z k H μ ( Δ R ) ,   with   b k = i = 1 m Γ ( γ i + k + μ β i + 1 ) Γ ( γ i + δ i + k + μ β i + 1 ) > 0 ,
with the same radius of convergence R > 0 and the same signs of the coefficients in their series expansions.

2.3. Some Special Cases of GFC Operators

We emphasize here only some operators of FC that are recently exploited very often in publications on FC operators of SF. In [6] and the author’s other papers as well as in works by other authors, there many other particular cases of linear integral and differential operators provided and used with applications in geometric (univalent) function theory, in differential and integral equations of integer and fractional order, operational calculus, transmutation theory, special functions theory, mathematical models of phenomena of fractional order, etc.
For m = 1 the kernel-functions of the generalized fractional integrals and derivatives (18) and (23) can be represented as
H 1 , 1 1 , 0 σ ( γ + δ , 1 / β ) ( γ , 1 / β ) = β σ β 1 G 1 , 1 1 , 0 σ β γ + δ γ = β σ β γ + β 1 ( 1 σ β ) δ 1 Γ ( δ ) ,
therefore we have the E-K and R-L ( γ = 0 , β = 1 ) operators of classical FC. Many other integration and differentiation operators introduced and used by different authors appear as special cases of I β , 1 γ , δ = I β γ , δ , D β , 1 γ , δ = D β γ , δ , R δ and D δ .
When m = 2 , the kernels H 2 , 2 2 , 0 and G 2 , 2 2 , 0 reduce to a Gauss hypergeometric lfunction:
H 2 , 2 2 , 0 σ ( γ 1 + δ 1 + 1 1 β , 1 β ) , ( γ 2 + δ 2 + 1 1 β , 1 β ) ( γ 1 + 1 1 β , 1 β ) , ( γ 2 + 1 1 β , 1 β ) = 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 σ β ) .
In this case, the generalized fractional integrals are known as hypergeometric fractional integrals, and some of them are introduced and studied by Love, Saxena, Kalla, Saigo (see in next Section 6), Hohlov, etc.
In the case m = 3 , a recently very popular example is with the Marichev–Saigo–Maeda (M-S-M) operators. These FC integration operators are introduced and studied by Marichev [39] and by Saigo et al. [40]. Their kernel-function, the Appel F 3 function (Horn function)
F 3 a , a , b , b , c , z , ξ = m , n = 0 ( a ) m ( a ) n ( b ) m ( b ) n ( c ) m + n z m ξ n m ! n ! ,   | z | < 1 ,   | ξ | < 1   ( see , e . g . ,   [ 7 , 15 ] ) ,
appears as a case of the G 3 , 3 3 , 0 -function and of H 3 , 3 3 , 0 -function. Indeed, according to [7], p.727, (2); and as observed in Kiryakova [6], p.21:
( 1 σ ) c 1 Γ ( c ) F 3 a , a , b , b , c , 1 1 σ , 1 σ = G 3 , 3 3 , 0 σ a + b , c a , c b a , b , c a b = H 3 , 3 3 , 0 σ ( a + b , 1 ) , ( c a , 1 ) , ( c b , 1 ) ( a , 1 ) , ( b , 1 ) , ( c a b , 1 ) ,   R e   c > 0 .
Therefore, our generalized fractional integrals reduce in this case to the M-S-M integral operators.
Let m 1 be an arbitrary integer, but all δ ’s are integers, say δ 1 = = δ m = 1 . Then we have the Bessel type integral and differential operators of arbitrary (higher) integer order, introduced by Dimovski [41] (see also [42]) and named as hyper-Bessel operators by Kiryakova ([6], Ch.3), as shown related to the hyper-Bessel functions of Delerue [43] as their eigenfunctions (see Example 16, in next Section 8). The studies on these operators gave rise to our GFC, since they appeared as “fractional” integrals and derivatives of integer multi-orders ( 1 , 1 , , 1 ) and for λ > 0 their fractional powers have multi-orders ( λ , λ , , λ ) . In Section 8, we will discuss also the Gelfond–Leontiev operators generated by the multi-index M-L functions, as more general operators of arbitrary multiplicity m > 1 and arbitrary fractional multi-order.
As mentioned, here we stress on only a few particular examples of GFC operators I ( β k ) , m ( γ k ) , ( δ k ) , D ( β k ) , m ( γ k ) , ( δ k ) , that are involved in results on to the topic of this survey. This is because many other authors’ works handle the evaluation of images of various elementary or special functions under the classical or some “generalized operators of FC”—such as the operators of R-L, E-K, Saigo, Marichev–Saigo–Maeda. Say, one takes first the cosine or Bessel function, later the generalized Bessel (Bessel-Maitland) function, then an M-L or generalized M-L function, etc., so as to produce new publications by same standard techniques. Very rarely observed, or mostly is ignored, the fact from relation (21) that these are 2-tuple ( m = 2 ), respectively 3-tuple ( m = 3 ), or m-tuple (arbitrary m > 1 ) compositions of Erdélyi–Kober operators. Therefore, the task can be done at once, if one knows how an E-K operator acts on such special functions, all being cases of Wright g.h.f. (4), and then applying the procedure a suitable number of times (2-, 3-, or m). Thus, the result can be predicted in advance, having in mind the general statements in the next sections.

3. Erdélyi-Kober and Riemann–Liouville Images of p Ψ q , p F q and Simpler Special Functions

Some basic classical results on the topic exist from the previous century that should not be forgotten and on which our approach was built. Namely, the image of a generalized hypergeometric function p F q , with p q + 1 , under the R-L fractional integral/ derivative is shown to be the same special function with indices p and q increased by 1:
R δ z ν 1 p F q ( a 1 , , a p ; b 1 , , b q ; λ z ) = Γ ( ν ) Γ ( δ + ν ) z δ + ν 1 p + 1 F q + 1 ( a 1 , , a p , ν ; b 1 , , b q , δ + ν ; λ z ) ,
with R e   δ > 0 , R e   ν > 0 , p q + 1 ; λ 0 , z C and if p = q + 1 : | λ z | < 1 is additionally required. For this, we can refer to Erdélyi et al. ([44], Vol.2), Ch. XIII, Equations (95)–(97); Askey ([45], p.19), and emphasize the survey by Lavoie-Osler-Tremblay [46], a table on p.261. Then, to make use of (32), in our older works on the topic since 1984-1985, we started from the R-L images of the p F q functions with the lowest possible indices: 1 F 0 , 1 F 1 and 0 F 1 , given by three basic elementary functions, see for example in Kiryakova ([6], Ch.4). Some particular illustrative cases in this direction were mentioned there, as
I 2 1 , ν + 1 2 cos   z = R z 2 ν + 1 2 cos   z z = π 2 ν z ν J ν ( z ) = π 2 ν 1 z ν + 1 0 F 1 ν + 1 ; z 2 4 ,   ν > 1 2 ,
R δ z ν 1 exp ( λ z ) = Γ ( ν ) Γ ( δ + ν ) z δ + ν 1 1 F 1 ( ν ; δ + ν ; λ z ) ,   R e   δ > 0 ,   R e   ν > 0 ,
R δ z ν 1 ( z λ ) μ = ( λ ) μ z δ + ν 1 Γ ( ν ) Γ ( δ + ν ) 2 F 1 ( μ , ν ; δ + ν , b 1 , , b q ; z / λ ) ,   R e   δ > 0 ,   R e   ν > 0 .
Note that (33) is an interpretation of the Poisson integral formula for the Bessel function, that has been generalized in [42] and ([6], Ch.4) to represent the hyper-Bessel functions J ν 1 , , ν m ( m ) with multi-indices ( ν 1 , , ν m ) (the Bessel function J ν is the case for m = q p = 1 with one index ν ), that is, to represent the 0 F q p -functions by means of “generalized cosine” cos m .
R-L integrals/ derivatives of the most general G- and H-functions are also well known in the literature (for example, from [44]), and these are the same type of functions but with increased orders and additional parameters.
Along with the mentioned old classical results, recently, new articles are published on the evaluation of classical (R-L, E-K) or generalized FC operators of classical SF or of SF of FC almost every day (e.g., in 2020: [47,48]), and also of their multivariate or matrix variants. Just as one example on fractional operators for the matrix Wright hypergeometric functions (5), is a 2020 paper [49].
The classical results (32) and (33)–(35) have been extended in our works (as in ([6], Ch.4), [22,24,50]) in terms of the Erdélyi–Kober operators (12) and (15) and for their counterparts of the GFC: I ( β k ) , m ( γ k ) , ( δ k ) , D ( β k ) , m ( γ k ) , ( δ k ) , not only for p F q but for p Ψ q as well. To reduce the Wright g.h.f. p Ψ q in the general case to three basic simplest functions with lowest indices p and q, we also apply modifications as the Wright–Erdélyi–Kober multiple operators with a Bessel–Maitland kernel-function and in general, GFC operators with H-functions like in (20) but with different parameters 1 / β k > 0 and 1 / λ k > 0 in the upper and low row. Details are in Kiryakova [24].
Now we provide some basic statements necessary for the topic of this survey, repeating in a few lines the ideas of the proofs, so as to clarify the approach used.
Lemma 1.
The image of a Wright g.h.f. p Ψ q under the Erdélyi–Kober fractional integral (12) is the same type of function in which the indices p- and q are increased by one, and so, has two additional parameters:
I β γ , δ z c   p Ψ q ( a 1 , A 1 ) , , ( a p , A p ) ( b 1 , B 1 ) , , ( b q , B q ) λ z μ = z c p + 1 Ψ q + 1 ( a i , A i ) 1 p , ( γ + 1 + c / β , μ / β ) ( b j , B j ) 1 q , ( γ + δ + 1 + c / β , μ / β ) λ z μ .
It is supposed that R e   δ > 0 ,   R e   γ > 1 ,   μ > 0 ,   λ 0 , c is arbitrary (real), and if p = q + 1 , then | λ z μ | < 1 is additionally required.
Proof. 
In a simpler case with c = 0 , this is Lemma 1 from Kiryakova [2]. There, a proof is based on the Formula (44) (Section 4) for the integral (the Mellin transform) of a product of two H-functions, since both the p Ψ q -function and the kernel of the E-K operator are cases of H-functions; compare (5) and Section 2.3. This approach will be discussed later for the more general case of GFC operators.
As a very standard technique, to prove (36), one can use term-by-term integration in series (4), similarly to that in Kilbas ([16], Th.2) for the particular case of R-L integral, with c : = ν 1 there:
R δ z c p Ψ q ( a 1 , A 1 ) , , ( a p , A p ) ( b 1 , B 1 ) , , ( b q , B q ) λ z μ = z c + δ p + 1 Ψ q + 1 ( a i , A i ) 1 p , ( c + 1 , μ ) ( b j , B j ) 1 q , ( c + δ + 1 , μ ) λ z μ .
Note that the known simpler formula (32), written by means of E-K integrals, appears as a special case of (36) if all A 1 = = A p = B 1 = = B q = 1 and β = 1 , see Equation (4.2.2′) in Kiryakova [6], for R e   δ > 0 ,   R e   γ > 1 :
I 1 γ , δ p F q ( a 1 , , a p ; b 1 , , b q ; λ z ) = Γ ( γ + 1 ) Γ ( γ + δ + 1 ) p + 1 F q + 1 ( a 1 , , a p , γ + 1 ; b 1 , , b q , γ + δ + 1 ; λ z ) .
Lemma 2.
The image of a Wright g.h.f. p Ψ q under the E-K fractional derivative (15) is the same kind of function but with indices p and q increased by 1, and with 2 the additional parameters:
D β γ , δ z c p Ψ q ( a 1 , A 1 ) , , ( a p , A p ) ( b 1 , B 1 ) , , ( b q , B q ) λ z μ = z c p + 1 Ψ q + 1 ( a i , A i ) 1 p , ( γ + δ + 1 + c / β , μ / β ) ( b j , B j ) 1 q , ( γ + 1 + c / β , μ / β ) λ z μ ,
provided R e   δ > 0 ,   R e   γ > 1 ,   μ > 0 ,   λ 0 , and if p = q + 1 , we require | λ z μ | < 1 .
Proof. 
For c = 0 this is Lemma 3 in Kiryakova [2], and for the case of γ = 0 , β = 1 we have the formula for the R-L fractional derivative from Kilbas ([16], Th.4) (where ν 1 : = c , R e   c > 1 , and the same other conditions):
D δ z c p Ψ q ( a 1 , A 1 ) , , ( a p , A p ) ( b 1 , B 1 ) , , ( b q , B q ) λ z μ = z c δ p + 1 Ψ q + 1 ( a i , A i ) 1 p , ( c + 1 , μ ) ( b j , B j ) 1 q , ( c + 1 δ , μ ) λ z μ .
The same standard term-by-term integration/differentiation technique can be used for the proof of (39).
Consider also the simplest case as an analog of (39), when β = 1 , c = 0 and all A 1 = = A p = B 1 = = B q = 1 . This is our Lemma 4.3.1 from [6] for the p F q -functions. Namely,
D 1 γ , δ p F q ( a 1 , , a p ; b 1 , , b q ; λ z ) = Γ ( γ + δ + 1 ) Γ ( γ + 1 ) p + 1 F q + 1 ( a 1 , , a p , γ + δ + 1 ; b 1 , , b q , γ + 1 ; λ z ) .
In the proof of (41) given in [2], we used the relation (16) between the E-K derivative (15) and the R-L derivative D δ f ( z ) = ( d d z ) n R n δ f ( z ) , n = [ δ ] + 1 , combined with the result (32). Then, employed a formula from ([7], Section 7.2.3, (51)) for differentiation of integer order n of a generalized hypergeometric function p F q with specific parameters as above.
Here we demonstrate a new proof of Lemma 2 for the more general case of Wright function p Ψ q . For simplicity, β = 1 and μ = 1 . Interpreting the E-K derivative (15) as in (16), we have subsequently:
The   L . H . S .   of   ( 41 ) = z γ D δ z γ + δ z c p Ψ q ( a i , A i ) 1 p ( b j , B j ) 1 q λ z
= z γ d d z n R n δ z γ + δ + c p Ψ q ( a i , A i ) 1 p ( b j , B j ) 1 q λ z
( and   due   to   ( 37 ) ) = z γ d d z n z γ + n + c p + 1 Ψ q + 1 ( a i , A i ) 1 p , ( γ + δ + c + 1 , 1 ( b j , B j ) 1 q , ( γ + n + c + 1 , 1 ) λ z .
Now, we use the representation (5) of the Wright g.h.f. as an H-function, and may apply a formula for differentiation of integer order n of the H-function, say Equation (1.69) from Mathai-Saxena-Haubold [51], to continue as follows:
= z γ d d z n z γ + n + c H p + 1 , q + 2 1 , p + 1 λ z ( 1 a i , A i ) 1 p , ( γ δ c , 1 ) ( 0 , 1 ) , ( 1 b j , B j ) 1 q , ( γ n c , 1 )
= z γ z γ + c H p + 2 , q + 3 1 , p + 2 λ z ( γ n c , 1 ) , ( 1 a i , A i ) 1 p , ( γ δ c , 1 ) ( 0 , 1 ) , ( 1 b j , B j ) 1 q , ( γ n c , 1 ) , ( γ c , 1 )
and because of the coincidence of the terms ( γ n c , 1 ) in upper and low parameters’ rows, according to the reduction order formula for the H-function: (1.56) in [51], see also (E.8) in [6], and [7,9]), we have
= z c H p + 1 , q + 2 1 , p + 1 λ z ( 1 a i , A i ) 1 p , ( γ δ c , 1 ) ( 0 , 1 ) , ( 1 b j , B j ) 1 q , ( γ c , 1 ) ,
which, by using again (5) to go back to a Wright g.h.f., gives the result (39). In case β 1 , substitution z z 1 / β is necessary, and same for μ 1 .
Yet another approach to check the validity of (39) is to use the identity D β γ , δ I β γ , δ f ( z ) = f ( z ) for f ( z ) : = z c p Ψ q ( a 1 , A 1 ) , , ( a p , A p ) ( b 1 , B 1 ) , , ( b q , B q ) λ z μ , the result from Lemma 1 and reduction of the intermediate result p + 2 Ψ q + 2 to a p + 1 Ψ q + 1 -function, since the last two equal parameters in the upper and low rows of its series eliminate each other. □
Remark 2.
The corresponding result for the Caputo-type E-K derivative * D β γ , δ for images of the Wright p Ψ q -functions, and in particular also for the p F q -functions and for the simpler case of operators with β 1 = = β m = β > 0 , will be presented in another separate work.
Example 1.
The classical result (34) can be extended to an Erdélyi–Kober integral if we use Lemma 1 for exp ( z ) = 0 F 0 ( ; ; z ) = 1 F 1 ( 0 ; 0 ; z ) = 1 Ψ 1 ( 0 , 1 ) ( 0 , 1 ) z (which is also a 0 Ψ 0 , G 0 , 1 1 , 0 , and H 0 , 1 1 , 0 -function):
I β γ , δ z c exp ( λ z ) = I β γ , δ z c 1 Ψ 1 ( 0 , 1 ) ( 0 , 1 ) λ z = z c 2 Ψ 2 ( 0 , 1 ) , ( γ + 1 + c / β , 1 / β ) ( 0 , 1 ) , ( γ + δ + 1 + c / β , 1 / β ) z = 1 Ψ 1 ( γ + 1 + c / β , 1 / β ) ( γ + δ + 1 + c / β , 1 / β ) λ z , r e d u c i b l e   t o   ( 34 )   f o r   γ = 0 ,   β = 1 ,   ν = c + 1 > 0 .

4. Results for the Generalized Fractional Calculus Operators of Special Functions

Here we present our results on evaluating operators of generalized fractional calculus (in the sense of [6] and of Riemann–Liouville type) of wide classes of special functions as the Wright generalized hypergeometric functions p Ψ q and even of the Fox H-functions (thus incorporating the SF of FC) and in particular, of the p F q - and Meijer G-functions (thus having general results also for the “classical” SF).
We start with the most general result, presented in Kiryakova ([2], Th.3) and mentioned in ([3], Th.4.3).
Theorem 2.
The generalized (m-tuple) fractional integral I ( β k ) , m ( γ k ) , ( δ k ) of a H-function is again an H-function:
I ( β k ) , m ( γ k ) , ( δ k ) H u , v s , t λ z ( c i , C i ) 1 u ( d j , D j ) 1 v = H u + m , v + m s , t + m λ z ( c i , C i ) 1 t , ( γ k ) 1 m , ( c i , C i ) t + 1 u ( d j , D j ) 1 s , ( γ k δ k ) 1 m , ( d j , D j ) s + 1 v .
Note that three of the orders of the H-function are increased by the multiplicity m, and additional m + m parameters appear depending on those of the operator.
Proof. 
The following known formula for integral (can be seen as a Mellin transform) of product of two Fox H-functions is very important for evaluating integrals of products of special functions of general nature, because almost all of them can be presented as H-functions (([9], Section 5.1, (5.1.1)), ([7], Section 2.25, (1)), see also (E.21′) in ([6], Appendix)):
0 σ α 1 H u , v s , t ϰ σ ( c i , C i ) 1 u ( d l , D l ) 1 v   H p , q m , n ω σ r ( a j , A j ) 1 p ( b k , B k ) 1 q d σ = ϰ α H p + v , q + u m + t , n + s ω ϰ r ( a j , A j ) 1 n , ( 1 d l α D l , r D l ) 1 v , ( a j , A j ) n + 1 p ( b k , B k ) 1 m , ( 1 c i α C i , r C i ) 1 u , ( b k , B k ) m + 1 q ,
under the conditions Δ > 1 , a * = Δ + 1 > 0 ) (in terms of (2)).
To prove (43) we use the definition (18) of I ( β k ) , m ( γ k ) , ( δ k ) , the fact that the kernel H m , m m , 0 -function vanishes for | σ | > 1 and so the limits ( 0 , 1 ) of the integral can be changed into ( 0 , ) , and the above Formula (44):
I ( β k ) , m ( γ k ) , ( δ k ) H u , v s , t λ z ( c i , C i ) 1 u ( d j , D j ) 1 v = 0 1 d σ = 0 d σ
= 0 H m , m m , 0 σ ( γ k + δ k + 1 1 β k , 1 β k ) 1 m ( γ k + 1 1 β k , 1 β k ) 1 m H u , v s , t λ z σ ( c i , C i ) 1 u ( d j , D j ) 1 v d σ
= H u + m , v + m s + 0 , t + m λ z ( c i , C i ) 1 t , ( 1 γ k 1 + 1 β k 1 β k , 1 β k ) 1 m , ( c i , C i ) t + 1 u ( d j , D j ) 1 s , ( 1 γ k δ k 1 + 1 β k 1 β k , 1 β k ) 1 m , ( d j , D j ) s + 1 v
= H u + m , v + m s , t + m λ z ( c i , C i ) 1 t , ( γ k , 1 β k ) 1 m , ( c i , C i ) t + 1 u ( d j , D j ) 1 s , ( γ k δ k , 1 β k ) 1 m , ( d j , D j ) s + 1 v .
 □
For β 1 = β 2 = = β m = β > 0 , the image of an arbitrary G-function under the simpler GFC integrals with Meijer’s G m , m m , 0 -kernels has been provided earlier.
Corollary 1
(Lemma 1.2.2, [6])). The I β , m ( γ k ) , ( δ k ) -image of a G-function is also a G-function in which the three orders are increased by the multiplicity m and has additional m + m parameters depending on those of the GFC operator:
I β , m ( γ k ) , ( δ k ) G u , v s , t λ z β ( c i ) 1 u ( d j ) 1 v = G u + m , v + m s , t + m λ z β ( c i ) 1 t , ( γ k ) 1 m , ( c i ) t + 1 u ( d j ) 1 s , ( γ k δ k ) 1 m , ( d j ) s + 1 v .
Proof. 
In this case one can use a formula for the integral of product of two arbitrary G-functions, simpler than (44) (see for example, ([7], Section 2.24, (1)]) and with a proof in ([6], App., (A.29))), and be reminded again that the G m , m m , 0 -function vanishes outside the unit disc. Thus the integral (45)
I = 0 1 d σ = 0 G m , m m , 0 σ ( γ k + δ k ) 1 m ( γ k ) 1 m G u , v s , t λ z β ( c i ) 1 u ( d j ) 1 v d σ ,
gives the required image G-function. Because the G u , v s , t -function is from the space of an analytic function in a disk centered at the origin, and has the following asymptotic behavior
G u , v s , t λ z β = O z d *   as   | z | 0 ,   with   d * = β min j d j > max k [ β ( γ k + 1 ) ] ,
the conditions for the used formula to hold on are satisfied. □
Formulas (43) and (45) can be used to evaluate practically all (classical and generalized) operators of FC of arbitrary SF ( which are representable either as a G- or as a more general H-function).
Now we present the main result on the topic of this survey paper, which comes from Kiryakova [2], Th.1.
Theorem 3.
Assume that the conditions δ k 0 , γ k > 1 , β k > 0 , k = 1 , , m and μ > 0 , λ 0 hold. The image of a Wright g.h.f. p Ψ q ( z ) by a generalized fractional integral (20) (multiple, m-tuple Erdélyi–Kober integral) is another Wright g.h.f. with indices p and q increased by the multiplicity m and with additional parameters coming from those of the GFC integral:
I ( β k ) 1 m , m ( γ k ) 1 m , ( δ k ) 1 m z c   p Ψ q ( a 1 , A 1 ) , , ( a p , A p ) ( b 1 , B 1 ) , , ( b q , B q ) λ z μ = z c   p + m Ψ q + m ( a i , A i ) 1 p , ( γ i + 1 + c β i , μ β i ) 1 m ( b j , B j ) 1 q , ( γ i + δ i + 1 + c β i , μ β i ) 1 m λ z μ .
Proof. 
Here we briefly repeat the proof from Kiryakova [2], Th.1, in order to exhibit the main ideas on which this survey paper is based.
As one approach to prove (46), the general integral formula (44) can be used. This theorem can be seen also as a consequence of Theorem 2. It is because the kernel-function of the operator is a H m , m m , 0 -function and the p Ψ q -function is a H p , q + 1 1 , p -function, see (5). Then, according to (45) the result will be a H p + m , q + 1 + m 1 , p + m -function that should be recognized as a p + 1 Ψ q + 1 -function, because it is reduced to a H p + 1 , q + 2 1 , p + 1 -function in view of the coincidence of ( m 1 ) parameters in the upper and low row (use “reduction order” property of the H-function, [7], Section 8.3, 6.; [6], App. (E.8), etc.).
However, to clarify our main idea it is more instructive to refer to the decomposition property (21) presenting the generalized fractional integral (18) as a product of commuting (classical) Erdélyi–Kober operators. In the simplest case, we use subsequently m-times (36) from Lemma 1, to get:
I ( β k ) 1 m , m ( γ k ) 1 m , ( δ k ) 1 m p Ψ q ( a i , A i ) 1 p ( b j , B j ) 1 q λ z = I ( β k ) 1 m 1 , m 1 ( γ k ) 1 m 1 , ( δ k ) 1 m 1 I β m γ m , δ m p Ψ q ( a i , A i ) 1 p ( b j , B j ) 1 q λ z = I ( β k ) 1 m 2 , m 2 ( γ k ) 1 m 2 , ( δ k ) 1 m 2 I β m 1 γ m 1 , δ m 1 p + 1 Ψ q + 1 ( a i , A i ) 1 p , ( γ m + 1 , 1 / β m ) ( b j , B j ) 1 q , ( γ m + δ m + 1 , 1 / β m ) λ z = = I β 1 γ 1 , δ 1 I β 2 γ 2 , δ 2 p + m 2 Ψ q + m 2 ( a i , A i ) 1 p , ( γ r + 1 , 1 / β r ) 1 m 2 ( b j , B j ) 1 q , ( γ r + δ r + 1 , 1 / β r ) 1 m 2 λ z = I β 1 γ 1 , δ 1 p + m 1 Ψ q + m 1 ( a i , A i ) 1 p , ( γ r + 1 , 1 / β r ) 1 m 1 ( b j , B j ) 1 q , ( γ r + δ r + 1 , 1 / β r ) 1 m 1 λ z = p + m Ψ q + m ( a i , A i ) 1 p , ( γ r + 1 , 1 / β r ) 1 m ( b j , B j ) 1 q , ( γ r + δ r + 1 , 1 / β r ) 1 m λ z .
To derive the general relation (46) we apply to the above result the property for “generalized commutation” from Kiryakova ([6], Ch.5, (5.1.28)), namely:
I ( β k ) , m ( γ k ) , ( δ k ) z c f ( z μ ) = z c I ( β k m ) , m ( γ k + c β k ) , ( δ k ) f ( z μ ) ,   with   μ > 0 .
 □
Corollary 2.
(Lemma 4.2.1., Equations (4.2.2)–(4.2.2’) in Kiryakova ([6], Ch.4)) The image of a p F q g.h.f. (6) under a generalized (m-tuple) fractional integral I 1 , m ( γ k , δ k ) : = I ( 1 , , 1 ) , m ( γ k , δ k ) , (19) with (for simplicity) all β k = β = 1 , k = 1 , , m , is another g.h.f. of the same kind with indices increased by the multiplicity m:
I 1 , m ( γ k ) 1 m , ( δ k ) 1 m p F q a 1 , , a p ; b 1 , , b q ; λ z = p + m F q + m a 1 , , a p , ( γ i + 1 ) 1 m ; ( b 1 , , b q , ( γ i + δ i + 1 ) 1 m ; λ z .
The above results (46) and (49) can be interpreted alternatively as the assertions stated in our earlier works ([6,50], Ch.4) (in the simpler case of Corollary 2), and later in [24] (in more general case of Theorem 3). That is, a p + m Ψ q + m -function (resp. a p + m F q + m -function) of the form below can be represented by means of a multiple (m-tuple) operator of GFC
I ˜ = I ( 1 / β i ) i = 1 m , m ( a p + i 1 ) i = 1 m , ( b q + i a p + i ) i = 1 m
of a p Ψ q -function (resp. a p F q -function), with orders reduced by m, namely:
p + m Ψ q + m ( a i , A i ) i = 1 p ; ( a p + i , 1 / β i ) i = 1 m ( b j , B j ) k = 1 q ; ( b q + i , 1 / β i ) i = 1 m λ z = I ˜ p Ψ q ( a j , A j ) j = 1 p ( b k , B k ) k = 1 q λ z .
In the case of Wright function with arbitrary parameters A p + i , B q + i , i = 1 , , m :
p + m Ψ q + m ( a i , A i ) i = 1 p ;   ( a p + i , A p + i ) i = 1 m ( b j , B j ) k = 1 q ;   ( b q + i , B q + i ) i = 1 m z , such kind of result is presented in [24] by means of more general operators I ˜ , the so-called Wright-Erdélyi–Kober operators. This means that using a suitable number of times of a procedure similar to that in proof of Theorem 3, from any p Ψ q -function (resp. p F q -function) we can go down to one of the three basic generalized hypergeometric functions, depending on if p < q , p = q or p = q + 1 : 0 Ψ q p , 1 Ψ 1 , 2 Ψ 1 ; resp. to: 0 F q p (hyper-Bessel f. and cos m -f.), 1 F 1 (confluent h.f. and exp-f.), 2 F 1 (Gauss f. and beta-distribution of form z α ( 1 z ) β ). This is the reason that we classified the g.h.f. to be of three basic types, as: “g.h.f. of Bessel/cosine type”, “g.h.f. of confluent/exp type” and “g.h.f. of Gauss/beta-distribution type”. Details on this approach and such a classification of the SF can be found in Kiryakova ([6,22,24,50], Ch.4).
Analogously to Theorem 3, we have also a relation (image) for the generalized fractional derivatives of g.h.f., presented as Theorem 2 in Kiryakova [2]. The more general formula (as below) is available in Kiryakova, ([3], Theorem 4.2.).
Theorem 4.
D ( β k ) 1 m , m ( γ k ) 1 m , ( δ k ) 1 m z c p Ψ q ( a 1 , A 1 ) , , ( a p , A p ) ( b 1 , B 1 ) , , ( b q , B q ) λ z μ = z c   p + m Ψ q + m ( a i , A i ) 1 p , ( γ k + δ k + 1 + c β k , μ β k ) 1 m ( b j , B j ) 1 q , ( γ k + 1 + c β k , μ β k ) 1 m λ z μ .
Proof. 
One possible approach to derive this, is to use a decomposition formula for the generalized (multiple) fractional derivatives (23),
D ( β k ) , m ( γ k ) , ( δ k ) f ( z ) = D β m γ m , δ m D β m 1 γ m 1 , δ m 1 D β 1 γ 1 , δ 1 f ( z ) ,
as sequential derivatives. Then, we apply m-times the result (39) of Lemma 2.
We can verify (51) also directly, in the same way as in the end of the proof of Lemma 2, using the basic relation D ( β k ) , m ( γ k ) , ( δ ) k I ( β ) k , m ( γ k ) , ( δ k ) f ( z ) = f ( z ) . □
The case of images under the Caputo type generalized fractional derivatives D ( β k ) 1 m , m ( γ k ) 1 m , ( δ k ) 1 m defined in (24) will be discussed in a separate work, see also Remark 2.

5. Examples of Erdélyi–Kober and Riemann–Liuoville Operators of Some Special Functions

In the beginning of Section 3 we already acknowledged the contributions by some classical authors, such as Erdélyi et al., Askey, Lavoie-Osler-Tremblay, to provide the images of some special and elementary functions under the Riemann–Liuoville fractional integral/derivative, see Formulas (32)–(35). We may refer also to works where detailed tables of images under Riemann–Liuoville operators are provided, for example the book Erdélyi et al. [44], some recent surveys, including in this Journal, such as by Rogosin [18] (as for M-L type functions), Garrappa-Kaslik-Popolizio [52] (images of elementary functions expressed by M-L functions).
As mentioned in Section 3, the proof of Lemma 1, a most general result for Riemann–Liuoville operators of special functions (in sense of g.h.f.) is formula (37) from Kilbas ([16], Th.2):
R δ z c p Ψ q ( a 1 , A 1 ) , , ( a p , A p ) ( b 1 , B 1 ) , , ( b q , B q ) λ z μ = z c + δ p + 1 Ψ q + 1 ( a i , A i ) 1 p , ( c + 1 , μ ) ( b j , B j ) 1 q , ( c + δ + 1 , μ ) λ z μ ,
and for the R-L derivative, the corresponding result is as in Equation (40).
It may be instructive to repeat (as from [2]) some very special cases of the images (37) and (40) under the R-L integral R δ f ( z ) = z δ I 1 0 , δ f ( z ) , that have been derived by the cited authors by the standard term-by-term integration/differentiation. Naturally, these come also as specifications of our results from Lemmas 1 and 2 for the Erdélyi–Kober operators (case m = 1 ).
Example 2.
The R-L fractional integral of the weighted Bessel function, for R e   δ > 0 ,   R e   ν > 1 ,   R e   ( γ + ν ) > 0 , is given by Kilbas-Sebastian ([53], Cor.1, (28)), in the form
R δ z γ 1 J ν ( z ) = z γ + ν + δ 1 2 ν 1 Ψ 2 ( γ + ν , 2 ) ( γ + ν + δ , 2 ) , ( ν + 1 , 1 ) 1 4 z 2 .
We can use J ν ( z ) = 0 Ψ 1 ( ν + 1 , 1 ) 1 4 z 2 , so to see (52) as an immediate corollary of (37), and also of our E-K result (36). Because a 1-tuple fractional calculus operator (R-L, or E-K) is applied, the preliminary expectation is confirmed to have as a result a 0 + 1 Ψ 1 + 1 -function.
Example 3.
The R-L fractional integral of the generalized Bessel function J ν κ (usually called Bessel–Maitland function, a name that should correctly be called the Bessel-Wright function) is derived in Kilbas ([16], Th.8, (26)), and extends the above formula (52): Note the representation J ν κ ( z ) = 0 Ψ 1 ( ν + 1 , κ ) z , then from our result in Lemma 1, and in particular, from (37), it is expected to have the result as a 1 Ψ 2 -function:
R δ z γ 1 J ν κ ( λ z μ ) = z γ + δ 1 1 Ψ 2 ( γ , μ ) ( γ + δ , μ ) , ( ν + 1 , κ ) λ z μ ,
for R e δ > 0 , R e ( γ 1 ) > 1 , κ > 1 , μ > 0 . The sign “minus” in the argument of RHS was missing in [16] due to a possible typo.
The same result, in terms of the (classical) Wright function is presented in the same paper, Kilbas ([16], Th.6, (18)), with the true argument sign (we slightly change the denotations to be similar as in the first row of our (11)),
R δ z γ 1 ϕ ( κ , ν + 1 ; λ z μ ) = z γ + δ 1 1 Ψ 2 ( γ , μ ) ( γ + δ , μ ) , ( ν + 1 , κ ) λ z μ .
A more useful result, in the sense that the R-L integral transforms a generalized Bessel function/resp. Wright function, into same kind of function but with increased index comes if we put γ 1 = ν , κ = μ in (53), see, for example, ([16], Cor.8.1, (28)), wth R e δ > 0 , R e ν > 1 , μ > 0 :
R δ z ν J ν μ ( λ z μ ) = z ν + δ J ν + δ + 1 μ ( λ z μ ) , R δ z ν ϕ ( μ , ν + 1 ; λ z μ ) = z ν + δ ϕ ( μ , ν + δ + 1 ; λ z μ ) ,   λ 0 .
Next, we mention an example with the so-called generalized M-series. Namely, (M.) Sharma and Jain [54] introduced the special function p M α , β q ( z ) , as an extension of both g.h.f. p F q ( z ) and the (2-parameters) M-L function E α , β ( z ) :
p M α , β q a 1 , , a p ; b 1 , s , b q ; z = k = 0 ( a 1 ) k ( a p ) k ( b 1 ) k ( b q ) k z k Γ ( α k + β ) = κ p + 1 Ψ q + 1 ( a 1 , 1 ) , , ( a p , 1 ) , ( 1 , 1 ) ( b 1 , 1 ) , , ( b q , 1 ) , ( β , α ) z .
Here z , α , β C , R e α > 0 , p q are the integer orders, and if p = q + 1 we require additionally that | z | < R = α α , and κ : = j = 1 q Γ ( b j ) / i = 1 p Γ ( a i ) . Usually the following particular cases are always mentioned: (1) β = 1 : this is the (simpler) M-series, introduced by M. Sharma (2008, in same journal as [54]); (2) p = q = 0 (that is, no upper and no lower parameters): this is the M-L function E α , β ( z ) ; (3) p = 0 , q = 1 , b 1 = 1 : one has the Wright function ϕ ( α , β , z ) , or the generalized Bessel–Maitland function; (4) p = q = 1 , a 1 = γ , b 1 = 1 : this is the Prabhakar M-L type function (8), (5) α = β = 1 : we have the g.h.f. p F q ( a 1 , , a p ; b 1 , , b q ; z ) , etc.
Since (56) is a p + 1 Ψ q + 1 -function, all FC operators of the form (20) (and their particular cases as R-L, E-K, Saigo, M-S-M) of the M-series can be evaluated using our formulas in Lemmas 1 and 2 and Theorems 3 and 4.
Example 4.
In [54], the images of the generalized M-series are derived in the case of R-L fractional integral and derivative of order δ > 0 :
R δ p M α , β q a 1 , , a p ; b 1 , , b q ; z = z δ Γ ( 1 + δ ) p + 1 M α , β q + 1 a 1 , , a p , 1 ; b 1 , , b q , 1 + δ ; z ,
D δ p M α , β q a 1 , , a p ; b 1 , , b q ; z = z δ Γ ( 1 δ ) p + 1 M α , β q + 1 a 1 , , a p , 1 ; b 1 , , b q , 1 δ ; z ,
using term-by-term integration/differentiation of the series (56). However, having in mind the representations in both sides as Wright g.h.f., one can get these results directly from the corollaries of Lemmas 1 and 2, the R-L integral (37) and derivative (40). See also in Lavault [55].
The formulas in Theorems 3 and 4 can easily be reduced to corresponding results for generalized fractional integrals and derivatives (20) and (23) of the M-series, to appear in terms of p + m M α , β q + m ( z ) , with additional parameters depending on ( γ k ) 1 m , ( δ k ) 1 m . Again in view of (56), evaluation of other particular FC operators, such as E-K, Saigo, M-S-M, of the M-series can be done. For example, the M-S-M images were evaluated by Kumar and Saxena [56].

6. Saigo Hypergeometric Operators of Various Special Functions

In a series of papers since 1978, such as [57] (for more references see in [6,58]), Saigo introduced a linear integral operator with Gauss function in the kernel, and applied it first for studying BVP for PDE as the Euler-Darboux equation. Later on, this operator was used by him and collaborators in geometric function theory (classes of univalent functions). It happens that, as a case of the hypergeometric integral operators, the Saigo operator has also a role as an FC operator and this has recently become a reason for great interest for researchers in FC, and mainly to authors whose job is to evaluate images of Saigo operator(s) of various special functions. A search in Google for the phrase “Saigo operator” + “function” returns now more than 1060 results (of course some of them may also concern the more general Marichev–Saigo–Maeda, discussed in next Section 7).
First, let us remind the definition and two basic properties of the Saigo operators.
For complex α , β , η and R e α > 0 , the Saigo fractional integration operator (the LHS version) is
I α , β , η f ( z ) = z α β Γ ( α ) 0 z ( z ξ ) α 1 2 F 1 ( α + β , η ; α ; 1 ξ z ) f ( ξ ) d ξ = z β Γ ( α ) 0 1 s ( 1 σ ) α 1 2 F 1 ( α + β , η ; α ; 1 σ ) f ( z σ ) d σ ,
and we skip the discussion on the RHS versions of the Saigo integrals, as similar. The Saigo fractional derivative is used as: D α , β , η f ( z ) = ( d / d z ) n I α + n , β n , η n f ( z ) with n = R e α + 1 . For its explicit differ-integral expression, see for example in ([6], Ch.1). More details can be found in Kiryakova ([6,58], Ch.1, Ch.5) and other our papers dealing with these operators in classes of univalent functions (some of which are joined with Professor Megumi Saigo). A basic formula (known from the original Saigo works) that all authors use (and sometimes derive again) is for the image of a power function:
I α , β , η z p = Γ ( p + 1 ) Γ ( p + η β + 1 ) / Γ ( p β + 1 ) Γ ( p + α + η + 1 ) z p β ,
for R e α > 0 , R e ( p + 1 ) > max [ 0 , R e ( β η ) ] .
As mentioned in Section 2.3, the Saigo operators are cases of the hypergeometric operators of FC, and of the GFC operators for m = 2 , simply because according to (30) the Gauss kernel function is representable as the kernel of (19) and (18) with m = 2 , β = 1 , γ 1 = η β , γ 2 = 0 , δ 1 = η , δ 2 = α + η :
( 1 σ ) α 1 Γ ( α ) 2 F 1 ( α + β , η ; α ; 1 σ ) = G 2 , 2 2 , 0 σ β , α + η η β , 0 .
Thus, the Saigo operator is a generalized (2-tuple fractional integral) of the form (20) and therefore in view of (21), it is also a commutable composition of two classical E-K fractional integrals, see for example ([6], Ch.1):
I α , β , η f ( z ) = z β I ( 1 , 1 ) , 2 ( η β , 0 ) , ( η , α + η ) f ( z ) = z β I 1 η β , η I 1 0 , α + η f ( z ) = I 1 η , η I 1 β , α + η z β f ( z ) = I ( 1 , 1 ) , 2 ( η , β ) , ( η , α + η ) z β f ( z ) = R η z α β R α + η f ( z ) .
The relation between the first and second lines follows by application of the “generalized commutation” between (multiple) Erdélyi–Kober operators and power functions (([6], Ch.1, (1.3.3)), ([34], Th.4), etc.). For particular parameters α , β , η , the Saigo operator can reduce to one E-K operator or an R-L operator, say for β = α , η = 0 it is an R-L integral; and for η = α , an E-K integral: I α , β , α = z β I 1 α β , α .
Therefore, the Saigo image of some special function, which can be represented as a Wright function p Ψ q , can be written as a particular case of the general formulas (46), resp. (51), or also, as a subsequent two-times application of classical E-K operators. Therefore, the Saigo image of a p Ψ q -function can always be predicted to result into a p + 2 Ψ q + 2 -function (unless some parameters in upper and lower rows eliminate each other, and so the indices can be reduced). Our result, as a corollary of Theorem 3 and Corollary 2 states as follows.
Lemma 3.
The images of the Wright g.h.f. p Ψ q , and in particular of the g.h.f. p F q , under the Saigo operator (59) are the same kind of functions with orders increased by 2:
I α , β , η z c p Ψ q ( a i , A i ) 1 p ( b j , B j ) 1 q λ z μ = z c β p + 2 Ψ q + 2 ( a i , A i ) 1 p , ( η β + 1 + c , μ ) , ( 1 + c , μ ) ( b j , B j ) 1 q , ( β + 1 + c , μ ) , ( α + η + 1 + c , μ ) λ z μ ,
(for c = 0 , μ = 1 , this is Cor. 3 in [2]) and
I α , β , η p F q a 1 , , a p ; b 1 , , b q ; λ z = z β p + 2 F q + 2 a 1 , , a p , η β + 1 , 1 ; b 1 , , b q , β + 1 , α + η + 1 ; λ z ,
under the mentioned conditions in the definition of (59).
The following examples for Saigo operators of particular functions from our previous papers [2,3] are repeated here as an illustration for the general result in Lemma 3.
Example 5.
The Saigo fractional integral (59) of a weighted Bessel function was evaluated in Kilbas-Sebastian ([53], Th.1), for R e α > 0 , R e ν > 1 , R e ( γ + ν ) > max 0 , R e ( β η ) :
I α , β , η z γ 1 J ν ( z ) = z γ + ν β 1 2 ν 2 Ψ 3 ( γ + ν , 2 ) , ( γ + η + ν β , 2 ) ( γ + ν β , 2 ) , ( γ + ν + α + η , 2 ) , ( ν + 1 , 1 ) 1 4 z 2 .
To apply (61) from Lemma 3, let us remind the reader again that z γ 1 J ν ( z ) = z γ 1 0 Ψ 1 ( ν , 1 ) 1 4 z 2 and so, the result should be expected to appear as a 0 + 2 Ψ 1 + 2 . Alternatively, to exhibit the use of decomposition of the Saigo operator in two R-L operators (the last relation in (60)) combined with (52) from Example 2, we may proceed as follows:
R η z α β R α + η z γ 1 0 Ψ 1 ( ν , 1 ) 1 4 z 2
= R η z α β 2 ν z γ + ν + α + η 1 1 Ψ 2 ( γ + ν , 2 ) ( γ + ν + α + η , 2 ) , ( ν + 1 , 1 ) 1 4 z 2
= 2 ν R η z γ β + η + ν 1 1 Ψ 2 ( γ + ν , 2 ) ( γ + ν + α + η , 2 ) , ( ν + 1 , 1 ) 1 4 z 2 = 2 Ψ 3 ,   a s   i n   ( 63 ) .
Example 6.
The more special case for Saigo fractional integral of a (weighted) cosine function is the formula from the same paper of Kilbas-Sebastian ([53], Th.5, (47)), for R e α > 0 , R e γ > max 0 , R e ( β η ) :
I α , β , η z γ 1 cos z = π z γ β 1 2 Ψ 3 ( γ , 2 ) , ( γ + η β , 2 ) ( γ β , 2 ) , ( γ + η + α , 2 ) , ( 1 2 , 1 ) 1 4 z 2 .
Note that cos z = π z / 2 J 1 / 2 ( z ) , and use the result (63) of Example 4 with ν = 1 / 2 . To use our general approach, we can present the cos-function as a 0 Ψ 1 -function, and predict the result to be a 0 + 2 Ψ 1 + 2 -function in (64).
Next, we consider a case with a more general special function, called generalized K-series. In [59] (K.) Sharma introduced an extension of both a g.h.f. p F q ( z ) and Prabhakar (three-parameter Mittag-Leffler) function E α , β γ ( z ) (see (8)):
p K α , β ; γ q a 1 , , a p ; b 1 , , b q ; z : = p K α , β ; γ q ( z ) = k = 0 ( a 1 ) k ( a p ) k ( b 1 ) k ( b q ) k ( γ ) k z k Γ ( α k + β ) ,
with z , α , β C , R e α > 0 , integers p q (and additional requirement | z | < R = α α if p = q + 1 ). When γ = 1 it reduces to the (generalized) M-series (56) by Sharma-Jain [54], Example 4.
Example 7.
Recently, Lavault [55] represented the above K-series in terms of a Wright g.h.f.:
p K α , β ; γ q a 1 , , a p ; b 1 , , b q ; z = j = 1 q Γ ( b j ) Γ ( γ ) i = 1 p Γ ( a i ) p + 2 Ψ q + 2 ( a 1 , 1 ) , ( a p , 1 ) , ( γ , 1 ) , ( 1 , 1 ) ( b 1 , 1 ) , , ( b q , 1 ) , ( 1 , 1 ) , ( β , α ) z ,
and calculated some of its FC operators, as the R-L, Saigo and M-S-M operators. As should be expected, the image of a p K α , β ; γ q -function under the R-L integral is a p + 1 K α , β ; γ q + 1 -function (Th. 4.1 there), similarly to Example 4 for the M-series.
The Saigo operator is also derived by Lavault in [55]: for the M-series—in Th. 4.2, and for the K-series—in Cor. 4.3. Namely, Equation (4.10), [55] reads as:
I α , β , γ t σ 1 p K ξ , η ; ν q ( c z μ ) = 1 q Γ ( b j ) 1 p Γ ( a i ) z σ β 1 Γ ( ν ) p + 3 Ψ q + 3 ( a i , 1 ) 1 p , ( σ , μ ) , ( β + γ + σ , μ ) , ( ν , 1 ) ( b j ) 1 q , ( β + σ , μ ) , ( α + γ + σ , μ ) , ( η , ξ ) c z μ .
Let us note that the K-series is a p + 2 Ψ q + 2 -function (66), and from our Lemma 3 the expected result should be a p + 4 Ψ q + 4 -function, with indices increased by two. However, pairs of upper and lower rows’ parameters appear the same and eliminate each other, therefore the result reduces to a p + 3 Ψ q + 3 , as above.

7. Marichev–Saigo–Maeda (M-S-M) Operators of Various Special Functions

As mentioned in Section 2.3, there is an interesting particular case of the GFC operators (20) and (23) for m = 3 , often abbreviated as M-S-M (MSM) operators. These operators have also become very popular in works dedicated to evaluate FC images of special functions. A search in Google for “Marichev–Saigo–Maeda” returns at least 2430 results, and for “MSM operator”—some 2670 results.
This operator appeared in a paper by Marichev of 1974, [39], see also in the book ([31], Section 8.4.51); and further was introduced and studied by Saigo, Saigo and Maeda in 1996, see [40], also by Saigo and Saxena (1996, 1998, 2001), details on references are in ([6,37,58,60], Ch.1), etc.
For complex parameters a , a , b , b , c , R e c > 0 , the Marichev–Saigo–Maeda (M-S-M) integral operator, of which the kernel is the Appel function, or Horn’s function F 3 (see ([15], Vol.1), also [7])
F 3 a , a , b , b , c , z , ξ = m , n = 0 ( a ) m ( a ) n ( b ) m ( b ) n ( c ) m + n z m ξ n m ! n ! , | z | < 1 , | ξ | < 1 ,
is defined as the linear integral operator
I a , a , b , b , c f ( z ) = z a Γ ( c ) 0 z ( z ξ ) c 1 ξ a F 3 ( a , a , b , b ; c ; 1 ξ z , 1 z ξ ) f ( ξ ) d ξ = z c a a 0 1 ( 1 σ ) c 1 Γ ( c ) σ a F 3 ( a , a , b , b ; c ; 1 σ , 1 1 σ ) f ( z σ ) d σ .
Observing the representation (31) of the kernel F 3 -function as a kernel of the generalized fractional integrals (20) (see Section 2.3), it is evident that the M-S-M operator is nothing but their special case for m = 3 . Then, in view of (21), it is also a composition of three commutable classical E-K integrals (see Kiryakova [6,37,58]). This fact seems to be unknown to the other authors (or is continuously ignored):
I a , a , b , b , c f ( z ) = z c a a 0 1 σ a G 3 , 3 3 , 0 σ a + b , c a , c b a , b , c a b f ( z σ ) d σ = z c a a 0 1 G 3 , 3 3 , 0 σ a a + b , c 2 a , c a b a a , b a , c 2 a b f ( z σ ) d σ = z c a a 0 1 H 3 , 3 3 , 0 σ ( a a + b , 1 ) , ( c 2 a , 1 ) , ( c a b , 1 ) ( a a , 1 ) , ( b a , 1 ) , ( c 2 a b , 1 ) f ( z σ ) d σ = z c a a I ( 1 , 1 , 1 ) , 3 ( a a , b a , c 2 a b ) , ( b , c a b , a ) f ( z ) = z c a a I 1 a a , b I 1 b a , c a b I 1 c 2 a b , a f ( z ) .
The relations (31) and (69) have been recently denied and argued in the Response of authors [61] to our critical Commentary [5] to their paper [62]. Then, I needed to support (by my footnote remark to [61]) the truth of (69) as appearing also in the basic FC book by Samko-Kilbas-Marichev [31], see there Equation (10.38) (for decomposition of Saigo operator) and Equation (10.46), p. 193 (for decomposition of the M-S-M operator).
For the above reasons, to evaluate M-S-M images of special functions, which are representable as Wright g.h.f., one can use the general result of Theorem 3. Thus we have:
Lemma 4.
The image of a Wright g.h.f. under the M-S-M fractional integral is given by the formula
I a , a , b , b , c z ν p Ψ q ( a i , A i ) 1 p ( b j , B j ) 1 q λ z μ = z c a a p + 3 Ψ q + 3 ( a i , A i ) 1 p , ( a a + 1 + ν , 1 ) , ( b a + 1 + ν , 1 ) , ( c 2 a b + 1 + ν , 1 ) ( b j , B j ) 1 q , ( a a + b + 1 + ν , 1 ) , ( c 2 a + 1 + ν , 1 ) , ( c a b + 1 + ν , 1 ) λ z μ .
The corresponding simpler result for ν = 0 is given by Corollary 4 in Kiryakova [2].
The M-S-M fractional derivatives D a , a , b , b , c , denoted also by I a , a , b , b , c with R e c 0 , are considered by Saigo and Maeda and by the next authors as originally defined by analogy with the Saigo derivatives D α , β , η . In view of (69), they can be considered also as special cases of the generalized fractional derivatives (23) with m = 3 , namely as D a , a , b , b , c = D ( 1 , 1 , 1 ) , 3 ( a a , b a , c 2 a b ) , ( b , c a b , a ) z c .
The authors after Saigo-Maeda use to derive first a formula for the M-S-M image of a power function z p , ignoring the fact that it exists in the original paper (1996) (and follows also as a particular case of our (27) in Section 2). Then, to find the M-S-M fractional integral or derivative of a particular special function, they use the standard techniques of term-by-term integration/or differentiation of the corresponding powers series. However, our general approach says that we know in advance the image of a p Ψ q -function expected as a p + 3 Ψ q + 3 , see (70).
We provide a few illustrative examples for other authors’ results, mentioned also in Kiryakova [2].
Example 8.
The formula for the M-S-M generalized fractional integral of a weighted Bessel function:
I a , a , b , b , c z γ 1 J ν ( z ) = z γ + ν a a + c 1 2 ν × 3 Ψ 4 ( γ + ν , 2 ) , ( γ + ν + c a a b , 2 ) , ( γ + ν + b a , 2 ) ( γ + ν + b , 2 ) , ( γ + ν + c a a , 2 ) , ( γ + ν + c a b , 2 ) , ( ν + 1 , 1 ) z 2 4 ,
can be found in Purohit-Suthar-Kalla ([60], Th.2.1, (10)). It is supposed that R e c > 0 , R e ν > 1 , R e ( γ + ν ) > max 0 , R e ( a + a + b c ) , R e ( a b ) . The same result, however, can be obtained in the way as discussed in Example 5, using the M-S-M operator’s representation (69) and the result from Lemma 4. Then, by analogy with Example 6, one can derive the particular result from ([60], Cor.3.1, (24)) for the M-S-M image of z γ 1 cos z , again in terms of 3 Ψ 4 ( z 2 4 ) .
Example 9.
Mondal-Nisar ([63], Th.3, (11)) evaluated the M-S-M integral (68) of the so-called generalized Bessel function
W p , β , γ ( z ) = k = 0 ( 1 ) k γ k Γ ( p + β 2 + 1 2 + k ) k ! z 2 2 k + p .
Evidently, it is a variant (up to variable substitution) of the Bessel–Maitland-Wright function J ν κ ( z ) and of the Wright function ϕ ( z ) , representable as 0 Ψ 1 -function of ( z 2 / 4 ) . Then as well expected, the result comes as a 3 Ψ 4 -function, since the indices are increased by 3.
Example 10.
Next, in Nisar-Mondal-Agrawal ([64], Th.1), the authors derive the M-S-M operator of the Bessel-Struve function, which is representable as a 2 × 2 -indices (multi-index, m = 2 ) Mittag-Leffler function (9), see Examples (11) mentioned in Section 2.1, as well as a Wright g.h.f. 1 Ψ 1 ,
S ν ( z ) = Γ ( ν ) π k = 0 Γ ( 1 2 + k 2 ) Γ ( ν + 1 + k 2 ) z k k ! = Γ ( ν ) π 1 Ψ 1 ( 1 2 , 1 2 ) ( ν + 1 , 1 2 ) z .
Then, the result for I a , a , b , b , c t γ 1 S ν ( λ z ) is expected, written in terms of a 1 + 3 Ψ 1 + 3 ( λ z ) = 4 Ψ 4 ( λ z ) , with parameters following from the general scheme.
Example 11.
The M-S-M operator (68) of a generalized multi-index Mittag-Leffler function
E ( α j , β j ) 1 m γ , κ ( z ) = k = 0 ( γ ) κ k z k j = 1 m Γ ( α j k + β j ) k ! = 1 Γ ( γ ) 1 Ψ m ( γ , κ ) ( β j , α j ) 1 m z ,
is handled in Agarwal-Rogosin-Trujillo [29]. When m = 1 it was studied also by Srivastava-Tomovski [65]. Note that for γ = κ = 1 the above function reduces to the ( 2 m ) multi-index Mittag-Leffler function (9). This appeared also in Saxena-Nishimoto [66] and was studied in Saxena-Pogany-Ram-Daiya [67]. The result from ([29], Th.3.1, (3.2)) is the following:
I a , a , b , b , c z ρ 1 E ( α j , β j ) 1 m γ , κ ( λ z μ ) = z ρ + c a a 1 Γ ( γ ) × 4 Ψ m + 3 ( γ , κ ) , ( ρ , μ ) , ( ρ + c a a b , μ ) , ( ρ + b a , μ ) ( α j , β j ) 1 m , ( ρ + b , μ ) , ( ρ + c a a , μ ) , ( ρ + c b a , μ ) λ z μ .
Using the representations of the M-S-M operator as three-tuple generalized fractional integral (69) and of this special function as a 1 Ψ m -function, the same formula can be evaluated by the general result in Theorem 3, that is the image is again a Wright g.h.f. but its indices are increased by three.
Example 12.
We were stuck on a paper by Kumar-Gupta-Rawat [68] (very fast accepted and published with a lot of typographical problems). The authors there aim to “establish certain generalized fractional differentiation involving M-L type function with four parameters, recently introduced by Garg et al. (2016)”. Namely, they have evaluated its image under the Marichev–Saigo–Maeda derivative D a , a , b , b , c , corresponding to the integral operator (68). Their result, Theorem 1 (p.205), reads as follows:
D a , a , b , b , c t ρ 1 ξ , γ E μ , ν ( λ z σ ) = z ρ + a + a c 1 Γ ( ξ )
× 5 Ψ 4 ( ξ , γ ) , ( ρ , σ ) , ( ρ + a + a + b c , σ ) , ( ρ + a b , σ ) , ( 1 , 1 ) ( ρ b , σ ) , ( ρ + a + a c , σ ) , ( ν , μ ) , ( ρ + a + b c , σ ) λ z σ .
The authors did not observe the fact that the considered M-L type function is a case of the generalized Wright function (4), the definition of which is also given in the mentioned paper, namely:
ξ , γ E μ , ν ( z ) = k = 0 ( ξ ) γ k Γ ( μ k + ν ) z k = 1 Γ ( ξ ) k = 0 Γ ( ξ + γ k ) Γ ( k + 1 ) Γ ( μ k + ν ) z k k ! = 1 Γ ( ξ ) 2 Ψ 1 ( ξ , γ ) , ( 1 , 1 ) ( ν , μ ) z .
Then, its image under the M-S-M differentiation, as a three-tuple generalized fractional derivative, is well expected to be the Wright function 2 + 3 Ψ 1 + 3 , in view of our general results as (51) and (70).
Example 13.
In [62], Agarwal-Jain-Baleanu considered the M-S-M images of the generalized Lommel-Wright function
J ω , θ φ , m ( z ) = ( z 2 ) ω + 2 θ k = 0 ( 1 ) k ( z 2 ) 2 k Γ ( θ + k + 1 ) m Γ ( ω + k φ + θ + 1 ) = ( z 2 ) ω + 2 θ 1 Ψ m + 1 ( 1 , 1 ) ; ( θ + 1 , 1 ) , , ( θ + 1 , 1 ) , ( ω + θ + 1 , φ ) ; z 2 / 4 ,
which is a Wright g.h.f. (see Equation (1.1) there). We can note that it is also example of the multi-index M-L function (9), namely J ω , θ φ , m ( z ) = ( z 2 ) ω + 2 θ ( z 2 ) ω + 2 θ E ( 1 , , 1 , φ ) , ( θ + 1 , , θ + 1 , ω + θ + 1 ) ( m + 1 ) ( z 2 ) 2 . Then the result, as calculated by the authors, follows directly from Theorem 3 and especially from Lemma 4 (below, A : = χ + ω + 2 θ , φ > 0 ):
I 0 + ξ , ξ , ρ , ρ , ϰ t χ 1 J ω , θ φ , m ( t z ) ( x ) = x A ξ ξ + ϰ 1 ( z 2 ) ω + 2 θ × 4 Ψ 4 + m ( A , 2 ) , ( A + ϰ ξ ξ ρ , 2 ) , ( A + ρ ξ , 2 ) , ( 1 , 1 ) ( A + ρ , 2 ) , ( A + ϰ ξ ξ , 2 ) , ( A + ϰ ξ ρ , 2 ) , ( ω + θ + 1 , φ ) , ( θ + 1 , 1 ) ( t z ) 2 4 ,
to be again a Wright g.h.f. but with indices increased by three, that is, a 4 Ψ m + 4 -function. In [62] also many special cases are derived, such as Beta-transform (that is E-K integral), Saigo operator, path integral, of the function (73) and of its particular cases. As in the Commentary [5] we discussed the possibilities to use our unified approach, the authors tried to argue with the facts in their Response [61]. The curious readers are recommended to read Kiryakova’s footnote comments at the bottom to this Response [61].

8. Multiple Gel’fond-Leontiev Operators of Multi-Index Mittag-Leffler Functions; Hyper-Bessel Operators and Functions

We consider now GFC images of the multi-index M-L functions (9).
Lemma 5.
Taking in general m n (m-tuple operators of GFC and 2 n -indexed M-L functions), we have
I ( β k ) , m ( γ k ) , ( δ k ) z c E ( α i ) , ( ν i ) ( n ) ( λ z μ ) = I ( β k ) , m ( γ k ) , ( δ k ) z c 1 Ψ n ( 1 , 1 ) ( ν i , α i ) 1 n ) λ z μ = z c 1 + m Ψ n + m ( 1 , 1 ) , ( γ k + 1 + c / β k , μ / β k ) 1 m ( ν i , α i ) 1 n , ( γ k + δ k + 1 + c / β k , μ / β k ) 1 m λ z m .
This is an easy corollary of Theorem 3. In particular, for c = 0 , μ = 1 , m = n , and for GFC parameters taken to be γ k = ν k 1 , β k = 1 / α k , k = 1 , 2 , , m , it happens that the parameters ( γ k + 1 , 1 / β k ) 1 m in the upper row and ( ν i , α i ) 1 m in bottom row appear equal and cancel each other, and then the W. g.h.f. 1 + m Ψ m + m reduces to 1 Ψ m , again a multi-index M-L function!
Then, as proved (in other direct way) in our previous papers, we have:
Example 14
(Kiryakova, [21,22]). For each fixed j = 1 , , m , a (classical) E-K integral of (9) reads as:
I 1 / α j β j 1 , δ j E ( α i ) , ( β i ) ( λ z ) = E ( α i ) , ( β 1 , , β j 1 , β j + δ j , β j + 1 , , β m ) ( λ z ) .
This is an extension of the result for an E-K integral (12) of an M-L function: I 1 / α β 1 , δ E α , β ( z ) = E α , β + δ ( z ) , where its second index is increased by the order of fractional integral. After m-times application of the above relation with respect to each j = 1 , , m , we obtain that a GFC integral (18) with suitably chosen parameters transforms a multi-index M-L function into the same kind of multi-index M-L function of which the indices of the second set are increased by the multi-order of fractional integration:
I ( 1 / α k ) , m ( β k 1 ) , ( δ k ) E ( α i ) 1 m , ( β i ) 1 m ( λ z ) = E ( α i ) 1 m , ( β i + δ i ) 1 m ( λ z ) ,   w i t h   R e δ k > 0 , γ k > 1 , α k > 0 , k = 1 , , m , λ 0 .
If we take δ k = α k , k = 1 , , m , formula (77) has the form (Kiryakova, [21,22])
I ( 1 / α k ) , m ( β k 1 ) , ( α k ) E ( α i ) , ( β i ) ( λ z ) = E ( α i ) , ( β i + α i ) ( λ z ) = ( λ z ) 1 E ( α i ) , ( β i ) ( λ z ) 1 i = 1 m Γ ( β i ) .
According to the operational rules of the GFC (([6], Ch.5)) and rewriting the above relation for the generalized fractional derivative D ( 1 / α k ) , m ( β k 1 α k ) , ( α k ) defined as in (23), we have
D ( 1 / α k ) , m ( β k 1 α k ) , ( α k ) E ( α i ) , ( β i ) ( λ z ) = ( λ z ) E ( α i ) , ( β i ) ( λ z ) + i = 1 m Γ ( β i α i ) 1 .
Here one can see an analogy with the results (10.6), (10.9) from Haubold-Mathai-Saxena [17] for the R-L operators (in the case m = 1 ).
Next, let us consider the special cases of GFC operators for which the multi-index M-L functions (9) appear as eigenfunctions, that is, these special functions are transformed into the same kind of functions with the same multi-indices.
Example 15.
The so-called Gelfond–Leontiev (G-L) operators are operators of generalized integration and differentiation, defined for functions f ( z ) = j = 0 a j z j analytic in a disk | z | < R , and are generated by the coefficients of a given entire function φ ( σ ) , used as a multipliers’ sequence. They were introduced in a paper by Gelfond and Leontiev of 1951 (for details and references see our works, such as [6,20,21]). In the case when φ ( σ ) = E ( α i ) , ( β i ) ( σ ) is the multi-index M-L function (9), these operators were considered by Kiryakova [20], see also [21], etc., and called (multiple) Dzrbashjan–Gelfond–Leontiev (D-G-L) operators. We defined them as follows:
D f ( z ) = j = 1 a j Γ ( β 1 + j α 1 ) Γ ( β m + j α m ) Γ ( β 1 + ( j 1 ) α 1 ) Γ ( β m + ( j 1 ) α m ) z j 1 , L f ( z ) = j = 0 a j Γ ( β 1 + j α 1 ) Γ ( β m + j α m ) Γ ( β 1 + ( j + 1 ) α 1 ) Γ ( β m + ( j + 1 ) α m ) z j + 1 ,
and noted that the image functions D f ( z ) , L f ( z ) are also analytic functions in the same disk | z | < R .
We have shown (e.g., in [21,22]) that the operators (78) can be analytically extended (outside a disk, to holomorphic functions in starlike domain) to operators of GFC, namely to generalized integrals and derivatives of fractional multi-order ( α 1 , , α m ) as follows:
L f ( z ) = z 1 I ( 1 / α k ) , m ( β k 1 ) , ( α k ) f ( z ) ,   D f ( z ) = z 1 D ( 1 / α k ) , m ( β k 1 α k ) , ( α k ) f ( z ) z 1 f ( 0 ) k = 1 m Γ ( β k ) Γ ( β k α k ) .
Then, as proved in [21,22], etc., and seen also in the end of Example 14, the multi-index M-L function (9) is a solution of the differential equation of fractional multi-order
D E ( α i ) , ( β i ) ( λ z ) = E ( α i ) , ( β i ) ( λ z ) ,   t h a t   i s ,   E ( α i ) , ( β i ) ( z )   i s   a n   e i g e n f u n c t i o n   o f   t h e   o p e r a t o r   D .
In view of this relation, the multi-index M-L function serves as an eigenfunction of the D-G-L differentiation generated by its own coefficients!
In the case of the 3 m -parametric M-L type functions (multi-index Prabhakar functions, [26])
E ( α i ) , ( β i ) ( γ i ) , m ( z ) = k = 0 ( γ 1 ) k ( γ m ) k Γ ( α 1 k + β 1 ) Γ ( α m k + β m ) z k ( k ! ) m ,   with   Pochhamer   symbols   ( γ i ) k : = Γ ( γ i + k ) Γ ( γ i ) ,
the R-L, classical E-K operators and some multiple E-K operators are evaluated in the works of Paneva-Konovska, for example, the book [19].
Example 16.
In 1966, and his later works, Dimovski [41] introduced a very general class of differential operators of arbitrary (integer) order generalizing the Bessel operators of second order. His aims were to develop operational calculus for these operators, both via a Laplace-type integral transform and by the Mikusinski algebraical approach. These operators have the alternative representations
B f ( t ) = t α 0 d d t t α 1 d d t t α m 1 d d t t α m f ( t ) = t β P m t d d t f ( t ) = t β k = 1 m t d d t + β γ k f ( t ) ,   t > 0 ,
with arbitrary parameters α 0 ; α k , γ k , k = 1 , , m ; β > 0 ; P m a polynomial of degree m, and their different cases were studied by many authors as appearing in various equations of mathematical physics, problems in analysis, etc., disciplines. The name “hyper-Bessel differential operator” for (80) was introduced by Kiryakova in the further studies on the topic, for example ([6], Ch.3), and next ones as [69]. For the linear right inverse operator denoted by L and called hyper-Bessel integral operator (such that B L f ( t ) = f ( t ) ), we have found a representation by an integral operator with Meijer’s G m , m m , 0 -function in the kernel, and later, the same kind of integral representation also for its fractional powers L λ , λ > 0 . These results were the hint of how to introduce the operators of GFC: the generalized integration and differentiation (20) and (26) of arbitrary fractional multi-order ( δ 1 , δ 2 , , δ m ) instead of the multi-order ( λ , λ , , λ ) for L λ . The story is explained in [69]. Due to the representation of the hyper-Bessel operators in the form: B = t 1 D ( β , β , , β ) , m ( γ k 1 ) , ( 1 , 1 , , 1 ) , L = t I ( β , β , , β ) , m ( γ k ) , ( 1 , 1 , , 1 ) , these operators are important examples of the generalized “fractional” derivatives and also of the Gelfond–Leontiev operators (78). Indeed, for simplicity we take β = 1 and γ m = 0 , then we have that B = t 1 D ( 1 ) 1 m , m ( γ k 1 ) 1 m , ( 1 ) 1 m is a particular case of the operator denoted by D in the previous Example 15, with modified denotations.
Consider the m-th order (that is, of multi-order ( 1 , 1 , , 1 ) ) hyper-Bessel differential equation
B   y ( t ) = λ   y ( t ) ,   λ 0 .
In ([6], Ch.3, Th.3.4.3) and the next corollaries (see also [37]), we proved that the functions, j = 1 , , m :
y j ( t ) = G 0 , m 1 , 0 λ t γ j , γ 1 , , γ j 1 , γ j + 1 , , γ m 1 , 0 = ( λ t ) γ j k = 1 m Γ ( γ k + 1 ) 0 F m 1 ( 1 + γ i γ j ) i j ; λ t : = J ( 1 + γ i γ j ) i j ( m 1 ) ( λ t ) ,   t h e   h y p e r - B e s s e l   f u n c t i o n s ,
form a fundamental system of solutions of equation (81) in a neighborhood of origin t = + 0 . Under the assumptions of Th.3.4.3 in [6], γ 1 < < γ m < γ 1 + 1 and γ m = 0 , we have that γ j ( 0 , 1 ) for all j = 1 , , m , and one of these solutions, for j = m , can be written as
y m ( t ) = k = 1 m ( γ k + 1 ) 1   0 F m 1 ( 1 + γ i ) 1 m 1 ; λ t = k = 1 m ( γ k + 1 ) 1 1 F m 1 ; ( 1 + γ i ) 1 m 1 , 1 ; λ t = 1 Ψ m ( 1 , 1 ) ( 1 + γ i ) 1 m 1 , 1 λ z = E ( 1 , 1 , , 1 ) , ( 1 + γ i ) 1 m ( λ z ) ,   a   c a s e   o f   t h e   m u l t i - i n d e x   M - L   f u n c t i o n s .
Therefore, Example 16 appears a special case of Example 15, and shows that the multi-index Mittag-Leffler functions (9) can be seen also as “fractional indices” analogs, extensions of the hyper-Bessel functions (82), which themselves are multi-index variants of the classical Bessel function.

9. Some “New” Special Functions and Their FC Images

Recently, some authors claimed to introduce and consider “new” SF. Among these are examples of the so-called k-analogs of the Bessel and Mittag-Leffler functions, some generalized multi-index Bessel and Mittag-Leffer functions, and some S-functions. The mentioned k-analogs are based on the use of the k-Γ-function, which, however, can be rewritten in terms of the “classical” Γ -function:
Γ k ( s ) = 0 exp ( t k k ) t s 1 d t = k s k 1 Γ ( s k ) ,   s C , Re ( s ) > 0 , k > 0 ;   Γ ( . )   the   Gamma - function .
Usually, the denotations include also the k-analogs of the Pochhamer symbol:
( η ) ν , κ : = Γ k ( η + ν κ ) / Γ k ( λ ) ,   η C \ { 0 } , ν C ,
and in view of (83) are representable again by means of classical Gamma-functions.
Then, one can easily observe that such “new SF” are just cases of the Wright generalized hypergeometric function p Ψ q . Therefore, all the results provided by the mentioned authors to evaluate FC operators of these special functions follow from our general ones, say from Theorems 3 and 4, or the special cases as Lemmas 1 and 2 (for E-K operators, incl. R-L ones), Lemma 3 (for Saigo operators), Lemma 4 (for M-S-M operators), and so on. As an illustration, we repeat some examples from Kiryakova [4].
Example 17.
A generalization of the Bessel function, called generalized k-Bessel function was introduced by Gehlot [70] and studied by Mondal [71], Shaktawat et al. [72], defined as
W ν , c k ( z ) = n = 0 ( c ) n Γ k ( n k + ν + k ) · ( z / 2 ) 2 n + ν k n ! ,   z C , k > 0 , R e ( ν ) > 1 , c C .
Lets us note that this function is practically a Wright g.h.f. 0 Ψ 1 , and even a simpler g.h.f. 0 F 1 of the same type as the classical Bessel function:
W ν , c k ( z ) = ( z / 2 ) ν / k n = 0 [ c ( z / 2 ) 2 ] n k n + 1 + ( ν / k ) Γ ( n + 1 + ( ν / k ) ) Γ ( n + 1 ) = ( z / 2 ) ν / k k 1 + ( ν / k ) n = 0 [ c ( z / 2 ) 2 ] n k n Γ ( 1 + ( ν / k ) + n . 1 ) Γ ( 1 + n . 1 ) = ( z / 2 ) ν / k k 1 + ( ν / k ) n = 0 [ ( c / k ) ( z / 2 ) 2 ] n Γ ( 1 + ( ν / k ) + n . 1 ) Γ ( 1 + n . 1 ) = ( z / 2 ) ν / k k 1 + ( ν / k ) 1 Ψ 2 ( 1 , 1 ) ( 1 + ν k , 1 ) , ( 1 , 1 ) c k z 2 2 = ( z / 2 ) ν / k k 1 + ( ν / k ) 0 Ψ 1 ( 1 + ν k , 1 ) c k z 2 2 = ( z / 2 ) ν / k k 1 + ( ν / k ) Γ ( 1 + ν ) 0 F 1 ; 1 + ν k ; c k z 2 4 .
Naturally, for k = 1 , c = 1 , (85) becomes the classical Bessel function:
W ν , 1 1 ( z ) = ( z / 2 ) ν Γ ( 1 + ν ) 0 F 1 ; 1 + ν ; z 2 4 .
In the case c = 1 , Gehlot [70] considered (85) as a solution of a k-Bessel differential equation. Mondal [71] studied its properties for complex c C . Shaktawat et al. [72] evaluated the M-S-M operators of FC of this function. In view of Lemma 4, the result there (Th.1, (18)) is well expected to appear in terms of the 3 Ψ 4 -function (because the 3-tuple FC integral increases by three the indices of the initial 0 Ψ 1 -function).
Example 18.
The simplest k-analogs of the M-L function are considered by Dorrego-Cerruti [73] and Gupta and Parihar [74], and very recently (2020/2021) studied also by Ali et al. [47]:
E k , α , β ( z ) = n = 0 z n Γ k ( α n + β ) ,   r e s p .   E k ; ν , ρ δ = n = 0 ( δ ) η , k z n Γ k ( ν n + ρ ) n ! .
Various further extensions appeared, as the generalized k-Mittag-Leffler function, studied by Gupta and Parihar [74] and Nisar-Eata-Dhaifalla-Choi [75] in the form (note that the index p was missing in these authors’ original denotation):
E κ , α , β η , δ , p , q ( z ) : = n = 0 ( η ) q n , κ Γ k ( α n + β ) ( δ ) p n , κ z n ,   κ , p , q R + ; α , β , η , δ C ,
with
m i n { R e   ( α ) , R e   ( β ) , R e   ( η ) , R e   ( δ ) } > 0 ;   q R e   ( α ) + p ;
the k-Pochhammer symbol as in (84). Again, the function (87) can be rewritten as a Wright g.h.f., now as 2 Ψ 2 . Using the representations for (83) and (84) we have, respectively:
( η ) q n , κ = = k q n κ / k Γ ( η + q n κ k ) / Γ ( η k ) ;   ( δ ) p n , κ = = k p n κ / k Γ ( δ + p n κ k ) / Γ ( δ k ) ;
a n d   Γ k ( α n + β ) = k α n + β / k · k 1 Γ ( α n + β k ) .
Then,
E κ , α , β η , δ , p , q ( z ) = k 1 β k Γ ( δ / k ) Γ ( η / k ) n = 0 Γ ( η k + n · q κ k ) Γ ( 1 + n . 1 ) Γ ( δ k + n · p κ k ) Γ ( β k + n · α k ) · k q κ p κ α k z n n !
= k 1 β k Γ ( δ / k ) Γ ( η / k ) 2 Ψ 2 ( η k , q κ k ) , ( 1 , 1 ) ( δ k , p κ k ) , ( β k , α k ) k ( q p ) κ α k z .
By the standard techniques, Nisar-Eata-Dhaifalla-Choi [75] evaluated FC operators of the functions (87). In view of our general results, as expected, the results are 5 Ψ 5 -functions (for the MSM operators, Ths. 1–2, 3–4) there), and in particular, 4 Ψ 4 -functions (for the Saigo operators, Cor. 3.1–3.2, there). Also, the pathway integrals (that are related to E-K integrals) are calculated.
Example 19.
The so-called multi-index Bessel function:
J ( β j ) m , κ , b ( α j ) m , γ , c ( z ) = k = 0 c k ( γ ) κ k j = 1 m Γ ( α j k + β j + b + 1 2 ) · z k k ! ,   m = 1 , 2 , 3 , ,
with the Pochhammer symbol ( γ ) κ k , were introduced and studied in a series of papers by Nisar at al., see, for example, Nisar-Purohit-Parmar [76]. The authors proposed a result for the R-L fractional integral of (88), unfortunately written wrongly in their Theorem 1, Equation (2.4) of [76] as a 2 Ψ 2 -function, although it is evidently a 2 Ψ m + 1 -function having ( m + 1 ) parameters in the low row. The true result should be
R λ t δ 1 J ( β j ) m , ϰ , b ( α j ) m , γ , c ( z ) = 1 Γ ( γ ) z δ + λ 1 2 Ψ m + 1 ( γ , ϰ ) , ( δ , 1 ) ( β j + b + 1 2 , α j ) 1 m , ( λ + δ , 1 ) ; c z .
However, it is easily seen that (88) is: J ( β j ) m , κ , b ( α j ) m , γ , c ( z ) = 1 Γ ( γ 1 Ψ m ( γ , ϰ ) ( β j + b + 1 2 , α j ) j = 1 m ; c z , so this result follows directly from our Lemma 1. Note that the function (88) is also a special case of the generalized multi-index M-L function in Example 11 with β j β j + ( b + 1 ) / 2 , and then the results for its images under FC operators follow from these in Agarwal-Rogosin-Trujillo [29].
Very recently (published 24 September 2020), in Mubeen etal. [48], the authors considered integral transforms, including FC operators, of yet more general SF called “extended generalized multi-index Bessel function" introduced by Kamarujjama-Khan-Khan (2019) with an additional member ( δ ) k in the denominator of the series, as:
J ( β j ) m , κ , b , δ ( α j ) m , γ , c ( z ) = k = 0 ( γ ) κ k ( c z ) k ( δ ) k j = 1 m Γ ( α j k + β j + b + 1 2 ) .
Following similar manipulations as we did in [4] (Section 5.3, Equation (48)) for the case of (88), one can show that
J ( β j ) m , κ , b , δ ( α j ) m , γ , c ( z ) = Γ ( δ ) Γ ( γ ) 2 Ψ m + 1 ( γ , κ ) , ( 1 , 1 ) ( δ , 1 ) , ( β j + b + 1 2 , α j ) 1 m c z ,
and is evidently reducible to (88) for δ = 1 . Therefore, the M-S-M fractional integral will be a 5 Ψ m + 4 -function—presented as an explicit SF, instead of the authors’ hardly visible result in form of some unknown complicated series, compared with Th.5.3, [48] for the extension (1.17) of the function (90).
Example 20.
The S-function was introduced in Saxena-Daiya [77] as a “new” special function extending the M-L function ( p = q = 0 , k = 1 ), the Prabhakar function (8), the M-series (56) of Sharma and Jain [54] ( γ = 1 , k = 1 ), etc., by
S [ z ] : = S ( p , q ) α , β , γ , τ , k ( a 1 , , a p ; b 1 , , b q ; z ) = n = 0 ( a 1 ) n ( a p ) n · ( γ ) n τ , k ( b 1 ) n ( b q ) n · Γ k ( n α + β ) · z n n ! ,
with k R , ; α , β , γ , τ C ; R e ( α ) > 0 ; R e ( α ) > k , R e ( τ ) , p < q + 1 . For p = q = 0 it reduces to the generalized k-Mittag-Leffler function E k , α , β γ , τ ( z ) , see in Example 18, the simplest case by Gupta and Parihar [74].
However, as shown in Kiryakova ([4], Section 5.4), this S-function (91) appears to be a Wright g.h.f. (4) of the form p + 1 Ψ q + 1 z k τ α k , namely:
S [ z ] = k 1 β k Γ ( b 1 ) Γ ( b q ) Γ ( a 1 ) Γ ( a p ) · Γ ( γ k ) p + 1 Ψ q + 1 ( a 1 , 1 ) , , ( a p , 1 ) , ( γ k , τ ) ( b 1 , 1 ) , , ( b q , 1 ) , ( β k , α k ) ; z k τ α k .
Unfortunately, this fact has not been observed neither by the authors of [77] introducing it, nor by their numerous followers. Then, all results for images of FC operators, such as R-L, E-K, Saigo, M-S-M, the Euler-transform (which is in fact E-K operator), Laplace transform, follow as images of the Wright function according to our general results, say Theorem 4.1. Then, as evaluated in [77], Th. 2.10, (32), the Euler transform is a function p + 2 Ψ q + 2 ( z k τ α k ) , because the indices are increased by one for the E-K operator; the Saigo operators will increase the indices by two; the M-S-M integral will be in terms of Wright function with indices increased by three, namely: ( p + 1 ) + 3 Ψ ( q + 1 ) + 3 , etc.
Special cases of (91) are the generalized K-series (65) and M-series (56), see Examples 4 and 7.
Example 21.
The generalized k-Wright function (multi-parametric k-M-L function) is introduced by Purohit and Badguzer [78] as a k-extension of the Wright g.h.f. (4):
p Ψ q k ( z ) = p Ψ q k ( a 1 , A 1 ) , , ( a p , A p ) ( b 1 , B 1 ) , , ( b q , B q ) z = n = 0 Γ k ( a 1 + n A 1 ) Γ k ( a p + n A p ) Γ k ( b 1 + n B 1 ) Γ k ( b q + n B q ) z n n ! , k > 0 .
However, from the representation (83), it is seen that this “new” function is again a Wright generalized hypergeometric function, namely:
c o n s t p + 1 Ψ q + 1 ( a i / k , A i / k ) i = 1 p ( b j / k , B j / k ) j = 1 q k ( A 1 + + A p B 1 B q ) / k · z .
Then the M-S-M operators evaluated for (92) by these authors can appear directly from our general results (say, Lemma 4) in terms of p + 4 Ψ q + 4 -functions.

10. Conclusions

10.1. The researchers on the topic can be advised to follow a procedure like this:
(1)
Check if the considered SF can be presented as a Wright g.h.f. p Ψ q or as simpler p F q -function; in more complicated cases, or if it is a Fox H-function or a Meijer G-function;
(2)
Check if the operator of FC to be evaluated is some special case of the GFC operators, that is, if it can be presented as a composition of classical R-L or E-K operators (also in the form (20) and (26));
(3)
Then, apply a general result like Theorem 3, Theorem 4 (or more generally, Theorem 2) and their special cases (Lemmas 1–4) and the examples, provided in this survey.
10.2. In Section 3 we first give the images of the SF (the generalized hypergeometric functions p Ψ q , p F q and their simplest cases) for the classical FC operators: E-K and R-L, and show that these are the same kinds of functions of which the indices p , q are increased by 1. Then, the images under the GFC operators are obtained by m-times application of these results, in Section 4. Our result states that the image of a p Ψ q -function (resp. p F q -function) under any (m-tuple) GFC operator can be predicted by Theorems 3 and 4 to be a p + m Ψ q + m -function (resp. p + m F q + m -function) with additional parameters depending on those of the FC operators. Using this general approach, one can avoid application of the standard term-by-term integration/differentiation of the power series for each particular special function the authors choose to treat.
10.3. For the proofs of Theorems 3 and 4 and their corollaries, see Kiryakova [2,3,4], and for their alternative interpretations—in other our works as ([6,22,24,50], Ch.4), [27]. The basic idea is that by means of a multiple (m-tuple) operator of GFC each p + m Ψ q + m -function (resp. a p + m F q + m -function) can be reduced to a p Ψ q -function (resp. a p F q -function), see comments and formula (50) before Theorem 4. Thus, by a suitable number of steps, from any p Ψ q -function (resp. p F q -function) we can reach to one of the three basic generalized hypergeometric functions, depending on either p < q , p = q or p = q + 1 : 0 Ψ q p , 1 Ψ 1 , 2 Ψ 1 (resp. 0 F q p , 1 F 1 , 2 F 1 ). Additionally, by an Erdélyi–Kober operator these are reducible to one of the three elementary functions: cos z , z α exp ( z ) or ( 1 z ) α z β . Details are given in Kiryakova [27], submitted to this Journal. As a conclusion, we have classified the g.h.f., that is the SF, in three basic classes: “g.h.f. of Bessel/cosine type”, “g.h.f. of confluent/exp type” and “g.h.f. of Gauss/beta-distribution type”, each of these classes with own specific properties. Thus, the title of Kiryakova [50] appeared as: “All the special functions are fractional differintegrals of elementary functions”.
10.4. In some papers, the authors evaluate an operator of FC of a particular special function in terms of another special function. Or even, the final result is written only as a series not recognized as some SF. However, for both theoretical reasons and possible applications, the results can be useful when a GFC operator transforms a special function from some class into a special function of the same class, although with changed (increased/decreased) indices and additional parameters. Such are our Theorems 3 and 4, showing that a p Ψ q -function (resp. a p F q -function) is transformed into a p + m Ψ q + m -function (resp. a p + m F q + m -function). We discussed similar results also for some particular cases of SF, such as for the (classical) Wright function ϕ ( α , β , z ) , the M-series p M α , β q ( z ) . Among the illustrative examples for SF having FC images of the same class, are our formulas: (32), (36), (37), (39), (40), (41), (46), (49), (51), and their corollaries like (55)–(58),(61), (62), (70), (77).
10.5. Next goal: the most useful results on the topic are when we can specify an operator of GFC corresponding to the considered special function, so that this function is to be its eigenfunction. In this survey we give such examples, say for the multi-index Mittag-Leffler functions E ( α i ) , ( β i ) ( z ) and the hyper-Bessel functions of Delerue 0 F m 1 ( 1 + γ i ) 1 m 1 ; z . These are the formulas (79), (81) and (82), etc. Another example, for the eigenfunction of the simplest fractional order differential equation, is the Rabotnov function z α 1 E α , α ( z α ) (called also fractional exponent), namely:
D α y α ( z ) = λ y α ( z ) ,   where   y α ( z ) = z α 1 E α , α ( λ z α ) ,   α > 0 ,   λ 0 .
10.6. Many authors are publishing results on the images of particular special functions under some integral transforms like the Laplace transform, Mellin transform, Euler (Beta) transform, Whittaker transform. Observe that the Euler transform (called so after the Euler integral formula for the Gauss function) is just a case of the Erdélyi–Kober fractional integral (12), as an extension of the Riemann Liuoville fractional integral (14). Therefore, there is no need to separately evaluate these two transforms (Euler transform and Riemann–Liouville operator), and what is more, to repeat such calculations for each particular special function. One can just apply the general result, as in Lemma 1. Note that the so-called pathway-transform is also closely related to the E-K integral. To evaluate a Laplace transform, say for any special function which is an H-function, one can use the general integral formula (44) and the representation of the kernel exponential function as a Wright g.h.f. (see (42), Example 1), then also as a H-function: exp ( z ) = H 0 , 1 1 , 0 z ( 0 , 1 ) . As already mentioned in Remark 1, a basic approach to evaluate integral transforms (also FC operators) of special functions relies on their images under the Mellin transform in terms of Gamma-functions, to which the fundamental book by Marichev [8] is devoted.

Author Contributions

The ideas and results in this paper survey and reflect the author’s (V.K.) sole contributions, resulting from more than 30 years of research on the topic. The author haves read and agreed to the published version of the manuscript.

Funding

This research received no financial funding.

Acknowledgments

This paper is done under the working programs on bilateral collaboration contracts of Bulgarian Academy of Sciences with Serbian and Macedonian Academies of Sciences and Arts, and under the COST program, COST Action CA15225 ‘Fractional’.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Machado, J.A.T.; Kiryakova, V. Recent history of the fractional calculus: Data and statistics. In Handbook of Fractional Calculus with Applications. Volume 1: Basic Theory; Kochubei, A., Luchko, Y., Eds.; De Gryuter: Berlin, Germany, 2019; pp. 1–12. [Google Scholar] [CrossRef]
  2. Kiryakova, V. Fractional calculus operators of special functions?—The result is well predictable! Chaos Solitons Fractals 2017, 102, 2–15. [Google Scholar] [CrossRef]
  3. Kiryakova, V. Use of fractional calculus to evaluate some improper integrals of special functions. AIP Conf. Proc. 2017, 1910, 050012. [Google Scholar] [CrossRef] [Green Version]
  4. Kiryakova, V. Fractional calculus of some “new” but not new special functions: k-, multi-index-, and S-analogues. AIP Conf. Proc. 2019, 2172, 0500088. [Google Scholar] [CrossRef]
  5. Kiryakova, V. Commentary: “A remark on the fractional integral operators and the image formulas of generalized Lommel-Wright function”. Front. Phys. 2019, 7, 145. [Google Scholar] [CrossRef]
  6. Kiryakova, V. Generalized Fractional Calculus and Applications; Longman—J. Wiley: Harlow, UK; New York, NY, USA, 1994. [Google Scholar]
  7. Prudnikov, A.P.; Brychkov, Y.; Marichev, O.I. Integrals and Series, Vol. 3: More Special Functions; Gordon and Breach Sci. Publ.: New York, NY, USA; London, UK; Paris, France; Tokyo, Japan, 1992. [Google Scholar]
  8. Marichev, O.I. Handbook of Integral Transforms of Higher Transcendental Functions, Theory and Algorithmic Tables; Ellis Horwood: Chichester, UK, 1983. Translated from Russian; Method of Evaluation of Integrals of Special Functions (In Russian); Nauka i Teknika, Minsk, Belarus, 1978 [Google Scholar]
  9. Srivastava, H.M.; Gupta, K.S.; Goyal, S.P. The H-Functions of One and Two Variables with Applications; South Asian Publs: New Delhi, India, 1982. [Google Scholar]
  10. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  11. Podlubny, I. Fractional Differential Equations; Acad. Press: Boston, MA, USA, 1999. [Google Scholar]
  12. Yakubovich, S.; Luchko, Y. The Hypergeometric Approach to Integral Transforms and Convolutions; Ser. Mathematics and Its Applications 287; Kluwer Acad. Publ.: Dordrecht, The Netherlands; Boston, MA, USA; London, UK, 1994. [Google Scholar]
  13. Mathai, A.M.; Haubold, H.J. Special Functions for Applied Scientists; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  14. Gorenflo, R.; Kilbas, A.; Mainardi, F.; Rogosin, S. Mittag-Leffler Functions, Related Topics and Applications, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2014. [Google Scholar] [CrossRef]
  15. Erdélyi, A. (Ed.) Higher Transcendental Functions; McGraw Hill: New York, NY, USA, 1953; Volume 1–3. [Google Scholar]
  16. Kilbas, A.A. Fractional calculus of the generalized Wright function. Fract. Calc. Appl. Anal. 2005, 8, 113–126. [Google Scholar]
  17. Haubold, H.J.; Mathai, A.M.; Saxena, R.K. Mittag-Leffler functions and their applications. J. Appl. Math. 2011, 51, 298628. [Google Scholar] [CrossRef] [Green Version]
  18. Rogosin, S. The role of the Mittag-Leffler function in fractional modeling. Mathematics 2015, 3, 368–381. [Google Scholar] [CrossRef]
  19. Paneva-Konovska, J. From Bessel to Multi-Index Mittag-Leffler Functions: Enumerable Families, Series in Them and Convergence; World Scientific Publishing: London, UK, 2016. [Google Scholar]
  20. Kiryakova, V. Multiple Dzrbashjan-Gelfond-Leontiev fractional differintegrals. Recent Adv. Appl. Math. 1996, 96, 281–294. [Google Scholar]
  21. Kiryakova, V. Multiple (multiindex) Mittag-Leffler functions and relations to generalized fractional calculus. J. Comput. Appl. Math. 2000, 118, 241–259. [Google Scholar] [CrossRef] [Green Version]
  22. Kiryakova, V. The multi-index Mittag-Leffler functions as important class of special functions of fractional calculus. Comput. Math. Appl. 2010, 59, 1885–1895. [Google Scholar] [CrossRef] [Green Version]
  23. Kilbas, A.A.; Koroleva, A.A.; Rogosin, S.V. Multi-parametric Mittag-Leffler functions and their extension. Fract. Calc. Appl. Anal. 2013, 16, 378–404. [Google Scholar] [CrossRef]
  24. Kiryakova, V. The special functions of fractional calculus as generalized fractional calculus operators of some basic functions. Comput. Math. Appl. 2010, 59, 1128–1141. [Google Scholar] [CrossRef] [Green Version]
  25. Dzrbashjan, M.M. On the integral transformations generated by the generalized Mittag-Leffler function. Izv. Arm. SSR 1960, 13, 21–63. (In Russian) [Google Scholar]
  26. Paneva-Konovska, J. Multi-index (3m-parametric) Mittag-Leffler functions and fractional calculus. Compt. Rend. Acad. Bulg. Sci. 2011, 64, 1089–1098. [Google Scholar]
  27. Kiryakova, V. A guide to special functions in fractional calculus. Math. Spec. Issue Spec. Funct. Math. Phys. Part II 2020. Submitted. [Google Scholar]
  28. Luchko, Y.; Kiryakova, V. The Mellin integral transform in fractional calculus. Fract. Calc. Appl. Anal. 2013, 16, 405–430. [Google Scholar] [CrossRef]
  29. Agarwal, P.; Rogosin, S.V.; Trujillo, J.J. Certain fractional integral operators and the generalized multi-index Mittag-Leffler functions. Proc. Indian Acad. Sci. (Math. Sci.) 2015, 125, 291–306. [Google Scholar] [CrossRef]
  30. Paneva-Konovska, J.; Kiryakova, V. On the multi-index Mittag-Leffler functions and their Mellin transforms. Int. J. Appl. Math. 2020, 33, 549–571. [Google Scholar] [CrossRef]
  31. Samko, S.; Kilbas, A.; Marichev, O. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach: Yverdon, Switzerland, 1993. [Google Scholar]
  32. Sandev, T.; Tomovski, Ž. Fractional Equations and Models (Theory and Applications); Springer: Berlin/Heidelberg, Germany, 2019. [Google Scholar] [CrossRef]
  33. Sneddon, I.N. The use in mathematical analysis of Erdélyi-Kober operators and some of their applications. In Fractional Calculus and Its Applications (Proc. Internat. Conf. Held in New Haven); Ross, B., Ed.; Lecture Notes in Math. 457; Springer: New York, NY, USA, 1975; pp. 37–79. [Google Scholar]
  34. Kiryakova, V. Generalized fractional calculus operators with special functions. In Handbook of Fractional Calculus with Applications. Volume 1: Basic Theory; Kochubei, A., Luchko, Y., Eds.; De Gryuter: Berlin, Germany, 2019; pp. 87–110. [Google Scholar] [CrossRef]
  35. Kiryakova, V.; Luchko, Y. Riemann-Liouville and Caputo type multiple Erdélyi-Kober operators. Cent. Eur. J. Phys. 2013, 11, 1314–1336. [Google Scholar] [CrossRef]
  36. Luchko, Y.; Trujillo, J.J. Caputo type modification of the Erdélyi-Kober fractional derivative. Fract. Calc. Appl. Anal. 2007, 10, 249–267. [Google Scholar]
  37. Kiryakova, V. A brief story about the operators of the generalized fractional calculus. Fract. Calc. Appl. Anal. 2008, 11, 203–220. [Google Scholar]
  38. Kiryakova, V. Gel’fond-Leont’ev integration operators of fractional (multi-)order generated by some special functions. AIP Conf. Proc. 2018, 2048, 050016. [Google Scholar] [CrossRef]
  39. Marichev, O.I. Volterra equation of Mellin convolutional type with a Horn function in the kernel. Izv. AN BSSR, Ser. Fiz.-Mat. Nauk 1974, 1, 128–129. (In Russian) [Google Scholar]
  40. Saigo, M.; Maeda, N. More generalization of fractional calculus. In Transform Methods & Special Functions, Varna’96 (Proc. Second In- ternat. Workshop); Rusev, P., Dimovski, I., Kiryakova, V., Eds.; Science Culture Technology Publishing: Singapore, 1998; pp. 386–400. [Google Scholar]
  41. Dimovski, I. Operational calculus for a class of differental operators. C. R. Acad. Bulg. Sci. 1966, 19, 1111–1114. [Google Scholar]
  42. Dimovski, I.; Kiryakova, V. Generalized Poisson transmutations and corresponding representations of hyper-Bessel functions. C. R. Acad. Bulg. Sci. 1986, 39, 29–32. [Google Scholar]
  43. Delerue, P. Sur le calcul symboloque à n variables et fonctions hyperbesseliennes (II). Annales Soc. Sci. Bruxelles, Ser. 1 1953, 3, 229–274. [Google Scholar]
  44. Erdélyi, A. (Ed.) Tables of Integral Transforms; McGraw Hill: New York, NY, USA, 1954; Volume 1–2. [Google Scholar]
  45. Askey, R. Orthogonal Polynomials and Special Functions; SIAM: Philadelphia, PA, USA, 1975. [Google Scholar]
  46. Lavoie, J.L.; Osler, T.J.; Tremblay, R. Fractional derivatives and special functions. SIAM Rev. 1976, 18, 240–268. [Google Scholar] [CrossRef]
  47. Ali, R.S.; Mubeen, S.; Ahmad, M.M. A class of fractional integral operators with multi-index Mittag-Leffler k-function and Bessel k-function of first kind. J. Math. Comput. Sci. 2020, 22, 266–281. [Google Scholar] [CrossRef]
  48. Mubeen, S.; Ali, R.S.; Nayab, I.; Rahman, G.; Abdeljavad, T.; Nisar, K.S. Integral transforms of an extended generalized multi-index Bessel function. AIMS Math. 2020, 5, 7531–7546. [Google Scholar] [CrossRef]
  49. Abdalla, M. Fractional operators for the Wright hypergeometric matrix functions. Adv. Differ. Equ. 2020, 2020, 246. [Google Scholar] [CrossRef]
  50. Kiryakova, V. All the special functions are fractional differintegrals of elementary functions. J. Phys. A Math. Gen. 1997, 30, 5085–5103. [Google Scholar] [CrossRef] [Green Version]
  51. Mathai, A.M.; Saxena, R.K.; Haubold, H.J. The H-Function; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  52. Garrappa, R.; Kaslik, E.; Popolizio, M. Evaluation of fractional integrals and derivatives of elementary functions: Overview and tutorial. Mathematics 2019, 7, 407. [Google Scholar] [CrossRef] [Green Version]
  53. Kilbas, A.A.; Sebastian, N. Generalized fractional integration of Bessel function of first kind. Integr. Transf. Spec. Funct. 2008, 19, 869–883. [Google Scholar] [CrossRef]
  54. Sharma, M.; Jain, R. A note on a generalized M-series as a special function of fractional calculus. Fract. Calc. Appl. Anal. 2009, 12, 449–452. [Google Scholar]
  55. Lavault, C. Fractional calculus and generalized Mittag-Leffler type functions. arXiv 2017, arXiv:1703.01912. [Google Scholar]
  56. Kumar, D.; Saxena, R.K. Generalized fractional calculus of the M-Series involving F3 hypergeometric function. Sohag. J. Math. 2015, 2, 17–22. [Google Scholar] [CrossRef]
  57. Saigo, M. A remark on integral operators involving the Gauss hypergeometric functions. Math. Rep. Coll. Gen. Ed. Kyushu Univ. 1978, 11, 135–143. [Google Scholar]
  58. 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]
  59. Sharma, K. An introduction to the generalized fractional integration. Bol. Soc. Paran. Math. 2012, 30, 85–90. [Google Scholar] [CrossRef]
  60. Purohit, S.D.; Suthar, D.L.; Kalla, S.L. Marichev-Saigo-Maeda fractional integration operators of the Bessel functions. Le Mat. 2012, LXVII, 21–32. [Google Scholar] [CrossRef]
  61. Agarwal, R.; Jain, S.; Agarwal, R.P.; Baleanu, D. Response: Commentary: A remark on the fractional integral operators and the image formulas of generalized Lommel-Wright function. Front. Phys. 2020, 8, 72. [Google Scholar] [CrossRef]
  62. Agarwal, R.; Jain, S.; Agarwal, R.P.; Baleanu, D. A remark on the fractional integral operators and the image formulas of generalized Lommel-Wright function. Front. Phys. 2018, 6, 79. [Google Scholar] [CrossRef]
  63. Mondal, S.R.; Nisar, K.S. Marichev-Saigo-Maeda fractional integration operators involving generalized Bessel functions. Math. Probl. Eng. 2014, 11, 274093. [Google Scholar] [CrossRef]
  64. Nisar, K.S.; Mondal, S.R.; Agarwal, P. Composition formulas of Bessel-Struve kernel function. arXiv 2016, arXiv:1602.00279v1. [Google Scholar] [CrossRef] [Green Version]
  65. Srivastava, H.M.; Tomovski, Ž. Fractional calculus with an integral operator containing generalized Mittag-Leffler function in the kernel. Appl. Math. Comput. 2009, 211, 198–210. [Google Scholar] [CrossRef]
  66. Saxena, R.K.; Nishimoto, K. N-fractional calculus of generalized Mittag-Leffler functions. J. Fract. Calc. 2010, 37, 43–52. [Google Scholar]
  67. Saxena, R.K.; Pogany, T.K.; Ram, J.; Daiya, J. Dirichlet averages of generalized multi-index Mittag-Leffler functions. Armen. J. Math. 2010, 3, 174–187. [Google Scholar]
  68. Kumar, D.; Kumar Gupta, R.; Singh Rawat, D. Marichev-Saigo-Maeda fractional differential operator involving mittag-Leffler type function with four parameters. J. Chem. Biol. Phys. Sci. Sect. C: Phys. Sci. 2017, 7, 201–210. [Google Scholar]
  69. Kiryakova, V. From the hyper-Bessel operators of Dimovski to the generalized fractional calculus. Fract. Calc. Appl. Anal. 2014, 17, 977–1000. [Google Scholar] [CrossRef]
  70. Gehlot, K.S. Differential equation of k-Bessel’s function and its properties. Nonl. Anal. Differ. Equ. 2014, 2, 61–67. [Google Scholar] [CrossRef] [Green Version]
  71. Mondal, S.R. Representation formulae and monotonicity of the generalized k-Bessel functions. arXiv 2016, arXiv:1611.07499. [Google Scholar]
  72. Shaktawat, B.S.; Rawat, D.S.; Gupta, R.K. On generalized fractional calculus of the generalized k-Bessel function. J. Rajasthan Acad. Phys. Sci. 2017, 16, 9–19. [Google Scholar]
  73. Dorrego, G.A.; Cerruti, R.A. The k-Mittag-Leffler function. Int. J. Contemp. Math. Sci. 2012, 7, 705–716. [Google Scholar]
  74. Gupta, A.; Parihar, C.L. k-New generalized Mittag-Leffler function. J. Fract. Calc. Appl. 2014, 5, 165–176. [Google Scholar]
  75. Nisar, K.S.; Eata, A.F.; Al-Dhaifallah, M.; Choi, J. Fractional calculus of generalized k-Mittag-Leffler function and its applications to statistical distribution. Adv. Differ. Equ. 2016, 304. [Google Scholar] [CrossRef] [Green Version]
  76. Nisar, K.S.; Purohit, S.D.; Parmar, R.K. Fractional calculus and certain integrals of generalized multiindex Bessel function. arXiv 2017, arXiv:1706.08039. [Google Scholar]
  77. Saxena, R.K.; Daiya, J. Integral transforms of S-functions. Le Mat. 2015, LXX, 147–159. [Google Scholar]
  78. Purohit, M.; Badguzer, A. MSM fractional integration and differentiation operators of multi-parametric K-Mittag Leffler function and generalized multi-index Bessel function. Intern. J. Stat. Appl. Math. 2018, 3, 156–161. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kiryakova, V. Unified Approach to Fractional Calculus Images of Special Functions—A Survey. Mathematics 2020, 8, 2260. https://0-doi-org.brum.beds.ac.uk/10.3390/math8122260

AMA Style

Kiryakova V. Unified Approach to Fractional Calculus Images of Special Functions—A Survey. Mathematics. 2020; 8(12):2260. https://0-doi-org.brum.beds.ac.uk/10.3390/math8122260

Chicago/Turabian Style

Kiryakova, Virginia. 2020. "Unified Approach to Fractional Calculus Images of Special Functions—A Survey" Mathematics 8, no. 12: 2260. https://0-doi-org.brum.beds.ac.uk/10.3390/math8122260

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