Integral Representations and Algebraic Decompositions of the Fox-Wright Type of Special Functions

The manuscript surveys the special functions of the Fox-Wright type. These functions are generalizations of the hypergeometric functions. Notable representatives of the type are the Mittag-Leffler functions and the Wright function. The integral representations of such functions are given and the conditions under which these function can be represented by simpler functions are demonstrated. The connection with generalized Erdélyi-Kober fractional differential and integral operators is demonstrated and discussed.


Introduction
This paper is concerned with integral representations of the Fox-Wright functions and their relationship to fractional calculus. The first characteristic exemplar of this function family has been introduced by E. M. Wright, who generalized the concept in a series of papers in 1930s. The Fox-Wright special functions are very general mathematical objects, which have broad applications in mathematical physics, notably in descriptions of wave phenomena, heat and mass transfer. They encompass the generalized hypergeometric functions p F q and are related to the family of the Bessel functions.
The conditions for existence of the generalized Wright function together with its representation in terms of the Mellin-Barnes integral and of the H-function can be found in [1]. Fox-Wright functions encompass other important families of functions, such as the the Mittag-Leffler functions (surveys in [2,3]). The Mittag-Leffler function in turn expresses the solution of fractional order integral or fractional order differential equations. It has applications in the theory of random walks, Levy flights, superdiffusive transport, among others. Another important example is the M-Wright function, which expresses the fundamental solution of the time-fractional diffusion-wave equation [4]. A recent survey about the properties of the function can be found in [5].
The objective of the present paper is to give a self-contained treatment of the generalized fractional calculus Erdélyi-Kober operators, which appear as re-parametrizations of the Euler integrals. The actions of the Erdélyi-Kober operators are thus expressed in a natural way as adding or removing parameters of multi-parameter Fox-Wright functions.
Many authors introduce the Fox-Wright functions from their representation as H-functions, which are in turn defined as Mellin transforms pairs. Such presentation tends to obfuscate the utility of Fox-Wright functions. In this paper, the Fox-Wright functions are represented as generalized hypergeoemtric series (GHG) and related to the theory of the Euler Gamma and Beta functions.
This provides some advantages, as for example Theorem 3, which to the present author's knowledge has not been stated in such form before.

Preliminaries and Notation
The generalized hypergeometric functions are defined by the infinite hypergeometric (HG) series The defining property fo HG series is that the coefficients are rational functions of the index variable (i.e., k).
In the following sections, we will use the parametric notation similar to the one adopted by Oldham and Spanier [6].
The classical hypergeometric series p F q obey the differential identities These relationships can be suitably generalized for fractional operators. The Fox-Wright functions are further generalizations of the hypergeometric (HG) functions of the form For this generalization, one can not expect that in general the coefficients are rational functions of the index variable.
For convergence of the series the condition will be assumed everywhere [7,8]. At this point we introduce some extended notation under the convention p+1Ψq (z) = a 1 , . . . , a p (A, a) z b 1 , . . . , b q − , p+1Ψq (0) = 1, In this notation, the hypergeometric indices of the function are written first while the non-simplified indices are left second. The non-simplified indices result in factors of the form or their reciprocals, respectively, and follow the usual convention established in literature. The order in the parametric convention for the arguments of the Gamma function follows the usual convention, used for example in [9,10]. This is unfortunately converse to the order of the more conventional Wright function and Mittag-Leffler type of functions.
The following simplifying convention will be used further: and This example shows different ways to write a hypergeometric function.

Algebraic Decomposition
The coefficients of the GHG series can be identified by means of the following Lemma: where q k is a linear expression in the index k or under the same convention. Then respectively, where a = q k+1 mod q k and A = q k mod k Proof. We prove the first case only since the second one follows identical reasoning. By hypothesis c k = Γ(A + ka) for some unknown A and a. Let's form the ratio Then (A + a + ka) mod (A + k a) = a and A + k a mod k = A.
The generalized hypergeometric series can be decomposed in symmetric (even) and anti-symmetric (odd) series as follows: Theorem 1 (GH Series Parity Decomposition). Suppose that the generalized hypergeometric series S = S e + S o is absolutely convergent at z. Denote S e as the even part while S o as the odd part.
If S is of the form Proof. Let the even part and odd series be S e and S o , respectively. We prove only the first statement because the second one can be proved in an identical way. For simplicity of calculations suppose that S is of the form − (A, a) z − − For the even part: so that the ratio of the coefficients is Therefore, For the odd part starting from k = 1 so the the ratio of the coefficients is Therefore, after an index shift m → m + 1 From this result, it can be seen that the hypergeometric (HG) series are not closed with regard to the parity decomposition. In contrast, the GHG series of the Wright type are closed with regard to this decomposition.
The simplest example illustrating the Theorem is given by the exponential series.
A non-trivial example of the present Theorem is the following For the first part by the Gamma duplication formula for the argument k + 1/2: For the second part by the Gamma duplication formula for the argument 2(k + 1): The negative multiplicative parameters can be raised to the numerator by the application of the following Theorem: Proof. Consider the monomial By the reflection formula This can be embedded in the complex monomial expression Assuming that z is real, the original expression B k is the imaginary part of C k . Further, C k has modulus so that the infinite series for C k converges and so does its imaginary part. Finally the GHG parameters can be read off from the arguments of the Gamma functions.

