Next Article in Journal
Compact Resolutions and Analyticity
Next Article in Special Issue
Gauss’ Second Theorem for F12(1/2)-Series and Novel Harmonic Series Identities
Previous Article in Journal
A Semantic Enhancement Framework for Multimodal Sarcasm Detection
Previous Article in Special Issue
Novel Formulas for B-Splines, Bernstein Basis Functions, and Special Numbers: Approach to Derivative and Functional Equations of Generating Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Going Next after “A Guide to Special Functions in Fractional Calculus”: A Discussion Survey

by
Virginia Kiryakova
*,† and
Jordanka Paneva-Konovska
Institute of Mathematics and Informatics, Bulgarian Academy of Sciences, 1113 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2024, 12(2), 319; https://doi.org/10.3390/math12020319
Submission received: 17 December 2023 / Revised: 8 January 2024 / Accepted: 16 January 2024 / Published: 18 January 2024
(This article belongs to the Special Issue Integral Transforms and Special Functions in Applied Mathematics)

Abstract

:
In the survey Kiryakova: “A Guide to Special Functions in Fractional Calculus” (published in this same journal in 2021) we proposed an overview of this huge class of special functions, including the Fox H-functions, the Fox–Wright generalized hypergeometric functions p Ψ q and a large number of their representatives. Among these, the Mittag-Leffler-type functions are the most popular and frequently used in fractional calculus. Naturally, these also include all “Classical Special Functions” of the class of the Meijer’s G- and p F q -functions, orthogonal polynomials and many elementary functions. However, it so happened that almost simultaneously with the appearance of the Mittag-Leffler function, another “fractionalized” variant of the exponential function was introduced by Le Roy, and in recent years, several authors have extended this special function and mentioned its applications. Then, we introduced a general class of so-called (multi-index) Le Roy-type functions, and observed that they fall in an “Extended Class of SF of FC”. This includes the I-functions of Rathie and, in particular, the H ¯ -functions of Inayat-Hussain, studied also by Buschman and Srivastava and by other authors. These functions initially arose in the theory of the Feynman integrals in statistical physics, but also include some important special functions that are well known in math, like the polylogarithms, Riemann Zeta functions, some famous polynomials and number sequences, etc. The I- and H ¯ -functions are introduced by Mellin–Barnes-type integral representations involving multi-valued fractional order powers of Γ -functions with a lot of singularities that are branch points. Here, we present briefly some preliminaries on the theory of these functions, and then our ideas and results as to how the considered Le Roy-type functions can be presented in their terms. Next, we also introduce Gelfond–Leontiev generalized operators of differentiation and integration for which the Le Roy-type functions are eigenfunctions. As shown, these “generalized integrations” can be extended as kinds of generalized operators of fractional integration, and are also compositions of “Le Roy type” Erdélyi–Kober integrals. A close analogy appears with the Generalized Fractional Calculus with H- and G-kernel functions, thus leading the way to its further development. Since the theory of the I- and H ¯ -functions still needs clarification of some details, we consider this work as a “Discussion Survey” and also provide a list of open problems.

1. Introduction

“Special functions are particular mathematical functions that have more or less established names and notations due to their importance in mathematical analysis, functional analysis, geometry, physics, or other applications. The term is defined by consensus, and thus lacks a general formal definition, but the list of mathematical functions contains functions that are commonly accepted as special…” (see more, in https://en.wikipedia.org/wiki/Special_functions, accessed on 16 December 2023). Almost by the end of the 20th century, under terms such as “Higher Transcendental Functions”, “Special Functions of Mathematical Physics” or “Named Special functions”, scientists had in mind these appearing as solutions of differential equations of integer order (usually 2nd order) or as integrals of elementary functions. One can briefly refer to the variety of all these, like the Bessel and cylindrical functions; the Gauss, Kummer, Tricomi, confluent and generalized hypergeometric functions; the classical orthogonal polynomials (such as Laguerre, Jacobi, Gegenbauer, Legendre, Tchebisheff, Hermite, etc.); the incomplete Gamma and Beta functions; the Error functions; and the Airy, Whittaker, etc., functions. Their definitions, properties and tables have appeared in many basic handbooks on the topic of that era, but we just give a quick mention here to the three-volume Bateman project [1] (planned by Bateman as “A Guide to the Mathematical Functions” but edited posthumously by his collaborators).
In the 1970–1980s, fractional calculus (FC: calculus where the operators of differentiation and integration can be of arbitrary, i.e., “fractional”, order) enjoyed a revival and fast development as not only an exotic mathematical theory but also due to its acknowledged useful applications in mathematical models of various processes and phenomena of the real physical and social world. Thus, the solutions of differential and integral equations of (arbitrary) fractional order became unavoidable tools. Incorporating the “classical” special functions but also much more general ones, the class of so-called “Special Functions of Fractional Calculus” (SF of FC) has been extensively investigated, and is a topic of many modern books (such as [2,3,4,5,6,7,8], and so on) and a huge number of surveys and articles. One of these works is our survey, Kiryakova: “A Guide to Special Functions in Fractional Calculus” (published in this same journal in 2021) whose title is borrowed by Bateman’s plans for the project, appearing as [1]. As the main representatives of such functions, we consider the Fox H-functions and the Fox–Wright generalized hypergeometric functions p Ψ q , which are extensions of Meijer’s G- and p F q -functions and of other “classical” ones. There is a very long list of such functions, among which those named after Mittag-Leffler are the most popular and often used in FC. In the mentioned survey [9], we present an attempt to overview the definitions, some properties and applications, and the works dedicated to these SF of FC.
However, many other special functions that have not up to now been mentioned under the scheme of “SF of FC” also appear to be important in the analysis, number theory, statistics and various branches of applied mathematics and physics. Here, we pay attention to the so-called I-functions of Rathie [10], and to their specifications, the H ¯ -functions of Inayat-Hussain [11], see also Buschman and Srivastava [12]. These functions were introduced long ago (at the end of the 20th century) but were somehow neglected since they looked too complicated or artificial. However, they have appeared in the theory of the Feynman integrals in statistical physics, the partition of the Gaussian model, etc., and it so happens that some well-known SF are their particular cases, like the polylogarithms, Riemann Zeta functions and others.
In recent years, a new direction in the theory of SF has appeared to extend the Le Roy functions. While studying a general class of the so-called Le Roy-type functions, the authors of this paper observed that these functions can fall in an “Extended Class of SF of FC”, and after our publications on the analytical properties of the Le Roy-type functions as important entire functions ([13,14,15,16]), here, we present some ideas about their relations to the I- and H ¯ -functions and their role for the further development of Generalized Fractional Calculus.
As the theory of the I- and H ¯ -functions still needs some clarifications and further development, and in view of the appearing multi-valued fractional order powers of Γ -functions with a lot of singularities that are branch points, we consider this work to be a “Discussion Survey” and also discuss a list of Open Problems.

2. Preliminary Definitions

From the Bateman Project [1] (Vol.1, p.49): “…Of all integrals which contain Gamma functions in their integrands the most important ones are the so-called Mellin-Barnes integrals. Such integrals were first introduced by S. Pincherle ([17], 1888); their theory has been developed in 1910 by H. Mellin …and they were used for a complete integration of the hypergeometric differential equation by E.W. Barnes, 1908.” See Mainardi–Pagnini [18].
Definition 1
(Ch. Fox [19], 1961; see books such as [3,4,6,7,8,20], etc.). The generalized hypergeometric function, defined by means of the Mellin–Barnes-type contour integral in the complex plane:
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 , z 0 ,
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 ) ,
is known as the Fox H-function. 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 are arbitrary complexes such that A i ( b j + l ) B j ( a i l 1 ) , l , l = 0 , 1 , 2 , ; i = 1 , , n ; j = 1 , , m , and so, the poles of the Gamma functions in the numerator do not coincide. The single-valued branch of z s is chosen as z s = exp [ s { l o g | z | + i a r g z } ] , | a r g z | < π . Note that L is an infinite contour that separates all the poles of Γ ( b j + B j s ) to the left and all the poles of Γ ( 1 a i A i s ) to the right of L , and can be one of the three following variants, as described, e.g., in the book by Kilbas–Saigo [3]:
( i ) L = L is a left loop placed in the horizontal strip starting at the point + i φ 1 and terminating at the point + i φ 2 with < φ 1 < φ 2 < + .
( i i ) L = L + is a right loop placed in the horizontal strip starting at the point + + i φ 1 and terminating at the point + + i φ 2 with < φ 1 < φ 2 < + .
( i i i ) L = L i c is a contour that starts at the point c i and terminates at the point c + i , with a suitable abscissa c ( , + ) . Often, it is also denoted as L = ( c i , c + i ) .
Note that when the Mellin transform of the H-function exists (see Th.2.2 in [3]), it is equal to the integrand H p , q m , n ( s ) . Moreover, in some books and other works, the term z s in (1) is taken to be z s (then the orientation of the contours is changed).
The following parameters are used to characterize the behavior and the properties of the H-function:
R = i = 1 p A i A i j = 1 q B j B j ; Δ = j = 1 q B j i = 1 p A i ; δ = m + n p + q 2 , μ = 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 .
Depending on these, the H-function can be a function analytic in the disks | z | < R or | z | > R , in some angular sectors or in the whole complex plane. For more details on the properties of the Fox H-functions, the reader can consult the above-mentioned and other contemporary handbooks on special functions and their applications.
It is important to know about the behavior of the H-function around the singular points. This matter is rather complicated and, especially when there is a third singular point on the circle of convergence ( | z | = R ), there are known results only in some particular cases. For more information, the reader can see the work of Karp [21], which comments on and revisits the old, basic results of Braaksma [22].
The so-called Meijer G-function (C.S. Meijer [23], 1936–1941, etc.) is a simpler and more popular case of the H-function (1) when A i = B j = 1 , i = 1 , , p ; j = 1 , , q , and details on its definition, properties and examples appeared in [1] (Vol.1, Ch.5), as well as in the 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 .
(Hereafter, we prefer denotations by z s instead of z s , as this is more conventional, which only yields a change in the orientation of contours.) The values of the characteristic parameters in (2) for the G-function reduce as follows: R = 1 , Δ = q p , δ = m + n p + q 2 , and its behavior depends on them. The G-function is simpler than the H-function but yet is general enough and incorporates most of the Classical Special Functions (known also as Named SF, SF of Mathematical Physics), many orthogonal polynomials and elementary functions. Long lists of examples can be found, e.g., in [1] (Vol.1), [20] (Appendix C) and other recent books.
The following generalized hypergeometric function was introduced and studied by Sir Edward Maitland (E.-M.) Wright in series of his works such as [24,25], etc., 1933–1940, see also Fox [26]. It is an example of a H-function that in general (unless A j , B k are integer or rational) does not reduce to a G-function.
Definition 2
(see, e.g., [2,7,20] (Appendix E), etc.). The Wright generalized hypergeometric function p Ψ q ( z ) , also called the Fox–Wright function, is defined by the power series:
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 ! ,
but is representable also as a H-function:
p Ψ q ( a i , A i ) 1 p ( b j , B j ) 1 q z = 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 ) .
According to the denotations for the corresponding parameters (2), the power series (4) defines an entire function of z if Δ > 1 ; it is absolutely convergent in the disk { | z | < R } for Δ = 1 , and the same holds for | z | = R if R e ( μ ) > 1 / 2 . Aside from the numerous handbooks on SF, the details are summarized in Gorenflo–Luchko–Mainardi [27].
The Wright g.h.f. reduces to the more popular generalized hypergeometric  p F q -function, and thus to a Meijer G-function (3), when A 1 = = A p = 1 , B 1 = = B q = 1 :
p Ψ q ( a 1 , 1 ) , , ( a p , 1 ) ( b 1 , 1 ) , , ( b q , 1 ) z = c k = 0 ( a 1 ) k ( a p ) k ( b 1 ) k ( b q ) k z k k !
= c · p F q ( a 1 , , a p ; b 1 , , b q ; z ) = 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 ) 0 = 1 , ( a ) k = a ( a + 1 ) ( a + k 1 ) = Γ ( a + k ) / Γ ( a ) . Here, it is pertinent to mention that the identity ( a ) k = Γ ( a + k ) / Γ ( a ) holds true only if ( a + k ) 1 , 2 , .
Details on the theory and examples of the p F q -functions can be found in the old sources such as [1] (Vol.1, Ch.4). Series (6) converges for all finite z if p q ; for | z | < 1 if p = q + 1 ; and diverges for z 0 if p > q + 1 . The simplest (lowest indices) particular cases are the Gauss hypergeometric function 2 F 1 , the Kummer (confluent hypergeometric) function 1 F 1 and the Bessel type functions 0 F 1 . Corresponding to these three possible cases for the orders p and q, in Kiryakova [20] (Ch.4), see also [9] (Sect.8, Th.7–Th8), we suggested a classification of the p F q and p Ψ q functions in three basic groups of SF.
Since in the next sections we often explore the Mittag-Leffler function E α (introduced in [28]) and its extensions, let us briefly give a reminder of some definitions.
The Mittag-Leffler (M-L) functions  E α , β (with two indices) are defined by the power series (details in [2,29,30,31], etc. )
E α , β ( z ) = k = 0 z k Γ ( α k + β ) , α > 0 , β > 0 ; f o r β = 1 : E α ( z ) = k = 0 z k Γ ( α k + 1 ) .
These are entire functions of order ρ = 1 / α and type 1, and are also representable as
E α , β ( z ) = 1 Ψ 1 ( 1 , 1 ) ( β , α ) z = H 1 , 2 1 , 1 z ( 0 , 1 ) ( 1 β , α ) .
A three-parameter variant, usually called a Prabhakar function ([32]; see detailed studies, e.g., in Garra–Garrappa [33], Giusti et al. [34]), has an additional parameter τ > 0 , namely:
E α , β τ ( z ) = k = 0 ( τ ) k Γ ( α k + β ) z k k ! = 1 Γ ( τ ) · 1 Ψ 1 ( τ , 1 ) ( β , α ) z
= 1 Γ ( τ ) H 1 , 2 1 , 1 z ( 1 τ , 1 ) ( 0 , 1 ) , ( 1 β , α ) ,
where ( τ ) 0 = 1 , ( τ ) k = Γ ( τ + k ) / Γ ( τ ) , k = 1 , 2 , 3 , . denotes the Pochhammer symbol. For τ = 1 , we have the M-L function E α , β , and if, additionally, β = 1 , it is E α .
The following multi-index extensions of the M-L function (7) (also called the vector index extension) were introduced by Luchko et al. [35,36] and Kiryakova [37], and studied extensively in our works such as [38,39,40], etc.
Let m 1 be an integer, then by means of the sets of the real “multi-indices”: ( α 1 > 0 , α 2 > 0 , , α m > 0 ) and ( β 1 , β 2 , , β m ) , we define the multi-index Mittag-Leffler function (multi-M-L f.) as follows:
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 .
Later, Kilbas—Koroleva–Rogosin [41] studied these functions without the restrictions for all of the α i (or their real parts) to be obligatorily non-negative.
In [42,43], etc., Paneva-Konovska introduced and studied the ( 3 m ) -parametric (multi-index) Mittag-Leffler–Prabhakar functions, similar to (9) but with an additional set of Prabhakar parameters ( τ 1 > 0 , , τ m > 0 ) :
E ( α i ) , ( β i ) ( τ i ) , m ( z ) = k = 0 ( τ 1 ) k ( τ m ) k Γ ( α 1 k + β 1 ) Γ ( α m k + β m ) z k ( k ! ) m
= T · m Ψ 2 m 1 ( τ 1 , 1 ) , , ( τ m , 1 ) ( β 1 , α 1 ) , , ( β m , α m ) ; ( 1 , 1 ) , , ( 1 , 1 ) ( m 1 ) t i m e s z
= T · H m , 2 m 1 , m z ( 1 τ i , 1 ) 1 m ( 0 , 1 ) 1 m ; ( 1 β i , α i ) 1 m , where T = i = 1 m Γ ( τ i ) 1 .
These are extensions of both the Prabhakar function (8) and the ( 2 m )-multi-index M-L functions (9). The Prabhakar function appears for m = 1 , while for τ 1 = = τ m = 1 , these are shown in (9). It is proved that the multi-index M-L-type functions (9) and (10) are entire functions, and their Mellin transforms are proposed in our joint paper [44].
More details and numerous particular cases of the Mittag-Leffler-type functions (7)–(10), appearing as very important tools in FC, are provided in our recent survey “A Guide to Special Functions of Fractional Calculus” [9]. A reminder of some of these special functions are given again in brief in Section 6.
However, it seems that many other SF, that have not been mentioned as of yet under the scheme of “SF of FC” (in the sense of Kiryakova [9]), also appear to be important in the analysis, number theory, statistics and various branches of applied mathematics. Such examples and other extensions of the Fox H-functions, Meijer G-functions, Fox–Wright p Ψ q - and p F q -functions were introduced at the end of the 20th century but were somehow neglected as they looked too complicated or artificial. Next, we pay attention to the so-called I-functions of Rathie [10], and to their simpler variants, the H ¯ -functions of Inayat-Hussain [11], see also Buschman and Srivastava [12]. These are further extensions of the classes of the special functions beyond the cases of the Wright generalized hypergeometric functions (4) and Fox H-functions (1).
Definition 3.
The I-function is defined by a kind of Mellin–Barnes-type contour integral (Rathie [10]):
I p , q m , n z ( a j , A j , α j ) 1 p ( b k , B k , β k ) 1 q = 1 2 π i L I p , q m , n ( s ) z s d s , z 0 ,
where I ( s ) stands for
I p , q m , n ( s ) = k = 1 m Γ β k ( b k B k s ) j = 1 n Γ α j ( 1 a j + A j s ) k = m + 1 q Γ β k ( 1 b k + B k s ) j = n + 1 p Γ α j ( a j A j s )
with the power exponents for the Gamma functions: α j , j = 1 , , p and β k , k = 1 , , q , that in general, are not positive integers, and the contours L are discussed below.
Clearly, for non-integer values of these “powers” parameters, the I-function (11) is not expressible as a H-function, nor as other “classical” special functions.
Most of the denotations and some conditions on the orders and parameters in (11) and (12) are similar to those used for the H-function. Here, a j , j = 1 , , p and b k , k = 1 , , q can also be complex numbers such that no singularity of Γ β k ( b k B k s ) , k = 1 , , m coincides with any singularity of Γ α j ( 1 a j + A j s ) , j = 1 , , n . Again, we prefer the definition with term z s so as to be in agreement with the denotations in the initial works on the topic, and suppose that z s = exp [ s { log | z | + i arg z } ] , 0 arg z < 2 π .
In general, for non-integer α j and β k , the mentioned singularities are no longer poles but are converted to branching points of the multi-valued fractional powers of the Gamma functions. According to Rathie [10], the corresponding branch cuts should be chosen in a way so that the path of integration can be one of the three type of contours L :
(a) L starts from c i and goes to c + i , where the real abscissa c is chosen so that the singularities of Γ β k ( b k B k s ) , k = 1 , , m , lie to the right of L , and all singularities of Γ α j ( 1 a j + A j s ) , j = 1 , , n , lie to the left of it;
(b) L is a loop that begins and ends at + and encircles all the singularities of Γ β k ( b k B k s ) , k = 1 , , m , once in the clockwise direction, but none of the singularities of Γ α j ( 1 a j + A j s ) , j = 1 , , n ;
(c) L is a loop beginning and ending at and encircling all the singularities of Γ α j ( 1 a j + A j s ) , j = 1 , , n , once in the anti-clockwise direction, but none of the singularities of Γ β k ( b k B k s ) , k = 1 , , m .
Rathie [10] noted that contour (a) can be considered as a particular case of contour (c). In the case of contour (b), the sign of the function is changed with respect to case (c). When more than one of the above definitions of L make sense, they lead to the same result.
As in the theory of the Fox H-function, the following values related to the parameters (analogues of these in (2)) are important to characterize the behavior of the I-function:
Δ = k = 1 m β k B k k = m + 1 q β k B k + j = 1 n α j A j j = n + 1 p α j A j ; μ = k = 1 q β k B k j = 1 p α j A j ; = j = 1 p α j [ R e ( a j ) 1 / 2 ] k = 1 q β k [ R e ( b k ) 1 / 2 ] ; R = k = 1 q B k β k B k / j = 1 p A j α j A j .
In terms of (13), it is proved (Rathie [10]) that: the I-function (11) for contour L defined by (a), converges when | arg z | < Δ π / 2 , if Δ > 0 ; if | arg z | = Δ π / 2 , Δ 0 , the integral (11) converges absolutely when: (i) μ = 0 if > 1 ; (ii) μ 0 , if with s = σ + i t , σ and t real, σ chosen so that for | t | , we have + σ μ > 1 . Otherwise, the integral (11) for contour L defined by (b) converges for q 1 and either μ > 0 , or μ = 0 inside the disk | z | < R . Alternatively, the integral (11) for contour L defined by (c) converges when p 1 and either μ < 0 , or μ = 0 if | z | > R .
To determine single-valued branches of the multi-valued Γ -functions to the fractional powers in (11), one needs to draw suitable branch cuts. The details of such a procedure were proposed in the recent paper by Rogosin and Dubatovskaya [45] (Sect.3), for a similar function denoted there as (3.1) by H G , with L taken as a contour of the so-called “Slater” type ([46], see Marichev [5] (Ch.4)): L + starting at + + i φ 2 and finishing at + + i φ 1 , with suitable φ 1 < φ 2 . However, in the case of the I-functions, the ideas like in the Slater theorem for the G-functions (see same place) would need to be expanded when we have ratios of Γ -functions to arbitrary (fractional) powers.
Definition 4.
Another generalization of the H-function, denoted by the symbol H ¯ and called the Inayat-Hussain function ([11]) or bar-H-function, appears as a slightly simpler case of the I-function:
H ¯ p , q m , n z ( a j , A j , α j ) 1 n , ( a j , A j , 1 ) n + 1 p ( b k , B k , 1 ) 1 m , ( b k , B k , β k ) m + 1 q
= 1 2 π i i i k = 1 m Γ ( b k B k s ) j = 1 n Γ α j ( 1 a j + A j s ) k = m + 1 q Γ β k ( 1 b k + B k s ) j = n + 1 p Γ ( a j A j s ) z s d s .
Compared to the definition of the I-function, part of the powers of the Γ -functions are taken to be equal to 1, namely: β k = 1 , k = 1 , , m and α j = 1 , j = n + 1 , , p , and the contour L is chosen to be specifically as of type (a) (the imaginary axis, or if necessary shifted suitably as ( c i , + i )). In general (again like the I-function), the H ¯ -function does not reduce to a H-function, unless the additional power parameters in the Gamma functions are positive integers.
For the bar-H-function (14), the values of the parameters in (13) have corresponding but similar forms.
Definition 5.
A particularly important case with m = 1 , n = p (see Inayat-Hussain [11] (27)) corresponds to a kind of generalized hypergeometric function, similar to the Fox–Wright function (4). This is called the extended Fox–Wright function:
H ¯ p , q + 1 1 , p z ( 1 a j , 1 , α j ) 1 p ( 0 , 1 ) , ( 1 b k , 1 , β k ) 1 q : = p Ψ ¯ q ( a j , 1 ; α j ) j = 1 p ( b k , 1 ; β k ) k = 1 q z = j = 1 p Γ α j ( a j ) k = 1 q Γ β k ( b k ) n = 0 j = 1 p ( a j ) n α j k = 1 q ( b k ) n β k z n n ! , w i t h ( a ) 0 = 1 , ( a ) n : = a ( a + 1 ) ( a + n 1 ) .
In this case, there are substantial simplifications for the parameters (13) characterizing the behavior of the function (15):
Δ = 1 k = 1 q β k + j = 1 p α j ; μ = k = 1 q β k j = 1 p α j ; R = 1 ,
that is, under some conditions (if μ = 0 ), the function (15) is analytic inside | z | < 1 , or outside it; compare with the discussions in Inayat-Hussain [11], Rathie [10] and Saxena [47].
However, we go further, and introduce an even more general extension of the Fox–Wright function that we do need and use in our studies.
Definition 6.
We define the Generalized Fox–Wright function:
p Ψ ˜ q ( a j , A j , α j ) j = 1 p ( b k , B k , β k ) k = 1 q z : = H ¯ p , q + 1 1 , p z ( 1 a j , A j , α j ) 1 p ( 0 , 1 ) , ( 1 b k , B k , β k ) 1 q ,
with arbitrary positive parameters A j , B k as in the “classical” case (4) (instead of these taken to be 1 in (15)) but again with “fractional” power parameters α j > 0 , β k > 0 .
Next, for this p Ψ ˜ q -function, the characteristic parameters (13) take the form:
μ = 1 + k = 1 q β k B k j = 1 p α j A j ; R = k = 1 q B k β k B k / j = 1 p A j α j A j ( the   radius of | analycity ,   if   μ = 0 ) .
We observe (see next Section 4) that these generalized Fox–Wright functions, under suitable conditions for the singularities of the Gamma functions, sound very similar to the Le Roy-type functions, including the general multi-index case of the F m -function (see Section 3, Definition 7), introduced and studied in our recent works [13,14,15], and of course, the original Le Roy function.
To this end, we need to compare the corresponding Mellin–Barnes (M-B)-type integrals (14) for the H ¯ -functions with the M-B-type representations for the Le Roy functions (from our above-mentioned works), but also to compare the corresponding power series and their behaviors.
First, we provide a reminder of a general result for the “computational“ representation of some specific cases of the I-function and H ¯ -function in the form of power series. Such a representation for an I-function with β k = 1 , k = 1 , , m (named bar-I-function), and when all the poles of Γ ( b k B k s ) , k = 1 , , m are simple, can be found in Rathie [10] (6.8), and, resp., in Saxena [47] (3.1), where some techniques resembling the residue theorem are used. In the simpler case of the bar-H-function, when it is supposed additionally that α j = 1 , j = n + 1 , , p (see definition (14)), we can cite the following result.
Theorem 1.
The function H ¯ p , q m , n ( z ) can be represented by means of the power series ([10,47]):
H ¯ p , q m , n ( z ) = r = 0 h = 1 m j = 1 n Γ α j ( 1 a j + A j ξ ) j = n + 1 p Γ 1 ( a j A j ξ ) · k = 1 , k h m Γ 1 ( b k B k ξ ) k = m + 1 q Γ β k ( 1 b k + B k ξ ) · ( 1 ) r z ξ B h r ! .
Here, for shortness, ξ : = b h + r / B h . Series (18) exists for: 0 < | z | < if q 1 and μ > 0 ; or for 0 < | z | < R if μ = 0 , supposing that β k ( B h + ν 1 ) β h ( B k + ν 2 ) for k h ; k , h = 1 , , m ; ν 1 , ν 2 = 0 , 1 , 2 , .
One can observe an evident analogy of (18) with the power series representation for the Fox H-function from the handbooks, such as: [7] (§8.3.2: 3. and 4.), [4], etc.
Then, as a particular case, we can derive a series representation for the generalized Fox–Wright function (16), introduced here. Below, we deliberately add “stars” to the denotations of parameters and change some summation indices ( r k , but k i ) so as to distinguish these from the other traditional notations we use in the next sections, as previously adopted in several of the oldest and some recent works on Mittag-Leffler and Le Roy-type functions.
Theorem 2.
For the generalized Fox–Wright function p Ψ ˜ q we have: m = 1 , n = p , and assume parameters b 0 * = 0 , B 0 * = 1 , β 0 * = 1 (taken as the first ones of the q + 1 parameters in the H ¯ p , q + 1 1 , p -function below). Then, ξ : = b 0 * + k / B 0 * = k , and ( 1 ) k z ξ = ( z ) k . Thus, from (18), we obtain the power series representation for (16):
p Ψ ˜ q ( a j * , A j * ; α j * ) j = 1 p ( b i * , B i * ; β i * ) i = 1 q z = H ¯ p , q + 1 1 , p z ( 1 a j * , A j * , α j * ) 1 p ( 0 , 1 ) , ( 1 b i * , B i * , β i * ) 1 q
= k = 0 j = 1 p Γ α j * ( A j * k + a j * ) i = 1 q Γ β i * ( B i * k + b i * ) · z k k ! .
The characteristic parameters (13) are written now as:
μ = 1 + i = 1 q β i * B i * j = 1 p α j * A j * , R = i = 1 q ( B i * ) β i * B i * / j = 1 p ( A j * ) α j * A j * .
Series (19) defines an entire function (that is, absolutely convergent for all 0 < | z | < ) if μ > 0 (the other condition for the H ¯ -function, here: q + 1 1 , is satisfied ad hoc), or analytical one in 0 < | z | < R if μ = 0 .
Observe the close analogy with the series representation (4) for the Fox–Wright function. Here, following the discussions after Definitions 3 and 4, suitable cuts are to be inserted in C so as to fix single-valued branches for z k and for the included multi-valued fractional powers of the Γ -functions.

