# .

In mathematics, a generalized hypergeometric series is a series in which the ratio of successive coefficients indexed by n is a rational function of n. The series, if convergent, defines a generalized hypergeometric function, which may then be defined over a wider domain of the argument by analytic continuation. The generalized hypergeometric series is sometimes just called the hypergeometric series, though this term also sometimes just refers to the Gaussian hypergeometric series. Generalized hypergeometric functions include the (Gaussian) hypergeometric function and the confluent hypergeometric function as special cases, which in turn have many particular special functions as special cases, such as elementary functions, Bessel functions, and the classical orthogonal polynomials.

Notation

A hypergeometric series is formally defined as a power series

$$\beta_0 + \beta_1 z + \beta_2 z^2 + \dots = \sum_{n \geqslant 0} \beta_n z^n \,\!$$

in which the ratio of successive coefficients is a rational function of n. That is,

$$\frac{\beta_{n+1}}{\beta_n} = \frac{A(n)}{B(n)}$$

where A(n) and B(n) are polynomials in n.

For example, in the case of the series for the exponential function,

$$1+\frac{z}{1!}+\frac{z^2}{2!}+\frac{z^3}{3!}+\dots,$$

βn = n!−1 and βn+1n = 1/(n+1). So this satisfies the definition with A(n) = 1 and B(n) = n+1.

It is customary to factor out the leading term, so β0 is assumed to be 1. The polynomials can be factored into linear factors of the form (aj+n) and (bk+n) respectively, where the aj and bk are complex numbers.

For historical reasons, it is assumed that (1+n) is a factor of B. If this is not already the case then both A and B can be multiplied by this factor; the factor cancels so the terms are unchanged and there is no loss of generality.

The ratio between consecutive coefficients now has the form

$$\frac{c(a_1+n)\dots(a_p+n)}{d(b_1+n)\dots(b_q+n)(1+n)},$$

where c and d are the leading coefficients of A and B. The series then has the form

$$1 + \frac{a_1\dots a_p}{b_1\dots b_q.1}\frac{cz}{d} + \frac{a_1\dots a_p}{b_1\dots b_q.1} \frac{(a_1+1)\dots(a_p+1)}{(b_1+1)\dots (b_q+1).2}\left(\frac{cz}{d}\right)^2+\dots,$$

or, by scaling z by the appropriate factor and rearranging,

$$1 + \frac{a_1\dots a_p}{b_1\dots b_q}\frac{z}{1!} + \frac{a_1(a_1+1)\dots a_p(a_p+1)}{b_1(b_1+1)\dots b_q(b_q+1)}\frac{z^2}{2!}+\dots.$$

This has the form of an exponential generating function. The standard notation for this series is

$$\,{}_pF_q(a_1,\ldots,a_p;b_1,\ldots,b_q;z) or \,{}_pF_q \left[\begin{matrix} a_1 & a_2 & \ldots & a_{p} \\ b_1 & b_2 & \ldots & b_q \end{matrix} ; z \right]$$

Using the rising factorial or Pochhammer symbol:

$$(a)_n=a(a+1)(a+2)...(a+n-1),\,(a)_0 = 1,$$

this can be written

$$\,{}_pF_q(a_1,\ldots,a_p;b_1,\ldots,b_q;z) = \sum_{n=0}^\infty \frac{(a_1)_n\dots(a_p)_n}{(b_1)_n\dots(b_q)_n} \, \frac {z^n} {n!}$$

(Note that this use of the Pochhammer symbol is not standard, however it is the standard usage in this context.)
Special cases

Many of the special functions in mathematics are special cases of the confluent hypergeometric function or the hypergeometric function; see the corresponding articles for examples.

Some of the functions related to more complicated hypergeometric functions include:

Dilogarithm:

$$\operatorname{Li}_2(x) = \sum_{n>0}\,{x^n}{n^{-2}} = x \; {}_3F_2(1,1,1;2,2;x)$$

