Meijer G-function

Meijer G-function

In mathematics, the G-function was introduced by Cornelis Simon Meijer (1936) as a very general function intended to include most of the known special functions as particular cases. This was not the only attempt of its kind: the generalized hypergeometric function and the MacRobert E-function had the same aim, but Meijer's G-function was able to include those as particular cases as well. The first definition was made by Meijer using a series; nowadays the accepted and more general definition is via a path integral in the complex plane, introduced firstly by Arthur Erdélyi in 1953. With the current definition, the majority of the special functions can be represented in terms of the G-function and of the gamma function.

A still more general function, which introduces additional parameters into Meijer's G-function is Fox's H-function.

Contents

Definition of the Meijer G-function

A general definition of the Meijer G-function is given by the following line integral in the complex plane:


G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1, \dots, a_p \\ b_1, \dots, b_q \end{matrix} \; \right| \, z \right) = \frac{1}{2 \pi i} \int_L \frac{\prod_{j=1}^m \Gamma(b_j - s) \prod_{j=1}^n \Gamma(1 - a_j +s)} {\prod_{j=m+1}^q \Gamma(1 - b_j + s) \prod_{j=n+1}^p \Gamma(a_j - s)} z^s \,ds.

This integral is of the so-called Mellin–Barnes type, and may be viewed as an inverse Mellin transform. The definition holds under the following assumptions:

  • 0 ≤ mq and 0 ≤ np, where m, n, p and q are integer numbers
  • akbj ≠ 1, 2, 3, ... for k = 1, 2, ..., n and j = 1, 2, ..., m, which implies that no pole of any Γ(bjs), j = 1, 2, ..., m, coincides with any pole of any Γ(1 − ak + s), k = 1, 2, ..., n
  • z ≠ 0

Note that for historical reasons the first lower and second upper index refer to the top parameter row, while the second lower and first upper index refer to the bottom parameter row. One often encounters the following more synthetic notation using vectors:


G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1, \dots, a_p \\ b_1, \dots, b_q \end{matrix} \; \right| \, z \right) =
G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \,  z \right) .

Implementions of the G-function in computer algebra systems typically employ separate vector arguments for the four (possibly empty) parameter groups a1 ... an, an+1 ... ap, b1 ... bm, and bm+1 ... bq, and thus can omit the orders p, q, n, and m as redundant.

The L in the integral represents the path to be followed while integrating. Three choices are possible for this path:

1. L runs from −i∞ to +i∞ such that all poles of Γ(bjs), j = 1, 2, ..., m, are on the right of the path, while all poles of Γ(1 − ak + s), k = 1, 2, ..., n, are on the left. The integral then converges for |arg z| < δ π, where

\delta = m + n - \tfrac{1}{2} (p+q) ;
an obvious prerequisite for this is δ > 0. The integral additionally converges for |arg z| = δ π ≥ 0 if (q − p) (σ + 12) > Re(ν) + 1, where σ represents Re(s) as the integration variable s approaches both +i∞ and −i∞, and where

\nu = \sum_{j = 1}^q b_j - \sum_{j = 1}^p a_j .
As a corollary, for |arg z| = δ π and p = q the integral converges independent of σ whenever Re(ν) < −1.
2. L is a loop beginning and ending at +∞, encircling all poles of Γ(bjs), j = 1, 2, ..., m, exactly once in the negative direction, but not encircling any pole of Γ(1 − ak + s), k = 1, 2, ..., n. Then the integral converges for all z if q > p ≥ 0; it also converges for q = p > 0 as long as |z| < 1. In the latter case, the integral additionally converges for |z| = 1 if Re(ν) < −1, where ν is defined as for the first path.
3. L is a loop beginning and ending at −∞ and encircling all poles of Γ(1 − ak + s), k = 1, 2, ..., n, exactly once in the positive direction, but not encircling any pole of Γ(bjs), j = 1, 2, ..., m. Now the integral converges for all z if p > q ≥ 0; it also converges for p = q > 0 as long as |z| > 1. As already stated for the second path, in the case of p = q the integral also converges for |z| = 1 when Re(ν) < −1.

The conditions for convergence are readily established by applying Stirling's asymptotic approximation to the gamma functions in the integrand. When the integral converges for more than one of these paths, the results of integration can be shown to agree; if it converges for only one path, then this is the only one to be considered. In fact, numerical path integration in the complex plane constitutes a practicable and sensible approach to the calculation of Meijer G-functions.

As a consequence of this definition, the Meijer G-function is an analytic function of z with possible exception of the origin z = 0 and of the unit circle |z| = 1.

Differential equation

The G-function satisfies the following linear differential equation of order max(p,q):


\left[ (-1)^{p - m - n} \;z \prod_{j = 1}^p \left( z \frac{d}{dz} - a_j + 1 \right) - \prod_{j = 1}^q \left( z \frac{d}{dz} - b_j \right) \right] G(z) = 0.

For a fundamental set of solutions of this equation in the case of pq one may take:


G_{p,q}^{\,1,p} \!\left( \left. \begin{matrix} a_1, \dots, a_p \\ b_h, b_1, \dots, b_{h-1}, b_{h+1}, \dots, b_q \end{matrix} \; \right| \, (-1)^{p-m-n+1} \;z \right), \quad h = 1,2,\dots,q,

