gamma function

aesthetics  →
being  →
complexity  →
database  →
enterprise  →
ethics  →
fiction  →
history  →
internet  →
knowledge  →
language  →
licensing  →
linux  →
logic  →
method  →
news  →
perception  →
philosophy  →
policy  →
purpose  →
religion  →
science  →
sociology  →
software  →
truth  →
unix  →
wiki  →
essay  →
feed  →
help  →
system  →
wiki  →
critical  →
discussion  →
forked  →
imported  →
original  →
gamma function
[ temporary import ]
please note:
- the content below is remote from Wikipedia
- it has been imported raw for GetWiki
{{pp-vandalism|small=yes}}{{about||the gamma function of ordinals|Veblen function|the gamma distribution in statistics|Gamma distribution|the function used in video and image color representations|Gamma correction}}{{use dmy dates|date=December 2016}}(File:Gamma plot.svg|thumb|right|325px|The gamma function along part of the real axis)In mathematics, the gamma function (represented by {{math|Γ}}, the capital Greek alphabet letter gamma) is one commonly used extension of the factorial function to complex numbers. The gamma function is defined for all complex numbers except the non-positive integers, and for any positive integer {{math|n}},
Gamma(n) = (n-1)! .
Derived by Daniel Bernoulli, for complex numbers with a positive real part the gamma function is defined via a convergent improper integral:
Gamma(z) = int_0^infty x^{z-1} e^{-x},dx, operatorname{Re}z > 0 .
The gamma function then is defined as the analytic continuation of this integral function to a meromorphic function that is holomorphic in the whole complex plane except the non-positive integers, where the function has simple poles.Other extensions of the factorial function do exist, but the gamma function is the most popular and useful. It is a component in various probability-distribution functions, and as such it is applicable in the fields of probability and statistics, as well as combinatorics.The gamma function has no zeroes, so the reciprocal gamma function {{math|1/Γ(z)}} is a holomorphic function. In fact the gamma function corresponds to the Mellin transform of the negative exponential function:
Gamma(z) = { mathcal M e^{-x} } (z)


(File:Factorial Interpolation.svg|thumb|250px|The gamma function interpolates the factorial function to non-integer values.)The gamma function can be seen as a solution to the following interpolation problem:
"Find a smooth curve that connects the points {{math|(x, y)}} given by {{math|y {{=}} (x − 1)!}} at the positive integer values for {{math|x}}."
A plot of the first few factorials makes clear that such a curve can be drawn, but it would be preferable to have a formula that precisely describes the curve, in which the number of operations does not depend on the size of {{math|x}}. The simple formula for the factorial, {{math|x! {{=}} 1 × 2 × … × x}}, cannot be used directly for fractional values of {{math|x}} since it is only valid when {{math|x}} is a natural number (i.e., a positive integer). There are, relatively speaking, no such simple solutions for factorials; no finite combination of sums, products, powers, exponential functions, or logarithms will suffice to express {{math|x!}}; but it is possible to find a general formula for factorials using tools such as integrals and limits from calculus. A good solution to this is the gamma function.There are infinitely many continuous extensions of the factorial to non-integers: infinitely many curves can be drawn through any set of isolated points. The gamma function is the most useful solution in practice, being analytic (except at the non-positive integers), and it can be characterized{{ambiguous|reason=What does it mean that it can be {{}}characterized{{}}?|date=September 2018}} in several ways. However, it is not the only analytic function which extends the factorial, as adding to it any analytic function which is zero on the positive integers, such as {{math|k sin m{{pi}}x}}, will give another function with that property.(File:Gamma plus sin pi z.png|thumb|250px|The gamma function, {{math| Γ(z)}} in blue, plotted along with {{math| Γ(z) + sin({{pi}}z)}} in green. Notice the intersection at positive integers, both are valid analytic continuations of the factorials to the non-integers)A more restrictive property than satisfying the above interpolation is to satisfy the recurrence relation defining a translated version of the factorial function,
f(1) = 1, f(x+1) = x f(x),
for {{math|x}} equal to any positive real number. But this would allow for multiplication by any periodic analytic function which evaluates to one on the positive integers, such as {{math|e'k sin m{{pi}}x}}. There's a final way to solve all this ambiguity: Bohr–Mollerup theorem states that when the condition that {{math|f}} be logarithmically convex (or "super-convex"JOURNAL, Kingman, J. F. C., A Convexity Property of Positive Matrices, The Quarterly Journal of Mathematics, 1961, 12, 1, 283–284, 10.1093/qmath/12.1.283, 1961QJMat..12..283K, ) is added, it uniquely determines {{math|f}} for positive, real inputs. From there, the gamma function can be extended to all real and complex values (except the negative integers and zero) by using the unique analytic continuation of {{math|f}}.{{MathWorld | urlname=Bohr-MollerupTheorem | title=Bohr–Mollerup Theorem}} Also see Euler's infinite product definition below where the properties {{math|1=f(1) = 1}} and {{math|1=f(x + 1) = x f(x)}} together with the asymptotic requirement that {{math|1=limn→+∞ (n−1)! n'x / f(n+x) = 1}} uniquely define the same function.


Main definition

