# .

In mathematics, the Gaussian or ordinary hypergeometric function 2F1(a,b;c;z) is a special function represented by the hypergeometric series, that includes many other special functions as specific or limiting cases. It is a solution of a second-order linear ordinary differential equation (ODE). Every second-order linear ODE with three regular singular points can be transformed into this equation.

For systematic lists of some of the many thousands of published identities involving the hypergeometric function, see the reference works by Arthur Erdélyi, Wilhelm Magnus, and Fritz Oberhettinger et al. (1953), Abramowitz & Stegun (1965), and Daalhuis (2010).

History

The term "hypergeometric series" was first used by John Wallis in his 1655 book Arithmetica Infinitorum. Hypergeometric series were studied by Euler, but the first full systematic treatment was given by Gauss (1813), Studies in the nineteenth century included those of Ernst Kummer (1836), and the fundamental characterisation by Bernhard Riemann of the hypergeometric function by means of the differential equation it satisfies. Riemann showed that the second-order differential equation (in z) for the 2F1, examined in the complex plane, could be characterised (on the Riemann sphere) by its three regular singularities.

The cases where the solutions are algebraic functions were found by H. A. Schwarz (Schwarz's list).
The hypergeometric series

The hypergeometric function is defined for |z| < 1 by the series

$$\,_2F_1(a,b;c;z) = \sum_{n=0}^\infty \frac{(a)_n(b)_n}{(c)_n} \, \frac {z^n} {n!}$$

provided that c is not 0, −1, −2, … Notice that the series terminates if either "a" or "b" is a negative integer. The Pochhammer symbol is defined by

$$(x)_n = \left\{ \begin{array}{ll} 1 & \mbox{if } n = 0 \\ x(x+1) \cdots (x+n-1) & \mbox{if } n > 0 \end{array} \right.$$

For other complex values of z it can be analytically continued along any path in the complex plane that avoids the branch points 0 and 1.
Special cases

Many of the common mathematical functions can be expressed in terms of the hypergeometric function, or as limiting cases of it. Some typical examples are

$$\ln(1+z)=z\,_2F_1(1,1;2;-z).$$
$$(1-z)^{-a} = \,_2F_1(a,b;b;z)$$
$$\arcsin z = z \,_2F_1\left(\tfrac{1}{2}, \tfrac{1}{2}; \tfrac{3}{2};z^2\right)$$

The confluent hypergeometric function (or Kummer's function) can be given as a limit of the hypergeometric function

$$M(a,c,z) = \lim_{b\rightarrow\infty}{}_2F_1(a,b;c;b^{-1}z)$$

so all functions that are essentially special cases of it, such as Bessel functions, can be expressed as limits of hypergeometric functions. These include most of the commonly used functions of mathematical physics.

Legendre functions are solutions of a second order differential equation with 3 regular singular points so can be expressed in terms of the hypergeometric function in many ways, for example

$${}_2F_1(a,1-a;c;z) = \Gamma(c)z^{\tfrac{1-c}{2}}(1-z)^{\tfrac{c-1}{2}}P_{-a}^{1-c}(1-2z)$$

Several orthogonal polynomials, including Jacobi polynomials P(α,β)
n and their special cases Legendre polynomials, Chebyshev polynomials, Gegenbauer polynomials can be written in terms of hypergeometric functions using

$${}_2F_1(-n,\alpha+1+\beta+n;\alpha+1;x) = \frac{n!}{(\alpha+1)_n}P^{(\alpha,\beta)}_n(1-2x)$$

Other polynomials that are special cases include Krawtchouk polynomials, Meixner polynomials, Meixner–Pollaczek polynomials.

Elliptic modular functions can sometimes be expressed as the inverse functions of ratios of hypergeometric functions whose arguments a, b, c are 1, 1/2, 1/3, ... or 0. For examples, if

$$\tau = {\rm{i}}\frac{{}_2F_1(\frac{1}{2},\frac{1}{2};1;1-z)}{{}_2F_1(\frac{1}{2},\frac{1}{2};1;z)}$$

then

$$z = \kappa^2(\tau) = \frac{\theta_2(\tau)^4}{\theta_3(\tau)^4}$$

is an elliptic modular function of τ.

Incomplete beta functions Bx(p,q) are related by

$$B_x(p,q) = \frac{x^p}{p}{}_2F_1(p,1-q;p+1;x)$$

The complete elliptic integrals K and E are given by

$$K(k) = \tfrac{\pi}{2}\, _2F_1\left(\tfrac{1}{2},\tfrac{1}{2};1;k^2\right)$$

$$E(k) = \tfrac{\pi}{2}\, _2F_1\left(-\tfrac{1}{2},\tfrac{1}{2};1;k^2\right)$$

The hypergeometric differential equation

The hypergeometric function is a solution of Euler's hypergeometric differential equation

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

which has three regular singular points: 0,1 and ∞. The generalization of this equation to three arbitrary regular singular points is given by Riemann's differential equation. Any second order differential equation with three regular singular points can be converted to the hypergeometric differential equation by a change of variables.
Solutions at the singular points

Solutions to the hypergeometric differential equation are built out of the hypergeometric series $$\;_2F_1(a,b;c;z).$$ . The equation has two linearly independent solutions. At each of the three singular points 0, 1, ∞, there are usually two special solutions of the form xs times a holomorphic function of x, where s is one of the two roots of the indicial equation and x is a local variable vanishing at the regular singular point. This gives 3 × 2 = 6 special solutions, as follows.

Around the point z = 0, two independent solutions are, if c is not an integer,

$$\;_2F_1(a,b;c;z)$$

and

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

Around z = 1, if c − a − b is not an integer, one has two independent solutions

$$\;_2F_1(a,b;1+a+b-c;1-z)$$

and

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

Around z = ∞, if a − b is not an integer, one has two independent solutions

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

and

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

Any 3 of these 6 solutions satisfy a linear relation as the space of solutions is 2-dimensional, giving (6
3) = 20 linear relations between them called connection formulas.
Kummer's 24 solutions

A second order Fuchsian equation with n singular points has a group of symmetries acting (projectively) on its solutions, isomorphic to the Coxeter group Dn of order 2n-1×n!. For the hypergeometric equation n=3, so the group is of order 24 and is isomorphic to the symmetric group on 4 points, and was first described by Kummer. The isomorphism with the symmetric group is accidental and has no analogue for more than 3 singular points, and it is sometimes better to think of the group as an extension of the symmetric group on 3 points (acting as permutations of the 3 singular points) by a Klein 4-group (whose elements change the signs of the differences of the exponents at an even number of singular points). Kummer's group of 24 transformations is generated by the three transformations taking a solution F(a, b; c; z) to one of

$$\,(1-z)^{-a} F(a,c-b;c;z/(z-1))$$
$$\,F(a,b;1+a+b-c;1-z)$$
$$(\,1-z)^{-b} F(c-a,b;c;z/(z-1))$$

which correspond to the transpositions (12), (23), and (34) under an isomorphism with the symmetric group on 4 points 1, 2, 3, 4.

Applying Kummer's 24=6×4 transformations to the hypergeometric function gives the 6 = 2×3 solutions above corresponding to each of the 2 possible exponents at each of the 3 singular points, each of which appears 4 times because of the identities

$$\;_2F_1(a,b;c;z)= (1-z)^{c-a-b} \;_2F_1(c-a,c-b;c;z)$$ (Euler transformation)
$$\;_2F_1(a,b;c;z)=(1-z)^{-a} \;_2F_1(a,c-b;c;z/(z-1))$$ (Pfaff transformation)
$$\;_2F_1(a,b;c;z)=(1-z)^{-b} \;_2F_1(c-a,b;c;z/(z-1))$$ (Pfaff transformation)

Q-form

The hypergeometric differential equation may be brought into the Q-form

$$\frac{d^2u}{dz^2}+Q(z)u(z) = 0$$

by making the substitution w = uv and eliminating the first-derivative term. One finds that

$$Q=\frac{z^2[1-(a-b)^2] +z[2c(a+b-1)-4ab] +c(2-c)}{4z^2(1-z)^2}$$

and v is given by the solution to

$$\frac{d}{dz}\log v(z) = - \frac {c-z(a+b+1)}{2z(1-z)}.$$

The Q-form is significant in its relation to the Schwarzian derivative.
Schwarz triangle maps

The Schwarz triangle maps or Schwarz s-functions are ratios of pairs of solutions.

$$s_k(z) = \frac{\phi_k^{(1)}(z)}{\phi_k^{(0)}(z)}$$

where k is one of the points 0, 1, ∞. The notation

$$\,D_k(\lambda,\mu,\nu;z)=s_k(z)$$

is also sometimes used. Note that the connection coefficients become Möbius transformations on the triangle maps.

Note that each triangle map is regular at z ∈ {0, 1, ∞} respectively, with

$$s_0(z)=z^\lambda (1+\mathcal{O}(z))$$
$$s_1(z)=(1-z)^\mu (1+\mathcal{O}(1-z))$$

and

$$s_\infty(z)=z^\nu (1+\mathcal{O}(1/z)). In the special case of λ, μ and ν real, with \( 0\le|\lambda|,|\mu|,|\nu|<1$$ then the s-maps are conformal maps of the upper half-plane H to triangles on the Riemann sphere, bounded by circular arcs. This mapping is a special case of a Schwarz–Christoffel mapping. The singular points 0,1 and ∞ are sent to the triangle vertices. The angles of the triangle are πλ, πμ and πν respectively.

Furthermore, in the case of λ=1/p, μ=1/q and ν=1/r for integers p, q, r, then the triangle tiles the sphere, and the s-maps are inverse functions of automorphic functions for the triangle group \langle p,q,r\rangle=\Delta (p,q,r).
Monodromy group

The monodromy of a hypergeometric equation describes how fundamental solutions change when analytically continued around paths in the z plane that return to the same point. That is, when the path winds around a singularity of $$\;_2F_1$$ , the value of the solutions at the endpoint will differ from the starting point.

Two fundamental solutions of the hypergeometric equation are related to each other by a linear transformation; thus the monodromy is a mapping (group homomorphism):

$$\pi_1(z_0,\mathbb{C}\setminus\{0,1\}) \to GL(2,\mathbb{C})$$

where \pi_1 is the fundamental group. In other words the monodromy is a two dimensional linear representation of the fundamental group. The monodromy group of the equation is the image of this map, i.e. the group generated by the monodromy matrices.
Integral formulas
Euler type

If B is the beta function then

$$\Beta(b,c-b)\,_2F_1(a,b;c;z) = \int_0^1 x^{b-1} (1-x)^{c-b-1}(1-zx)^{-a} \, dx \quad \real(c) > \real(b) > 0,$$

provided |z| < 1 or |z| = 1 and both sides converge, and can be proved by expanding (1 − zx)−a using the binomial theorem and then integrating term by term. This was given by Euler in 1748 and implies Euler's and Pfaff's hypergeometric transformations.

Other representations, corresponding to other branches, are given by taking the same integrand, but taking the path of integration to be a closed Pochhammer cycle enclosing the singularities in various orders. Such paths correspond to the monodromy action.
Barnes integral

Barnes used the theory of residues to evaluate the Barnes integral

$$\frac{1}{2\pi i}\int_{-i\infty}^{i\infty} \frac{\Gamma(a+s)\Gamma(b+s)\Gamma(-s)}{\Gamma(c+s)} (-z)^s \, ds$$

as

$$\frac{\Gamma(a)\Gamma(b)}{\Gamma(c)}\,_2F_1(a,b;c;z),$$

where the contour is drawn to separate the poles 0, 1, 2... from the poles −a, −a - 1, ..., −b, −b − 1, ... .
John transform

The Gauss hypergeometric function can be written as a John transform (Gelfand, Gindikin & Graev 2003, 2.1.2).
Gauss' contiguous relations

The six functions $$\displaystyle{}_2F_1 (a\pm 1,b;c;z), \displaystyle{}_2F_1 (a,b\pm 1;c;z)$$ , and $$\displaystyle{}_2F_1 (a,b;c\pm 1;z)$$ are called contiguous to $$\displaystyle{}_2F_1 (a,b;c;z)$$ . Gauss showed that \displaystyle{}_2F_1 (a,b;c;z) can be written as a linear combination of any two of its contiguous functions, with rational coefficients in terms of a,b,c, and z. This gives (6
2)=15 relations, given by identifying any two lines on the right hand side of

\begin{align} z\frac{dF}{dz} = z\frac{ab}{c}F(a+,b+,c+) &=a(F(a+)-F)\\ &=b(F(b+)-F)\\ &=(c-1)(F(c-)-F)\\ &=\frac{(c-a)F(a-)+(a-c+bz)F}{1-z}\\ &=\frac{(c-b)F(b-)+(b-c+az)F}{1-z}\\ &=z\frac{(c-a)(c-b)F(c+)+c(a+b-c)F}{c(1-z)}\\ \end{align}

In the notation above, $$\displaystyle{}F= {}_2F_1(a,b;c;z), F(a+)={}_2F_1(a+1,b;c;z)$$ and so on.

Repeatedly applying these relations gives a linear relation over C(z) between any three functions of the form $$\displaystyle{}_2F_1 (a+m,b+n;c+l;z)$$ , where m, n, and l are integers.
Gauss's continued fraction
Main article: Gauss continued fraction

Gauss used the contiguous relations to give several ways to write a quotient of two hypergeometric functions as a continued fraction, for example:

$$\frac{{}_2F_1(a+1,b;c+1;z)}{{}_2F_1(a,b;c;z)} = \cfrac{1}{1 + \cfrac{\frac{(a-c)b}{c(c+1)} z}{1 + \cfrac{\frac{(b-c-1)(a+1)}{(c+1)(c+2)} z}{1 + \cfrac{\frac{(a-c-1)(b+1)}{(c+2)(c+3)} z}{1 + \cfrac{\frac{(b-c-2)(a+2)}{(c+3)(c+4)} z}{1 + {}\ddots}}}}}$$

Transformation formulas

Transformation formulas relate two hypergeometric functions at different values of the argument z.
Fractional linear transformations

Euler's transformation is

$$\displaystyle{}_2F_1 (a,b;c;z) = (1-z)^{c-a-b}{}_2F_1 (c-a, c-b;c ; z).$$

It follows by combining the two Pfaff transformations

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

which in turn follow from Euler's integral representation.

If two of the numbers 1 − c, c − 1, a − b, b − a, a + b − c, c − a − b are equal or one of them is 1/2 then there is a quadratic transformation of the hypergeometric function, connecting it to a different value of z related by a quadratic equation. The first examples were given by Kummer (1836), and a complete list was given by Goursat (1881). A typical example is

$$F(a,b;2b;z) = (1-z)^{-(1/2)a}F(\tfrac{1}{2}a, b-\tfrac{1}{2}a; b+\tfrac{1}{2}; \frac{z^2}{4z-4})\,$$

Higher order transformations

If 1−c, a−b, a+b−c differ by signs or two of them are 1/3 or −1/3 then there is a cubic transformation of the hypergeometric function, connecting it to a different value of z related by a cubic equation. The first examples were given by Goursat (1881). A typical example is

$$F(\tfrac{3}{2}a,\tfrac{1}{2}(3a-1);a+\tfrac{1}{2};-z^{2}/3) = (1+z)^{1-3a}F(a-\tfrac{1}{3}, a, 2a, 2z(3+z^2)(1+z)^{-3})\,$$

There are also some transformations of degree 4 and 6. Transformations of other degrees only exist if a, b, and c are certain rational numbers.
Values at special points z

See (Slater 1966, Appendix III) for a list of summation formulas at special points, most of which also appear in (Bailey 1935). (Gessel & Stanton 1982) gives further evaluations at more points. (Koepf 1995) shows how most of these identities can be verifed by computer algorithms.
Special values at z = 1

Gauss's theorem, named for Carl Friedrich Gauss, is the identity

$$\;_2F_1 (a,b;c;1)= \frac{\Gamma(c)\Gamma(c-a-b)}{\Gamma(c-a)\Gamma(c-b)}, Re(c)>Re(a+b)$$

which follows from Euler's integral formula by putting z = 1. It includes the Vandermonde identity, first found by Zhu Shijie (= Chu Shi-Chieh), as a special case.

Dougall's formula generalizes this to the bilateral hypergeometric series at z = 1.
Kummer's theorem (z = −1)

There are many cases where hypergeometric functions can be evaluated at z = −1 by using a quadratic transformation to change z = −1 to z = 1 and then using Gauss's theorem to evaluate the result. A typical example is Kummer's theorem, named for Ernst Kummer:

$$\;_2F_1 (a,b;1+a-b;-1)= \frac{\Gamma(1+a-b)\Gamma(1+\tfrac12a)}{\Gamma(1+a)\Gamma(1+\tfrac12a-b)}$$

which follows from Kummer's quadratic transformations

\begin{align}\;_2F_1(a,b;1+a-b;z)&= (1-z)^{-a} \;_2F_1 \left(\frac a 2, \frac{1+a}2-b; 1+a-b; -\frac{4z}{(1-z)^2}\right)\\ &=(1+z)^{-a} \;_2F_1\left(\frac a 2, \frac{a+1}2; 1+a-b; \frac{4z}{(1+z)^2}\right)\end{align}

and Gauss's theorem by putting z = −1 in the first identity.
Values at z = 1/2

Gauss's second summation theorem is

$$\;_2F_1 \left(a,b;\tfrac12\left(1+a+b\right);\tfrac12\right) = \frac{\Gamma(\tfrac12)\Gamma(\tfrac12\left(1+a+b\right))}{\Gamma(\tfrac12\left(1+a)\right)\Gamma(\tfrac12\left(1+b\right))}.$$

Bailey's theorem is

$$\;_2F_1 \left(a,1-a;c;\tfrac12\right)= \frac{\Gamma(\tfrac12c)\Gamma(\tfrac12\left(1+c\right))}{\Gamma(\tfrac12\left(c+a\right))\Gamma(\tfrac12\left(1+c-a\right))}.$$

Other points

There are many other formulas giving the hypergeometric function as an algebraic number at special rational values of the parameters, some of which are listed in (Gessel & Stanton 1982) and (Koepf 1995). Some typical examples are given by

$${}_2F_1(a,-a;1/2;x^2/4(x-1)) = ((1-x)^a+(1-x)^{-a})/2\,$$

Generalizations

Generalizations of the hypergeometric function include:

Appell series, a 2-variable generalization of hypergeometric series
Basic hypergeometric series where the ratio of terms is a periodic function of the index
Bilateral hypergeometric series pHp are similar to generalized hypergeometric series, but summed over all integers
Elliptic hypergeometric series where the ratio of terms is an elliptic function of the index
Fox H-function, an extension of the Meijer G-function
Fox–Wright function, a generalization of the generalized hypergeometric function
Generalized hypergeometric series pFq where the ratio of terms is a rational function of the index
Heun function, solutions of second order ODE's with four regular singular points
Horn function, 34 distinct convergent hypergeometric series in two variables
Hypergeometric function of a matrix argument, the multivariate generalization of the hypergeometric series
Lauricella hypergeometric series, hypergeometric series of three variables
MacRobert E-function, an extension of the generalized hypergeometric series pFq to the case p>q+1.
Meijer G-function, an extension of the generalized hypergeometric series pFq to the case p>q+1.

References

Abramowitz, Milton; Stegun, Irene A., eds. (1965), "Chapter 15", Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, New York: Dover, pp. 555, ISBN 978-0486612720, MR0167642.
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.
W. N. Bailey Generalized Hypergeometric Series Cambridge 1935
Frits Beukers, Gauss' hypergeometric function (2002) (Lecture notes reviewing basics, as well as triangle maps and monodromy)
Daalhuis, Adri B. Olde (2010), "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
Erdélyi, Arthur; Magnus, Wilhelm; Oberhettinger, Fritz & Tricomi, Francesco G. (1953). Higher transcendental functions. Vol. I. New York–Toronto–London: McGraw–Hill Book Company, Inc.. ISBN 978-0-89874-206-0. MR0058756.
George Gasper and Mizan Rahman, Basic Hypergeometric Series, 2nd Edition, (2004), Encyclopedia of Mathematics and Its Applications, 96, Cambridge University Press, Cambridge. ISBN 0-521-83357-4.
Gauss, Carl Friedrich (1813). "Disquisitiones generales circa seriam infinitam 1 + \tfrac {\alpha \beta} {1 \cdot \gamma} ~x + \tfrac {\alpha (\alpha+1) \beta (\beta+1)} {1 \cdot 2 \cdot \gamma (\gamma+1)} ~x~x + \mbox{etc.} " (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)
Gelfand, I. M.; Gindikin, S. G.; Graev, M. I. (2003) [2000], Selected topics in integral geometry, Translations of Mathematical Monographs, 220, Providence, R.I.: American Mathematical Society, ISBN 978-0-8218-2932-5, MR2000133
Gessel, Ira & Stanton, Dennis (1982). "Strange evaluations of hypergeometric series". SIAM Journal on Mathematical Analysis 13 (2): 295–308. doi:10.1137/0513021. ISSN 0036-1410. MR647127.
Goursat, Édouard (1881). "Sur l'équation différentielle linéaire, qui admet pour intégrale la série hypergéométrique" (in French). Annales Scientifiques de l'École Normale Supérieure 10: 3–142. Retrieved 2008-10-16.
Klein, Felix (1981) (in German). Vorlesungen über die hypergeometrische Funktion. Grundlehren der Mathematischen Wissenschaften. 39. Berlin, New York: Springer-Verlag. ISBN 978-3-540-10455-1. MR668700.
Koepf, Wolfram (1995), "Algorithms for m-fold hypergeometric summation", Journal of Symbolic Computation 20 (4): 399–417, doi:10.1006/jsco.1995.1056, ISSN 0747-7171, MR1384455
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.)
Kummer, Ernst Eduard (1836). "Über die hypergeometrische Reihe 1 + \tfrac {\alpha\cdot \beta} {1 \cdot \gamma} ~x + \tfrac {\alpha (\alpha+1) \beta (\beta+1)} {1 \cdot 2 \cdot \gamma (\gamma+1)} x^2 + \tfrac {\alpha (\alpha+1)(\alpha+2) \beta (\beta+1)(\beta+2)} {1 \cdot 2\cdot 3 \cdot \gamma (\gamma+1)(\gamma+2)} x^3 +\mbox{etc.} " (in German). Journal für die reine und angewandte Mathematik 15: 39–83, 127–172. ISSN 0075-4102.
Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007), "Section 6.13. Hypergeometric Functions", Numerical Recipes: The Art of Scientific Computing (3rd ed.), New York: Cambridge University Press, ISBN 978-0-521-88068-8
Slater, Lucy Joan (1960). Confluent hypergeometric functions. Cambridge, UK: Cambridge University Press. MR0107026.
Slater, Lucy Joan (1966). Generalized hypergeometric functions. Cambridge, UK: Cambridge University Press. ISBN 0-521-06483-X. MR0201688. (there is a 2008 paperback with ISBN 978-0-521-09061-2)
Wall, H. S. (1948). Analytic Theory of Continued Fractions. D. Van Nostrand Company, Inc..
Whittaker, E. T. & Watson, G. N. (1927). A Course of Modern Analysis. Cambridge, UK: Cambridge University Press.
Yoshida, Masaaki (1997). Hypergeometric Functions, My Love: Modular Interpretations of Configuration Spaces. Braunschweig/Wiesbaden: Friedr. Vieweg & Sohn. ISBN 3-528-06925-2. MR1453580.