Hahn polynomials:

$$Q_n(x;a,b,N)= {}_3F_2(-n,-x,n+a+b+1;a+1,-N+1;1).\$$

Wilson polynomials:

$$p_n(t^2)=(a+b)_n(a+c)_n(a+d)_n \; {}_4F_3\left( \begin{matrix} -n&a+b+c+d+n-1&a-t&a+t \\ a+b&a+c&a+d \end{matrix} ;1\right).$$

Terminology

When all the terms of the series are defined and it has a non-zero radius of convergence, then the series defines an analytic function. Such a function, and its analytic continuations, is called the hypergeometric function.

The case when the radius of convergence is 0 yields many interesting series in mathematics, for example the incomplete gamma function has the asymptotic expansion

$$\Gamma(a,z) \sim z^{a-1}e^{-z}\left(1+\frac{a-1}{z}+\frac{(a-1)(a-2)}{z^2}\dots\right)$$

which could be written za−1e−z 2F0(1−a,1;;−z−1). However, the use of the term hypergeometric series is usually restricted to the case where the series defines an actual analytic function.

The ordinary hypergeometric series should not be confused with the basic hypergeometric series, which, despite its name, is a rather more complicated and recondite series. The "basic" series is the q-analog of the ordinary hypergeometric series. There are several such generalizations of the ordinary hypergeometric series, including the ones coming from zonal spherical functions on Riemannian symmetric spaces.

The series without the factor of n! in the denominator (summed over all integers n, including negative) is called the bilateral hypergeometric series.
Convergence conditions

There are certain values of the aj and bk for which the numerator or the denominator of the coefficients is 0.

If any aj is a non-positive integer (0, −1, −2, etc.) then the series only has a finite number of terms and is, in fact, a polynomial of degree −aj.
If any bk is a non-positive integer (excepting the previous case with −bk < aj) then the denominators become 0 and the series is undefined.

Excluding these cases, the ratio test can be applied to determine the radius of convergence.

If p ≤ q then the ratio of coefficients tends to zero. This implies that the series converges for any finite value of z. An example is the power series for the exponential function.
If p = q+1 then the ratio of coefficients tends to one. This implies that the series converges for |z|<1 and diverges for |z|>1. Analytic continuation can be employed for larger values of z.
If p > q+1 then the ratio of coefficients grows without bound. This implies that, besides z=0, the series diverges. This is then a divergent or asymptotic series, or it can be interpreted as a symbolic shorthand for a differential equation that the sum satisfies.

The question of convergence for p=q+1 when z is on the unit circle is more difficult. It can be shown that the series converges absolutely at z=1 if

$$\Re\left(\sum b_k - \sum a_j\right)>0.$$

Basic properties

It is immediate from the definition that the order of the parameters a_j, or the order of the parameters b_k can be changed without changing the value of the function. Also, if any of the parameters a_j is equal to any of the parameters b_k, then the matching parameters can be "cancelled out", with certain exceptions when the parameters are non-positive integers. For example,

$$\,{}_2F_1(3,1;1;z) = \,{}_2F_1(1,3;1;z) = \,{}_1F_0(3;;z).$$

Euler's integral transform

The following basic identity is very useful as it relates the higher-order hypergeometric functions in terms of integrals over the lower order ones [1]:

$$_{A+1}F_{B+1}\left[ \begin{array}{c} a_{1},\ldots ,a_{A},c \\ b_{1},\ldots ,b_{B},d \end{array} ;z\right] =\frac{\Gamma (d)}{\Gamma (c)\Gamma (d-c)} \int_{0}^{1}t^{c-1}(1-t)_{{}}^{d-c-1}\ {}_{A}F_{B}\left[ \begin{array}{c} a_{1},\ldots ,a_{A} \\ b_{1},\ldots ,b_{B} \end{array} ;tz\right] dt$$