3. Multi-Index M-L-P Functions of Le Roy Type as I - and H ¯ -Functions

The special function F ( γ ) , defined in the whole complex plane C by the power series
F ( γ ) ( z ) = k = 0 z k k ! γ = k = 0 z k Γ ( k + 1 ) γ , z C , γ > 0 ,
is known as the Le Roy function. It was introduced by the French mathematician Le Roy [48,49] who used it to study the asymptotics of the analytic continuation of the sum of power series. His idea and goals were similar to those of the Swedish mathematician Mittag-Leffler when introducing the function (7) in [28]. The Mittag-Leffler function is now very popular, and often called the Queen function of Fractional Calculus, see [2], also [50],
E α , β ( z ) = k = 0 z k Γ ( α k + β ) , z C , α > 0 , β C .
It is also known as the “fractional” exponential function due to the “fractional” parameter α in the Γ -function replacing the factorial k ! in the exponential series, while, for (20), it emphasizes that the “fractional” index γ > 0 appears as the power of the factorial k ! (further extended to a Γ -function or to several Γ -functions).
In recent years, several studies appeared to revitalize the interest in (20) via its extensions and useful applications in various areas. It is important to mention that even in the simple case γ = 1 / 2 , the Le Roy function (20) appeared as R ( z ) in a study by Kolokoltsov [51] in quantum stochastics and particle systems. As he himself commented, “the function R ( z ) plays the same role for stochastic equations as the exponential function and Mittag-Leffler functions for deterministic equations”.
Among its generalizations, as a starting point, we mention the Mittag-Leffler function of Le Roy-type (MLR-function) as an extension of both the Le Roy and Mittag-Leffler functions:
F α , β ( γ ) ( z ) = k = 0 z k Γ ( α k + β ) γ , α > 0 , β > 0 , γ > 0 .
It was introduced by Gerhold [52] and Garra and Polito [53], and studied further by Garrappa–Rogosin–Mainardi [54] and Garra–Orsingher–Polito [55], and then also by Gorska–Horzela [56], Simon [57], Mehrez–Das [58] and Mehrez [59]. In [55], it is mentioned that the Le Roy-type functions (21) are used in probability in the context of the studies of COM-Poisson distributions (in the sense of Conway and Maxwell), which are special classes of weighted Poisson distributions. It is shown that as a first probabilistic application, these functions can also be useful in the construction of a new generalization of the COM-Poisson distribution that can be interesting for statistical applications and in physics in the case of generalized coherent states.
Next, a multi-index analogue of (21), inspired by the multi-index Mittag-Leffler functions (9) (Kiryakova [38,39,40]; Luchko et al. [35,60]), was considered by Rogosin and Dubatovskaya [45,61] with α i > 0 , β i > 0 , γ i > 0 (and also under more general conditions):
F ( α , β ) m ( γ ) m ( z ) = k = 0 z k Γ ( α 1 k + β 1 ) γ 1 Γ ( α m k + β m ) γ m .
The following Le Roy-type analogue of the Prabhakar (three-parameter Mittag-Leffler) function (8), that we call the Prabhakar function of Le Roy type, was introduced by Tomovski and Mehrez [62], and studied by Paneva-Konovska [15]:
F α , β ; τ ( γ ) ( z ) = k = 0 ( τ ) k Γ ( α k + β ) γ z k k ! , where ( τ ) k = Γ ( τ + k ) Γ ( τ ) , τ > 0 .
Further, in our recent works [13,14,16], we introduced and studied a very general Le Roy-type special function, as a multi-index analogue of the previously mentioned functions.
Definition 7.
The function, called the multi-index Mittag-Leffler–Prabhakar function of Le Roy-type (abbrev. as multi-MLPR), is defined by taking 4 m parameters α i , β i , τ i , γ i , i = 1 , , m :
F m ( z ) : = F α i ; β i ; τ i γ i ; m ( z )
= k = 0 ( τ 1 ) k ( τ m ) k Γ ( α 1 k + β 1 ) γ 1 Γ ( α m k + β m ) γ m · z k ( k ! ) m
= 1 i = 1 m Γ ( τ i ) · k = 0 i = 1 m Γ ( 1 . k + τ i ) Γ m 1 ( 1 . k + 1 ) i = 1 m Γ γ i ( α i . k + β i ) · z k k !
= k = 0 c k z k , w i t h c k = i = 1 m Γ ( k + τ i ) Γ ( k + 1 ) · 1 Γ ( τ i ) · 1 Γ ( α i k + β i ) γ i .
Note that in (25), we replaced part of the m factorials from (24) with ( k ! ) m 1 = Γ m 1 ( 1 . k + 1 ) .
In view of the applications of such kinds of functions and of their particular cases, here we consider only real positive parameters:
α i > 0 , β i > 0 , γ i > 0 , τ i > 0 , i = 1 , , m c o n d i t i o n i = 1 m α i γ i > 0 .
In [13] (Th.2), we proved that under these conditions, the F m -function (24) is an entire function of order ρ and type σ , with, respectively:
ρ = 1 α 1 γ 1 + + α m γ m , σ = 1 ρ i = 1 m ( α i ) α i γ i ρ .
The values of ρ and σ were also specified for more general requirements on the parameters, like
α i , β i , γ i , τ i C , R e ( α i ) > 0 , R e ( γ i ) > 0 , and i = 1 m R e ( α i γ i ) > 0 .
Results for the analytical properties of the multi-MLPR functions (24) are published in our recent works [13,14,16], including Mellin–Barnes-type contour integral representations, images under Laplace transform and Erdélyi–Kober operators of fractional integration, etc. For the case of the functions (22), see the corresponding results of Rogosin–Dubatovskaya [45,61].
Next, we develop the idea that the F m -function, and the other mentioned Le Roy-type functions, can be represented in terms of some special functions generalizing the Fox H-function and the Fox–Wright function p Ψ q , namely, as some I-, H ¯ - and p Ψ ˜ q -functions. The following result is derived with a short sketch of the details. For shortness, we further denote T : = i = 1 m 1 / Γ ( τ i ) .
Theorem 3.
Under the assumptions (26), the multi-index Mittag-Leffler–Prabhakar function of Le Roy type (24) can be represented in terms of the generalized Fox–Wright function (16), and thus, also as I-function (11) and H ¯ -function (14). Namely,
F m ( z ) : = F α i ; β i ; τ i γ i ; m ( z )
= T · m Ψ ˜ 2 m 1 ( τ i , 1 , 1 ) 1 m ( 1 , 1 , 1 ) ( m 1 ) t i m e s , ( β i , α i , γ i ) 1 m z
= T · H ¯ m , 2 m 1 , m z ( 1 τ i , 1 , 1 ) 1 m ( 0 , 1 ) m t i m e s , ( 1 β i , α i , γ i ) 1 m
= T · I m , 2 m 1 , m z ( 1 τ i , 1 , 1 ) 1 m ( 0 , 1 , 1 ) m t i m e s , ( 1 β i , α i , γ i ) 1 m .
Proof. 
For the power series representation (24) of the function F m , we apply Theorem 2. Here, p = m , q = 2 m 1 and we use the sets of parameters denoted as follows:
a j * = τ j , A j * = 1 , α j * = 1 , j = 1 , , m ; b i * = 1 , B i * = 1 , β i * = 1 , i = 1 , , m 1 ; b i + m 1 * = β i , B i + m 1 * = α i , β i + m 1 * = γ i , i = 1 , , m .
In the case of the F m -function, for the related H ¯ - and I-functions, we have the conditions of existence (13) satisfied, namely: q + 1 = 2 m 1 , and the other characteristic parameter is: μ = 1 + ( m 1 ) + i = 1 m α i γ i m = i = 1 m α i γ i > 0 , according to the assumed conditions (26). This allows us to conclude from Theorem 2 that (28)–(30) are entire functions of z in the complex plane.
For the definition of the H ¯ -function (29), and, resp., of the I-function (30), as was mentioned before, the path of integration can be either a shifted imaginary axis L = ( c i , c + i ) or a “Slater”-type contour L + starting at + + i φ 2 and finishing at + + i φ 1 , with suitable φ 1 < φ 2 (following the ideas from Marichev [5] (Ch.4), Rogosin–Dubatovskaya [45]). It is supposed that these contours are chosen not to cross the branch cuts and to separate the singularities of the Gamma functions in the numerator of the corresponding integrals.
Indeed, let us analyze the singularities of the involved Gamma functions in the representation of (29):
F m ( z ) = T · H ¯ m , 2 m 1 , m z ( 1 τ i , 1 , 1 ) 1 m ( 0 , 1 ) m t i m e s , ( 1 β i , α i , γ i ) 1 m
= T 2 π i L Γ 1 ( 1 s ) · i = 1 m Γ 1 ( τ i + s ) Γ m 1 ( 1 + s ) · i = 1 m Γ γ i ( β i + α i s ) z s d s .
The poles of Γ ( 1 s ) are s l = 1 + l , l = 0 , 1 , 2 , ; therefore, it has no singularities for s < 1 . The functions Γ ( τ i + s ) have poles at s i k = τ i k , k = 0 , 1 , 2 , , and have no singularities for s > m 0 , where m 0 : = min { τ 1 , , τ m } . Then, the Gamma functions in the numerator have no singularities for s ( m 0 , 1 ) .
The poles of the ( m 1 ) functions Γ ( 1 + s ) are s l = 1 l , l = 0 , 1 , 2 , , i.e., s = 1 , 2 , 3 , . Finally, for the Γ γ i ( β i + α i s ) , i = 1 , , m in the denominator, the singularities appear at s i n = β i α i n α i β i α i , n = 0 , 1 , 2 , and so, if we denote m 0 ˜ : = min { β 1 / α 1 , , β m / α m } , there will be no their singularities for s > m 0 ˜ . When the γ i are not integers, to avoid multi-valueness, we need to make a branch cut to the left of s = m 0 ˜ < 0 .
It is now evident that the vertical strip ( c i , c + ) with c min ( m 0 , m 0 ˜ ) , 1 is free of any singularities of the involved Gamma functions. However, to have a contour that does not intersect the branch cuts to the left of s = m 0 ˜ < 0 and also to the left of s < 0 (for single values of z s ), we can take a contour L = ( c * i , c * + i ) with c * ( 0 , 1 ) to ensure the existence of the bar-H-function in (29), and therefore also for (30). □