Integral Representations by Beta Integrals
The main result of this section is given by the theorem below. The result allows for the representation of a GHG function of order (p + 1, q + 1) in terms of an integral of a GHG function of order (p, q + 1) or in special cases (p, q). In this section everywhere the argument of a GHG function will be considered real-valued.

Theorem 3 (Beta integral representation).
For B > A > 0 and b ≥ a the following representation holds By change of variables t = 1/(1 + u) Proof. The proof follows from the hypothesis of absolute convergence of the series. Therefore, the order of integration and summation can be switched. Let w = Γ(B)/Γ(A) and suppose that a = b. Observe that by Equation (A3) Therefore, by absolute convergence of the series Therefore, Furthermore, let now a = b = 1. It can be further observed that for the monomial term This representation step reduces a (p + 1, q + 1) GHG series into a (p, q + 1) GHG series. It can be seen that the reduction via Beta integral is not complete except if a = b. Therefore, it can be instructive to distinguish homogeneous GHG series with indices a i = b i and different multiplicies. This is the subject of the following results: where Hankel path Ha + starts at infinity on the real axis, encircling 0 in a positive sense, and returns to infinity along the real axis, respecting the cut along the positive real axis, while Ha − is its reflection.
Proof. From the Heine's formula for the reciprocal Gamma function representation Ha + Moreover, the complex integral can be converted to a real integral for suitable functions [11].

Theorem 5 (Real GH Series Representation).
Suppose that all indices a i and b i are real. Then for some real A > 0 and a > 0 and z < 1 provided that all parameters are real.
In summary, the section shows that a (p, q) GHG series can be reduced to a p + q multiple integrals of the Euler type.

Mittag-Leffler Functions
The 2 parameter Mittag-Leffler function [13,14] under the present convention will be denoted as This immediately gives the complex integral representation according to Corollary 4 However, in this case the contour encloses the curve |1 − z/τ a | = 1. Another example is the 3-parameter Mittag-Leffler function generalization, that is the Prabhakar function [15] defined as which leads to an integral involving the Wright function.
An interesting special case is the function E γ a,1 (z) which is a confluent Kummer ( 1 F 1 ) hypergeometric function. In this case for a > γ

The Kummer-Wright Function
In particular, the following representation can be stated for the basic GH function (the Kummer-Wright function) And also This Wright function in turn can be represented as (see Equation (A4)): For the homogeneous case then according to the representation of the reciprocal Gamma function. Therefore, according to the Homogeneous Reduction Corollary: which is a generalization of the Kummer function.

Generalized Fractional Erdélyi-Kober Operations
The theory of GHG series has an interesting relationship with the generalized fractional calculus. The right Erdélyi-Kober (E-K) fractional integrals are defined as [10]: while another equivalent form is [16] [Ch. 18]: The two forms are related by the change of variables u = z τ 1/β . Therefore, from the preceding presentation it follows that This corresponds to the findings of Kiryakova [17]. The E-K operator reduces to the Riemann-Liouville fractional integral for β = 1 and γ = 0 as It is interesting to note further that the EK operators map the functions of the Dimovski space D α into one another [18]: I γ,δ β : D α → D α , α > −β(γ + 1) where the function space is given by the following: Definition 1 (Dimovski Space [19]). The space of functions D α consists of all functions f (x), x > 0, that can be represented in the form f (x) = x p f 1 (x) with p > α and f 1 ∈ C([0, ∞)).
The corresponding generalized Erdélyi-Kober fractional derivative is defined by a composition product as where < δ > is the fractional part and [δ] is the integral part of the number. For β = 1 the E-K fractional derivative operator reduces to the Riemann-Liouville fractional derivative of order δ as The last equation can be used for index reduction.

Discussion
These results demonstrate that the homogeneous class of GHG series are closed with respect to the (generalized) fractional calculus operations. One can, therefore, expect that the solutions of fractional differential equations could be expressed in terms of Fox-Wright functions. Therefore, the Fox-Wright functions are the most general functions of the mathematical physics.
Furthermore, the results demonstrate a strong link between the special function theory and the theory of the Euler Beta and Gamma functions. It appears that the E-K operators can be thought of Euler's Beta integrals in disguise. Moreover, the Gamma integral representations can be interpreted as Laplace or Inverse Laplace transforms.
Finally, the main consequence of so-presented results is that all GHG functions of the Fox-Wright type can be represented as multiple (complex) integrals of three primitive special functions of orders (1, 0), (0, 1) and (1, 1) respectively. This corroborates the findings of Kiryakova [9,10,17]. These multiple integrals can be denoted as generalized fractional differ-integrals; however, this line of representation is superfluous to the necessities of the numerical and physical modeling.

Appendix A. Euler Integrals
The Gamma integral, i.e., the Euler integral of second kind is defined as while for all z / ∈ Z − Γ(z) = 1 2i sin πz Ha − e τ τ z−1 dτ, τ ∈ C The complex representation for the reciprocal Gamma function is given by the Heine's integral as The contour is depicted in Figure A1. For non-integral arguments the branch cut is selected as the negative real axis.
The Beta integral (i.e., the Euler integral of first kind) is given as The Beta function can be continued analytically along the self-intersecting Pochhammer contour as DIMITER PRODANOV  Kiryakova [5]. These multiple integrals can be denoted as generalized fractional differeitegrals [6], however this line of representation is superfluous to the necessities of the numerical (i.e. physical) modeling.

Appendix B. The Wright Function
The function W ( λ, µ| z), named after E. M. Wright, is defined as the infinite series