and similarly in the case of pq:


G_{p,q}^{\,q,1} \!\left( \left. \begin{matrix} a_h, a_1, \dots, a_{h-1}, a_{h+1}, \dots, a_p \\ b_1, \dots, b_q \end{matrix} \; \right| \, (-1)^{q-m-n+1} \;z \right), \quad h = 1,2,\dots,p.

These particular solutions are analytic except for a possible singularity at z = 0 (as well as a possible singularity at z = ∞), and in the case of p = q also an inevitable singularity at z = (−1)pmn. As will be seen presently, they can be identified with generalized hypergeometric functions pFq−1 of argument (−1)pmn z that are multiplied by a power zbh, and with generalized hypergeometric functions qFp−1 of argument (−1)qmn 1z that are multiplied by a power zah−1, respectively.

Relationship between the G-function and the generalized hypergeometric function

If the integral converges when evaluated along the second path introduced above, and if no confluent poles appear among the Γ(bjs), j = 1, 2, ..., m, then the Meijer G-function can be expressed as a sum of residues in terms of generalized hypergeometric functions pFq−1 (Slater's theorem):


G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) = \sum_{h=1}^m \frac{\prod_{j=1}^m \Gamma(b_j - b_h)^* \prod_{j=1}^n \Gamma(1+b_h - a_j) \; z^{b_h}} {\prod_{j=m+1}^q \Gamma(1+b_h - b_j) \prod_{j=n+1}^p \Gamma(a_j - b_h)} \times

\times \; _{p}F_{q-1} \!\left( \left. \begin{matrix} 1+b_h - \mathbf{a_p} \\ (1+b_h - \mathbf{b_q})^* \end{matrix} \; \right| \, (-1)^{p-m-n} \; z \right) .

For the integral to converge along the second path one must have either p < q, or p = q and |z| < 1, and for the poles to be distinct no pair among the bj, j = 1, 2, ..., m, may differ by an integer or zero. The asterisks in the relation remind us to ignore the contribution with index j = h as follows: In the product this amounts to replacing Γ(0) with 1, and in the argument of the hypergeometric function, if we recall the meaning of the vector notation,


1 + b_h - \mathbf{b_q} = (1 + b_h - b_1), \,\dots, \,(1 + b_h - b_j), \,\dots, \,(1 + b_h - b_q),

this amounts to shortening the vector length from q to q−1.

Note that when m = 0, the second path does not contain any pole, and so the integral must vanish identically,


G_{p,q}^{\,0,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) = 0,

if either p < q, or p = q and |z| < 1.

Similarly, if the integral converges when evaluated along the third path above, and if no confluent poles appear among the Γ(1 − ak + s), k = 1, 2, ..., n, then the G-function can be expressed as:


G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) = \sum_{h=1}^n \frac{\prod_{j=1}^n \Gamma(a_h - a_j)^* \prod_{j=1}^m \Gamma(1-a_h + b_j) \; z^{a_h-1}} {\prod_{j=n+1}^p \Gamma(1-a_h + a_j) \prod_{j=m+1}^q \Gamma(a_h - b_j)} \times

\times \; _{q}F_{p-1} \!\left( \left. \begin{matrix} 1-a_h + \mathbf{b_q} \\ (1-a_h + \mathbf{a_p})^* \end{matrix} \; \right| \, \frac{(-1)^{q-m-n}}{z} \right) .

For this, either p > q, or p = q and |z| > 1 are required, and no pair among the ak, k = 1, 2, ..., n, may differ by an integer or zero. For n = 0 one consequently has:


G_{p,q}^{\,m,0} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) = 0,

if either p > q, or p = q and |z| > 1.

On the other hand, any generalized hypergeometric function can readily be expressed in terms of the Meijer G-function:


\; _{p}F_{q} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right)
= \frac  {\Gamma(\mathbf{b_q})}{\Gamma(\mathbf{a_p})} \; G_{p,\,q+1}^{\,1,\,p} \!\left( \left. \begin{matrix} 1-\mathbf{a_p} \\ 0,1 - \mathbf{b_q} \end{matrix} \; \right| \, -z \right)
= \frac {\Gamma(\mathbf{b_q})}{\Gamma(\mathbf{a_p})}  \; G_{q+1,\,p}^{\,p,\,1} \!\left( \left. \begin{matrix} 1,\mathbf{b_q} \\ \mathbf{a_p} \end{matrix} \; \right| \, \frac{-1}{z} \right),

where we have made use of the vector notation:


\Gamma(\mathbf{a_p}) = \prod_{j=1}^p \Gamma(a_j).

This holds unless a nonpositive integer value of at least one of its parameters ap reduces the hypergeometric function to a finite polynomial, in which case the gamma prefactor of either G-function diverges and the parameter sets of the G-functions violate the requirement akbj ≠ 1, 2, 3, ... for k = 1, 2, ..., n and j = 1, 2, ..., m from the definition above. Apart from this restriction, the relationship is valid whenever the generalized hypergeometric series pFq(z) converges, i. e. for any finite z when pq, and for |z| < 1 when p = q + 1. In the latter case, the relation with the G-function automatically provides the analytic continuation of pFq(z) to |z| ≥ 1 with a branch cut from 1 to ∞ along the real axis. Finally, the relation furnishes a natural extension of the definition of the hypergeometric function to orders p > q + 1. By means of the G-function we can thus solve the generalized hypergeometric differential equation for p > q + 1 as well.