The notation {{math|Γ(z)}} is due to Legendre. If the real part of the complex number {{math|z}} is positive ({{math|Re(z) > 0}}), then the integral
Gamma(z) = int_0^infty x^{z-1} e^{-x}, dx
converges absolutely, and is known as the Euler integral of the second kind (the Euler integral of the first kind defines the beta function). Using integration by parts, one sees that:
Gamma(z+1) & = int_0^infty x^{z} e^{-x} , dx [4pt]
& = Big[-x^z e^{-x}Big]_0^infty + int_0^infty z x^{z-1} e^{-x}, dx [4pt]
& = lim_{xto infty}(-x^z e^{-x}) - (0 e^{-0}) + zint_0^infty x^{z-1} e^{-x}, dx
Recognizing that -x^z e^{-x}to 0 as xto infty,
Gamma(z+1) & = zint_0^infty x^{z-1} e^{-x}, dx [6pt]
& = zGamma(z)
We can calculate Gamma(1)text{:}
Gamma(1) & = int_0^infty x^{1-1} e^{-x},dx [6pt]
& = Big[-e^{-x}Big]_0^infty [6pt]
& = lim_{xto infty}(-e^{-x}) - (-e^{-0}) [6pt]
& = 0 - (-1) [6pt]
& = 1
Given that Gamma(1) = 1 and Gamma(n+1) = nGamma(n),
Gamma(n) = 1 cdot 2 cdot 3 cdots (n-1) = (n-1)!
for all positive integers {{mvar|n}}. This can be seen as an example of proof by induction.The identity Gamma(z) = frac {Gamma(z + 1)} {z} can be used (or, yielding the same result, analytic continuation can be used) to uniquely extend the integral formulation for {{math|Γ(z)}} to a meromorphic function defined for all complex numbers {{mvar|z}}, except integers less than or equal to zero. It is this extended version that is commonly referred to as the gamma function.

Alternative definitions

Euler's definition as an infinite product

When seeking to approximate {{math|z!}} for a complex number {{math|z}}, it turns out that it is effective to first compute {{math|n!}} for some large integer {{math|n}}, then use that to approximate a value for {{math|(n+z)!}}, and then use the recursion relation {{math|1={{math|m!}} = m {{math|(m−1)!}}}} backwards {{math|n}} times, to unwind it to an approximation for {{math|z!}}. Furthermore, this approximation is exact in the limit as {{math|n}} goes to infinity.Specifically, for a fixed integer {{mvar|m}}, it is the case that
lim_{n to infty} frac{n! ; (n+1)^m}{(n+m)!} = 1,,
and we can ask that the same formula is obeyed when the arbitrary integer {{mvar|m}} is replaced by an arbitrary complex number {{mvar|z}}
lim_{n to infty} frac{n! ; (n+1)^z}{(n+z)!} = 1,.
Multiplying both sides by {{math|z!}} gives
z! &= lim_{n to infty} n!frac{z!}{(n+z)!} (n+1)^z [8pt]
&= lim_{n to infty} frac{1 cdots n}{(1+z) cdots (n+z)} (n+1)^z [8pt]
&= lim_{n to infty} frac{1 cdots n}{(1+z) cdots (n+z)} left[left(1+frac11right)left(1+frac12right)cdots left(1+frac1nright)right]^z [8pt]
&= prod_{n=1}^infty left[ frac{1}{1+frac{z}{n}} left(1 + frac{1}{n}right)^z right].
end{align}This infinite product formula converges for all complex numbers {{mvar|z}} except the negative integers, which fail because trying to use the recursion relation {{math|1={{math|m!}} = m {{math|(m âˆ’ 1)!}}}} backwards through the value {{math|1=m = 0}} involves a division by zero.Similarly for the gamma function, the definition as an infinite product due to Euler is valid for all complex numbers z except the non-positive integers:
Gamma(z) = frac{1}{z} prod_{n=1}^{infty} frac{left(1 + frac{1}{n}right)^z}{1+frac{z}{n}} ,.
By this construction, the gamma function is the unique function that simultaneously satisfies Gamma(1) = 1, Gamma(z+1) = z Gamma(z) for all complex numbers z except the non-positive integers, and lim_{n to infty} frac{Gamma(n+z)}{(n-1)!;n^z} = 1 for all complex numbers z.

Weierstrass's definition

The definition for the gamma function due to Weierstrass is also valid for all complex numbers {{math|z}} except the non-positive integers:
Gamma(z) = frac{e^{-gamma z}} z prod_{n=1}^infty left(1 + frac z n right)^{-1} e^{z/n}
where gamma approx 0.577216 is the Euler–Mascheroni constant.

In terms of generalized Laguerre polynomials

A representation of the incomplete gamma function in terms of generalized Laguerre polynomials is
Gamma(z,x)=x^z e^{-x} sum_{n=0}^infty frac{L_n ^{(z)}(x)}{n+1},
which converges for operatorname{Re}(z) > -1 and x>0.{{dlmf|authorlink=Richard Askey|first=R. A.|last= Askey|first2= R.|last2= Roy |id=8.7 |title=Series Expansions|ref=none}}A somewhat unusual representation of the gamma function in terms of Laguerre polynomials is given by
Gamma(z)=t^z sum_{n=0}^{infty} frac{L_n^{(z)}(t)}{z+n},, which converges for Re(z)  Gamma'(x).
The last of these statements is, essentially by definition, the same as the statement that psi^{(1)}(x) > 0, where psi^{(1)} is the polygamma function of order 1. To prove the logarithmic convexity of the gamma function, it therefore suffices to observe that psi^{(1)} has a series representation which, for positive real {{math|x}}, consists of only positive terms.Logarithmic convexity and Jensen's inequality together imply, for any positive real numbers x_1, ldots, x_n and a_1, ldots, a_n,
Gammaleft(frac{a_1x_1 + cdots + a_nx_n}{a_1 + cdots + a_n}right) le bigl(Gamma(x_1)^{a_1} cdots Gamma(x_n)^{a_n}bigr)^{frac{1}{a_1 + cdots + a_n}}.
There are also bounds on ratios of gamma functions. The best-known is Gautschi's inequality, which says that for any positive real number {{math|x}} and any {{math|s ∈ (0, 1)}},
x^{1 - s} < frac{Gamma(x + 1)}{Gamma(x + s)} < (x + 1)^{1 - s}.

Stirling's formula