Differentiation

The generalized hypergeometric function satisfies

$$(z\frac{{\rm{d}}}{{\rm{d}}z} + a_j) \; {}_pF_q(a_1,\dots,a_j,\dots,a_p;b_1,\dots,b_q;z) = a_j \; {}_pF_q(a_1,\dots,a_j+1,\dots,a_p;b_1,\dots,b_q;z).$$

and

$$(z\frac{{\rm{d}}}{{\rm{d}}z} + b_k - 1) \; {}_pF_q(a_1,\dots,a_p;b_1,\dots,b_k,\dots,b_q;z) = (b_k - 1) \; {}_pF_q(a_1,\dots,a_p;b_1,\dots,b_k-1,\dots,b_q;z).$$

Combining these gives a differential equation satisfied by w = pFq:

$$z\prod_{n=1}^{p}(z\frac{{\rm{d}}}{{\rm{d}}z} + a_n)w = z\frac{{\rm{d}}}{{\rm{d}}z}\prod_{n=1}^{q}(z\frac{{\rm{d}}}{{\rm{d}}z} + b_n-1)w.$$

Contiguous function and related identities

Let \vartheta be the operator $$z\frac{d}{dz}$$. From the differentiation formulas given above, the linear space spanned by $$\; {}_pF_q(a_1,\dots,a_p;b_1,\dots,b_q;z)$$ and $$\vartheta\; {}_pF_q(a_1,\dots,a_p;b_1,\dots,b_q;z)$$contains each of

$$\; {}_pF_q(a_1,\dots,a_j+1,\dots,a_p;b_1,\dots,b_q;z), \; {}_pF_q(a_1,\dots,a_p;b_1,\dots,b_k-1,\dots,b_q;z),$$
$$z\; {}_pF_q(a_1+1,\dots,a_p+1;b_1+1,\dots,b_q+1;z), \; {}_pF_q(a_1,\dots,a_p;b_1,\dots,b_q;z).$$

Since the space has dimension 2, any three of these p+q+2 functions are linearly dependent. These dependencies can be written out to generate a large number of identities involving $${}_pF_q$$.

For example, in the simplest non-trivial case,

$$\; {}_0F_1(;a;z) = (1) \; {}_0F_1(;a;z),$$
$$\; {}_0F_1(;a-1;z) = (\frac{\vartheta}{a-1}+1) \; {}_0F_1(;a;z),$$
$$z \; {}_0F_1(;a+1;z) = (a\vartheta) \; {}_0F_1(;a;z),$$

So

$$\; {}_0F_1(;a-1;z)- \; {}_0F_1(;a;z) = \frac{z}{a(a-1)} \; {}_0F_1(;a+1;z).$$

This, and other important examples,

$$\; {}_1F_1(a+1;b;z)- \, {}_1F_1(a;b;z) = \frac{z}{b} \; {}_1F_1(a+1;b+1;z),$$
$$\; {}_1F_1(a;b-1;z)- \, {}_1F_1(a;b;z) = \frac{az}{b(b-1)} \; {}_1F_1(a+1;b+1;z),$$
$$\; {}_1F_1(a;b-1;z)- \, {}_1F_1(a+1;b;z) = \frac{(a-b+1)z}{b(b-1)} \; {}_1F_1(a+1;b+1;z)$$

$$\; {}_2F_1(a+1,b;c;z)- \, {}_2F_1(a,b;c;z) = \frac{bz}{c} \; {}_2F_1(a+1,b+1;c+1;z),$$
$$\; {}_2F_1(a+1,b;c;z)- \, {}_2F_1(a,b+1;c;z) = \frac{(b-a)z}{c} \; {}_2F_1(a+1,b+1;c+1;z),$$
$$\; {}_2F_1(a,b;c-1;z)- \, {}_2F_1(a+1,b;c;z) = \frac{(a-c+1)bz}{c(c-1)} \; {}_2F_1(a+1,b+1;c+1;z),$$