Polynomial cases

To express polynomial cases of generalized hypergeometric functions in terms of Meijer G-functions, a linear combination of two G-functions is needed:


\; _{p+1}F_{q} \!\left( \left. \begin{matrix} -h, \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) = h! \; \frac{\prod_{j=n+1}^p \Gamma(1 - a_j) \prod_{j=m+1}^q \Gamma(b_j)} {\prod_{j=1}^n \Gamma(a_j) \prod_{j=1}^m \Gamma(1 - b_j)} \times

\times\left[ G_{p+1,\,q+1}^{\,m+1,\,n} \!\left( \left. \begin{matrix} 1-\mathbf{a_p}, h+1 \\ 0, 1-\mathbf{b_q} \end{matrix} \; \right| \, (-1)^{p-m-n} \; z \right) + (-1)^h \; G_{p+1,\,q+1}^{\,m,\,n+1} \!\left( \left. \begin{matrix} h+1, 1-\mathbf{a_p} \\ 1-\mathbf{b_q}, 0 \end{matrix} \; \right| \, (-1)^{p-m-n} \; z \right) \right] ,

where h = 0, 1, 2, ... equals the degree of the polynomial p+1Fq(z). The orders m and n can be chosen freely in the ranges 0 ≤ mq and 0 ≤ np, which allows to avoid that specific integer values or integer differences among the parameters ap and bq of the polynomial give rise to divergent gamma prefactors or to a conflict with the definition of the G-function. Again, the formula can be verified by expressing the two G-functions as sums of residues; no cases of confluent poles permitted by the definition of the G-function need be excluded here.

Basic properties of the G-function

As can be seen from the definition of the G-function, if equal parameters appear among the ap and bq determining the factors in the numerator and the denominator of the integrand, the fraction can be simplified, and the order of the function thereby be reduced. Whether the order m or n will decrease depends of the particular position of the parameters in question. Thus, if one of the ak, k = 1, 2, ..., n, equals one of the bj, j = m + 1, ..., q, the G-function lowers its orders p, q and n:


G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1, a_2, \dots, a_p \\ b_1, \dots, b_{q-1}, a_1 \end{matrix} \; \right| \, z \right) =
G_{p-1,\,q-1}^{\,m,\,n-1} \!\left( \left. \begin{matrix} a_2, \dots, a_p \\ b_1, \dots, b_{q-1} \end{matrix} \; \right| \, z \right), \quad n,p,q \geq 1.

For the same reason, if one of the ak, k = n + 1, ..., p, equals one of the bj, j = 1, 2, ..., m, then the G-function lowers its orders p, q and m:


G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1, \dots, a_{p-1}, b_1 \\ b_1, b_2, \dots, b_q \end{matrix} \; \right| \, z \right) =
G_{p-1,\,q-1}^{\,m-1,\,n} \!\left( \left. \begin{matrix} a_1, \dots, a_{p-1} \\ b_2, \dots, b_q \end{matrix} \; \right| \, z \right), \quad m,p,q \geq 1.

Starting from the definition, it is also possible to derive the following properties:


z^{\rho} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) =
G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} + \rho \\ \mathbf{b_q} + \rho \end{matrix} \; \right| \, z \right),

G_{p+2,\,q}^{\,m,\,n+1} \!\left( \left. \begin{matrix} \alpha, \mathbf{a_p}, \alpha' \\ \mathbf{b_q} \end{matrix} \; \right| \, z  \right) =
(-1)^{\alpha'-\alpha} \; G_{p+2,\,q}^{\,m,\,n+1} \!\left( \left. \begin{matrix} \alpha', \mathbf{a_p}, \alpha \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right), \quad n \leq p, \; \alpha'-\alpha = 0,\pm1,\pm2,\dots,

G_{p,\,q+2}^{\,m+1,\,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \beta, \mathbf{b_q}, \beta' \end{matrix} \; \right| \, z \right) =
(-1)^{\beta'-\beta} \; G_{p,\,q+2}^{\,m+1,\,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \beta', \mathbf{b_q}, \beta \end{matrix} \; \right| \, z \right), \quad m \leq q, \; \beta'-\beta = 0,\pm1,\pm2,\dots,

G_{p+1,\,q+1}^{\,m,\,n+1} \!\left( \left. \begin{matrix} \alpha, \mathbf{a_p} \\ \mathbf{b_q}, \beta \end{matrix} \; \right| \, z \right) =
(-1)^{\beta-\alpha} \; G_{p+1,\,q+1}^{\,m+1,\,n} \!\left( \left. \begin{matrix} \mathbf{a_p}, \alpha \\ \beta, \mathbf{b_q} \end{matrix} \; \right| \, z \right), \quad m \leq q, \; \beta-\alpha = 0,1,2,\dots,
 
G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) =
G_{q,p}^{\,n,m} \!\left( \left. \begin{matrix} 1-\mathbf{b_q} \\ 1-\mathbf{a_p} \end{matrix} \; \right| \, \frac{1}{z} \right),

G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) =
\frac{k^{1+\nu+(p-q)/2}} {(2 \pi)^{(k-1) \delta}} \; G_{k p, \, k q}^{\, k m, \, k n} \!\left( \left. \begin{matrix} a_1/k, \dots, (a_1+k-1)/k, \dots, a_p/k, \dots, (a_p+k-1)/k \\ b_1/k, \dots, (b_1+k-1)/k, \dots, b_q/k, \dots, (b_q+k-1)/k \end{matrix} \; \right| \, \frac{z^k} {k^{k(q-p)}} \right), \quad k = 1,2,3,\dots.