(File:Gamma1.png|thumb|Representation of the gamma function in the complex plane. Each point {{nowrap|z}} is colored according to the argument of {{nowrap|Gamma(z).}} The contour plot of the modulus {{nowrap||Gamma(z)|}} is also displayed.)(File:Gamma abs 3D.png|thumb|3-dimensional plot of the absolute value of the complex gamma function)The behavior of Gamma(z) for an increasing positive variable is simple: it grows quickly — faster than an exponential function. Asymptotically as zto infty , the magnitude of the gamma function is given by Stirling's formula
Gamma(z+1)simsqrt{2pi z}left(frac{z}{e}right)^z,
where the symbol sim means that the ratio of the two sides converges to 1 or asymptotically converges.


The behavior for non-positive z is more intricate. Euler's integral does not converge for z le 0, but the function it defines in the positive complex half-plane has a unique analytic continuation to the negative half-plane. One way to find that analytic continuation is to use Euler's integral for positive arguments and extend the domain to negative numbers by repeated application of the recurrence formula,
choosing n such that z + n is positive. The product in the denominator is zero when z equals any of the integers 0, -1, -2,cdots. Thus, the gamma function must be undefined at those points to avoid division by zero; it is a meromorphic function with simple poles at the non-positive integers.We can rewrite this definition as:
(z+n) Gamma(z)=frac{Gamma(z+n+1)}{z(z+1)cdots(z+n-1)}.
For a function f of a complex variable z, at a simple pole c, the residue of f is given by:
operatorname{Res}(f,c)=lim_{zto c}(z-c)f(z).
When z = -n,
Gamma(z+n+1) = Gamma(1) = 1
z(z+1)cdots(z+n-1) = (-1)^n n!
So the residues of the gamma function at those points are:
operatorname{Res}(Gamma,-n)=frac{(-1)^n}{n!}.{{mathworld|urlname=GammaFunction |title=Gamma Function}}
The gamma function is non-zero everywhere along the real line, although it comes arbitrarily close to zero as {{math|z → −∞}}. There is in fact no complex number z for which Gamma (z) = 0, and hence the reciprocal gamma function frac {1}{Gamma (z)} is an entire function, with zeros at z = 0, -1, -2, cdots.


The gamma function has a local minimum at z_{min} approx 1.46163 where it attains the value Gammaleft(z_{min}right) approx 0.885603. The gamma function must alternate sign between the poles because the product in the forward recurrence contains an odd number of negative factors if the number of poles between z and z + n is odd, and an even number if the number of poles is even.

Integral representations

There are many formulas, besides the Euler integral of the second kind, that express the gamma function as an integral. For instance, when the real part of {{math|z}} is positive,Whittaker and Watson, 12.2 example 1.
Gamma(z) = int_0^1 left(log frac{1}right)^{z-1},dt.
Binet's first integral formula for the gamma function states that, when the real part of {{math|z}} is positive, then:Whittaker and Watson, 12.31.
log Gamma(z) = left(z - frac{1}{2}right)log z - z + frac{1}{2}log(2pi) + int_0^infty left(frac{1}{2} - frac{1} + frac{1}{e^t - 1}right)frac{e^{-tz}},dt.
The integral on the right-hand side may be interpreted as a Laplace transform. That is,
logleft(Gamma(z)left(frac{e}{z}right)^zsqrt{2pi z}right) = mathcal{L}left(frac{1}{2t} - frac{1}{t^2} + frac{1}{t(e^t - 1)}right)(z).
Binet's second integral formula states that, again when the real part of {{math|z}} is positive, then:Whittaker and Watson, 12.32.
log Gamma(z) = left(z - frac{1}{2}right)log z - z + frac{1}{2}log(2pi) + 2int_0^infty frac{arctan(t/z)}{e^{2pi t} - 1},dt.
Let {{math|C}} be a Hankel contour, meaning a path that begins and ends at the point {{math|∞}} on the Riemann sphere, whose unit tangent vector converges to {{math|−1}} at the start of the path and to {{math|1}} at the end, which has winding number 1 around {{math|0}}, and which does not cross {{math|[0, ∞)}}. Fix a branch of log(-t) by taking a branch cut along {{math|[0, ∞)}} and by taking log(-t) to be real when {{math|t}} is on the negative real axis. Assume {{math|z}} is not an integer. Then Hankel's formula for the gamma function is:Whittaker and Watson, 12.22.
Gamma(z) = -frac{1}{2isin pi z}int_C (-t)^{z-1}e^{-t},dt,
where (-t)^{z-1} is interpreted as exp((z-1)log(-t)). The reflection formula leads to the closely related expression
frac{1}{Gamma(z)} = frac{i}{2pi}int_C (-t)^{-z}e^{-t},dt,
again valid whenever {{math|z}} is not an integer.

Fourier series expansion

The logarithm of the gamma function has the following Fourier series expansion for 0 < z < 1:
lnGamma(z) = left(frac{1}{2} - zright)(gamma + ln 2) + (1 - z)lnpi - frac{1}{2}lnsin(pi z) + frac{1}{pi}sum_{n=1}^infty frac{ln n}{n}sin (2pi n z),
which was for a long time attributed to Ernst Kummer, who derived it in 1847.BOOK, Harry, Bateman, Arthur, Erdélyi, Higher Transcendental Functions, McGraw-Hill, 1955, BOOK, H. M., Srivastava, J., Choi, Series Associated with the Zeta and Related Functions, Kluwer Academic, The Netherlands, 2001, However, Iaroslav Blagouchine discovered that Carl Johan Malmsten first derived this series in 1842.JOURNAL, 10.1007/s11139-013-9528-5, Iaroslav V., Blagouchine, Rediscovery of Malmsten's integrals, their evaluation by contour integration methods and some related results, Ramanujan J., 35, 1, 21–110, 2014,weblink JOURNAL, 10.1007/s11139-015-9763-z, Iaroslav V., Blagouchine, Erratum and Addendum to "Rediscovery of Malmsten's integrals, their evaluation by contour integration methods and some related results", Ramanujan J., 42, 3, 777–781, 2016,

Raabe's formula