can be used to generate continued fraction expressions known as Gauss's continued fraction.

Similarly, by applying the differentiation formulas twice, there are $$\binom{p+q+3}{2}$$ such functions contained in $$\{1, \vartheta, \vartheta^2\}\; {}_pF_q(a_1,\dots,a_p;b_1,\dots,b_q;z)$$, which has dimension three so any four are linearly dependent. This generates more identities and the process can be continued. The identities thus generated can be combined with each other to produce new ones in a different way.

A function obtained by adding \pm 1 to exactly one of the parameters \)$$a_j, b_k in \;$$ $${}_pF_q(a_1,\dots,a_p;b_1,\dots,b_q;z)$$ is called contiguous to \; $${}_pF_q(a_1,\dots,a_p;b_1,\dots,b_q;z)$$. Using the technique outlined above, an identity relating \; $${}_0F_1(;a;z)$$ and its two contiguous functions can be given, six identities relating $$\; {}_1F_1(a;b;z)$$ and any two of its four contiguous functions, and fifteen identities relating $$\; {}_2F_1(a,b;c;z)$$and any two of its six contiguous functions have been found. (The first one was derived in the previous paragraph. The last fifteen were given by Gauss in his 1812 paper.)
Identities
For identities involving the Gauss hypergeometric function 2F1, see hypergeometric function.

A number of other hypergeometric function identities were discovered in the nineteenth and twentieth centuries.
Saalschütz's theorem

Saalschütz's theorem[2] (Saalschütz 1890) is

$${}_3F_2 (a,b, -n;c, 1+a+b-c-n;1)= \frac{(c-a)_n(c-b)_n}{(c)_n(c-a-b)_n}$$

for n a positive integer, where

$$\;(x)_n = x(x+1) \cdots (x+n-1), \, \, (x)_0 = 1$$

is the Pochhammer symbol.
Dixon's identity
Main article: Dixon's identity

Dixon's identity,[3] first proved by Dixon (1902), gives the sum of a well-poised 3F2 at 1:

$${}_3F_2 (a,b,c;1+a-b,1+a-c;1)= \frac{\Gamma(1+\frac{a}{2})\Gamma(1+\frac{a}{2}-b-c)\Gamma(1+a-b)\Gamma(1+a-c)} {\Gamma(1+a)\Gamma(1+a-b-c)\Gamma(1+\frac{a}{2}-b)\Gamma(1+\frac{a}{2}-c)}.$$

Dougall's formula

Dougall's formula (Dougall 1907) gives the sum of a terminating well-poised series:

$${}_7F_6\left(\begin{matrix}a&1+\frac{a}{2}&b&c&d&e&-m\\&\frac{a}{2}&1+a-b&1+a-c&1+a-d&1+a-e&1+a+m\\ \end{matrix};1\right)$$

$$=\frac{(1+a)_m(1+a-b-c)_m(1+a-c-d)_m(1+a-b-d)_m}{(1+a-b)_m(1+a-c)_m(1+a-d)_m(1+a-b-c-d)_m}$$

provided that m is a non-negative integer (so that the series terminates) and

$$1+2a=b+c+d+e-m\,.$$

Many of the other formulas for special values of hypergeometric functions can be derived from this as special or limiting cases.
Generalization of Kummer's transformations and identities for _2F_2

$$e^{-x} \; {}_2F_2(a,1+d;c,d;x)= \; {}_2F_2(c-a-1,f+1;c,f;-x)$$

where $$f=\frac{d(a-c+1)}{a-d};$$

$$e^{-\frac x 2} \, {}_2F_2 \left(a, 1+b; 2a+1, b; x\right)= \, {}_0F_1 \left(;a+\frac 1 2; \frac {x^2} {16}\right) - \frac{x\left(1-\frac{2a} b\right)}{2(2a+1)}\; {}_0F_1 \left(;a+\frac 3 2; \frac {x^2} {16}\right),$$