The abbreviations ν and δ were introduced in the definition of the G-function above.

Derivatives and antiderivatives

Concerning derivatives of the G-function, there are these relationships:


\frac{d}{dz} \left[ z^{1-a_1} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) \right] =
z^{-a_1} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1 - 1, a_2, \dots, a_p \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right), \quad n \geq 1,

\frac{d}{dz} \left[ z^{1-a_p} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) \right] =
- z^{-a_p} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1, \dots, a_{p-1}, a_p - 1 \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right), \quad n < p.

\frac{d}{dz} \left[ z^{-b_1} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) \right] =
- z^{-1-b_1} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ b_1 + 1, b_2, \dots, b_q \end{matrix} \; \right| \, z \right), \quad m \geq 1,

\frac{d}{dz} \left[ z^{-b_q} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) \right] =
z^{-1-b_q} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ b_1, \dots, b_{q-1}, b_q + 1 \end{matrix} \; \right| \, z \right), \quad m < q,

From these four, equivalent relations can be deduced by simply calculating the derivative on the left-hand side and manipulating a bit. One obtains for example:


z \frac{d}{dz} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) =
G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1 -1, a_2, \dots, a_p \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) +
(a_1 - 1) \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right), \quad n \geq 1.

Moreover, for derivatives of arbitrary order k, one has


z^k \frac{d^k}{dz^k} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) =
G_{p+1,\,q+1}^{\,m,\,n+1} \!\left( \left. \begin{matrix} 0, \mathbf{a_p} \\ \mathbf{b_q}, k \end{matrix} \; \right| \, z \right) =
(-1)^k \; G_{p+1,\,q+1}^{\,m+1,\,n} \!\left( \left. \begin{matrix} \mathbf{a_p}, 0 \\ k, \mathbf{b_q} \end{matrix} \; \right| \, z \right),

z^k \frac{d^k}{dz^k} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, \frac{1}{z} \right) =
G_{p+1,\,q+1}^{\,m+1,\,n} \!\left( \left. \begin{matrix} \mathbf{a_p}, 1-k \\ 1, \mathbf{b_q} \end{matrix} \; \right| \, \frac{1}{z} \right) =
(-1)^k \; G_{p+1,\,q+1}^{\,m,\,n+1} \!\left( \left. \begin{matrix} 1-k, \mathbf{a_p} \\ \mathbf{b_q}, 1 \end{matrix} \; \right| \, \frac{1}{z} \right),

which hold for k < 0 as well, thus allowing to obtain the antiderivative of any G-function as easily as the derivative. By choosing one or the other of the two results provided in either formula, one can always prevent the set of parameters in the result from violating the condition akbj ≠ 1, 2, 3, ... for k = 1, 2, ..., n and j = 1, 2, ..., m that is imposed by the definition of the G-function. Note that each pair of results becomes unequal in the case of k < 0.

From these relationships, corresponding properties of the Gauss hypergeometric function and of other special functions can be derived.

Recurrence relations

By equating different expressions for the first-order derivatives, one arrives at the following 3-term recurrence relations among contiguous G-functions:


(a_p - a_1) \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) =
G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1-1, a_2, \dots, a_p \\ b_1, \dots, b_q \end{matrix} \; \right| \, z \right) +
G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1, \dots, a_{p-1}, a_p-1 \\ b_1, \dots, b_q \end{matrix} \; \right| \, z \right), \quad 1 \leq n < p,

(b_1 - b_q) \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) =
G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1, \dots, a_p \\ b_1+1, b_2, \dots, b_q \end{matrix} \; \right| \, z \right) +
G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1, \dots, a_p \\ b_1, \dots, b_{q-1}, b_q+1 \end{matrix} \; \right| \, z \right), \quad 1 \leq m < q,

(b_1 - a_1 + 1) \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) =
G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1-1, a_2, \dots, a_p \\ b_1, \dots, b_q \end{matrix} \; \right| \, z \right) +
G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1, \dots, a_p \\ b_1+1, b_2, \dots, b_q \end{matrix} \; \right| \, z \right), \quad n \geq 1, \; m \geq 1,

(a_p - b_q - 1) \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right) =
G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1, \dots, a_{p-1}, a_p-1 \\ b_1, \dots, b_q \end{matrix} \; \right| \, z \right) +
G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1, \dots, a_p \\ b_1, \dots, b_{q-1}, b_q+1 \end{matrix} \; \right| \, z \right), \quad n < p, \; m < q.