In 1840 Joseph Ludwig Raabe proved that
int_a^{a+1}lnGamma(z), dz = tfrac12ln2pi + aln a - a,quad a>0.
In particular, if {{math|a {{=}} 0}} then
int_0^1lnGamma(z), dz = tfrac12ln2pi.
The latter can be derived taking the logarithm in the above multiplication formula, which gives an expression for the Riemann sum of the integrand. Taking the limit for 'm'→∞ gives the formula.

Pi function

An alternative notation which was originally introduced by Gauss and which was sometimes used is the pi function, which in terms of the gamma function is
Pi(z) = Gamma(z+1) = z Gamma(z) = int_0^infty e^{-t} t^z, dt,
so that {{math|Π(n) {{=}} n!}} for every non-negative integer {{math|n}}.Using the pi function the reflection formula takes on the form
Pi(z) Pi(-z) = frac{pi z}{sin( pi z)} = frac{1}{operatorname{sinc}(z)}
where {{math|sinc}} is the normalized sinc function, while the multiplication theorem takes on the form
Pileft(frac{z}{m}right) , Pileft(frac{z-1}{m}right) cdots Pileft(frac{z-m+1}{m}right) = (2 pi)^{frac{m-1}{2}} m^{-z-frac12} Pi(z).
We also sometimes find
pi(z) = frac{1}{Pi(z)},
which is an entire function, defined for every complex number, just like the reciprocal gamma function. That pileft(zright) is entire entails it has no poles, so Pileft(zright), like Gammaleft(zright), has no zeros.The volume of an {{math|n}}-ellipsoid with radii {{math|r{{sub|1}}, …, r{{sub|n}}}} can be expressed as
V_n(r_1,dotsc,r_n)=frac{pi^{frac{n}{2}}}{Pileft(frac{n}{2} right)} prod_{k=1}^n r_k.

Relation to other functions

  • In the first integral above, which defines the gamma function, the limits of integration are fixed. The upper and lower incomplete gamma functions are the functions obtained by allowing the lower or upper (respectively) limit of integration to vary.
  • The gamma function is related to the beta function by the formula

Beta(x,y)=frac{Gamma(x) ; Gamma(y)}{Gamma(x+y)}.

pi^{-frac{z}{2}} ; Gammaleft(frac{z}{2}right) zeta(z) = pi^{-frac{1-z}{2}} ; Gammaleft(frac{1-z}{2}right) ; zeta(1-z).
It also appears in the following formula:
zeta(z) Gamma(z) = int_0^infty frac{u^{z}}{e^u - 1} , frac{du}{u},
which is valid only for {{math|Re(z) > 1}}. The logarithm of the gamma function satisfies the following formula due to Lerch:
logGamma(x) = zeta_{H}'(0,x) - zeta'(0),
where {{math|ζH}} is the Hurwitz zeta function, {{math|ζ}} is the Riemann zeta function and the prime ({{math|′}}) denotes differentiation in the first variable.

langletau^nrangle equiv int_0^infty dt, t^{n-1}, e^{ - left( frac{tau} right)^beta} = frac{tau^n}{beta}Gamma left({n over beta }right).

Particular values

Some particular values of the gamma function are:
Gammaleft(-tfrac{3}{2}right) &=& tfrac{4}{3} sqrt{pi} &approx& 2.363,271,801,207 Gammaleft(-tfrac{1}{2}right) &=& -2sqrt{pi} &approx& -3.544,907,701,811 Gammaleft(tfrac{1}{2}right) &=& sqrt{pi} &approx& 1.772,453,850,906 Gamma(1) &=& 0! &=& 1 Gammaleft(tfrac{3}{2}right) &=& tfrac{1}{2}sqrt{pi} &approx& 0.886,226,925,453 Gamma(2) &=& 1! &=& 1 Gammaleft(tfrac{5}{2}right) &=& tfrac{3}{4}sqrt{pi} &approx& 1.329,340,388,179 Gamma(3) &=& 2! &=& 2 Gammaleft(tfrac{7}{2}right) &=& tfrac{15}{8}sqrt{pi} &approx& 3.323,350,970,448 Gamma(4) &=& 3! &=& 6end{array}The complex-valued gamma function is undefined for non-positive integers, but in these cases the value can be defined in the Riemann sphere as {{math|∞}}. The reciprocal gamma function is well defined and analytic at these values (and in the entire complex plane):
frac{1}{Gamma(-3)} = frac{1}{Gamma(-2)} = frac{1}{Gamma(-1)} = frac{1}{Gamma(0)} = 0.

The log-gamma function