which links Bessel functions to {}_2F_2; this reduces to Kummer's second formula for b= 2 a:

$$e^{-\frac x 2} \, {}_1F_1(a,2a,x)= \, {}_0F_1 (;a+\frac 1 2; \frac{x^2}{16}).$$

\begin{align}_2F_2(a,b;c,d;x)=& \sum_{i=0} \frac{{b-d \choose i}{a+i-1 \choose i}}{{c+i-1 \choose i}{d+i-1 \choose i}} \; {}_1F_1(a+i;c+i;x)\frac{x^i}{i!} \\ =& e^x \sum_{i=0} \frac{{b-d \choose i}{a+i-1 \choose i}}{{c+i-1 \choose i}{d+i-1 \choose i}} \; {}_1F_1(c-a;c+i;-x)\frac{x^i}{i!},\end{align}

which is a finite sum if b-d is a non-negative integer.
Kummer's relation

Kummer's relation is

$$_2F_1\left(2a,2b;a+b+\tfrac 1 2;x\right)=\,_2F_1\left(a,b; a+b+\tfrac 1 2; 4x(1-x)\right).$$

Clausen's formula
Main article: Clausen's formula

Clausen's formula

$${}_3F_2(2c-2s-1, 2s, c-\frac 1 2; 2c-1, c; x)=\, {}_2F_1(c-s-\frac 1 2,s; c; x)^2$$

was used by de Branges to prove the Bieberbach conjecture.
Special cases
The series 0F0
Main article: Exponential function

As noted earlier, $$\,{}_0F_0(;;z) = e^z$$. The differential equation for this function is $$\frac{d}{dz}w = w$$, which has solutions $$w = ke^z$$ where k is a constant.
The series 1F0

Also as noted earlier, $$\,{}_1F_0(a;;z) = (1-z)^{-a$$}. The differential equation for this function is $$\frac{d}{dz}w = (z\frac{d}{dz}+a)w, or (1-z)\frac{dw}{dz} = aw$$, which has solutions $$w=k(1-z)^{-a}$$where k is a constant.
The series 0F1