4. Related “Eigen”-Operators for Some Classes of SF: Gelfond–Leontiev Operators and Operators of FC

In analysis, linear algebra, physics, etc., the notions related to the prefix “Eigen” (coming from the German word meaning “self” or “own”) play important roles. In short, an eigenfunction is a function that, when acted on by an operator, yields a scalar multiple of the function itself. The scalar value is called the eigenvalue. The eigenfunctions and eigenvalues are also important, for example, in quantum mechanics, because they allow us to describe the behavior of operators in a way that is easy to understand and calculate.
Aside from many other analytical properties of the special functions, and the evaluation of their images under integral transforms and operators of fractional calculus, it is an important but often still open problem to determine corresponding linear integral L and differential operators D that transform a function f into itself just multiplied by a scalar, e.g., D f = λ f . For shortness, here, we call such operators “eigen“-operators (eigenoperators) for the function f.
It so happens that a useful tool to resolve this (generally open) problem for some classes of special functions is the notion of Gelfond–Leontiev operators (G-L operators) for generalized integration and differentiation, introduced in these authors’ work [63] of 1951. Up to this point, we have developed the theory of the G-L operators to propose corresponding integral and differential operators for which the Mittag-Leffler functions (7) and their multi-index analogues (9) are eigenfunctions. It appears that these are also operators of Generalized Fractional Calculus (in the sense of Kiryakova [20]). See, for example, [20] (Ch.5), [37,38], etc., and survey [64].
Next, our goal is to propose eigenoperators corresponding to the considered Le Roy-type functions.

4.1. Preliminaries for the Gelfond–Leontiev Operators

We first provide a reminder of some of the basics of the theory of generalized integration and differentiation, and specifically, the notion of such Gelfond–Leontiev operators.
Consider a function
f ( z ) = k = 0 a k z k ,
analytic in a disk Δ R = { | z | < R } in the complex plane (in particular, it can be the unit disk with R = 1 ). Let b k k = 0 be an arbitrary sequence (multipliers’ sequence) satisfying suitable conditions and defining another analytic function in the same disk Δ R : b ( z ) = k = 0 b k z k . Then, let us define the operator
D b ; f ( z ) = ( b f ) ( z ) : = k = 0 a k b k z k ,
where ∘ denotes the Hadamard product (convolution). If b k for k , we can consider the operator D b ; f as a generalized differentiation. For b k 0 , k = 1 , 2 , , the inverse operation, or the convolution with the “reciprocal” function b * ( z ) = k = 0 z k / b k ,
I b ; f ( z ) = k = 0 a k b k z k = D b * ; f ( z )
can be called a generalized integration, see, e.g., Samko–Kilbas–Marichev [65] (§22.3).
Following this idea, the so-called Gelfond–Leontiev operators of generalized integration and differentiation can be considered, as introduced in [63]. There, the multipliers’ sequences b k k = 0 and 1 / b k k = 0 are constructed by means of the coefficients of an entire function.
Definition 8.
Let us take an entire function
ϕ ( λ ) = k = 0 ϕ k λ k
with a growth (order ρ > 0 and type σ 0 ) such that lim k k 1 ρ | ϕ k | k = ( σ e ρ ) 1 ρ . In [63], the following operator of generalized differentiation is introduced:
f ( z ) = k = 0 a k z k D ϕ D ϕ f ( z ) = k = 1 a k ϕ k 1 ϕ k z k 1 ,
together with its n-th powers D ϕ n , n = 0 , 1 , 2 , . We call the operation (32) a Gelfond–Leontiev (G-L) operator of generalized differentiation with respect to the function ϕ ( λ ) , and then can introduce the corresponding G-L operator of generalized integration, defined by
L ϕ f ( z ) = k = 0 a k ϕ k + 1 ϕ k z k + 1 .
Evidently, D ϕ L ϕ f ( z ) = f ( z ) .
That is, we take b k = ϕ k 1 ϕ k and 1 / b k + 1 = ϕ k + 1 ϕ k 0 . The conditions required for ϕ ( λ ) always hold for lim sup k . We suppose that there exists lim k ϕ k 1 ϕ k k = 1 , and therefore, based on the Cauchy–Hadamard formula, the above-defined series (31)–(33) for f ( z ) , D ϕ f ( z ) , L ϕ f ( z ) have the same radius of convergence (in particular, this can be R = 1 ).
Example 1.
For the simplest example we take ϕ ( λ ) = exp λ , with ϕ k = 1 / k ! = 1 / Γ ( k + 1 ) , that is b k = k , k = 0 , 1 , 2 , . Then, the operators (32) and (33) are the differentiation and integration of 1st order:
D f ( z ) = k = 1 k a k z k 1 = d d z f ( z ) , L f ( z ) = D 1 f ( z ) = k = 0 a k k + 1 z k + 1 = 0 z f ( ξ ) d ξ .

4.2. G-L Operators Related to M-L Functions and to FC Operators

In classical fractional calculus (FC), the most often-used definition of the integration of an arbitrary (fractional) order δ 0 is the Riemann–Liouville (R-L) operator of fractional integration (left-hand sided), defined by the formula
I δ f ( z ) = I 0 + δ f ( z ) = 1 Γ ( δ ) 0 z ( z ξ ) δ 1 f ( ξ ) d ξ
= z δ Γ ( δ ) 0 1 ( 1 σ ) δ 1 f ( z σ ) d σ , I 0 f ( z ) : = f ( z ) .
Then, the corresponding R-L fractional derivative, and, resp., the Caputo fractional derivative of order δ > 0 , are defined by compositions of the n-th order casual derivative ( n 1 < δ n , n N ) and fractional order integral of form (34), taken in the two cases in reverse order:
D δ f ( z ) = R L D δ f ( z )
: = D n I n δ f ( z ) = d d z n 1 Γ ( n δ ) 0 z ( z ξ ) n δ 1 f ( ξ ) d ξ ,
C D δ f ( z ) : = I n δ D n f ( z ) = 1 Γ ( n δ ) 0 z ( z ξ ) n δ 1 f ( n ) ( ξ ) d ξ , D 0 f ( z ) : = f ( z ) .
The theory of fractional calculus (FC) based on these (classical) operators is widely presented in many basic handbooks such as [4,65,66], etc. For its historical development, see survey [67].
Under these operators of FC, the images of analytic functions f ( z ) of the form (31) are as follows:
I δ f ( z ) = k = 0 a k Γ ( k + 1 ) Γ ( k + δ + 1 ) z k + δ , D δ f ( z ) = k = 0 a k Γ ( k + 1 ) Γ ( k δ + 1 ) z k δ .
This suggests that by replacing the multiplier sequences from Example 1 with more general ones, one can define G-L-type operators corresponding to some special functions.
Example 2.
The G-L operator of generalized integration generated by the Mittag-Leffler entire function ϕ ( λ ) : = E α , 1 ( λ ) = E α ( λ ) has the form
J α f ( z ) = k = 0 a k Γ ( α k + 1 ) Γ ( α k + α + 1 ) z k + 1 .
The following result for (38) is given in Samko–Kilbas–Marichev [65] (§22.3):
The G-L operator of generalized integration J α generated by the entire function E α ( λ ) is the following modification of the R-L integral of order α:
J α f ( z ) = 1 Γ ( α ) 0 1 ( 1 t ) α 1 f ( z t α ) d t = z 1 Γ ( α ) 0 z z 1 / α ξ 1 / α α 1 f ( ξ ) d ( ξ 1 / α ) .
The proof uses the known definition of the Beta-function
Γ ( α k + 1 ) Γ ( α k + α + 1 ) = B ( α , α k + 1 ) Γ ( α ) = 1 Γ ( α ) 0 1 ( 1 t ) α 1 t α k d t .
This integral representation allows us to extend the definition of the integration operator (38) from analytic functions f ( z ) in a disk to continuous functions (and even to integrable ones), given in a (complex) domain starlike with respect to the origin z = 0 .
A more general “classical” operator of fractional integration is the Erdélyi–Kober (E-K) fractional integral of order δ > 0 , with two additional parameters: real “weight” γ and some β > 0 :
I β γ , δ f ( z ) = z β ( γ + δ ) Γ ( δ ) 0 z ( z β ξ β ) δ 1 ξ β γ f ( ξ ) d ( ξ β )
= 1 Γ ( δ ) 0 1 ( 1 t ) δ 1 t γ f ( z t 1 / β ) d t = 0 1 β ( 1 σ β ) δ 1 σ β γ + β 1 Γ ( δ ) f ( z σ ) d σ
= 0 1 H 1 , 1 1 , 0 σ ( γ + δ + 1 1 β , 1 β ) ( γ + 1 1 β , 1 β ) f ( z σ ) d σ .
For γ = 0 , β = 1 , it reduces to the R-L fractional integral, namely: I 1 0 , δ f ( z ) = z δ I δ f ( z ) .
The involved additional parameters give more freedom and thus, much wider application, see, for example, Sneddon [68] and Kiryakova [20] (Ch.2). The corresponding Erdélyi–Kober (E-K) fractional derivative is introduced in [20] (Ch.2), (1.6.9) in the form
D β γ , δ f ( z ) = D η I β γ + δ , η δ f ( z ) ,
with η 1 < δ η , η N , and D η : = j = 1 η ( 1 β z d d z + γ + j ) , a   polynomial P ( d d z ) .
Let Ω be a complex domain, starlike with respect to the origin z = 0 (in particular, it can be a disk | z | < R ). We can consider the R-L and E-K fractional integrals, the corresponding R-L and E-K fractional derivative, and the operators of Generalized FC that are defined in Section 4.3 in functional spaces of the form
H μ ( Ω ) = f ( z ) = z μ f ˜ ( z ) ; μ 0 ; f ˜ H ( Ω ) the   space of   functions   analytic   in Ω , where 0 arg z < 2 π , i . e . ,   with   a   cut   along   the   positive   half - line { Re z 0 , Im z = 0 } , or π arg z < π , i . e . ,   with   a   cut   along   the   negative   half - line { Re z 0 , Im z = 0 } ,
so as to avoid multiplicities in members like z μ + k in f ( z ) .
In Kiryakova [20] (Ch.4), [69], etc., it is shown that the E-K operators transform the p F q - and p Ψ q -functions into the same kind of functions with increased orders p and q: resp., into p + 1 F q + 1 - and p + 1 Ψ q + 1 -functions. For the Le Roy-type functions (24), the images under the Erdélyi–Kober (E-K) fractional integration operator (40) are provided in our recent paper [14], and in particular cases, for the Riemann–Liouville fractional integrals (34) of the function (22), by Rogosin–Dubatovskaya [45]. Note also that R-L fractional integrals of the H ¯ -function are presented in Srivastava et al. [70].

4.3. G-L Operators for Multi-Index M-L Functions and Generalized FC