Similar relations for the diagonal parameter pairs a1, bq and b1, ap follow by suitable combination of the above. Again, corresponding properties of hypergeometric and other special functions can be derived from these recurrence relations.

Multiplication theorems

Provided that z ≠ 0, the following relationships hold:


G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, w z \right) =
w^{b_1} \sum_{k=0}^{\infty} \frac{(1 - w)^k}{k!} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ b_1+k, b_2, \dots, b_q \end{matrix} \; \right| \, z \right), \quad m \geq 1,

G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, w z \right) =
w^{b_q} \sum_{k=0}^{\infty} \frac{(w - 1)^k}{k!} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ b_1, \dots, b_{q-1}, b_q+k \end{matrix} \; \right| \, z \right), \quad m < q,

G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, \frac{z}{w} \right) =
w^{1-a_1} \sum_{k=0}^{\infty} \frac{(1 - w)^k}{k!} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1-k, a_2, \dots, a_p \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right), \quad n \geq 1,

G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, \frac{z}{w} \right) =
w^{1-a_p} \sum_{k=0}^{\infty} \frac{(w - 1)^k}{k!} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} a_1, \dots, a_{p-1}, a_p-k \\ \mathbf{b_q} \end{matrix} \; \right| \, z \right), \quad n < p.

These follow by Taylor expansion about w = 1, with the help of the basic properties discussed above. The radii of convergence will be dependent on the value of z and on the G-function that is expanded. The expansions can be regarded as generalizations of similar theorems for Bessel, hypergeometric and confluent hypergeometric functions.

Definite integrals involving the G-function

Among definite integrals involving an arbitrary G-function one has:


\int_0^{\infty} x^{s - 1} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, \eta x \right) dx =
\frac{\eta^{-s} \prod_{j = 1}^{m} \Gamma (b_j + s) \prod_{j = 1}^{n} \Gamma (1 - a_j - s)} {\prod_{j = m + 1}^{q} \Gamma (1 - b_j - s) \prod_{j = n + 1}^{p} \Gamma (a_j + s)}.

Note that the restrictions under which this integral exists have been omitted here. It is, of course, no surprise that the Mellin transform of a G-function should lead back to the integrand appearing in the definition above.

Euler-type integrals for the G-function are given by:


\int_0^1 x^{-\alpha} \; (1-x)^{\alpha - \beta - 1} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z x \right) dx =
\Gamma (\alpha - \beta) \; G_{p+1 ,\, q+1}^{\,m ,\, n+1} \!\left( \left. \begin{matrix} \alpha, \mathbf{a_p} \\ \mathbf{b_q}, \beta \end{matrix} \; \right| \, z \right),

\int_1^\infty x^{-\alpha} \; (x-1)^{\alpha - \beta - 1} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, z x \right) dx =
\Gamma (\alpha - \beta) \; G_{p+1 ,\, q+1}^{\,m+1 ,\, n} \!\left( \left. \begin{matrix} \mathbf{a_p}, \alpha \\ \beta, \mathbf{b_q} \end{matrix} \; \right| \, z \right).

Here too, the restrictions under which the integrals exist have been omitted. Note that, in view of their effect on the G-function, these integrals can be used to define the operations of fractional differentiation and fractional integration for a fairly large class of functions (Erdélyi–Kober operators).

A result of fundamental importance is that the definite integral over a product of two arbitrary G-functions can be represented by just another G-function (convolution theorem):


\int_0^{\infty} G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, \eta x \right)
G_{\sigma, \tau}^{\,\mu, \nu} \!\left( \left. \begin{matrix} \mathbf{c_{\sigma}} \\ \mathbf{d_\tau} \end{matrix} \; \right| \, \omega x \right) dx =

= \frac{1}{\eta} \; G_{q + \sigma ,\, p + \tau}^{\,n + \mu ,\, m + \nu} \!\left( \left. \begin{matrix} - b_1, \dots, - b_m, \mathbf{c_{\sigma}}, - b_{m+1}, \dots, - b_q \\ - a_1, \dots, -a_n, \mathbf{d_\tau} , - a_{n+1}, \dots, - a_p \end{matrix} \; \right| \, \frac{\omega}{\eta} \right) =
 = \frac{1}{\omega} \; G_{p + \tau ,\, q + \sigma}^{\,m + \nu ,\, n + \mu} \!\left( \left. \begin{matrix} a_1, \dots, a_n, -\mathbf{d_\tau} , a_{n+1}, \dots, a_p \\ b_1, \dots, b_m, -\mathbf{c_{\sigma}}, b_{m+1}, \dots, b_q \end{matrix} \; \right| \, \frac{\eta}{\omega} \right) .

Again, the restrictions under which the integral exists have been omitted here. Note how the Mellin transform of the result merely assembles the gamma factors from the Mellin transforms of the two functions in the integrand. Many of the amazing definite integrals listed in tables or produced by computer algebra systems are nothing but special cases of this formula.

The convolution formula can be derived by substituting the defining Mellin–Barnes integral for one of the G-functions, reversing the order of integration, and evaluating the inner Mellin-transform integral. The preceding Euler-type integrals follow analogously.

Laplace transform

Using the above convolution integral and basic properties one can show that:


\int_0^{\infty} e^{- \omega x} \; x^{- \alpha} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, \eta x \right) dx =
\omega^{\alpha - 1} \; G_{p + 1,\,q}^{\,m,\,n+1} \!\left( \left. \begin{matrix} \alpha, \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, \frac{\eta}{\omega} \right) ,

where Re(ω) > 0. This is the Laplace transform of a function G(ηx) multiplied by a power xα; if we put α = 0 we get the Laplace transform of the G-function. As usual, the inverse transform is then given by:


x^{- \alpha} \; G_{p,\,q+1}^{\,m,\,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q}, \alpha \end{matrix} \; \right| \, \eta x \right) = 
\frac{1}{2 \pi i} \int_{c - i \infty}^{c + i \infty} e^{\omega x} \; \omega^{\alpha - 1} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, \frac{\eta}{\omega} \right) d\omega,

where c is a real positive constant that places the integration path to the right of any pole in the integrand.

Another formula for the Laplace transform of a G-function is:


\int_{0}^{\infty} e^{- \omega x} \; G_{p,q}^{\,m,n} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, \eta x^2 \right) dx =
\frac{1}{\sqrt{\pi} \omega} \; G_{p+2,\,q}^{\,m,\,n+2} \!\left( \left. \begin{matrix} 0, \frac{1}{2}, \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \, \frac{4 \eta}{\omega^2} \right) ,

where again Re(ω) > 0. Details of the restrictions under which the integrals exist have been omitted in both cases.

Integral transforms using the G-function

In general, two functions k(z,y) and h(z,y) are called transform kernels if, for any suitable function f(z) or any suitable function g(z), the following two relationships hold simultaneously:


g(z) = \int_{0}^{\infty} k(z,y) \, f(y) \; dy, \quad
f(z) = \int_{0}^{\infty} h(z,y) \, g(y) \; dy.

The two kernels are said to be symmetric if k(z,y) = h(z,y).

Narain transform

Roop Narain (1962, 1963) showed that the functions:


k(z,y) = 2 \gamma \; (zy)^{\gamma - 1/2} \; G_{p+q,\,m+n}^{\,m,\,p} \!\left( \left. \begin{matrix} \mathbf{a_p}, \mathbf{b_q} \\ \mathbf{c_m}, \mathbf{d_n} \end{matrix} \; \right| \, (zy)^{2 \gamma} \right),

h(z,y) = 2 \gamma \; (zy)^{\gamma - 1/2} \; G_{p+q,\,m+n}^{\,n,\,q} \!\left( \left. \begin{matrix} -\mathbf{b_q}, -\mathbf{a_p} \\ -\mathbf{d_n}, -\mathbf{c_m} \end{matrix} \; \right| \, (zy)^{2 \gamma} \right)

are two asymmetric transform kernels, where γ > 0, np = mq > 0, and:


\sum_{j=1}^p a_j + \sum_{j=1}^q b_j = \sum_{j=1}^m c_j + \sum_{j=1}^n d_j,

along with further convergence conditions. In particular, if p = q, m = n, aj + bj = 0 for j = 1, 2, ..., p and cj + dj = 0 for j = 1, 2, ..., m, then the two kernels become symmetric. The well-known Hankel transform is a symmetric special case of the Narain transform (γ = 1, p = q = 0, m = n = 1, c1 = −d1 = ν2).

Wimp transform

Jet Wimp (1964) showed that these two functions are asymmetric transform kernels:


k(z,y) = G_{p+2,\,q}^{\,m,\,n+2} \!\left( \left. \begin{matrix} 1 - \nu + i z, 1 - \nu - i z, \mathbf{a_p} \\ \mathbf{b_q} \end{matrix} \; \right| \; y \right),

h(z,y) = \frac{i}{\pi} y e^{- \nu \pi i} \left[ e^{\pi y} A(\nu + i y, \nu - i y \,|\, z e^{i \pi} ) - e^{- \pi y} A(\nu - i y, \nu + i y \,|\, z e^{i \pi} ) \right],

where the function A(·) is defined as:


A(\alpha, \beta \,|\, z) = G_{p+2,\,q}^{\,q-m,\,p-n+1} \!\left( \left. \begin{matrix} -a_{n+1}, -a_{n+2}, \dots, -a_p, \alpha, -a_1, -a_2, \dots, -a_n, \beta \\ -b_{m+1}, -b_{m+2}, \dots, -b_q, -b_1, -b_2, \dots, -b_m \end{matrix} \; \right| \, z \right).

Generalized Laplace transform

The Laplace transform can be generalized in close analogy with Narain's generalization of the Hankel transform:


g(s) = 2 \gamma \int_0^{\infty} (st)^{\gamma + \rho - 1/2} \; G_{p,\,q+1}^{\,q+1,\,0} \!\left( \left. \begin{matrix} \mathbf{a_p} \\ 0, \mathbf{b_q} \end{matrix} \; \right| \, (st)^{2 \gamma} \right) f(t) \; dt,

f(t) = \frac {\gamma} {\pi i} \int_{c - i \infty}^{c + i \infty} (ts)^{\gamma - \rho - 1/2} \; G_{p,\,q+1}^{\,1,\,p} \!\left( \left. \begin{matrix} -\mathbf{a_p} \\ 0, -\mathbf{b_q} \end{matrix} \; \right| \, -(ts)^{2 \gamma} \right) g(s) \; ds,