The functions of the form \ $$,{}_0F_1(;a;z)$$ are called confluent hypergeometric limit functions and are closely related to Bessel functions. The relationship is: $$J_\alpha(x)=\frac{(x/2)^\alpha}{\Gamma(\alpha+1)} \;_0F_1 (;\alpha+1; -\tfrac{1}{4}x^2$$). The differential equation for this function is $$w = (z\frac{d}{dz}+a)\frac{dw}{dz} or z \( \frac{d^2w}{dz^2}+a\frac{dw}{dz}-w = 0$$. When a is not a positive integer, the substitution $$w = z^{1-a}u$$, gives a linearly independent solution $$z^{1-a}\;{}_0F_1(;2-a;z)$$, so the general solution is $$k\;{}_0F_1(;a;z)+l z^{1-a}\;{}_0F_1(;2-a;z)$$ where k, l are constants. (If a is a positive integer, the independent solution is given by the appropriate Bessel function of the second kind.)
The series 1F1
Main article: Confluent hypergeometric function

The functions of the form $$\,{}_1F_1(a;b;z)$$ are called confluent hypergeometric functions of the first kind, also written $$M(a;b;z)$$. The incomplete gamma function $$\gamma(a,z)$$ is a special case.

The differential equation for this function is ( $$z\frac{d}{dz}+a)w = (z\frac{d}{dz}+b)\frac{dw}{dz}$$ or $$z\frac{d^2w}{dz^2}+(b-z)\frac{dw}{dz}-aw = 0$$. When b is not a positive integer, the substitution $$w = z^{1-b}u$$, gives a linearly independent solution $$z^{1-b}\;{}_1F_1(1+a-b;2-b;z)$$, so the general solution is $$k\;{}_1F_1(a;b;z)+l z^{1-b}\;{}_1F_1(1+a-b;2-b;z)$$where k, l are constants.

When a is a non-positive integer, $$−n, \,{}_1F_1(-n;b;z)$$ is a polynomial. Up to constant factors, these are the Laguerre polynomials. This implies Hermite polynomials can be expressed in terms of $${}_1F_1$$ as well.
The series 2F0

This occurs in connection with the exponential integral function Ei(z).
The series 2F1
Main article: Hypergeometric function

Historically, the most important are the functions of the form $$\,{}_2F_1(a,b;c;z)$$. These are sometimes called Gauss' hypergeometric functions, classical standard hypergeometric or often simply hypergeometric functions. The term Generalized hypergeometric function is used for the functions ( \,{}_pF_q \)if there is risk of confusion. This function was first studied in detail by Carl Friedrich Gauss, who explored the conditions for its convergence.

The differential equation for this function is

$$(z\frac{d}{dz}+a)(z\frac{d}{dz}+b)w = (z\frac{d}{dz}+c)\frac{dw}{dz}$$

or

$$z(1-z)\frac{d^2w}{dz^2}+(c-(a+b+1)z)\frac{dw}{dz}-abw = 0.$$

It is known as the hypergeometric differential equation. When c is not a positive integer, the substitution

$$w = z^{1-c}u$$

gives a linearly independent solution $$z^{1-c}\; {}_2F_1(1+a-c,1+b-c;2-c;z)$$, so the general solution for |z|<1 is

$$k\; {}_2F_1(a,b;c;z)+l z^{1-c}\; {}_2F_1(1+a-c,1+b-c;2-c;z)$$

where k, l are constants. Different solutions can be derived for other values of z. In fact there are 24 solutions, known as the Kummer solutions, derivable using various identities, valid in different regions of the complex plane.

When a is a non-positive integer, $$−n, \,{}_2F_1(-n,b;c;z)$$ is a polynomial. Up to constant factors and scaling, these are the Jacobi polynomials. Several other classes of orthogonal polynomials, up to constant factors, are special cases of Jacobi polynomials, so these can be expressed using $$_2F_1$$as well. This includes Legendre polynomials and Chebyshev polynomials.

A wide range of integrals of elementary functions can be expressed using the hypergeometric function, e.g.:

$$\int_0^x\sqrt{1+y^\alpha}\,\mathrm{d}y\ =\ \left(2+\alpha\right)^{-1}x\left(\alpha\;{}_2F_1\left(\frac{1}{\alpha},\frac{1}{2};1+\frac{1}{\alpha};-x^\alpha\right)+2\sqrt{x^\alpha+1}\,\right)\ ,\ \ \ \alpha\neq0.$$

The series 3F1

This occurs in the theory of Bessel functions. It provides a way to compute Bessel functions of large arguments.

Generalizations

The generalized hypergeometric function is linked to the Meijer G-function and the MacRobert E-function. Hypergeometric series were generalised to several variables, for example by Paul Emile Appell; but a comparable general theory took long to emerge. Many identities were found, some quite remarkable. A generalization, the q-series analogues, called the basic hypergeometric series, were given by Eduard Heine in the late nineteenth century. Here, the ratios considered of successive terms, instead of a rational function of n, are a rational function of $$q^n$$. Another generalization, the elliptic hypergeometric series, are those series where the ratio of terms is an elliptic function (a doubly periodic meromorphic function) of n.

During the twentieth century this was a fruitful area of combinatorial mathematics, with numerous connections to other fields. There are a number of new definitions of general hypergeometric functions, by Aomoto, Israel Gelfand and others; and applications for example to the combinatorics of arranging a number of hyperplanes in complex N-space (see arrangement of hyperplanes).

Special hypergeometric functions occur as zonal spherical functions on Riemannian symmetric spaces and semi-simple Lie groups. Their importance and role can be understood through the following example: the hypergeometric series 2F1 has the Legendre polynomials as a special case, and when considered in the form of spherical harmonics, these polynomials reflect, in a certain sense, the symmetry properties of the two-sphere or, equivalently, the rotations given by the Lie group SO(3). In tensor product decompositions of concrete representations of this group Clebsch-Gordan coefficients are met, which can be written as 3F2 hypergeometric series.

Bilateral hypergeometric series are a generalization of hypergeometric functions where one sums over all integers, not just the positive ones.

Fox–Wright functions are a generalization of generalized hypergeometric functions where the Pochhammer symbols in the series expression are generalised to gamma functions of linear expressions in the index n.
Citations

^ (Slater 1966, Equation (4.1.2))
^ See (Slater 1966, Section 2.3.1) or (Bailey 1935, Section 2.2) for a proof.
^ See (Bailey 1935, Section 3.1) for a detailed proof. An alternative proof is in (Slater 1966, Section 2.3.3)

References

Askey, R. A.; Daalhuis, Adri B. Olde (2010), "Generalized hypergeometric function", in Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F. et al., NIST Handbook of Mathematical Functions, Cambridge University Press, ISBN 978-0521192255, MR2723248
Andrews, George E.; Askey, Richard & Roy, Ranjan (1999). Special functions. Encyclopedia of Mathematics and its Applications. 71. Cambridge University Press. ISBN 978-0-521-62321-6; 978-0-521-78988-2. MR1688958.
Bailey, W.N. (1935). Generalized Hypergeometric Series. Cambridge Tracts in Mathematics and Mathematical Physics. 32. London: Cambridge University Press. Zbl 0011.02303.
Dixon, A.C. (1902). "Summation of a certain series". Proc. London Math. Soc. 35 (1): 284–291. doi:10.1112/plms/s1-35.1.284. JFM 34.0490.02.
Dougall, J. (1907). "On Vandermonde's theorem and some more general expansions". Proc. Edinburgh Math. Soc. 25: 114–132. doi:10.1017/S0013091500033642.
Gasper, George; Rahman, Mizan (2004). Basic Hypergeometric Series. Encyclopedia of Mathematics and Its Applications. 96 (2nd ed.). Cambridge, UK: Cambridge University Press. ISBN 0-521-83357-4. MR2128719. Zbl 1129.33005. First edition with ISBN 0-521-35049-2.
Gauss, Carl Friedrich (1813). "Disquisitiones generales circa seriam infinitam" (in Latin). Commentationes societatis regiae scientarum Gottingensis recentiores (Göttingen) 2. (Gauss's original paper can be found in Carl Friedrich Gauss Werke, p. 125)
Heckman, Gerrit & Schlichtkrull, Henrik (1994). Harmonic Analysis and Special Functions on Symmetric Spaces. San Diego: Academic Press. ISBN 0-12-336170-2. (Part 1 treats hypergeometric functions on Lie groups.)
Saalschütz, L. (1890). "Eine Summationsformel" (in German). Zeitschrift für Mathematik und Physik 35: 186–188. JFM 22.0262.03.
Slater, Lucy Joan (1966). Generalized Hypergeometric Functions. Cambridge, UK: Cambridge University Press. ISBN 0-521-06483-X. MR0201688. Zbl 0135.28101. (there is a 2008 paperback with ISBN 978-0-521-09061-2)
Yoshida, Masaaki (1997). Hypergeometric Functions, My Love: Modular Interpretations of Configuration Spaces. Braunschweig/Wiesbaden: Friedr. Vieweg & Sohn. ISBN 3-528-06925-2. MR1453580.
Rathie, Arjun K.; Pogány, Tibor K. (2008). "New summation formula for 3F2(1/2) and a Kummer-type II transformation of 2F2(x)". Mathematical Communications 13: 63–66. MR2422088. Zbl 1146.33002.