(File:LogGamma Analytic Function.png|thumb|The analytic function {{math|log Γ(z)}})Because the gamma and factorial functions grow so rapidly for moderately large arguments, many computing environments include a function that returns the natural logarithm of the gamma function (often given the name lgamma or lngamma in programming environments or gammaln in spreadsheets); this grows much more slowly, and for combinatorial calculations allows adding and subtracting logs instead of multiplying and dividing very large values. It is often defined asWEB, Log Gamma Function,weblink Wolfram MathWorld, 3 January 2019,
ln Gamma ( z ) = - gamma z - ln z + sum _ { k = 1 } ^ { infty } left[ frac { z } { k } - ln left( 1 + frac { z } { k } right) right]
The digamma function, which is the derivative of this function, is also commonly seen.In the context of technical and physical applications, e.g. with wave propagation, the functional equation
ln Gamma(z) = ln Gamma(z+1) - ln z
is often used since it allows one to determine function values in one strip of width 1 in {{math|z}} from the neighbouring strip. In particular, starting with a good approximation for a {{math|z}} with large real part one may go step by step down to the desired {{math|z}}. Following an indication of Carl Friedrich Gauss, Rocktaeschel (1922) proposed for {{math|ln(Γ(z))}} an approximation for large {{math|Re(z)}}:
ln Gamma(z) approx (z - tfrac{1}{2}) ln z - z + tfrac{1}{2}ln(2pi).
This can be used to accurately approximate {{math|ln(Γ(z))}} for {{math|z}} with a smaller {{math|Re(z)}} via (P.E.Böhmer, 1939)
lnGamma(z-m) = lnGamma(z) - sum_{k=1}^{m} ln(z-k).
A more accurate approximation can be obtained by using more terms from the asymptotic expansions of {{math|ln(Γ(z))}} and {{math|Γ(z)}}, which are based on Stirling's approximation.
Gamma(z)sim z^{z - frac12} e^{-z} sqrt{2pi} left( 1 + frac{1}{12z} + frac{1}{288z^2} - frac{139}{51,840 z^3} - frac{571}{2,488,320 z^4}

as {{math|{{abs|z}} → ∞}} at constant {{math|{{abs|arg(z)}} < π}}.
In a more "natural" presentation:
ln Gamma(z) sim z ln (z) - z - tfrac12 ln left (frac{z}{2pi} right ) + frac{1}{12z} - frac{1}{360z^3} +frac{1}{1260 z^5} as {{math|{{abs|z}} → ∞}} at constant {{math|{{abs|arg(z)}} < π}}.
The coefficients of the terms with {{math|k > 1}} of {{math|z−k + 1}} in the last expansion are simply
where the {{math|Bk}} are the Bernoulli numbers.


The Bohr–Mollerup theorem states that among all functions extending the factorial functions to the positive real numbers, only the gamma function is log-convex, that is, its natural logarithm is convex on the positive real axis.In a certain sense, the {{math|ln(Γ)}} function is the more natural form; it makes some intrinsic attributes of the function clearer. A striking example is the Taylor series of {{math|ln(Γ)}} around 1:
ln Gamma(1+z)= -gamma z +sum_{k=2}^infty frac{zeta(k)}{k} , (-z)^k qquad forall; |z| < 1
with {{math|ζ(k)}} denoting the Riemann zeta function at {{math|k}}.So, using the following property:
zeta(s) Gamma(s) = int_0^infty frac{t^s}{e^t-1} , frac{dt}
we can find an integral representation for the {{math|ln(Γ)}} function:
ln Gamma(1+z)= -gamma z + int_0^infty frac{e^{-zt}-1+zt}{t(e^t -1)} , dt
or, setting {{math|z {{=}} 1}} and calculating {{math|γ}}:
ln Gamma(1+z)= int_0^infty frac{e^{-zt}-ze^{-t}-1+z}{t(e^t -1)} , dt
There also exist special formulas for the logarithm of the gamma function for rational {{math|z}}. For instance, if k and n are integers with k −1}}.It can also be written in terms of the Hurwitz zeta function:JOURNAL, Victor S., Adamchik, Polygamma functions of negative order, J. Comput. Appl. Math., 100, 2, 1998, 191–199, 10.1016/S0377-0427(98)00192-7, JOURNAL, R. W., Gosper, textstyle int_{n/4}^{m/6} log F(z) ,dz in special functions, q-series and related topics, J. Am. Math. Soc., 14, 1997,
int_0^z ln Gamma(x) , dx = frac{z}{2} ln(2 pi) + frac{z(1-z)}{2} - zeta'(-1) + zeta'(-1,z) .