where γ > 0 and pq, and where the constant c > 0 places the second integration path to the right of any pole in the integrand. For γ = 12, ρ = 0 and p = q = 0, this corresponds to the familiar Laplace transform.

Meijer transform

Two particular cases of this generalization were given by C.S. Meijer in 1940 and 1941. The case resulting for γ = 1, ρ = −ν, p = 0, q = 1 and b1 = ν may be written (Meijer 1940a, 1940b):


g(s) = \sqrt {2 / \pi} \int_0^{\infty} (st)^{1/2} \, K_{\nu}(st) \, f(t) \; dt,

f(t) = \frac {1} {\sqrt {2 \pi} \,i} \int_{c - i \infty}^{c + i \infty} (ts)^{1/2} \, I_{\nu}(ts) \, g(s) \; ds,

and the case obtained for γ = 12, ρ = −mk, p = q = 1, a1 = mk and b1 = 2m may be written (Meijer 1941a, 1941b):


g(s) = \int_0^{\infty} (st)^{-k-1/2} \, e^{-st/2} \, W_{k+1/2,\,m}(st) \, f(t) \; dt,

f(t) = \frac {\Gamma(1-k+m)} {2 \pi i \, \Gamma(1+2m)} \int_{c - i \infty}^{c + i \infty} (ts)^{k-1/2} \, e^{ts/2} \, M_{k-1/2,\,m}(ts) \, g(s) \; ds.

Here Iν and Kν are the modified Bessel functions of the first and second kind, respectively, Mk,m and Wk,m are the Whittaker functions, and constant scale factors have been applied to the functions f and g and their arguments s and t in the first case.

Representation of other functions in terms of the G-function

The following list shows how the familiar elementary functions result as special cases of the Meijer G-function:

 H(1-|x|) = G_{1,1}^{\,1,0} \!\left( \left. \begin{matrix} 1 \\ 0 \end{matrix} \; \right| \, x \right), \qquad \forall x
 H(|x|-1) = G_{1,1}^{\,0,1} \!\left( \left. \begin{matrix} 1 \\ 0 \end{matrix} \; \right| \, x \right), \qquad \forall x
 e^x = G_{0,1}^{\,1,0} \!\left( \left. \begin{matrix} - \\ 0 \end{matrix} \; \right| \, -x \right), \qquad \forall x
 \cos x = \sqrt{\pi} \; G_{0,2}^{\,1,0} \!\left( \left. \begin{matrix} - \\ 0,\frac{1}{2} \end{matrix} \; \right| \, \frac{x^2}{4} \right), \qquad \forall x
 \sin x = \sqrt{\pi} \; G_{0,2}^{\,1,0} \!\left( \left. \begin{matrix} - \\ \frac{1}{2},0 \end{matrix} \; \right| \, \frac{x^2}{4} \right), \qquad \frac{-\pi}{2} < \arg x \leq \frac{\pi}{2}
 \cosh x = \sqrt{\pi} \; G_{0,2}^{\,1,0} \!\left( \left. \begin{matrix} - \\ 0,\frac{1}{2} \end{matrix} \; \right| \, -\frac{x^2}{4} \right), \qquad \forall x
 \sinh x = -\sqrt{\pi}i \; G_{0,2}^{\,1,0} \!\left( \left. \begin{matrix} - \\ \frac{1}{2},0 \end{matrix} \; \right| \, -\frac{x^2}{4} \right), \qquad -\pi < \arg x \leq 0
 \arcsin x = \frac{-i}{2\sqrt{\pi}} \; G_{2,2}^{\,1,2} \!\left( \left. \begin{matrix} 1,1 \\ \frac{1}{2},0 \end{matrix} \; \right| \, -x^2 \right), \qquad -\pi < \arg x \leq 0
 \arctan x = \frac{1}{2} \; G_{2,2}^{\,1,2} \!\left( \left. \begin{matrix} \frac{1}{2},1 \\ \frac{1}{2},0 \end{matrix} \; \right| \, x^2 \right), \qquad \frac{-\pi}{2} < \arg x \leq \frac{\pi}{2}
 \operatorname{arccot}\, x = \frac{1}{2} \; G_{2,2}^{\,2,1} \!\left( \left. \begin{matrix} \frac{1}{2},1 \\ \frac{1}{2},0 \end{matrix} \; \right| \, x^2 \right), \qquad \frac{-\pi}{2} < \arg x \leq \frac{\pi}{2}
 \ln (1+x) = G_{2,2}^{\,1,2} \!\left( \left. \begin{matrix} 1,1 \\ 1,0 \end{matrix} \; \right| \, x \right), \qquad \forall x

Here, H denotes the Heaviside step function.