To explain the construction of the G-L operators generated by the multi-index M-L functions (9), we need to introduce briefly the notions of Generalized Fractional Calculus (GFC) from Kiryakova [20] and other works. Fractional calculus, as an exotic theory of the differentiation and integration of arbitrary (not-integer) order, arose in the 19th century, and has nowadays become a hot topic with a very wide scope of application. To read about its development, the reader can consult, e.g., the survey by Machado–Kiryakova [67]. We need to pay honor to Kalla [71] who introduced the notion of generalized operators of fractional integration by using an arbitrary special function as a kernel. Then, he and many other authors explored kernels such as the Gauss hypergeometric function, the Bessel functions, the Appell F 3 -function, and in general, the Fox H-function; for details, see [72]. However, the choice of suitable kernels as G m , m m , 0 - and H m , m m , 0 -functions was crucial for the possibility to develop a detailed and comparatively full theory of GFC [20].
Definition 9.
Take an integer (multiplicity) m 1 . By means of the sets of parameters: δ = ( δ 1 0 , , δ m 0 ) as a multi-order of fractional integration, γ = ( γ 1 , , γ m ) as a real multi-weight and additional multi-parameter β = ( β 1 > 0 , , β m > 0 ) , we define the integral operator:
I ( β k ) , m ( γ k ) , ( δ k ) f ( z ) = 0 1 H m , m m , 0 t ( γ i + δ i + 1 1 β i , 1 β i ) i = 1 m ( γ i + 1 1 β i , 1 β i ) i = 1 m f ( z t ) d t , i f i = 1 m δ i > 0 ,
and I ( β k ) , m ( γ k ) , ( δ k ) f ( z ) = f ( z ) if δ 1 = δ 2 = = δ m = 0 . We call it a multiple (m-tuple) Erdélyi–Kober fractional integration operator. More generally, all operators of the form
I f ( z ) = z δ 0 I ( β k ) , m ( γ k ) , ( δ k ) f ( z ) with δ 0 0 ,
are our generalized ( m -tuple) fractional integrals of multi-order ( δ 1 0 , , δ m 0 ) .
The kernel function in (42) is a suitably chosen case H m , m m , 0 of the Fox H-function (1), and when β i = β > 0 , i = 1 , , m , we have slightly simpler generalized fractional integrals with the kernel as Meijer G m , m m , 0 -function (3), see [20] (Ch.1).
The generalized fractional derivatives D ( β k ) , m ( γ k ) , ( δ k ) of multi-order ( δ 1 , , δ m ) , corresponding to (42) and (43), are also introduced explicitly, in a way similar to that for D δ and D β γ , δ but with more complicated differ-integral representations. See [20] and also [73,74].
Here, the operators of GFC are considered in the spaces H μ ( Ω ) of analytic functions with power weights of the form (41). In [20] (Ch.1, Ch.5) and in subsequent works, see, e.g., survey [73], the properties of these operators of GFC in various other functional spaces are studied and a full chain of operational rules is proposed, with a wide range of applications and many particular examples. Here, we mention only the following basic result for functions analytic in a disk that are useful for our next considerations.
Let the following conditions on the parameters be satisfied: β k ( γ k + 1 ) > μ , δ k 0 , k = 1 , , m . Then, the multiple E-K integral (42) maps the class H μ ( Ω ) of weighted analytic functions f ( z ) into itself. In particular, the images of the functions f ( z ) analytic in a disk Δ R = | z | < R :
f ( z ) = z μ k = 0 a k z k = z μ ( a 0 + a 1 z + ) H μ ( Δ R ) ,
have the same form:
I ( β k ) , m ( γ k ) , ( δ k ) f ( z ) = z μ k = 0 a k i = 1 m Γ ( γ i + k + μ β i + 1 ) Γ ( γ i + δ i + k + μ β i + 1 ) z k H μ ( Δ R ) ,
with the same radius of convergence R = lim sup k | a k | k 1 > 0 .
For m = 1 , (42) is the “classical” E-K integral (40). A very important property of the single-integral operators (42) with H-functions kernels (or with Meijer’s G-functions in the simpler case of equal β i = β > 0 , i = 1 , , m ) is the so-called composition/decomposition rule: they can also be represented by means of commutative compositions of classical E-K integrals ( m = 1 ):
I ( β k ) , m ( γ k ) , ( δ k ) f ( z ) = i = 1 m I β i γ i , δ i f ( z )
= 0 1 0 1 i = 1 m ( 1 t i ) δ i 1 t i γ i Γ ( δ i ) f z t 1 1 β 1 σ m 1 β m d t 1 d t m ,
that is, by iterated integrals without special functions involved in the kernel. This explains the wide and efficient use of the operators of the GFC from [20] because compositions like (46) often appear in problems related to applications, but their study and evaluation is based on the simple but effective tools of the theory of the kernel special functions (H- and G-functions) in definition (42).
Next, let us consider the G-L-type integration and differentiation operators we introduced as generated by the M-L function (7) and by the multi-index M-L function (9). The series representations for these operators can be analytically continued in the functional spaces H μ ( Ω ) as particular cases of the generalized fractional integrals I ( β k ) , m ( γ k ) , ( δ k ) and derivatives D ( β k ) , m ( γ k ) , ( δ k ) of fractional multi-order. We called these operators multiple Dzrbashjan–Gelfond–Leontiev operators, so as to honor the Armenian mathematician Dzrbashjan, author of the book [75], one of the most detailed studies of the 20th century on the M-L function (7) in the complex plane. In [29], he also introduced a 2 × 2 ( m = 2 ) M-L-type function as the first representative of the multi-index M-L functions (9). Now, we refer to these operators as G-L operators for the multi-index M-L functions, defined as follows.
Let f ( z ) be an analytic function in a disk Δ R = { | z | < R } and α i > 0 , β i R , i = 1 , , m be arbitrary parameters. The operators:
f ( z ) = k = 0 a k z k D ( α i ) , ( β i ) f ( z ) and L ( α i ) , ( β i ) f ( z ) ,
constructed as
D ( α i ) , ( β i ) f ( z ) = k = 1 a k Γ ( α 1 k + β 1 ) Γ ( α m k + β m ) Γ α 1 ( k 1 ) + β 1 ) Γ ( α m ( k 1 ) + β m z k 1 , L ( α i ) , ( β i ) f ( z ) = k = 0 a k Γ ( α 1 k + β 1 ) Γ ( α m k + β m ) Γ α 1 ( k + 1 ) + β 1 ) Γ ( α m ( k + 1 ) + β m z k + 1 ,
are called G-L operators generated by the multi-index M-L function (9).
Evidently, D ( α i ) , ( β i ) L ( α i ) , ( β i ) f ( z ) = f ( z ) for f ( z ) analytic in Δ R , and the coincidence of the radii of convergence of f ( z ) and of the series in (47) easily follows by the Cauchy–Hadamard formula and the asymptotic estimation of the Γ -function multipliers (details are, e.g., in [20] (Th.5.5.2)). In Kiryakova [20] (§5.4.ii) and subsequent works such as [37,39], etc., we provided the following integral representation of the G-L integration L ( α i ) , ( β i ) as its analytical extension in terms of the GFC operators of multi-order ( α 1 , , α m ) .
Theorem 4.
Let Ω Δ R be a domain in ℂ starlike with respect to the origin z = 0 . Then, the multi-index G-L integration operator in (47) can be analytically continued from H ( Δ R ) into H μ ( Ω ) by means of the single-integral operator
L ( α i ) , ( β i ) f ( z ) = z 0 1 H m , m m , 0 σ ( β i , α i ) 1 m ( β i α i , α i ) 1 m f ( z σ ) d σ = z I ( 1 / α i ) , m ( β i 1 ) , ( α i ) f ( z ) ,
that is, by a generalized fractional integral of form (43).
Here, we concentrate on the G-L operators of generalized integration and their representations by generalized fractional integrals, and skip the details on the corresponding differentiation analogues. However, let us mention that the multi-index G-L derivative in (47) can be represented by generalized fractional derivatives D ( β k ) , m ( γ k ) , ( δ k ) ([20,73]), and for analytic functions in H ( Ω ) H ( Δ R ) , and under restrictions β i α i 0 , 1 , 2 , ; z 0 ,
D ( α i ) , ( β i ) f ( z ) = z 1 D ( 1 / α i ) , m ( β i α i 1 ) , ( α i ) f ( z ) i = 1 m Γ ( β i ) Γ ( β i α i ) f ( 0 ) z .
By the rules of the GFC one can check that D ( α i ) , ( β i ) L ( α i ) , ( β i ) f ( z ) = f ( z ) .
The following result from [38,39,40] shows that the multi-index M-L functions (9) are “eigenfunctions” of the multi-index G-L operators (47)–(49). Namely:
The multi-index Mittag-Leffler functions (9) satisfy the integral relation
L ( α i ) , ( β i ) E ( α i ) , ( β i ) ( λ z ) = 1 λ E ( α i ) , ( β i ) ( λ z ) 1 λ i Γ ( β i ) , λ 0 ,
and the differential equation of fractional multi-order ( α 1 , , α m ) :
D ( α i ) , ( β i ) E ( α i ) , ( β i ) ( λ z ) = λ E ( α i ) , ( β i ) ( λ z ) , λ 0 .
Then, we may say that the G-L operators (47)–(49) are “eigen”-operators for the multi-index Mittag-Leffler functions (9).
Example 3.
For m = 1 , the operators (47) are the Dzrbashjan–Gelfond–Leontiev differentiation and integration operators, generated by the Mittag-Leffler function (7), as introduced in Kiryakova [20] (Ch.2):
D α , β f ( z ) = k = 1 a k Γ ( α k + β ) Γ ( α ( k 1 ) + β ) z k 1 , L α , β f ( z ) = k = 0 a k Γ ( α k + β ) Γ ( α ( k + 1 ) + β ) z k + 1 .
And their analytical continuations appear as (single , m = 1 ) E-K fractional integrals (40) and corresponding E-K derivatives:
L α , β f ( z ) = z 1 I 1 / α β 1 , α f ( z ) = z Γ ( α ) 0 1 ( 1 t ) α 1 t β 1 f ( z t α ) d t , D α , β f ( z ) = z 1 D 1 / α β α 1 , α f ( z ) f ( 0 ) Γ ( β ) Γ ( β α ) z 1 , z 0 .
For β = 1 , one has the simpler G-L operators (38) and (39) from Samko–Kilbas–Marichev [65], as in the previous Example 2.

4.4. G-L Operators Generated by the Le Roy-Type Functions F ( α , β ) m ( γ ) m

We consider only the case of the multi-index Mittag-Leffler functions (22) of Le Roy type F ( α , β ) m ( γ ) m , as the simpler cases introduced in Rogosin and Dubatovskaya’s articles [45,61]. It seems that the more general case of special functions (24) with additional Prabhakar parameters τ i 1 , i = 1 , , m , still cannot be treated with respect to corresponding eigenoperators (in the sense of Gelfond–Leontiev operators). At this point, we consider it to be an open problem.
First, let us present the corresponding (slightly simpler) representations of these Le Roy-type functions in terms of the generalizations of the Fox H-function and Fox–Wright function p Ψ q . As a corollary of Theorem 3, bearing in mind that τ i = 1 , T = 1 , and that some equal sets of parameters in the upper and lower rows in the Ψ ˜ - and I-functions cancel each other out (Property (7.1) in Rathie [10], quite similar to that for the H-functions, e.g., [7], etc.), we have:
F ( α , β ) m ( γ ) m ( z ) = k = 0 z k i = 1 m Γ γ i ( α i k + β i )
= m Ψ ˜ 2 m 1 ( 1 , 1 , 1 ) 1 m ( 1 , 1 , 1 ) ( m 1 ) t i m e s , ( β i , α i , γ i ) 1 m z
= 1 Ψ ˜ m ( 1 , 1 , 1 ) ( β i , α i , γ i ) 1 m z
= H ¯ m , 2 m 1 , m z ( 0 , 1 , 1 ) 1 m ( 0 , 1 , 1 ) m t i m e s , ( 1 β i , α i , γ i ) 1 m
= H ¯ 1 , m + 1 1 , 1 z ( 0 , 1 , 1 ) ( 0 , 1 , 1 ) , ( 1 β i , α i , γ i ) 1 m
= I m , 2 m 1 , m z ( 0 , 1 , 1 ) 1 m ( 0 , 1 , 1 ) m t i m e s , ( 1 β i , α i , γ i ) 1 m
= I 1 , m + 1 1 , 1 z ( 0 , 1 , 1 ) ( 0 , 1 , 1 ) , ( 1 β i , α i , γ i ) 1 m .
Now, let us construct the Gelfond–Leontiev generalized operators of integration and differentiation (32) and (33), generated by the coefficients of the mMLR-function F ( α , β ) m ( γ ) m :
ϕ k = 1 / i = 1 m Γ γ i ( α i k + β i ) .
For an entire function f ( z ) = k = 0 a k z k (or analytic in a disk), these have the following form of power series, to be a generalized G-L differentiation:
D f ( z ) : = D m M L R f ( z ) = k = 1 a k z k 1 · i = 1 m Γ γ i ( α i k + β i ) Γ γ i ( α i k + β i α i ) ,
and, respectively, a generalized G-L integration:
L f ( z ) : = L m M L R f ( z ) = k = 0 a k z k + 1 · i = 1 m Γ γ i ( α i k + β i ) Γ γ i ( α i k + β i + α i ) .
Then, D L f ( z ) = f ( z ) .
Theorem 5.
The generalized G-L operator of differentiation (32) is an eigenoperator for the function F ( α , β ) m ( γ ) m , that is, F ( α , β ) m ( γ ) m is an eigenfunction for this “differentiation” operator:
D F ( α , β ) m ( γ ) m ( λ z ) = D m M L R F ( α , β ) m ( γ ) m ( λ z ) = λ F ( α , β ) m ( γ ) m ( λ z ) , λ 0 .
For the G-L integration, we obtain the corresponding relation
L F ( α , β ) m ( γ ) m ( λ z ) = L m M L R F ( α , β ) m ( γ ) m ( λ z ) = 1 λ F ( α , β ) m ( γ ) m ( λ z ) 1 λ i = 1 m Γ γ i ( β i ) , λ 0 .
Proof. 
Indeed, using the power series representation (22),
D m M L R F ( α , β ) m ( γ ) m ( λ z ) = D m M L R k = 0 λ k i = 1 m Γ γ i ( α i k + β i ) z k
= k = 1 λ k i = 1 m Γ γ i ( α i k + β i ) z k 1 i = 1 m Γ γ i ( α i k + β i ) i = 1 m Γ γ i ( α i k + β i α i )
(where the groups of the Gamma functions Γ γ i ( α i k + β i ) cancel each other out, and we change the index of summation as k * : = k 1 , to have)
= λ k * = 0 ( λ z ) k * i = 1 m 1 Γ γ i ( α i k * + β i ) = λ F ( α , β ) m ( γ ) m ( λ z ) .
Analogously,
L F ( α , β ) m ( γ ) m ( λ z ) = L m M L R k = 0 λ k i = 1 m Γ γ i ( α i k + β 1 ) z k
= k = 0 λ k i = 1 m Γ γ i ( α i k + β i ) z k + 1 i = 1 m Γ γ i ( α i k + β i ) i = 1 m Γ γ i ( α i k + β i + α i )
= 1 λ k * = 1 ( λ z ) k * i = 1 m Γ γ i ( α i k * + β i ) = 1 λ F ( α , β ) m ( γ ) m ( λ z ) 1 λ i = 1 m Γ γ i ( β i ) .
Both relations (59) and (60) are exactly the expected analogues of those for the multi-index M-L functions, (51) and (50).
Now, we aim to find an integral representation of the G-L integral (58), as an analogue of the generalized fractional integration operator (48) in Theorem 4. Instead of the H-function H m , m m , 0 , we have now an I-function as a kernel. Then, we can admit that such an integral operator can also be somehow thought of as an operator of “fractional multi-order”.
Theorem 6.
The G-L integration operator (58) of an entire function f ( z ) , generated by means of the Le Roy-type function (22), can also be represented by means of the integral operator
I f ( z ) = L m M L R f ( z ) = z 0 1 I m , m m , 0 σ ( β i , α i , γ i ) 1 m ( β i α i , α i , γ i ) 1 m f ( z σ ) d σ .
This can be interpreted as a kind of a generalized fractional integration of multi-order ( α 1 , , α m ) .
We call the operators (58) and (61) the Gelfond–Leontiev–Le Roy (G-L-Le Roy) operators of generalized integration. To emphasize when these concern the general case of m = 1 , 2 , 3 , , we sometimes also use the notation I m and, for example, when m = 1 , the notation is I 1 . If necessary, sub-indices i are also used, I i 1 , so as to specify the relation to the particular parameters α i > 0 , β i > 0 , γ i > 0 .
Proof. 
First, let us analyse the nature of the kernel function I m , m m , 0 with the parameters as in (61). As an I-function, its definition (11) and (12) reads as
I m , m m , 0 ( σ ) = 1 2 π i L i = 1 m Γ γ i ( β i α i α i s ) i = 1 m Γ γ i ( β i α i s ) σ s d s .
The singular points for Γ γ i ( β i α i α i s ) are: s i k = β i α i + k α i 1 , k = 0 , 1 , 2 , , that is i = 1 , , m : s i k > 1 , and there are no singularities for s < 1 . For Γ γ i ( β i α i s ) , the singular points appear at s i l = β i α i + l α i , l = 0 , 1 , 2 , and all are s i l > 0 . Then, it seems that for s < 1 , none of the Gamma functions have singularities. The cut to ensure single values of the involved members can now be taken along the real half-line s 0 (see assumptions (41)). Therefore, we can have a contour L = ( c i , c + i ) with some c < 1 for which all singularities are to the right and which does not intersect the branch cut.
The parameters (13) for the I m , m m , 0 -function in (61) are as follows:
μ = i = 1 m α i γ i i = 1 m α i γ i = 0 ; Δ = i = 1 m α i γ i i = 1 m α i γ i = 0 ;
= i = 1 m γ i ( β i 1 / 2 ) i = 1 m γ i ( β i α i 1 / 2 ) = i = 1 m α i γ i > 0 ;
R = i = 1 m α i α i γ i / i = 1 m α i α i γ i = 1 .
This means that the kernel function I m , m m , 0 is an analytic function in the unit disk 0 < | z | < 1 and it vanishes for | z | > 1 . Then, the integral in (61) from 0 to 1, when necessary (for the evaluations below), can also be taken with limits from 0 to , having the same value.
Let us evaluate this integral operator of an entire function with power series f ( z ) = k = 0 a k z k :
I f ( z ) = z 0 1 I m , m m , 0 σ ( β i , α i , γ i ) 1 m ( β i α i , α i , γ i ) 1 m k = 0 a k z k σ k d σ
= z k = 0 a k z k 0 1 σ k I m , m m , 0 σ ( β i , α i , γ i ) 1 m ( β i α i , α i , γ i ) 1 m d σ
= z k = 0 a k z k 0 1 I m , m m , 0 σ ( β i + α i k , α i , γ i ) 1 m ( β i α i + α i k , α i , γ i ) 1 m d σ ,
where the exchanged order of integration and summation is admissible, and we also use the shifting property from Rathie [10] (7.3) for σ k I m , m m , 0 ( σ ) .
Now, we need to use an auxiliary result that in general will read as follows:
0 1 I m , m m , 0 σ ( a i , A i , α i ) 1 m ( b i , B i , β i ) 1 m d σ = 0 I m , m m , 0 σ ( a i , A i , α i ) 1 m ( b i , B i , β i ) 1 m d σ
= i = 1 m Γ β i ( b i + B i ) Γ α i ( a i + A i ) , for a i > b i > 0 .
This is an analogue of our “Auxiliary integral“ (Kiryakova [20] (App.E, (E.21))) for the similar H m , m m , 0 -function with a i > b i > 0 , i = 1 , , m :
0 1 H m , m m , 0 σ ( a i , A i ) 1 m ( b i , B i ) 1 m d σ = i = 1 m Γ ( b i + B i ) Γ ( a i + A i ) .
The simpler analogous formula for the case of G m , m m , 0 is proved in Kiryakova [20] (Lemma B.2), and the same one is repeated in Karp–Lopez [76].
For the H m , m m , 0 -function, Formula (64) follows by also taking s = 1 in the result for the Mellin transform of H p , m m , 0 -function (it can also be p = m ) from Karp–Prilepkina [77] (Th.6):
0 R ( o r ) σ s 1 H p , m m , 0 σ ( a i , A i ) ( b j , B j ) d σ = j = 1 m Γ ( b j + B j s ) i = 1 p Γ ( a i + A i s ) .
Analogously, the result for I m , m m , 0 follows from the Mellin transform image of the I-function ([78] (3.9)) taken for s = 1 . As in the case of the H-functions, the same result can be derived as a corollary by a more general integral formula, see Vellaisamy–Kataria [78] (Prop.3.1), and Lemma 1 in Section 4.5, involving a weighted product of two I-functions. For the corresponding formula in the case of the product of two H-functions, we refer to [3] (2.8.4) and [20] (E.21’).
For the particular parameters of the I m , m m , 0 -function here, from (63), we obtain
0 1 I m , m m , 0 σ ( β i + α i k , α i , γ i ) 1 m ( β i α i + α i k , α i , γ i ) 1 m d σ ( = 0 )
= i = 1 m Γ γ i ( β i α i + α i k + α i ) Γ γ i ( β i + α i k + α i ) = i = 1 m Γ γ i ( α i k + β i ) Γ γ i ( α i k + α i + β i ) ,
which we need to show that
I f ( z ) = I k = 0 a k z k = k = 0 a k z k + 1 i = 1 m Γ γ i ( α i k + β i ) Γ γ i ( α i k + α i + β i ) = L m M L R f ( z ) .
Remark 1.
Along with the check above that the two operators L m M L R f ( z ) and I f ( z ) are equivalent for entire functions f ( z ) given by power series, we can discuss the behavior of the kernel I m , m m , 0 around the singularities 0 and 1 in the improper integral (61). First, it is known (see, e.g., Rathie [10] (6.9)) that for small values of σ, i.e., for σ + 0 :
I m , m m , 0 σ ( β i , α i , γ i ) 1 m ( β i α i , α i , γ i ) 1 m σ c , w h e r e
c = min 1 i m { ( β i α i ) / α i } = min 1 i m { β i α i 1 } > 1 , f o r α i > 0 , β i > 0 .
The asymptotic behavior of the I-functions near the singular point σ = 1 (or in more general cases, near σ = R on the circle of convergence) is not yet well studied. The same situation stands for the Fox H-functions’ behavior at the third singular point when these are analytic in disks with finite radius | z | < R , or outside, in | z | > R . A particular result for the case of H p , m m , 0 -functions ( n = 0 , q = m ) with μ = j = 1 m B j i = 1 p A i = 0 is contained in Karp [21] (12.24). This can be paraphrased briefly as
H p , m m , 0 R · σ ( a i , A i ) 1 p ( b j , B j ) 1 m = ( 1 σ ) η 1 k = 0 h k ( 1 σ ) k Γ ( η + k ) ,
where in our case ( m = p = q ):
R = 1 a n d η = i = 1 p a i j = 1 m b j m p 2 = i = 1 p a i j = 1 m b j > 0 , t h a t i s , η 1 > 1 .
In the case of μ = j = 1 m B j i = 1 p A i = 0 (which is also true, in particular, if A i = B j ), the main term in the representation above depends only on the value of η given above.
The result for the simpler case of the Meijer G m , m m , 0 function (known also as the Meijer–Nørlund function) by Marichev (1981) is analogous, and was used by Kiryakova [20] (1.1.14):
G m , m m , 0 σ ( a k ) 1 m ( b k ) 1 m ( 1 σ ) η m * Γ ( η m * + 1 ) a s σ 1 , σ < 1 ,
with η m * = k = 1 m ( a k b k ) 1 > 1 in the case when (see [20] (Ch.1)) a k b k = ( γ k + δ k ) γ k = δ k > 0 , k = 1 , , m and k = 1 m δ k 1 , 2 , 3 , (to avoid the logarithmic case, which separately considered leads to similar asymptotics, ensuring η m * > 1 ). The same behavior follows from the result found by Karp–López [76] (8), for G m , m m , 0 ( σ ) in the disk | 1 σ | < 1 under the same conditions, showing that the behavior near the singular point σ = 1 depends only on this η > 1 .
In the case of I-function (62) of orders n = 0 , m = p = q , with parameters a i : = β i > β i α i : = b i , one has quite a similar result because of all of the equal second ( α i > 0 ) and third ( γ i > 0 ) parameters, i = 1 , , m . The main term of the asymptotics depends only on η:
I m , m m , 0 σ ( β i , α i , γ i ) ( β i α i , α i , γ i ) ( 1 σ ) η 1 , w i t h η = i = 1 m [ β i ( β i α i ) ] γ i = i = 1 m α i γ i > 0 η 1 > 1 , a s σ 1 , σ < 1 .
We clarify this situation in more detail for m = 1 in Section 4.5.
Alternatively, for the I m , m m , 0 -functions (62), one may interpret and specify the results (6.4) and (6.5) from Rathie [10] (Procedure 1) when μ = 0 and = i = 1 m α i γ i > 0 but is not integer. The technique used employs a known expansion for large s, involving the Bernoulli polynomials. Then, I ( σ ) can be expressed in a series of Beta functions in the unit disk that leads to similar asymptotics as those in (65). Otherwise, if = i = 1 m α i γ i is a positive integer, one falls in a logarithmic case, a result that is considered in [10] (6.6). In our situation (for shortness, we mention the result for m = 1 , = α γ = 1 , 2 , 3 , ), one can derive it in a form:
I m , m m , 0 ( σ ) ( 1 σ 1 / α ) α γ 1 ( α γ 1 ) ! + f o r σ 1 , σ < 1 w i t h ( α γ 1 ) > 1 .

4.5. Composition/Decomposition Property of the Gelfond–Leontiev–Le Roy Integrations

For the generalized fractional integrals (42) with H m , m m , 0 - (resp. G m , m m , 0 -) kernels in the Generalized Fractional Calculus (Kiryakova [20]), we have already mentioned the important composition/decomposition property (46). In the case of the G-L generalized integration in (53) for m = 1 , generated by the M-L function E α , β , we provide a reminder that its integral representation has the form
L α , β f ( z ) = z I 1 / α β 1 , α f ( z ) = z Γ ( α ) 0 1 ( 1 t ) α 1 t β 1 f ( z t α ) d t
= z 0 1 ( 1 σ 1 / α ) α 1 σ β / α 1 α Γ ( α ) f ( z σ ) d σ = z 0 1 H 1 , 1 1 , 0 σ ( β , α ) ( β α , α ) f ( z σ ) d σ ,
where the kernel H 1 , 1 1 , 0 analytic in | σ | < 1 can also be represented by the generalized binomial series
( 1 σ 1 / α ) α 1 σ β / α 1 α Γ ( α ) = Γ ( α ) σ β / α 1 α Γ ( α ) k = 0 α 1 k σ 1 / α k = σ β / α 1 α k = 0 σ 1 / α k Γ ( α k ) k ! .
Note that if α 1 = n , a non-negative integer, the k = n + 1 -term and all the next ones in the series are 0, i.e., the series is finite.
Next, we prove a similar result in that the Gelfond–Leontiev–Le Roy integrals (61) of “multi-order” ( α 1 , , α m ) can be represented as commutable compositions of m operators I 1 of the same kind for m = 1 , each one with different parameters. As mentioned before in the proof of Theorem 6, we use the notation I m for the case of arbitrary (multiplicity) m = 1 , 2 , 3 , , while a particular operator with m = 1 is denoted by I 1 , and we add a sub-index i for I i 1 to keep in mind that it is the case with parameters ( α i , β i , γ i ) . These operators I i 1 play the role of “Le Roy type” extensions of the Erdélyi–Kober operators (40) for the integration of “order” α i . In short, these can be called Erdélyi–Kober–Le Roy integrations (E-K-Le Roy), with each of them having a representation as for (61) with m = 1 :
I i 1 f ( z ) = z 0 1 I 1 , 1 1 , 0 σ ( β i , α i , γ i ) ( β i α i , α i , γ i ) f ( z σ ) d σ
= 0 z I 1 , 1 1 , 0 ζ z ( β i , α i , γ i ) ( β i α i , α i , γ i ) f ( ζ ) d ζ .
It is interesting to represent this kernel function I 1 , 1 1 , 0 in the way we did above for the H 1 , 1 1 , 0 , as a kind of similar “generalized” binomial series. Namely, for | σ | < 1 , γ > 0 , one has
I 1 , 1 1 , 0 σ ( β , α , γ ) ( β α , α , γ ) = σ β / α 1 α k = 0 σ 1 / α k Γ γ ( α γ k ) k ! .
The main term in such series gives an asymptotics as
I 1 , 1 1 , 0 σ ( β , α , γ ) ( β α , α , γ ) σ β / α 1 ( 1 σ 1 / α ) α γ 1 α Γ γ ( α γ ) ,
with β / α 1 > 1 for σ 0 , σ > 0 , and α γ 1 > 1 for σ 1 , σ < 1 .
This is to confirm again, for the simplest case of m = 1 , the described behavior of the kernel I m , m m , 0 near the singular points 0 and 1 in the improper integral (66).
Before proving the decomposition of the operator I m , we need the following auxiliary result that appears analogous to the formula for the Mellin transform, and thus, for an integral from 0 to of the product of two different H-functions. This can be found, for example, in some handbooks such as [3] (2.8.4) and (E.21’) in Kiryakova [20].
Here, it is derived from Proposition 3.1 in Vellaisamy–Kataria [78] for the Mellin transform of the product of two different I-functions. For our purposes, we reproduce it as an integral formula in a simpler form.
Lemma 1.
For λ 0 , ν 0 , s > 0 ,
0 σ s 1 · I p , q m , n λ σ ( a i , A i , α i ) 1 p ( b j , B j , β j ) 1 q · I u , v k , l ν σ ( c i , C i , φ i ) 1 u ( d j , D j , ψ j ) 1 v d σ
= 1 s · I q + u , p + v n + k , m + l ν λ ( c i , C i , φ i ) 1 l , ( 1 b j B j , B j , β j ) 1 q , ( c i , C i , φ i ) l + 1 u ( d j , D i , ψ j ) 1 k , ( 1 a i A i , A i , α i ) 1 p , ( d j , D j , ψ i ) k + 1 v .
A set of necessary operational properties for the I-functions that are rather similar to those for the Fox H-functions can be found in Rathie [10] (Sect. 7). We mention only a few of them that are needed below, such as: symmetry in some sets of parameters; shifting property (7.3) for σ c I p , q m , n ( σ ) ; and (7.5) for I p , q m , n ( σ ) I q , p n , m ( 1 / σ ) .
Now, we are ready to state and prove the composition/decomposition property of the Gelfond–Leontiev–Le Roy generalized integration.
Theorem 7.
For entire functions f ( z ) , the equivalent representations hold:
I m f ( z ) = i = 1 m I i 1 f ( z ) = I m 1 I m 1 1 I 1 1 f ( z ) .
This composition of operators is commutative.
Proof. 
Let us check the statement first for m = 2 . We consider the composition of two Le Roy-type Erdélyi–Kober operators of the form (66):
I 2 1 I 1 1 f ( z ) = 0 z I 1 , 1 1 , 0 ζ z ( β 2 , α 2 , γ 2 ) ( β 2 α 2 , α 2 , γ 2 ) I 1 1 f ( z ) d ζ
= 0 z I 1 , 1 1 , 0 ζ z ( β 2 , α 2 , γ 2 ) ( β 2 α 2 , α 2 , γ 2 ) 0 ζ I 1 , 1 1 , 0 τ ζ ( β 1 , α 1 , γ 1 ) ( β 1 α 1 , α 1 , γ 1 ) f ( τ ) d τ d ζ
= 0 z f ( τ ) d τ τ z I 1 , 1 1 , 0 ζ z ( β 2 , α 2 , γ 2 ) ( β 2 α 2 , α 2 , γ 2 ) · I 1 , 1 1 , 0 τ ζ ( β 1 , α 1 , γ 1 ) ( β 1 α 1 , α 1 , γ 1 ) d ζ
: = 0 z f ( τ ) I ˜ d τ .
The interchange of the order of integrals, as above, is well admissible. We have to note that both of the involved I 1 , 1 1 , 0 -functions in the above expression I ˜ vanish outside the unit disk. That is, the one with ( i = 1 ) parameters is 0 for ζ < τ , and the other one for ( i = 2 ) is 0 for ζ > z . Therefore, the above inner integral denoted by I ˜ will have the same value if its limits from τ to z are extended to an interval from 0 to . Moreover, we replace I 1 , 1 1 , 0 with the ( i = 1 ) parameters by I 1 , 1 0 , 1 of the reciprocal argument, according to the property (7.5) from Rathie [10] (quite similar to one for the H-functions) and then use Formula (67), taking λ : = z 1 , ν : = τ 1 and s = 1 . We have:
I 2 1 I 1 1 f ( z ) = 0 z I 2 , 2 0 , 2 z τ ( 1 β 1 + α 1 , α 1 , γ 1 ) , ( 1 β 2 + α 2 , α 2 , γ 2 ) ( 1 β 2 , α 1 , γ 2 ) , ( 1 β 1 , α 1 , γ 1 ) f ( τ ) d τ ,
and using again Rathie’s property (7.5) for I 2 , 2 0 , 2 I 2 , 2 2 , 0 , and symmetry of the I 2 , 2 2 , 0 -function with respect to the parameters in the upper and lower rows, this gives
I 2 1 I 1 1 f ( z ) = 0 z I 2 , 2 2 , 0 τ z ( β 1 , α 1 , γ 1 ) , ( β 2 , α 2 , γ 2 ) ( β 1 α 1 , α 1 , γ 1 ) , ( β 2 α 2 , α 2 , γ 2 ) f ( τ ) d τ , i . e . ,
I 2 1 I 1 1 f ( z ) = z 0 1 I 2 , 2 2 , 0 σ ( β 1 , α 1 , γ 1 ) , ( β 2 , α 2 , γ 2 ) ( β 1 α 1 , α 1 , γ 1 ) , ( β 2 α 2 , α 2 , γ 2 ) f ( z σ ) d σ = I 2 , 1 2 f ( z ) .
The statement for the composition of the arbitrary number m 1 of operators (66) can then be derived by mathematical induction, following the same procedure as above. In the final stage, we repeat the above evaluation for
I m 1 I m 1 f ( z ) = I m f ( z ) .
Note that due to symmetry of the I m , m m , 0 -function with respect to the parameters in the upper and lower rows, the composition is commutative, i.e., does not depend on the order of the operators I i 1 , i = 1 , 2 , , m in (66). □
Remark 2.
An alternative way to prove Theorem 7 is to show that the Mellin transform images M of the operators I m and of the composition i = 1 m I i 1 coincide for entire functions f ( z ) . This kind of proof is similar to that we used for the composition/decomposition property of the generalized fractional integrals (43) with H m , m m , 0 -kernels, see Th. 5.1.5 and Th. 5.2.1 in Kiryakova [20] (Ch.5). Using (67) again and the auxiliary formula (63), we have:
M I m f ( z ) ; s = i = 1 m Γ γ i ( β i s α i ) Γ γ i ( β i + α i s α i ) M f ( z ) ; s = i = 1 m M I i 1 f ( z ) ; s .
In a future study, we can look for a semigroup property for the new generalized operators of the fractional integration of multi-order ( α 1 , , α m ) of a form more general than that appearing in (61).

5. Illustration of Results from Section 3 and Section 4 for Particular Cases of the Le Roy-Type Functions

Here, we give examples as to how the discussed results work for the previously introduced and studied Le Roy-type functions F ( γ ) , F α , β ( τ ) , F ( α , β ) m ( γ ) m , and for some of their particular cases.
For m = 1 , the simplest case is the original Le Roy function (20), γ > 0 . According to Theorem 3, (28), etc., for α = β = 1 ,
F ( γ ) ( z ) = k = 0 z k ( k ! ) γ = k = 0 z k Γ γ ( 1 . k + 1 ) = 1 Ψ ˜ 1 ( 1 , 1 , 1 ) ( 1 , 1 , γ ) z
= H ¯ 1 , 2 1 , 1 z ( 0 , 1 , 1 ) ( 0 , 1 ) , ( 0 , 1 , γ ) = I 1 , 2 1 , 1 z ( 0 , 1 , 1 ) ( 0 , 1 , 1 ) , ( 0 , 1 , γ ) .
Naturally, for γ = 1 ,
F ( 1 ) ( z ) = k = 0 z k k ! = 1 Ψ 1 ( 1 , 1 ) ( 1 , 1 ) z = 1 F 1 ( 1 ; 1 ; z ) = exp ( z ) ,
as a “classical” generalized hypergeometric function, and in particular, the exponential one. Based on Theorem 4, the “eigen” integral operator generated by the Le Roy function for entire functions f ( z ) = k = 0 a k z k has the form
L 1 , 1 1 f ( z ) = z k = 0 a k Γ γ ( k + 1 ) Γ γ ( k + 2 ) z k = z 0 1 I 1 , 1 1 , 0 σ ( 1 , 1 , γ ) ( 0 , 1 , γ ) f ( z σ ) d σ ,
with L 1 , 1 1 F γ ( λ z ) = ( 1 / λ ) F γ ( λ z ) 1 / λ , and the corresponding G-L differentiation (57): D 1 , 1 1 F γ ( λ z ) = λ F γ ( λ z ) , λ 0 .
For the Mittag-Leffler-type Le Roy function (21) (Gerhold, Garra–Polito, Garrappa–Rogosin–Mainardi, etc.) with α > 0 , β > 0 , γ > 0 , the same results read as follows:
F α , β ( γ ) ( z ) = k = 0 z k Γ γ ( α k + β ) = 1 Ψ ˜ 1 ( 1 , 1 , 1 ) ( β , α , γ ) z
= H ¯ 1 , 2 1 , 1 z ( 0 , 1 , 1 ) ( 0 , 1 ) , ( 1 β , α , γ ) = I 1 , 2 1 , 1 z ( 0 , 1 , 1 ) ( 0 , 1 , 1 ) , ( 1 β , α , γ ) .
The corresponding eigenoperators L α , β γ , resp. B α , β γ satisfy the eigenfunction relations:
L α , β γ F α , β ( γ ) ( λ z ) = 1 λ F α , β ( γ ) ( λ z ) 1 λ Γ γ ( β ) , D α , β γ F α , β ( γ ) ( λ z ) = λ F α , β ( γ ) ( λ z ) ,
where we have the Le Roy-type analogue of the Erdélyi–Kober fractional integral:
L α , β γ f ( z ) = z k = 0 a k Γ γ ( α k + β ) Γ γ ( α k + α + β ) z k = z 0 1 I 1 , 1 1 , 0 σ ( β , α , γ ) ( β α , α , γ ) f ( z σ ) d σ .
For m = 2 , in our previous papers [13,14,16], we discussed the Le Roy-type function mentioned in Pogány [79], and kept the original denotations from there:
F ( p , q ; r , s ) α , β ( z ) = k = 0 z k Γ α ( p k + q ) Γ β ( r k + s ) .
Then, we provided the values for the order and type of this entire function. In view of the results of Section 3 and Section 4, the following new representations follow:
F ( p , q ; r , s ) α , β ( z ) = 2 Ψ ˜ 3 ( 1 , 1 , 1 ) ( 1 , 1 , 1 ) , ( q , p , α ) , ( s , r , β ) z
= H ¯ 2 , 4 1 , 2 ( z ) = I 2 , 4 1 , 2 z ( 0 , 1 , 1 ) ( 0 , 1 , 1 ) , ( 0 , 1 , 1 ) , ( 1 q , p , α ) , ( 1 s , r , β ) ,
and the corresponding Gelfond–Leontiev–Le Roy integration operator is represented by (with the adapted denotations)
L ( p , q ; r , s ) α , β f ( z ) = z 0 1 I 2 , 2 2 , 0 σ ( q , p ) , ( s , r ) ( q p , p ) , ( s r , r ) f ( z σ ) d σ .
The above kernel function can be considered as an analogue of the Gauss hypergeometric function related to G 2 , 2 2 , 0 , so such an integral operator can be thought as analogous to the popular Saigo fractional integrals. It is important to mention that a case of G-L-Le Roy-type integration with m = 3 and a kernel I 3 , 3 3 , 0 would resemble an extension of the so-called Marichev–Saigo–Maeda operators of FC, see more details in Kiryakova [80].
In the case of arbitrary m 1 , for the general M-L-P function of Le Roy type (24), the representations in terms of the generalized Fox–Wright function, H ¯ - and I-functions are given as in Theorem 3. However, the problem of defining corresponding Gelfond–Leontiev operators for which F m is an eigenfunction seems to still be open, in the general case.
In the case of Prabhakar parameters τ i = 1 , we have the following representations:
F ( α , β ) m ( γ ) m ( z ) = m Ψ ˜ 2 m 1 ( 1 , 1 , 1 ) 1 m ( 1 , 1 , 1 ) ( m 1 ) t i m e s , ( β i , α i , γ i ) 1 m z
= H ¯ m , 2 m 1 , m z ( 0 , 1 , 1 ) 1 m ( 0 , 1 ) m t i m e s , ( 1 β i , α i , γ i ) 1 m = I m , 2 m 1 , m z ( 0 , 1 , 1 ) 1 m ( 0 , 1 , 1 ) m t i m e s , ( 1 β i , α i , γ i ) 1 m ,
where some of the orders of the above Ψ ˜ , H ¯ - and I-functions can be suitably reduced because of the coincidence of some parameters ( 1 , 1 , 1 ) in the upper row with part of the same in the lower row. Then, these are analogous to the representations for the multi-index Mittag-Leffler functions (9) with γ i = 1 :
F ( α , β ) m ( 1 ) m : = E ( α i , β i ) ( m ) ( z ) = 1 Ψ m ( 1 , 1 ) ( β i , α i ) 1 m z = H 1 , m + 1 1 , 1 z ( 0 , 1 ) ( 0 , 1 ) , ( 1 β i , α i ) .
For the F ( α , β ) m ( γ ) m -functions, in Section 4 we developed the Gelfond–Leontiev generalized integrations and differentiation (58) and (57) as their eigenoperators, called Gelfond–Leontiev–Le Roy operators.
It is necessary to emphasize the evident analogy of the representations of (58) and (57), with the results for the operators (47): L ( α i , β i ) and D ( α i , β i ) generated by the multi-index Mittag-Leffler functions, in terms of the Generalized Fractional Calculus (GFC), see Kiryakova [20]. Compare the integral representation (48) with the new one of Le Roy-type, (61), when the kernel Fox H m , m m , 0 -function is replaced by a more general I m , m m , 0 -function. For m = 1 , we have L ( α , β ) in the form of an Erdélyi–Kober integral (40), and, respectively, the operator I 1 in (66) as a kind of “Erdélyi-Kober-Le Roy integral” (66).

6. Some Short Reminders from “Guide to SF of FC”

Here, we provide only a brief reminder of some of the SF of FC that are considered in detail in the survey of Kiryakova [9], and partly in the monographs by Kiryakova [20] and Paneva-Konovska [43]. The reason is because these special functions are to be compared with related extensions to the Le Roy type, I- and H ¯ -functions.

6.1. Classes of G- and H-Functions as Kernels of Laplace-Type Integral Transforms and of Operators of Generalized Fractional Calculus

In survey [9] (Sect.3), we attracted readers’ attention to the use of two basic classes of G- and H-functions with specific orders: (i)  G 0 , m m , 0 , resp., H 0 , m m , 0 with m = q , n = p = 0 ; and (ii)  G m , m m , 0 , resp., H m , m m , 0 with m = p = q , n = 0 .
Concerning the G- and H-functions of type (i), let us mention that the analogues of the Laplace transform, the so-called G- and H-transforms
G { f ( t ) ; s } = 0 G p , q m , n s t ( a j ) 1 p ( b k ) 1 q f ( t ) d t ,
and, respectively,
H { f ( t ) ; s } = 0 H p , q m , n s t ( a j , A j ) 1 p ( b k , B k ) 1 q f ( t ) d t ,
are said to be generalized integral transforms of Laplace type when
δ = m + n p + q 2 > 0 , resp . a * = j = 1 n A j j = n + 1 p A j + k = 1 m B k k = m + 1 q B k > 0 .
As their simplest examples, we mention the Laplace and Meijer transforms
L { f ( t ) ; s } = 0 exp ( s t ) f ( t ) d t , K ν { f ( t ) ; s } = 0 s t K ν ( s t ) f ( t ) d t ,
and also the Borel transform (considered by Dzrbasjan [75], see also Kiryakova [20] (Ch.2))
B ( ρ ) , ( μ ) f ( t ) ; s = 0 exp ( s ρ t ρ ) t μ ρ 1 f ( t ) d t , ρ > 0 , μ C .
In 1958, the Bulgarian mathematician Obrechkoff introduced a far-reaching generalization of both Laplace and Meijer transforms, particular cases of which were studied years later by many authors (Ditkin, Prudnikov, Krätzel, etc.). We mention, for example, the Krätzel transform ([81,82]),
L ν ( m ) f ( t ) ; s : = 0 Λ ( s , t ) f ( t ) d t , where Λ ( s , t ) is   the   K r ä t z e l   f u n c t i o n   as   below .
For details and studies on the latter function, see [81,82,83].
Note that the kernels of the above-mentioned integral transforms appear, resp., as
exp ( s ) = G 0 , 1 1 , 0 s 0 , K ν ( s ) = 1 2 G 0 , 2 2 , 0 s 2 4 ν 2 , ν 2 ,
exp ( s ρ t ρ ) = H 0 , 1 1 , 0 s t ( μ 1 ρ , 1 ρ ) , Λ ( s , t ) = s ν 1 + 2 m G 0 , m m , 0 s t 0 , ( ν + k 2 m ) k = 2 k = m .
These integral transforms have been used mainly for the purposes of operational calculi for different classes of differential operators, like
D = d d t , B ν = d d t t 1 ν d d t t ν , D 1 / ρ , B ν ( m ) = d d t t 1 m ν t 1 1 m d d t m 1 t ν + 1 2 m , etc .
The initial aims of Obrechkoff were to extend the theorem of S.N. Bernstein for absolutely monotonic functions representable by means of Laplace–Stieltjes transforms when the conditions for n-th derivatives are replaced by similar for differential operators, more general than those above. The Obrechkoff transform, as modified by Dimovski and studied in detail by Kiryakova (see [84], also [20,85] (Ch.3)), was defined originally in the form
O f ( t ) ; s = β 0 t β ( γ m + 1 ) 1 K ( s t ) β f ( t ) d t = β 0 λ ( t , s ) f ( t ) d t ,
with the kernel function
K ( s ) = 0 0 exp u 1 u m 1 s u 1 u m 1 k = 1 m u k γ m γ k 1 d u 1 d u m 1 .
Later, in Kiryakova [20] (Ch.3) (and in other works, see e.g., [84]), we proved that the kernel function (70) of the Obrechkoff transform is representable by a Meijer’s G 0 , m m , 0 -function, namely
λ ( t , s ) = s β ( γ m + 1 ) + 1 G m , m m , 0 ( s t ) β ( γ k 1 β + 1 ) 1 m .
Thus, the Obrechkoff transform appears to be a G-transform of Laplace type, related to the hyper-Bessel differential operators (Dimovski [86], Kiryakova [85], Dimovski–Kiryakova [84])
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 , α 1 , , α m , β : = m ( α 0 + α 1 + + α m ) > 0 , γ k : = ( α k + α k + 1 + + α m ) / β , k = 1 , , m , and P m a polynomial of degree m. Evidently, for particular choices of these parameters, (71) gives the above-mentioned simpler differential operators and many others frequently appearing in problems of mathematical physics.
Next, the generalized Obrechkoff transform (its fractionalized analogue) was introduced and studied by Kiryakova [20] (Ch.5) (and next in papers by Luchko, Al Mussalam and V.K. Tuan, and also in Yakubovich–Luchko [36]) with a Fox H m , m m , 0 -function as a kernel:
GO ( s ) : = B ( ρ i ) , ( μ i ) f ( t ) ; s = 0 H 0 , m m , 0 s t ( μ i 1 ρ i , 1 ρ i ) 1 m f ( t ) d t .
We called it multi-index Borel–Dzrbashjan transform because m = 1 is reduced to the mentioned Borel transform. The “fractional” Obrechkoff transform (72) can be considered not only as a H-transform but as a tool in operational calculus for the fractional multi-order analogues of hyper-Bessel differential operators (71), formally written as
D ( ρ i ) , ( μ i ) f ( t ) = t 1 i = 1 m t 1 + ( 1 μ i ) ρ i D t ρ i 1 / ρ i t ( μ i 1 ) ρ i f ( t ) ,
in the same way as the Laplace transform, Obrechkoff transform and its particular cases serve for the classical differentiation, for the hyper-Bessel operators (71), etc.
All the above details on the use of the functions G 0 , m m , 0 and H 0 , m m , 0 serve to pose an Open Problem (see 8.4. in the concluding Section 8) for the possibility to introduce and study a Laplace-type integral transform with I 0 , m m , 0 and H ¯ 0 , m m , 0 as kernel functions.
We also have to emphasize in terms of the other class of special functions the G m , m m , 0 and H m , m m , 0 -functions of type (ii), as discussed in survey [9], because here in Section 4, the classical FC operators of integration (such as R-L and E-K, (40)) are mentioned as having such kernels for m = 1 , and also the generalized fractional integrals (42) have these kernels for arbitrary order m 1 . In addition, in Section 4.4 (Theorem 6), their I m , m m , 0 -analogue appears as a kernel function of the integral operator (61), being an alternative representation of the series (58) for the introduced Gelfond–Leontiev generalized integration generated by the Le Roy-type functions (22). Thus, the need for the analogy is evident.

6.2. Some Important Cases of Mittag-Leffler Functions (7) and Multi-Index Mittag-Leffler Functions (9)

  • We mention a few cases of M-L functions ( m = 1 ) from §5.1 of [9]: Generalized trigonometric functions of higher integer orders m 1 , cos m and, resp., sine-functions; and next mentioned their “fractalized” analogues; also the Lorenzo–Hartley functions, the Rabotnov function, etc. Extensive literature is nowadays available on the theory of the M-L functions and their cases, with a few to mention such as: [2,30,31,39,50,87], etc.
  • For m = 2 , a not very popular function introduced by Dzrbashjan [29], in his and, respectively, in “our” denotations, is:
Φ ρ 1 , ρ 2 ( z ; μ 1 , μ 2 ) = k = 0 z k Γ ( μ 1 + k ρ 1 ) Γ ( μ 2 + k ρ 2 )
: = E ( 1 ρ 1 , 1 ρ 2 ) , ( μ 1 , μ 2 ) ( z ) = E ( α 1 , α 2 ) , ( β 1 , β 2 ) ( z ) , 1 / ρ i : = α i > 0 , μ i : = β i .
In addition to the classical M-L function, the geometric series and the Bessel functions J ν , we observed ([9]) that (73) also includes: the Struve and Lommel functions s μ , ν and H ν , and the (classical) Wright function [24,25] (details in [27] and in many books)
ϕ ( α , β ; z ) = k = 0 1 Γ ( α k + β ) z k k ! = 0 Ψ 1 ( β , α ) z = E ( α , 1 ) , ( β , 1 ) ( 2 ) ( z ) ,
and thus, also the Mainardi function  M ( z ; α / 2 ) = ϕ ( α / 2 , 1 α / 2 ; z ) , and the Airy function M ( z ; 1 / 3 ) = 3 2 / 3 A i ( z / 3 1 / 3 ) , etc.
In other denotations, the Wright function ϕ ( α , β ; z ) is also known as the generalized Bessel function, the Wright–Bessel function or misnamed as the Bessel–Maitland function:
J ν μ ( z ) = ϕ ( μ , ν + 1 ; z ) = 0 Ψ 1 ( ν + 1 , μ ) z
= k = 0 ( z ) k Γ ( ν + k μ + 1 ) k ! = E ( 1 / μ , 1 ) , ( ν + 1 , 1 ) ( 2 ) ( z ) ,
again as an example of the Dzrbashjan function. We discuss its extension as an H ¯ -function (81) in Section 7. In the same survey [9] (§5.2), further extensions of the Bessel function with 3, 4, etc., parameters are mentioned.
  • For the case m 2 , we mention the hyper-Bessel functions of Delerue [88] (for details and essential use of these functions, see [20] (Ch.3)):
J γ i , , γ m 1 ( m 1 ) ( z ) = z m i = 1 m 1 γ i E ( 1 , 1 , , 1 ) , ( γ 1 + 1 , γ 2 + 1 , , γ m 1 + 1 , 1 ) ( m ) ( z m ) m
= i = 1 m 1 Γ ( γ i + 1 ) 1 z m i = 1 m 1 γ i 0 F m 1 γ 1 + 1 , γ 2 + 1 , , γ m 1 + 1 ; ( z m ) m .
This representation suggests that the multi-index M-L functions (9) with arbitrary ( α 1 , , α m ) ( 1 , , 1 ) can be interpreted as fractional-indices analogues of the hyper-Bessel functions (75), which themselves are multi-index analogues of the Bessel function. Functions (75) are closely related to the hyper-Bessel differential operators (71) of Dimovski [86] and form a fundamental system of solutions of the differential equations of the form B y ( z ) = λ y ( z ) , as proved in Kiryakova [20] (Ch.3, Th.3.4.3).
In relation to the multi-index analogues of the Mittag-Leffler functions, we can also mention two other variants that appear as solutions of some fractional multi-order (multi-term) differential equations. In Gorenflo–Kilbas–Rogosin [89], the function
E α , μ , l ( z ) = k = 0 c k z k , w h e r e c k = j = 0 k 1 Γ [ α ( j μ + l ) + 1 ] Γ [ α ( j μ + l + 1 ) + 1 ] ,
is introduced and studied, with indices R e ( α ) > 0 , l C , μ R . One can be curious about a case when Le Roy-type fractional power indices are attached to the Gamma functions in the coefficient c k . Another variant is the function introduced and studied recently by Droghei [90]
W α ¯ , ν ¯ , n ( z ) = k = 0 i = 1 k j = 1 n Γ ( α n + 1 i + a j ) Γ ( α n + 1 i + b j ) z k Γ ( α n + 1 k + b n + 1 ) ,
with parameters α ¯ = α 1 , , α n + 1 ; ν ¯ = ν 1 , , ν n and a j = 1 + m = 1 j ν m 1 α m ; and b j = 1 + m = 1 j ν m 1 α m 1 , such that a j = b j α j with j = 1 , , n + 1 . In some particular cases, say for n = 1 , its three-parameter variant:
W α , β , ν ( z β ) = k = 0 i = 1 k Γ ( β i + 1 α ) Γ ( β i + 1 ) z β k Γ ( β k + 1 α + ν )
is shown to be a solution of the fractional order variant of the Bessel equation, and reducible to the Bessel function; to the Wright function ( α = 1 ); to the Mittag-Leffler function ( α = 0 : W 0 , α , β 1 ( z ) = E α , β ( z ) ); and to some cases of the hyper-Bessel and 2 × 2 Mittag-Leffler function of Dzrbashjan: W ν , ν , ν ( z ν ) = k = 0 z ν k Γ 2 ( ν k + 1 ) . Therefore, this can also be considered as a Le Roy-type function but with integer index γ = 2 .
Other cases of the p Ψ q and of multi-index M-L functions are mentioned in [9] (Ch.6), such as Virchenko and Ricci generalized hypergeometric functions, Mainardi–Massina and Paris generalized exponential integrals, generalized K-series and M-series, etc., as well as critics on the numerous so-called k-analogues defined by replacement of the Γ -functions by Γ k -functions that bring no essential novelties at all (details in Kiryakova [91]).
A multi-variable Mittag-Leffler multi-index function (called also multinomial M-L function) was introduced by Luchko et al. (e.g., [36,60]) as well. Further, a multinomial Prabhakar function was introduced by Bazhlekova–Bazhlekov, see Bazhlekova [92] (2.2), where an additional Prabhakar parameter δ is involved. The complete monotonicity of this new special function is studied via the Laplace transform image, and the appearance of the above-mentioned multinomial M-L type functions in resolving multi-term time-fractional evolution equations is discussed.

7. Other Important SF in the Scheme of the I - and H ¯ -Functions

Next, it is time to mention some other well-known or not so well-known special functions, but all with important significance, that also fall in the scheme of the H ¯ -functions of Inayat-Hussain [11] and of the more general I-functions of Rathie [10]. This class of generalized hypergeometric functions is introduced and studied as an extension of the Fox–Wright functions p Ψ q and of the Fox H-functions, and the functions are defined by kinds of “Mellin-Barnes”-type integrals (or more exactly, their analogues). They involve quotients of Gamma functions on arbitrary (fractional) powers, and in general, present multi-valued functions with branch points, thus needing precision of the meanings and conditions of existence. Several authors have studied the details of the I- and H ¯ -functions, along with the above-mentioned initiators, such as: Buschman–Srivastava [12], Saxena [47], Gupta–Soni [93], Srivastava–Lyn–Wang [70], Mathai–Saxena–Haubold [6], Vellaisamy–Kataria [78], Jolly [94], Srivastava [95], etc. Nevertheless, this theory has not yet been deeply developed and this has raised some open problems. However, we need to mention that the I- and H ¯ -functions arose from the needs for mathematical tools, such as for the evaluation of certain Feynman integrals (arising in perturbation calculations of the equilibrium properties of a magnetic model of phase transitions), partition function of the Gauss model from statistical mechanics, etc.
We skip the particular examples of I- and H ¯ -functions from the initial works of Inayat-Hussain and Rathie, invoking their appearance; these are also mentioned in Appendix A.5.1 in Mathai–Saxena–Haubold [6], as (A.60)–(A.62).
Here, we provide some examples of popular “mathematical” functions, starting from the simplest ones.
  • The polylogarithm function.
Consider (see for example, [1] (Vol.1), [5], etc.)
L i α ( z ) = k = 1 z k k α = z + z 2 2 α + z 3 3 α + , | z | < 1 , α C ,
which for particular choices of the index α appears as: ordinary logarithm L i 1 ( z ) = l n ( 1 z ) ; L i 0 ( z ) = z 1 z ; L i 1 ( z ) = z ( 1 z ) 2 ; L i n ( z ) = z n + 1 F n ( 1 , 1 , ; 2 , 2 , ; z ) , n = 0 , 1 , 2 , ; L i n ( z ) = z n F n 1 ( 2 , 2 , ; 1 , 1 , ; z ) , n = 1 , 2 , ; etc. Note also that for (the singular value) z = 1 when α > 1 , it gives the famous Riemann Zeta function, L i α ( 1 ) = ζ ( α ) = k = 1 1 / ( k α ) . The function (76) is also referred to as the Jonquière function. As particular examples, for special choice of argument and parameters, one can also mention some popular polynomials such as the Jonquière and Bernoulli polynomials, and number sequences such as the Stirling and Eulerian numbers, etc.
For the purposes here, we consider the real positive index  α > 0 . Contour-type integral representations of “Mellin-Barnes”-type have been considered in many of the above-mentioned works, like, e.g., Jolly [94] (1.1.24). We appeal to the representation as presented by Gerhold–Tomovski [96], (1.2)–(1.3):
L i α ( z ) ( z ) = 1 2 π i 1 / 2 i 1 / 2 + i π s α sin π s ( z ) s d s ,
and using [1] (Vol.1, Ch.1),
π sin π s = Γ ( s ) Γ ( 1 s ) ; s α = Γ ( 1 + s ) Γ ( s ) α , we   have :
L i α ( z ) = 1 2 π i 1 / 2 i 1 / 2 + i Γ α + 1 ( s ) Γ ( 1 s ) Γ α ( 1 + s ) ( z ) s d s
= H ¯ 1 , 2 1 , 1 z ( 1 , 1 , α + 1 ) ( 1 , 1 , 1 ) , ( 0 , 1 , α ) ,
as a H ¯ 1 , 2 1 , 1 -function of the same orders as the H ¯ -function (29) for our Mittag-Leffler–Prabhakar function of Le Roy type, F 1 for m = 1 (but with different parameters). Note that the singularities of the Gamma functions in the numerator of the above contour integral for L i α are all lying to the left of s = 0 and to the right of s = 1 , and so there are none in the interval c ( 0 , 1 ) , which explains the choice c = 1 / 2 for the “vertical” contour L = ( c i , c + i ) . Meanwhile, the series (76) can be rewritten as
L i α ( z ) = k = 1 Γ α ( k ) Γ α ( k + 1 ) z k ,
that is, in the form of an “extended” Le Roy-type function, as discussed in Section 4 of our joint work [14].
  • The generalized Riemann Zeta function (Hurwitz–Lerch Zeta function) ([1] eq.(1), §1.11).
As mentioned in [94] (1.1.27), see also [97] (1.1),
Φ ( z , α , b ) = n = 0 z n ( n + b ) α
= H ¯ 2 , 2 1 , 2 z ( 0 , 1 , 1 ) , ( 1 b , 1 , α ) ( 0 , 1 , 1 ) , ( b , 1 , α ) , b 0 , 1 , 2 , , | z | 1 ,
and more precisely, for α C when | z | < 1 , or R e ( α ) > 1 when | z | = 1 ; with a branch cut along with the positive semi-axis from 1 to + . Then, the function is analytic in z in the so-cut z-plane.
This contains not only the Riemann Zeta function ζ ( α ) ( z = 1 , b = 0 ) and the Hurwitz Zeta function ζ ( α , b ) = Φ ( 1 , α , b ) with R e ( α ) > 1 but also the above-mentioned polylogarithm function (76), as L i α ( z ) = z Φ ( z , α , 1 ) , α C when | z | < 1 , or R e ( α ) > 1 when | z | = 1 .
The more general Hurwitz–Lerch Zeta function is:
Φ μ , ν ρ , σ ( z , α , b ) = n = 0 ( μ ) ρ n ( ν ) σ n · z n ( n + b ) α , for | z | < δ = ρ ρ σ σ ,
involving Pochhamer symbols, and in some cases as:
Φ ν , ν σ , σ ( z , α , b ) = Φ ( z , α , b ) , and Φ μ , 1 1 , 1 ( z , α , b ) = n = 0 ( μ ) n n ! · z n ( n + b ) α ·
Details of the parameters and contoursfor the more general case are explained below.
  • The generalized Hurwitz–Lerch Zeta function (Srivastava–Saxena–Pogány–Saxena [97]).
It is assumed that λ , μ C ; b , ν = 0 , 1 , 2 , ; real ρ , σ , κ > 0 ; κ ρ σ > 1 when α , z C , or κ ρ σ = 1 and α C when | z | < δ * = ρ ρ σ σ κ κ , while κ ρ σ = 1 and Re ( α + ν λ μ ) > 1 when | z | = δ * . Then, the following generalized (extended) Hurwitz–Lerch Zeta functions is considered:
Φ λ , ν , μ ( ρ , σ , κ ) ( z , α , b ) = n = 0 ( λ ) ρ n ( μ ) σ n ( ν ) κ n n ! · z n ( n + b ) α ·
In Srivastava–Saxena–Pogány–Saxena [97], the Riemann–Liouville fractional derivatives of the above Hurwitz–Lerch Zeta-type functions are evaluated, as expected in terms of similar kind of functions, in full agreement with the results as in Kiryakova [69,80], and the ideas for the SF of FC ([9]). Functional relations are also presented of (78) with the generalized Wright–Fox functions p Ψ ˜ q of the form (16) from Section 2.
According to [97] (Th.4), the generalized Hurwitz–Lerch Zeta function (78) has the following contour integral representation:
Φ λ , ν , μ ( ρ , σ , κ ) ( z , α , b ) = Γ ( ν ) Γ ( λ ) Γ ( μ )
× L Γ ( s ) Γ ( λ + ρ s ) Γ ( μ + σ s ) Γ α ( s + b ) Γ ( ν + κ s ) Γ α ( s + b + 1 ) ( z ) s d s ,
for | a r g ( z ) | < π , and path of integration L = ( c i , c + i ) in the complex s-plane that separates the singularities of Γ ( s ) , Γ ( λ + ρ s ) , Γ ( μ + σ s ) and Γ α ( s + b ) . Then, the relation with the H ¯ -function is obtained, as in [97] (Th.5):
Φ λ , ν , μ ( ρ , σ , κ ) ( z , α , b ) = Γ ( ν ) Γ ( λ ) Γ ( μ )
× H ¯ 3 , 3 1 , 3 z ( 1 λ , ρ , 1 ) , ( 1 μ , σ , 1 ) , ( 1 b , 1 , α ) ( 0 , 1 ) , ( 1 ν , κ , 1 ) , ( b , 1 , α ) .
Several important special cases are then considered, such as asymptotic estimations, etc.
  • Generalized Wright–Bessel function (see Jolly [94], (1.1.27), p.13)
J ¯ λ ν , μ ( z ) = k = 0 ( z ) k Γ μ ( λ + k ν + 1 ) k ! = H ¯ 0 , 2 1 , 0 z ( 0 , 1 ) , ( λ , ν , μ ) .
Let us compare this with the popular Wright–Bessel/Bessel–Maitland function (e.g., Marichev [5]) considered often by Kiryakova at al., for example, (57) in [9], see in Section 6:
J λ ν ( z ) = 0 Ψ 1 ( λ + 1 , ν ) z = k = 0 ( z ) k Γ ( λ + k ν + 1 ) k ! = H 0 , 2 1 , 0 z ( 0 , 1 ) , ( λ , ν ) ,
which appears as a case of (81) if μ = 1 . See also Srivastava et al. [70] (3.2).
One can prolong this list with some more complicated or exotic-looking special functions that can fall in the scheme of the I-functions. We limit ourselves to the case below.
  • In Bhatter et al. [98], the authors introduce the so-called E-function:
τ E k h z ( ρ , a ) ; ( γ i , q i , α i ) 1 h ( α , β ) ; ( δ j , p j , r j ) 1 k = n = 0 c n ( 1 ) ρ n z a n + τ Γ ( α n + β )
with coefficients including the Pochhamer symbols c n = ( γ 1 ) q 1 n α 1 ( γ h ) q h n α h ( δ 1 ) p 1 n r 1 ( δ k ) p k n r k .
These lead to Γ -functions to (arbitrary) fractional powers both in the numerator and denominator of the power series (82), which is convergent for all finite values of
i = 1 h ( q i ) q i α i α α j = 1 k ( p j ) p j r j 1 | z a | .
We omit all assumptions on parameters and on the contour that are provided in the cited work. As shown in [98] (Th.5), this E-function allows a kind of Mellin–Barnes-type integral
τ E k h ( z ) = v = 1 k Γ r v ( δ v ) u = 1 h Γ α u ( γ u ) · z τ 2 π i L Γ ( s ) Γ ( 1 s ) i = 1 h Γ α i ( γ i q i s ) Γ ( β α s ) j = 1 k Γ r j ( δ j p j s ) ( 1 ) ρ ( z a ) s d s .
Then, it can be considered as a case of a H ¯ 2 , k 1 , 1 + h -function, and also has a form close to the “extended” Le Roy-type functions, discussed in Sect. 4 of our joint work [14].

8. Some Remarks and Open Problems

We now mention in brief some open or related problems that can be further discussed for the Le Roy-type functions (in the sense as discussed here), and more generally, for the I-functions.
8.1. 
Open problem to determine G-L-type operators of generalized integration (33) in the case of the Prabhakar-type Le Roy functions (24), and the corresponding “eigenoperators” as the generalized differentiation (32), for which these special functions can appear as eigenfunctions. As mentioned in Kiryakova [64], even for the simplest case m = 1 and γ = 1 , the problem stays open if τ 1 . In the particular case of Prabhakar parameters τ i = 1 , we proposed such operators in Section 4 here, see (58) as well as (61).
8.2. 
Open problem to find conditions of parameters for which the Le Roy-type functions (24), or the simpler (22), or the multi-index Mittag-Leffler functions (9) and (10), are Completely Monotone (CM). Over many years, several attempts have been made to study the CM of some simpler classes of special functions. As a reminder, a function f ( z ) is called CM if it is infinitely differentiable and ( 1 ) n f ( n ) ( z ) 0 , n = 0 , 1 , 2 , . According to the Bernstein theorem ([99]), a function f is CM if and only if it can be uniquely represented as a Laplace transform of a non-negative (weight) function.
The results for the Mittag-Leffler function E α , β : 0 α 1 , β α are known, see, e.g., [2] (after the earlier results by Pollard for β = 1 , and by Miller–Samko), for the classical Prabhakar function by Giusti at al. [34] and for the multinomial Prabhakar function by Bazhlekova [92]. See also Karp–Prilepkina [77], Berg–Çetinkaya–Karp [100] (Sect.3), etc., for some conditions related to the CM of some Fox H- and similar functions.
For the Le Roy-type function F ( α , β ) ( γ ) , there are results given by Gorska and Horzela [56] and Simon [57].
8.3. 
Open problem for the behavior of the I- and H ¯ -functions near the third singular points on the circle of convergence when these functions are analytic inside/outside disks with final radius R. Even for the Fox H-functions, in general, it is still an open problem, but some particular results are available in Karp [21] (12.24).
8.4. 
Open problem (as mentioned in Section 5) to define and study the properties of integral transforms of Laplace type with kernels I 0 , m m , 0 . These should be analogues of the Laplace; Borel-Dzrbashjan; Meijer; Obrechkoff; and “fractional” Obrechkoff transforms, where the kernel functions are, resp.: exp ( z ) ; z ρ 1 exp ( z ρ ) = G 0 , 1 1 , 0 ( z ρ ) ; K ν ( z ) = 1 2 G 0 , 2 2 , 0 z 2 4 ν 2 ν 2 ; G 0 , m m , 0 ( z ) ; H 0 , m m , 0 ( z ) , etc. See some details in Kiryakova [9] (Sect.3), and our works such as [20] (Ch.2, Ch.3, Ch.5); [37,38,84].
8.5. 
Open problem related to the notion of the non-holomonicity for the case of I-functions and Le Roy-type functions. A sequence { u ( n ) } is called holonomic (p-recursive) if it satisfies a homogeneous linear recurrence p 0 ( n ) u ( n ) + p 1 ( n ) u ( n + 1 ) + + p d ( n ) u ( n + d ) = 0 , n 0 with polynomials p k , and p d not identically zero. A formal power series f ( z ) = 0 u ( n ) z n is holonomic (d-finite) if it satisfies a homogenous linear ordinary differential equation p 0 ( z ) f ( z ) + p 1 ( z ) f ( z ) + + p d ( z ) f ( d ) ( z ) = 0 with polynomial coefficients. It is well known that such a power series is holonomic if and only if its coefficient sequence is.
For definitions of the above-mentioned notions, and non-holomonicity of the hypergeometric sequences  Γ α 1 ( n u 1 ) Γ α m ( n u m ) with non-integer α i , related to functions f ( z ) = Γ α 1 ( z u 1 ) Γ α m ( z u m ) , see Gerhold [101], Bell–Gerhold et al. [102] (Theorem 8), Flajolet–Gerhold–Salvy [103], Gerhold [52] (Sect.3), etc. Important examples of non-holonomic functions that are I-functions have been mentioned, such as the polylogarithms, Riemann Zeta functions, etc. More generally, I- and H ¯ -functions that involve such non-holonomic sequences are also expected to fall in the case of the non-holonomic functions, and in particular, the considered Le Roy-type functions.

Author Contributions

The ideas and results in this paper resulted from the authors’ recent joint research on the Le Roy-type functions. In an earlier survey paper in this journal [9], we confined the class of the special functions of fractional calculus only to the level of the Fox–Wright generalized hypergeometric functions, Fox H-functions, and especially, the Mittag-Leffler-type functions. Now, inspired by some works from the end of the 20th century by A.A. Inayat-Hussain, R.G. Buschman-H.M. Srivastava and A.K. Rathie, and by more recent ones, we observed that the Le Roy type functions as well as other important special functions, can be represented by the more general I-functions and H ¯ -functions. Thus, it sounds suitable to extend the class of functions related to fractional calculus. Since the theory of the latter-mentioned functions still puts forward several open problems, we title our work as a “Discussion Survey”. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no financial funding.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge the fruitful collaboration with Sergei Rogosin and Maryna Dubatovskaya for the joint research on the considered multi-index Mittag-Leffler–Prabhakar functions of Le Roy-type, as published in the very recent articles [14,16]. We are thankful also to Dmitrii Karp for the useful discussions on some subtle details in the theory of the Mellin–Barnes-type integral representations of the involved multi-valued functions. Thanks are due to the reviewers and editors who carefully checked the original version and made useful suggestions that improved this work.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Erdélyi, A.; Magnus, W.; Oberhettinger, F.T. (Eds.) Higher Transcendental Functions; McGraw Hill: New York, NY, USA, 1953–1955; Volume 1–3. [Google Scholar]
  2. Gorenflo, R.; Kilbas, A.; Mainardi, F.; Rogosin, S. Mittag-Leffler Functions, Related Topics and Applications, 2nd ed.; Springer: Berlin/Heidelberg, Germany, 2020. [Google Scholar] [CrossRef]
  3. Kilbas, A.A.; Saigo, M. H-Transforms: Theory and Applications; Series on Analytic Methods and Special Functions, 9; CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar]
  4. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  5. Marichev, O.I. Handbook of Integral Transforms of Higher Transcendental Functions, Theory and Algorithmic Tables; Ellis Horwood: Chichester, UK, 1983. Transl. from Russian Ed., Method of Evaluation of Integrals of Special Functions; Nauka i Teknika: Minsk, Russia, 1978. (In Russian) [Google Scholar]
  6. Mathai, A.M.; Saxena, R.K.; Haubold, H.J. The H-Function. Theory and Applications; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  7. Prudnikov, A.P.; Brychkov, Y.; Marichev, O.I. Integrals and Series, Vol. 3: More Special Functions; Gordon and Breach Science Publishers: New York, NY, USA; London, UK; Paris, France; Tokyo, Japan, 1992. [Google Scholar]
  8. Srivastava, H.M.; Gupta, K.S.; Goyal, S.P. The H-Functions of One and Two Variables with Applications; South Asian Publications: New Delhi, India, 1982. [Google Scholar]
  9. Kiryakova, V. A guide to special functions in fractional calculus. Mathematics 2021, 9, 106. [Google Scholar] [CrossRef]
  10. Rathie, A. A new generalization of the generalized hypergeometric functions. Le Matematiche 1997, LII, 297–310. [Google Scholar]
  11. Inayat-Hussain, A.A. New properties of hypergeometric series derivable from Feynman integrals: II. A generalization of the H-function. J. Phys. A Math. Gen. 1987, 20, 4119–4128. [Google Scholar]
  12. Buschman, R.G.; Srivastava, H.M. The H¯ functions associated with a certain class of Feynman integrals. J. Phys. A Math. Gen. 1990, 23, 4707–4710. [Google Scholar] [CrossRef]
  13. Kiryakova, V.; Paneva-Konovska, J. Multi-index Le Roy functions of Mittag-Leffler-Prabhakar type. Int. J. Appl. Math. 2022, 35, 743–766. [Google Scholar] [CrossRef]
  14. Kiryakova, V.; Paneva-Konovska, J.; Rogosin, S.; Dubatovskaya, M. Erdélyi-Kober fractional integrals (Part 2) of the multi-index Mittag-Leffler-Prabhakar functions of Le Roy type. Int. J. Appl. Math. 2023, 36, 605–623. [Google Scholar] [CrossRef]
  15. Paneva-Konovska, J. Prabhakar functions of Le Roy type: Inequalities and asymptotic formulae. Mathematics 2023, 11, 3768. [Google Scholar] [CrossRef]
  16. Paneva-Konovska, J.; Kiryakova, V.; Rogosin, S.; Dubatovskaya, M. Laplace transform (Part 1) of the multi-index Mittag-Leffler-Prabhakar functions of Le Roy type. Int. J. Appl. Math. 2023, 36, 455–474. [Google Scholar] [CrossRef]
  17. Pincherle, S. Sulle funzioni ipergeometriche generalizzate. Atti R. Accad. Lincei Rend. Cl. Sci. Fis. Mat. Natur. 1888, 4, 694–700, 792–799, (Reprinted in Salvatore Pincherle: Opere Scelte, UMI (Unione Matematica Italiana) Cremonese: Roma, Italy, 1954; Volume 1, pp. 223–239). [Google Scholar]
  18. Mainardi, F.; Pagnini, G. Salvatore Pincherle: The pioneer of the Mellin-Barnes integrals. J. Comput. Appl. Math. 2003, 153, 331–341. [Google Scholar] [CrossRef]
  19. Fox, C. The G and H-functions as symmetric Fourier kernels. Trans. Am. Math. Soc. 1961, 98, 395–429. [Google Scholar]
  20. Kiryakova, V. Generalized Fractional Calculus and Applications; Longman: Harlow, UK; J. Wiley: New York, NY, USA, 1994. [Google Scholar]
  21. Karp, D. Chapter 12—A note on Fox’s H-function in the light of Braaksma’s results. In Special Functions and Analysis of Differential Equations; Agarwal, P., Agarwal, R.P., Ruzhansky, M., Eds.; Chapman and Hall/CRC: New York, NY, USA, 2020; 12p, Available online: http://arxiv.org/abs/1904.10651v1 (accessed on 16 December 2023).
  22. Braaksma, B.L.J. Asymptotic expansions and analytic continuation for a class of Barnes integrals. Compos. Math. 1962–1964, 15, 239–341. [Google Scholar]
  23. Meijer, C.S. On the G-function. Indag. Math. 1946, 8, 124–134, 213–225, 312–324, 391–400, 468–475, 595–602, 661–670, 713–723. [Google Scholar]
  24. Wright, E.M. On the coefficients of power series having exponential singularities. J. Lond. Math. Soc. 1933, 8, 71–79. [Google Scholar] [CrossRef]
  25. Wright, E.M. The generalized Bessel function of order greater than one. Quart. J. Math. Oxf. Ser. 1940, 11, 36–48. [Google Scholar] [CrossRef]
  26. Fox, C. The asymptotic expansion of generalized hypergeometric functons. Proc. Lond. Math. Soc. Ser. 2 1928, 27, 389–400. [Google Scholar] [CrossRef]
  27. Gorenflo, R.; Luchko, Y.; Mainardi, F. Analytical properties and applications of the Wright function. Fract. Calc. Appl. Anal. 1999, 2, 383–414. [Google Scholar]
  28. Mittag-Leffler, G.M. Sur la nouvelle fonction Eα(x). C. R. de l’Acad. Sci. 1903, 137, 554–558. [Google Scholar]
  29. 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]
  30. Haubold, H.J.; Mathai, A.M.; Saxena, R.K. Mittag-Leffler functions and their applications. J. Appl. Math. 2011, 2011, 298628. [Google Scholar] [CrossRef]
  31. Rogosin, S. The role of the Mittag-Leffler function in fractional modeling. Mathematics 2015, 3, 368–381. [Google Scholar] [CrossRef]
  32. Prabhakar, T.R. A singular integral equation with a generalized Mittag-Leffler function in the kernel. Yokohama Math. J. 1971, 19, 7–15. [Google Scholar]
  33. Garra, R.; Garrappa, R. The Prabhakar or three parameter Mittag–Leffler function: Theory and application. Commun. Nonlinear Sci. Numer. Simul. 2018, 56, 314–319. [Google Scholar] [CrossRef]
  34. Giusti, A.; Colombaro, I.; Garra, R.; Garrappa, R.; Polito, F.; Popolizio, M.; Mainardi, F. A practical guide to Prabhakar fractional calculus. Fract. Calc. Appl. Anal. 2020, 23, 88–111. [Google Scholar] [CrossRef]
  35. Luchko, Y.F.; Srivastava, H.M. The exact solution of certain differential equations of fractional order by using operational calculus. Comput. Math. Appl. 1995, 29, 73–85. [Google Scholar] [CrossRef]
  36. Yakubovich, S.; Luchko, Y. The Hypergeometric Approach to Integral Transforms and Convolutions; Series Mathematics and Its Applications 287; Kluwer Academic Publishers: Dordrecht, The Netherlands; Boston, MA, USA; London, UK, 1994. [Google Scholar]
  37. Kiryakova, V. Multiple Dzrbashjan-Gelfond-Leontiev fractional differintegrals. In Recent Advances in Applied Mathematics’96 (Proceedings of International Workshop, Kuwait University), 1996; pp. 281–294. Available online: https://www.researchgate.net/publication/307122608_Multiple_Dzrbashjan-Gelfond-Leontiev_Fractional_Differintegrals_1 (accessed on 16 December 2023).
  38. Kiryakova, V. Multiindex Mittag-Leffler functions, related Gelfond-Leontiev operators and Laplace type integral transforms. Fract. Calc. Appl. Anal. 1999, 2, 445–462. [Google Scholar]
  39. Kiryakova, V. Multiple (multiindex) Mittag-Leffler functions and relations to generalized fractional calculus. J. Comput. Appl. Math. 2000, 118, 241–259. [Google Scholar] [CrossRef]
  40. 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]
  41. 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]
  42. Paneva-Konovska, J. Multi-index (3m-parametric) Mittag-Leffler functions and fractional calculus. C. R. Acad. Bulg. Sci. 2011, 64, 1089–1098. [Google Scholar]
  43. 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]
  44. 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]
  45. Rogosin, S.; Dubatovskaya, M. Multi-parametric Le Roy function revisited. Fract. Calc. Appl. Anal. 2023. Published Online First. [Google Scholar] [CrossRef]
  46. Slater, L.J. Generalized Hypergeomtric Functions; Cambridge University Press: London, UK; New York, NY, USA, 1966. [Google Scholar]
  47. Saxena, R.K. Functional relations involving generalized H-function. Le Matematiche 1998, LIII, 123–131. [Google Scholar]
  48. Le Roy, É. Sur les séries divergentes et les fonctions définies par un développement de Taylor. Ann. De La Fac. Des Sci. De Touluse 2e Sér. 1900, 2, 385–430. (In French) [Google Scholar]
  49. Le Roy, É. Valéurs asymptotiques de certaines séries procédant suivant les puissances entères et positives d’une variable réelle. Bulletin des Sci. Mathématiques, 2eme sér. 1900, 24, 245–268. (In French) [Google Scholar]
  50. Mainardi, F. Why the Mittag-Leffler function can be considered the Queen function of the fractional calculus? Entropy 2020, 22, 1359. [Google Scholar] [CrossRef]
  51. Kolokoltsov, V. The law of large numbers for quantum stochastic filtering and control of many particle systems. Theor. Math. Phys. 2021, 208, 937–957. [Google Scholar] [CrossRef]
  52. Gerhold, S. Asymptotics for a variant of the Mittag-Leffler function. Integral Transform. Spec. Funct. 2012, 23, 397–403. [Google Scholar] [CrossRef]
  53. Garra, R.; Polito, F. On some operators involving Hadamard derivatives. Integral Transform. Spec. Funct. 2013, 24, 773–782. [Google Scholar] [CrossRef]
  54. Garrappa, R.; Rogosin, S.; Mainardi, F. On a generalized three-parameter Wright function of le Roy type. Fract. Calc. Appl. Anal. 2017, 206, 1196–1215. [Google Scholar] [CrossRef]
  55. Garra, R.; Orsingher, E.; Polito, F. A note on Hadamard fractional differential equations with varying coefficients and their applications in probability. Mathematics 2018, 6, 4. [Google Scholar] [CrossRef]
  56. Gorska, K.; Horzela, A.; Garrappa, R. Some results on the complete monotonicity of Mittag-Leffler functions of Le Roy type. Fract. Calc. Appl. Anal. 2010, 22, 1284–1306. [Google Scholar] [CrossRef]
  57. Simon, T. Remark on a Mittag-Leffler function of Le Roy type. Integral Transform. Spec. Funct. 2022, 33, 108–114. [Google Scholar] [CrossRef]
  58. Mehrez, K.; Das, S. On some geometric properties of the Le Roy-type Mittag-Leffler functions. Hacet. J. Math. Stat. 2022, 51, 1085–1103. [Google Scholar] [CrossRef]
  59. Mehrez, K. Study of the analytic function related to the Le-Roy-type Mittag-Leffler function. Ukr. Math. J. 2023, 75, 719–743. [Google Scholar] [CrossRef]
  60. Luchko, Y. Operational method in fractonal calculus. Fract. Calc. Appl. Anal. 1999, 2, 463–488. [Google Scholar]
  61. Rogosin, S.; Dubatovskaya, M. Multi-parametric Le Roy function. Fract. Calc. Appl. Anal. 2023, 26, 54–69. [Google Scholar] [CrossRef]
  62. Tomovski, Ž.; Mehrez, K. Some families of generalized Mathieu–type power series, associated probability distributions and related inequalities involving complete monotonicity and log–convexity. Math. Inequal. Appl. 2017, 20, 973–986. [Google Scholar] [CrossRef]
  63. Gelfond, A.O.; Leontiev, A.F. On a generalization of the Fourier series. Mat. Sbornik 1951, 29, 477–500. (In Russian) [Google Scholar]
  64. 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]
  65. Samko, S.; Kilbas, A.; Marichev, O. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach: Yverdon, Switzerland, 1993. [Google Scholar]
  66. Podlubny, I. Fractional Differential Equations; Academic Press: Boston, MA, USA, 1999. [Google Scholar]
  67. 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; Chapter 1; pp. 1–21. [Google Scholar] [CrossRef]
  68. Sneddon, I.N. The use in mathematical analysis of Erdélyi-Kober operators and some of their applications. In Fractional Calculus and Its Applications, Proceedings of the International Conference, New Haven, CT, USA, June 1974; Ross, B., Ed.; Lecture Notes in Mathematics; Springer: New York, NY, USA, 1975; Volume 457, pp. 37–79. [Google Scholar]
  69. Kiryakova, V. Unified approach to fractional calculus images of special functions—A survey. Mathematics 2020, 8, 2260. [Google Scholar] [CrossRef]
  70. Srivastava, H.M.; Lyn, S.-D.; Wang, P.-Y. Some fractional-calculus results for the H¯-function associated with a class of Feynman integrals. Russ. J. Math. Phys. 2006, 13, 94–100. [Google Scholar]
  71. Kalla, S.L. Operators of Fractional Integration; Lecture Notes in Mathematics; Springer: Berlin/Heidelberg, Germany, 1980; Volume 798, pp. 258–280. [Google Scholar]
  72. Kiryakova, V. A brief story about the operators of the generalized fractional calculus. Fract. Calc. Appl. Anal. 2008, 11, 203–220. [Google Scholar]
  73. 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; Chapter 4; pp. 87–110. [Google Scholar] [CrossRef]
  74. 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]
  75. Dzrbashjan, M.M. Integral Transforms and Representations in the Complex Domain; Nauka: Moscow, Russia, 1966. (In Russian) [Google Scholar]
  76. Karp, D.; López, J.L. On a particular class of Meijer’s G functions appearing in fractional calculus. Int. J. Appl. Math. 2018, 31, 521–543. [Google Scholar] [CrossRef]
  77. Karp, D.; Prilepkina, E. Completely monotonic gamma ratio and infinitely divisible H-function of Fox. Comput. Methods Funct. Theory 2016, 16, 135–153. [Google Scholar] [CrossRef]
  78. Vellaisamy, P.; Kataria, K.K. The I-function distribution and its extensions. Teoria Veroyatnostej i ee Primenenia (Russ. Ed.) 2018, 63, 284–305. (In Russian) [Google Scholar] [CrossRef]
  79. Pogány, T. Integral form of Le Roy-type hypergeometric function. Integral Transform. Spec. Funct. 2018, 29, 580–584. [Google Scholar] [CrossRef]
  80. Kiryakova, V. Fractional calculus operators of special functions?—The result is well predictable! Chaos Solitons Fractals 2017, 102, 2–15. [Google Scholar] [CrossRef]
  81. Krätzel, E. Differentiationssätze der L-Transformation under Differentiagleichungen nach dem Operator. Math. Machrichten 1967, 35, 105–114. [Google Scholar] [CrossRef]
  82. Krätzel, E. Integral transformations of Bessel type. In Generalized Functions and Operational Calculus (Proc. Conf. Varna 1975); Bulgarian Academy of Sciences: Sofia, Bulgaria, 1979; pp. 148–155. [Google Scholar]
  83. Kilbas, A.A.; Saxena, R.K.; Trujillo, J.J. Krätzel function as a function of hypergeometric type. Fract. Calc. Appl. Anal. 2006, 9, 109–131. [Google Scholar]
  84. Dimovski, I.; Kiryakova, V. The Obrechkoff integral transform: Properties and relation to a generalized fractional calculus. Numer. Funct. Anal. Optimiz. 2000, 21, 121–144. [Google Scholar] [CrossRef]
  85. 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]
  86. Dimovski, I. Operational calculus for a class of differental operators. C. R. Acad. Bulg. Sci. 1966, 19, 1111–1114. [Google Scholar]
  87. Mainardi, F. A tutorial on the basic special functions of fractional calculus. WSEAS Trans. Math. 2020, 19, 74–98. [Google Scholar] [CrossRef]
  88. Delerue, P. Sur le calcul symboloque à n variables et fonctions hyperbesseliennes (II). Ann. Soc. Sci. Brux. Ser. 1 1953, 3, 229–274. [Google Scholar]
  89. Gorenflo, R.; Kilbas, A.A.; Rogosin, S. On the generalized Mittag-Leffler type function. Integral Transform. Spec. Funct. 1998, 7, 215–224. [Google Scholar] [CrossRef]
  90. Droghei, R. Properties of the multi-index special function Wα¯,ν¯(z). Fract. Calc. Appl. Anal. 2023, 26, 2057–2068. [Google Scholar] [CrossRef]
  91. Kiryakova, V. Fractional calculus of some “new” but not new special functions: k-, multi-index-, and S-analogues. AIP Conf. Proc. 2019, 2172, 050008. [Google Scholar] [CrossRef]
  92. Bazhlekova, E. Completely monotone multinomial Mittag-Leffler type functions and diffusion equations with multiple time-derivatives. Fract. Calc. Appl. Anal. 2021, 24, 88–111. [Google Scholar] [CrossRef]
  93. Gupta, K.C.; Soni, R.C. New properties of the hypergeometric series associated with Feynman integrals. Kyungpook Math. J. 2001, 41, 97–104. [Google Scholar]
  94. Jolly, N. New Investigations in Integral Transforms and Fractional Integral Operators Involving Generalized Extended Mittag-Leffer Function and Extended Hurwitz Lerch Zeta Function with Applications to the Solution of Fractional Differential Equations. Ph.D. Thesis, Malaviya National Institute of Technology, Jaipur, India, 2019. [Google Scholar]
  95. Srivastava, H.M. An introductory overwiew of fractional-calculus operators based upon the Fox-Wright and related higher transcendental functions. J. Adv. Eng. Comput. 2021, 5, 135–166. [Google Scholar] [CrossRef]
  96. Gerhold, S.; Tomovski, Ž. Asymptotic expansion of Mathieu power series and trigonometric Mathieu series. J. Math. Anal. Appl. 2019, 479, 1882–1892. [Google Scholar] [CrossRef]
  97. Stivastava, H.M.; Saxena, R.K.; Pogány, T.; Saxena, R. Integral and computational representations of the extended Hurwiz-Lerch zeta function. Integral Transform. Spec. Funct. 2011, 22, 487–506. [Google Scholar] [CrossRef]
  98. Bhatter, S.; Faisal, S.M.; Qureshi, M.I. A family of Mittag-Leffelr type functions and their properties. Palest. J. Math. 2015, 4, 367–373. [Google Scholar]
  99. Schilling, R.L.; Song, R.; Vondraček, Z. Bernstein Functions: Theory and Applications, 2nd ed.; De Gruyter: Berlin, Germany; Boston, MA, USA, 2012. [Google Scholar]
  100. Berg, C.; Çetinkaya, A.; Karp, D. Completely monotonic ratios of basic and ordinary gamma functions. Aequat. Math. 2021, 95, 569–588. [Google Scholar] [CrossRef]
  101. Gerhold, S. On some non-holonomic equences. Electr. J. Comb. 2004, 11, R87. [Google Scholar] [CrossRef]
  102. Bell, J.P.; Gerhold, S.; Klazar, M.; Luca, F. Non-holonomicity of sequences defined via elementary functions. arXiv 2006, arXiv:math/060514v1. [Google Scholar] [CrossRef]
  103. Flajolet, P.; Gerhold, S.; Salvy, B. Lindelöf representations and (non)-holonomic sequences. arXiv 2009, arXiv:0906.1957v2. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kiryakova, V.; Paneva-Konovska, J. Going Next after “A Guide to Special Functions in Fractional Calculus”: A Discussion Survey. Mathematics 2024, 12, 319. https://doi.org/10.3390/math12020319

AMA Style

Kiryakova V, Paneva-Konovska J. Going Next after “A Guide to Special Functions in Fractional Calculus”: A Discussion Survey. Mathematics. 2024; 12(2):319. https://doi.org/10.3390/math12020319

Chicago/Turabian Style

Kiryakova, Virginia, and Jordanka Paneva-Konovska. 2024. "Going Next after “A Guide to Special Functions in Fractional Calculus”: A Discussion Survey" Mathematics 12, no. 2: 319. https://doi.org/10.3390/math12020319

APA Style

Kiryakova, V., & Paneva-Konovska, J. (2024). Going Next after “A Guide to Special Functions in Fractional Calculus”: A Discussion Survey. Mathematics, 12(2), 319. https://doi.org/10.3390/math12020319

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