(File:Mplwp factorial gamma stirling.svg|thumb|right|upright=1.35|Comparison gamma (blue line) with the factorial (blue dots) and Stirling's approximation (red line))Complex values of the gamma function can be computed numerically with arbitrary precision using Stirling's approximation or the Lanczos approximation.The gamma function can be computed to fixed precision for {{math|Re(z) ∈ [1,2]}} by applying integration by parts to Euler's integral. For any positive number {{math|x}} the gamma function can be written
Gamma(z) &= int_0^x e^{-t} t^z , frac{dt} + int_x^infty e^{-t} t^z, frac{dt} &= x^z e^{-x} sum_{n=0}^infty frac{x^n}{z(z+1) cdots (z+n)} + int_x^infty e^{-t} t^z , frac{dt}.end{align}When {{math|Re(z) ∈ [1,2]}} and {{math|x ≥ 1}}, the absolute value of the last integral is smaller than {{math|(x + 1)e−x}}. By choosing a large enough {{math|x}}, this last expression can be made smaller than {{math|2−N}} for any desired value {{math|N}}. Thus, the gamma function can be evaluated to {{math|N}} bits of precision with the above series.A fast algorithm for calculation of the Euler gamma function for any algebraic argument (including rational) was constructed by E.A. Karatsuba,E.A. Karatsuba, Fast evaluation of transcendental functions. Probl. Inf. Transm. Vol.27, No.4, pp. 339–360 (1991).E.A. Karatsuba, On a new method for fast evaluation of transcendental functions. Russ. Math. Surv. Vol.46, No.2, pp. 246–247 (1991).E.A. Karatsuba "Fast Algorithms and the FEE Method".For arguments that are integer multiples of {{math|{{sfrac|1|24}}}}, the gamma function can also be evaluated quickly using arithmetic–geometric mean iterations (see particular values of the gamma function and {{harvtxt|Borwein|Zucker|1992}}).


One author describes the gamma function as "Arguably, the most common special function, or the least 'special' of them. The other transcendental functions […] are called 'special' because you could conceivably avoid some of them by staying away from many specialized mathematical topics. On the other hand, the gamma function {{math|y {{=}} Γ(x)}} is most difficult to avoid."Michon, G. P. "Trigonometry and Basic Functions". Numericana. Retrieved May 5, 2007.

Integration problems

The gamma function finds application in such diverse areas as quantum physics, astrophysics and fluid dynamics.Chaudry, M. A. & Zubair, S. M. (2001). On A Class of Incomplete Gamma Functions with Applications. p. 37 The gamma distribution, which is formulated in terms of the gamma function, is used in statistics to model a wide range of processes; for example, the time between occurrences of earthquakes.Rice, J. A. (1995). Mathematical Statistics and Data Analysis (Second Edition). p. 52–53The primary reason for the gamma function's usefulness in such contexts is the prevalence of expressions of the type {{math|f(t) e{{sup|−g(t)}}}} which describe processes that decay exponentially in time or space. Integrals of such expressions can occasionally be solved in terms of the gamma function when no elementary solution exists. For example, if {{math|f}} is a power function and {{math|g}} is a linear function, a simple change of variables gives the evaluation
int_0^infty t^b e^{-at} ,dt = frac{Gamma(b+1)}{a^{b+1}}.
The fact that the integration is performed along the entire positive real line might signify that the gamma function describes the cumulation of a time-dependent process that continues indefinitely, or the value might be the total of a distribution in an infinite space.It is of course frequently useful to take limits of integration other than 0 and {{math|∞}} to describe the cumulation of a finite process, in which case the ordinary gamma function is no longer a solution; the solution is then called an incomplete gamma function. (The ordinary gamma function, obtained by integrating across the entire positive real line, is sometimes called the complete gamma function for contrast.)An important category of exponentially decaying functions is that of Gaussian functions
and integrals thereof, such as the error function. There are many interrelations between these functions and the gamma function; notably, {{math|{{sqrt|π}}}} obtained by evaluating {{math|Γ({{sfrac|1|2}})}} is the "same" as that found in the normalizing factor of the error function and the normal distribution.The integrals we have discussed so far involve transcendental functions, but the gamma function also arises from integrals of purely algebraic functions. In particular, the arc lengths of ellipses and of the lemniscate, which are curves defined by algebraic equations, are given by elliptic integrals that in special cases can be evaluated in terms of the gamma function. The gamma function can also be used to calculate "volume" and "area" of {{math|n}}-dimensional hyperspheres.Another important special case is that of the beta function
Beta(x,y) = int_0^1 t^{x-1}(1-t)^{y-1},dt = frac{Gamma(x),Gamma(y)}{Gamma(x+y)}.

Calculating products

The gamma function's ability to generalize factorial products immediately leads to applications in many areas of mathematics; in combinatorics, and by extension in areas such as probability theory and the calculation of power series. Many expressions involving products of successive integers can be written as some combination of factorials, the most important example perhaps being that of the binomial coefficient
binom n k = frac{n!}{k!(n-k)!}.
The example of binomial coefficients motivates why the properties of the gamma function when extended to negative numbers are natural. A binomial coefficient gives the number of ways to choose {{math|k}} elements from a set of {{math|n}} elements; if {{math|k > n}}, there are of course no ways. If {{math|k > n}}, {{math|(n − k)!}} is the factorial of a negative integer and hence infinite if we use the gamma function definition of factorials—dividing by infinity gives the expected value of 0.We can replace the factorial by a gamma function to extend any such formula to the complex numbers. Generally, this works for any product wherein each factor is a rational function of the index variable, by factoring the rational function into linear expressions. If {{math|P}} and {{math|Q}} are monic polynomials of degree {{math|m}} and {{math|n}} with respective roots {{math|p{{sub|1}}, …, p{{sub|m}}}} and {{math|q{{sub|1}}, …, q{{sub|n}}}}, we have
prod_{i=a}^b frac{P(i)}{Q(i)} = left( prod_{j=1}^m frac{Gamma(b-p_j+1)}{Gamma(a-p_j)} right) left( prod_{k=1}^n frac{Gamma(a-q_k)}{Gamma(b-q_k+1)} right).
If we have a way to calculate the gamma function numerically, it is a breeze to calculate numerical values of such products. The number of gamma functions in the right-hand side depends only on the degree of the polynomials, so it does not matter whether {{math|b − a}} equals 5 or 105. By taking the appropriate limits, the equation can also be made to hold even when the left-hand product contains zeros or poles.By taking limits, certain rational products with infinitely many factors can be evaluated in terms of the gamma function as well. Due to the Weierstrass factorization theorem, analytic functions can be written as infinite products, and these can sometimes be represented as finite products or quotients of the gamma function. We have already seen one striking example: the reflection formula essentially represents the sine function as the product of two gamma functions. Starting from this formula, the exponential function as well as all the trigonometric and hyperbolic functions can be expressed in terms of the gamma function.More functions yet, including the hypergeometric function and special cases thereof, can be represented by means of complex contour integrals of products and quotients of the gamma function, called Mellin–Barnes integrals.

Analytic number theory

An elegant and deep application of the gamma function is in the study of the Riemann zeta function. A fundamental property of the Riemann zeta function is its functional equation:
Gammaleft(frac{s}{2}right)zeta(s)pi^{-frac{s}{2}} = Gammaleft(frac{1-s}{2}right)zeta(1-s)pi^{-frac{1-s}{2}}.
Among other things, this provides an explicit form for the analytic continuation of the zeta function to a meromorphic function in the complex plane and leads to an immediate proof that the zeta function has infinitely many so-called "trivial" zeros on the real line. Borwein et al. call this formula "one of the most beautiful findings in mathematics".BOOK, Borwein, J., Bailey, D. H. & Girgensohn, R., 2003, Experimentation in Mathematics, A. K. Peters, 133, 978-1-56881-136-9, Another champion for that title might be
zeta(s) ; Gamma(s) = int_0^infty frac{t^s}{e^t-1} , frac{dt}.
Both formulas were derived by Bernhard Riemann in his seminal 1859 paper "Über die Anzahl der Primzahlen unter einer gegebenen Größe" ("On the Number of Prime Numbers less than a Given Quantity"), one of the milestones in the development of analytic number theory—the branch of mathematics that studies prime numbers using the tools of mathematical analysis. Factorial numbers, considered as discrete objects, are an important concept in classical number theory because they contain many prime factors, but Riemann found a use for their continuous extension that arguably turned out to be even more important.


The gamma function has caught the interest of some of the most prominent mathematicians of all time. Its history, notably documented by Philip J. Davis in an article that won him the 1963 Chauvenet Prize, reflects many of the major developments within mathematics since the 18th century. In the words of Davis, "each generation has found something of interest to say about the gamma function. Perhaps the next generation will also."JOURNAL, Davis, P. J., 1959, Leonhard Euler's Integral: A Historical Profile of the Gamma Function, American Mathematical Monthly, 66, 10, 849–869,weblink 3 December 2016, 10.2307/2309786, 2309786,

18th century: Euler and Stirling

File:DanielBernoulliLettreAGoldbach-1729-10-06.jpg|thumb|Daniel Bernoulli's letter to Christian GoldbachChristian GoldbachThe problem of extending the factorial to non-integer arguments was apparently first considered by Daniel Bernoulli and Christian Goldbach in the 1720s, and was solved at the end of the same decade by Leonhard Euler. Euler gave two different definitions: the first was not his integral but an infinite product,
n! = prod_{k=1}^infty frac{left(1+frac{1}{k}right)^n}{1+frac{n}{k}},,
of which he informed Goldbach in a letter dated October 13, 1729. He wrote to Goldbach again on January 8, 1730, to announce his discovery of the integral representation
n!=int_0^1 (-ln s)^n, ds,,
which is valid for {{math|n > 0}}. By the change of variables {{math|t {{=}} −ln s}}, this becomes the familiar Euler integral. Euler published his results in the paper "De progressionibus transcendentibus seu quarum termini generales algebraice dari nequeunt" ("On transcendental progressions, that is, those whose general terms cannot be given algebraically"), submitted to the St. Petersburg Academy on November 28, 1729.Euler's paper was published in Commentarii academiae scientiarum Petropolitanae 5, 1738, 36–57. See E19 -- De progressionibus transcendentibus seu quarum termini generales algebraice dari nequeunt, from The Euler Archive, which includes a scanned copy of the original article. Euler further discovered some of the gamma function's important functional properties, including the reflection formula.James Stirling, a contemporary of Euler, also attempted to find a continuous expression for the factorial and came up with what is now known as Stirling's formula. Although Stirling's formula gives a good estimate of {{math|n!}}, also for non-integers, it does not provide the exact value. Extensions of his formula that correct the error were given by Stirling himself and by Jacques Philippe Marie Binet.

19th century: Gauss, Weierstrass and Legendre

(File:Euler factorial paper.png|thumb|250px|alt=De progressionibus transcendentibus, seu quarum termini generales algebraicae dari nequeunt|The first page of Euler's paper)Carl Friedrich Gauss rewrote Euler's product as
Gamma(z) = lim_{mtoinfty}frac{m^z m!}{z(z+1)(z+2)cdots(z+m)}
and used this formula to discover new properties of the gamma function. Although Euler was a pioneer in the theory of complex variables, he does not appear to have considered the factorial of a complex number, as instead Gauss first did.BOOK, Remmert, R., Kay, L. D., Classical Topics in Complex Function Theory, Springer, 2006, 978-0-387-98221-2, Gauss also proved the multiplication theorem of the gamma function and investigated the connection between the gamma function and elliptic integrals.Karl Weierstrass further established the role of the gamma function in complex analysis, starting from yet another product representation,
Gamma(z) = frac{e^{-gamma z}}{z} prod_{k=1}^infty left(1 + frac{z}{k}right)^{-1} e^frac{z}{k},
where {{math|γ}} is the Euler–Mascheroni constant. Weierstrass originally wrote his product as one for {{math|{{sfrac|1|Γ}}}}, in which case it is taken over the function's zeros rather than its poles. Inspired by this result, he proved what is known as the Weierstrass factorization theorem—that any entire function can be written as a product over its zeros in the complex plane; a generalization of the fundamental theorem of algebra.The name gamma function and the symbol {{math|Γ}} were introduced by Adrien-Marie Legendre around 1811; Legendre also rewrote Euler's integral definition in its modern form. Although the symbol is an upper-case Greek gamma, there is no accepted standard for whether the function name should be written "gamma function" or "Gamma function" (some authors simply write "{{math|Γ}}-function"). The alternative "pi function" notation {{math|Π(z) {{=}} z!}} due to Gauss is sometimes encountered in older literature, but Legendre's notation is dominant in modern works.It is justified to ask why we distinguish between the "ordinary factorial" and the gamma function by using distinct symbols, and particularly why the gamma function should be normalized to {{math|Γ(n + 1) {{=}} n!}} instead of simply using "{{math|Γ(n) {{=}} n!}}". Consider that the notation for exponents, {{math|xn}}, has been generalized from integers to complex numbers {{math|xz}} without any change. Legendre's motivation for the normalization does not appear to be known, and has been criticized as cumbersome by some (the 20th-century mathematician Cornelius Lanczos, for example, called it "void of any rationality" and would instead use {{math|z!}}).JOURNAL, Lanczos, C., 1964, A precision approximation of the gamma function, J. SIAM Numer. Anal. Ser. B, 1, Legendre's normalization does simplify a few formulae, but complicates most others. From a modern point of view, the Legendre normalization of the Gamma function is the integral of the additive character {{math|e−x}} against the multiplicative character {{math|xz}} with respect to the Haar measure tfrac{dx}{x} on the Lie group {{math|R+}}. Thus this normalization makes it clearer that the gamma function is a continuous analogue of a Gauss sum.

19th–20th centuries: characterizing the gamma function

It is somewhat problematic that a large number of definitions have been given for the gamma function. Although they describe the same function, it is not entirely straightforward to prove the equivalence. Stirling never proved that his extended formula corresponds exactly to Euler's gamma function; a proof was first given by Charles Hermite in 1900.BOOK, Knuth, D. E., The Art of Computer Programming, Volume 1 (Fundamental Algorithms), Addison-Wesley, 1997, Instead of finding a specialized proof for each formula, it would be desirable to have a general method of identifying the gamma function.One way to prove would be to find a differential equation that characterizes the gamma function. Most special functions in applied mathematics arise as solutions to differential equations, whose solutions are unique. However, the gamma function does not appear to satisfy any simple differential equation. Otto Hölder proved in 1887 that the gamma function at least does not satisfy any algebraic differential equation by showing that a solution to such an equation could not satisfy the gamma function's recurrence formula, making it a transcendentally transcendental function. This result is known as Hölder's theorem.A definite and generally applicable characterization of the gamma function was not given until 1922. Harald Bohr and Johannes Mollerup then proved what is known as the Bohr–Mollerup theorem: that the gamma function is the unique solution to the factorial recurrence relation that is positive and logarithmically convex for positive {{math|z}} and whose value at 1 is 1 (a function is logarithmically convex if its logarithm is convex).The Bohr–Mollerup theorem is useful because it is relatively easy to prove logarithmic convexity for any of the different formulas used to define the gamma function. Taking things further, instead of defining the gamma function by any particular formula, we can choose the conditions of the Bohr–Mollerup theorem as the definition, and then pick any formula we like that satisfies the conditions as a starting point for studying the gamma function. This approach was used by the Bourbaki group.Borwein & CorlessJOURNAL, Borwein, Jonathan M., Jonathan Borwein, Corless, Robert M., Gamma and Factorial in the Monthly, American Mathematical Monthly, en, Mathematical Association of America, 2017, 125, 5, 400–24, 1703.05349, 2017arXiv170305349B, 10.1080/00029890.2018.1420983, review three centuries of work on the gamma function.

Reference tables and software

Although the gamma function can be calculated virtually as easily as any mathematically simpler function with a modern computer—even with a programmable pocket calculator—this was of course not always the case. Until the mid-20th century, mathematicians relied on hand-made tables; in the case of the gamma function, notably a table computed by Gauss in 1813 and one computed by Legendre in 1825.File:Jahnke gamma function.png|thumb|300px|A hand-drawn graph of the absolute value of the complex gamma function, from Tables of Higher Functions by JahnkeJahnkeTables of complex values of the gamma function, as well as hand-drawn graphs, were given in Tables of Higher Functions by Jahnke and {{ill|Fritz Emde|de|lt=Emde}}, first published in Germany in 1909. According to Michael Berry, "the publication in J&E of a three-dimensional graph showing the poles of the gamma function in the complex plane acquired an almost iconic status."NEWS, Berry, M.,weblink Why are special functions special?, Physics Today, April 2001, There was in fact little practical need for anything but real values of the gamma function until the 1930s, when applications for the complex gamma function were discovered in theoretical physics. As electronic computers became available for the production of tables in the 1950s, several extensive tables for the complex gamma function were published to meet the demand, including a table accurate to 12 decimal places from the U.S. National Bureau of Standards.Abramowitz and Stegun became the standard reference for this and many other special functions after its publication in 1964.Double-precision floating-point implementations of the gamma function and its logarithm are now available in most scientific computing software and special functions libraries, for example TK Solver, Matlab, GNU Octave, and the GNU Scientific Library. The gamma function was also added to the C standard library (math.h). Arbitrary-precision implementations are available in most computer algebra systems, such as Mathematica and Maple. PARI/GP, MPFR and MPFUN contain free arbitrary-precision implementations. A little-known feature of the calculator app included with the Android operating system is that it will accept fractional values as input to the factorial function and return the equivalent gamma function value. The same is true for Windows Calculator (in scientific mode).

See also

{{div col}} {{div col end}}


  • {{Citizendium}}

Further reading

  • BOOK, Milton, Abramowitz, Irene A., Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, New York, Dover, 1972,weblink Chapter 6, Abramowitz and Stegun,
  • BOOK, G. E., Andrews, R., Askey, Richard Askey, R., Roy, Special Functions, Cambridge University Press, New York, 1999, 978-0-521-78988-2, Chapter 1 (Gamma and Beta functions),
  • BOOK, Artin, Emil, Emil Artin, The Gamma Function, Rosen, Michael, Exposition by Emil Artin: a selection, History of Mathematics, 30, Providence, RI, American Mathematical Society, 2006,
  • {{dlmf|authorlink=Richard Askey|first=R.|last= Askey|first2= R.|last2= Roy |id=5 |ref=none }}
  • JOURNAL, Birkhoff, George D., George David Birkhoff, Note on the gamma function, Bull. Amer. Math. Soc., 1913, 20, 1, 1–10, 1559418, 10.1090/s0002-9904-1913-02429-7,
  • BOOK, P. E., Böhmer, Differenzengleichungen und bestimmte Integrale, Differential Equations and Definite Integrals, Köhler Verlag, Leipzig, 1939,
  • JOURNAL, Philip J., Davis, Leonhard Euler's Integral: A Historical Profile of the Gamma Function, American Mathematical Monthly, 66, 10, 849–869, 1959, 10.2307/2309786, 2309786,
  • BOOK, Press, W. H., Teukolsky, S. A., Vetterling, W. T., Flannery, B. P., 2007, Numerical Recipes: The Art of Scientific Computing, 3rd, Cambridge University Press, New York, 978-0-521-88068-8, Section 6.1. Gamma Function,weblink
  • BOOK, O. R., Rocktäschel, Methoden zur Berechnung der Gammafunktion für komplexes Argument, Methods for Calculating the Gamma Function for Complex Arguments, Technische Universität Dresden, Technical University of Dresden, Dresden, 1922,
  • BOOK, Nico M., Temme, Special Functions: An Introduction to the Classical Functions of Mathematical Physics, John Wiley & Sons, New York, 978-0-471-11313-3, 1996,
  • BOOK, E. T. Whittaker, E. T., Whittaker, G. N., Watson, A Course of Modern Analysis, Cambridge University Press, 1927, 978-0-521-58807-2,

External links

{{Commons category|Gamma and related functions}} {{Authority control}}

- content above as imported from Wikipedia
- "gamma function" does not exist on GetWiki (yet)
- time: 5:48pm EDT - Wed, May 22 2019
[ this remote article is provided by Wikipedia ]
LATEST EDITS [ see all ]
M.R.M. Parrott