The subsequent list shows how some higher functions can be expressed in terms of the G-function:

 \gamma (\alpha,x) = G_{1,2}^{\,1,1} \!\left( \left. \begin{matrix} 1 \\ \alpha,0 \end{matrix} \; \right| \, x \right), \qquad \forall x
 \Gamma (\alpha,x) = G_{1,2}^{\,2,0} \!\left( \left. \begin{matrix} 1 \\ \alpha,0 \end{matrix} \; \right| \, x \right), \qquad \forall x
 J_\nu (x) = G_{0,2}^{\,1,0} \!\left( \left. \begin{matrix} - \\ \frac{\nu}{2}, \frac{-\nu}{2} \end{matrix} \; \right| \, \frac{x^2}{4} \right), \qquad \frac{-\pi}{2} < \arg x \leq \frac{\pi}{2}
 Y_\nu (x) = G_{1,3}^{\,2,0} \!\left( \left. \begin{matrix} \frac{- \nu - 1}{2} \\ \frac{\nu}{2}, \frac{-\nu}{2}, \frac{- \nu - 1}{2} \end{matrix} \; \right| \, \frac{x^2}{4} \right), \qquad \frac{-\pi}{2} < \arg x \leq \frac{\pi}{2}
 I_\nu (x) = i^{-\nu} \; G_{0,2}^{\,1,0} \!\left( \left. \begin{matrix} - \\ \frac{\nu}{2}, \frac{-\nu}{2} \end{matrix} \; \right| \, -\frac{x^2}{4} \right), \qquad -\pi < \arg x \leq 0
 K_\nu (x) = \frac{1}{2} \; G_{0,2}^{\,2,0} \!\left( \left. \begin{matrix} - \\ \frac{\nu}{2}, \frac{-\nu}{2} \end{matrix} \; \right| \, \frac{x^2}{4} \right), \qquad \frac{-\pi}{2} < \arg x \leq \frac{\pi}{2}
 \Phi (x,n,a) = G_{n+1,\,n+1}^{\,1,\,n+1} \!\left( \left. \begin{matrix} 0, 1-a, \dots, 1-a \\ 0, -a, \dots, -a \end{matrix} \; \right| \, -x \right), \qquad \forall x, \; n = 0,1,2,\dots
 \Phi (x,-n,a) = G_{n+1,\,n+1}^{\,1,\,n+1} \!\left( \left. \begin{matrix} 0, -a, \dots, -a \\ 0, 1-a, \dots, 1-a \end{matrix} \; \right| \, -x \right), \qquad \forall x, \; n = 0,1,2,\dots

Even the derivatives of γ(α,x) and Γ(α,x) with respect to α can be expressed in terms of the Meijer G-function. Here, γ and Γ are the lower and upper incomplete gamma functions, Jν and Yν are the Bessel functions of the first and second kind, respectively, Iν and Kν are the corresponding modified Bessel functions, and Φ is the Lerch transcendent.

References

External links


Wikimedia Foundation. 2010.

Игры ⚽ Нужна курсовая?

Look at other dictionaries:

  • Meijer G-Function — The G function was defined for the first time by the Dutch mathematician Cornelis Simon Meijer (1904 1974) in 1936 as an attempt to introduce a very general function that includes most of the known special functions as particular cases. This was… …   Wikipedia

  • Meijer (disambiguation) — Meijer is an American hypermarket chain headquartered in Grand Rapids, Michigan. Meijer may also refer to: Meijer G function a mathematical function on the complex plane Meijer (surname), a Dutch surname This disambiguation page lists articles… …   Wikipedia

  • Cornelis Simon Meijer — (August 17, 1904, Pieterburen – April 12, 1974) was a Dutch mathematician at the university of Groningen who introduced the Meijer G function, a very general function that includes most of the elementary and higher mathematical functions as… …   Wikipedia

  • Lambert W function — The graph of W(x) for W > −4 and x < 6. The upper branch with W ≥ −1 is the function W0 (principal branch), the lower branch with W ≤ −1 is the function W−1. In mathematics, the Lambert W function, also called the Omega function or product… …   Wikipedia

  • Meier function — In mathematics, Meier function might refer to: Kaplan–Meier estimator Meijer G function This disambiguation page lists mathematics articles associated with the same title. If an internal link led you here, you may …   Wikipedia

  • Erik Meijer (computer scientist) — Erik Meijer is a Dutch computer scientist who is currently a software architect at Microsoft Research. He was previously a professor at Utrecht University. He received his Ph.D from Nijmegen University in 1992.Meijer s research has included the… …   Wikipedia

  • Weibull fading — The Weibull fading can be used as a simple statistical model of fading (named after Waloddi Weibull). In wireless communications, the Weibull fading distribution seems to exhibit good fit to experimental fading channel measurements for both… …   Wikipedia

  • Интеграл Меллина — Барнса (Mellin Barnes integral) или интеграл Барнса (Barnes integral) в математике контурный интеграл от функции, содержащей произведение гамма функций. Интегралы такого типа тесно связаны с обобщёнными гипергеометрическими функциями. Они были… …   Википедия

  • G-функция — название нескольких функций в математике: G функция Барнса (Barnes G function) связана с гамма функцией, расширяет понятие cуперфакториала на поле комплексных чисел. G функция Зигеля (Siegel G function) класс функций в теории трансцендентности. G …   Википедия

  • List of mathematical functions — In mathematics, several functions or groups of functions are important enough to deserve their own names. This is a listing of pointers to those articles which explain these functions in more detail. There is a large theory of special functions… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”