Particular values of the gamma function

From testwiki
Jump to navigation Jump to search

Template:Short description The gamma function is an important special function in mathematics. Its particular values can be expressed in closed form for integer and half-integer arguments, but no simple expressions are known for the values at rational points in general. Other fractional arguments can be approximated through efficient infinite products, infinite series, and recurrence relations.

Integers and half-integers

For positive integer arguments, the gamma function coincides with the factorial. That is,

Γ(n)=(n1)!,

and hence

Γ(1)=1,Γ(2)=1,Γ(3)=2,Γ(4)=6,Γ(5)=24,

and so on. For non-positive integers, the gamma function is not defined.

For positive half-integers, the function values are given exactly by

Γ(n2)=π(n2)!!2n12,

or equivalently, for non-negative integer values of Template:Mvar:

Γ(12+n)=(2n1)!!2nπ=(2n)!4nn!πΓ(12n)=(2)n(2n1)!!π=(4)nn!(2n)!π

where Template:Math denotes the double factorial. In particular,

Γ(12) =π 1.7724538509055160273, Template:OEIS2C
Γ(32) =12π 0.8862269254527580137, Template:OEIS2C
Γ(52) =34π 1.3293403881791370205, Template:OEIS2C
Γ(72) =158π 3.3233509704478425512, Template:OEIS2C

and by means of the reflection formula,

Γ(12) =2π 3.5449077018110320546, Template:OEIS2C
Γ(32) =43π 2.3632718012073547031, Template:OEIS2C
Γ(52) =815π 0.9453087204829418812, Template:OEIS2C

General rational argument

In analogy with the half-integer formula,

Γ(n+13)=Γ(13)(3n2)!!!3nΓ(n+14)=Γ(14)(4n3)!!!!4nΓ(n+1q)=Γ(1q)(qn(q1))!(q)qnΓ(n+pq)=Γ(pq)1qnk=1n(kq+pq)

where Template:Math denotes the Template:Mvarth multifactorial of Template:Mvar. Numerically,

Γ(13)2.6789385347077476337 Template:OEIS2C
Γ(14)3.6256099082219083119 Template:OEIS2C
Γ(15)4.5908437119988030532 Template:OEIS2C
Γ(16)5.5663160017802352043 Template:OEIS2C
Γ(17)6.5480629402478244377 Template:OEIS2C
Γ(18)7.5339415987976119047 Template:OEIS2C.

As n tends to infinity,

Γ(1n)nγ

where γ is the Euler–Mascheroni constant and denotes asymptotic equivalence.

It is unknown whether these constants are transcendental in general, but Template:Math and Template:Math were shown to be transcendental by G. V. Chudnovsky. Template:Math has also long been known to be transcendental, and Yuri Nesterenko proved in 1996 that Template:Math, Template:Math, and Template:Math are algebraically independent.

For n2 at least one of the two numbers Γ(1n) and Γ(2n) is transcendental.[1]

The number Γ(14) is related to the lemniscate constant Template:Mvar by

Γ(14)=2ϖ2π

Borwein and Zucker have found that Template:Math can be expressed algebraically in terms of Template:Mvar, Template:Math, Template:Math, Template:Math, and Template:Math where Template:Math is a complete elliptic integral of the first kind. This permits efficiently approximating the gamma function of rational arguments to high precision using quadratically convergent arithmetic–geometric mean iterations. For example:

Γ(16)=3πΓ(13)223Γ(14)=2K(12)πΓ(13)=27/9πK(14(23))3312Γ(18)Γ(38)=824(21)πK(322)Γ(18)Γ(38)=2(1+2)K(12)π4

No similar relations are known for Template:Math or other denominators.

In particular, where AGM() is the arithmetic–geometric mean, we have[2]

Γ(13)=279π233112AGM(2,2+3)13
Γ(14)=(2π)32AGM(2,1)
Γ(16)=2149313π56AGM(1+3,8)23.

Other formulas include the infinite products

Γ(14)=(2π)34k=1tanh(πk2)

and

Γ(14)=A3eGππ216k=1(112k)k(1)k

where Template:Mvar is the Glaisher–Kinkelin constant and Template:Mvar is Catalan's constant.

The following two representations for Template:Math were given by I. Mező[3]

πeπ21Γ(34)2=ik=eπ(k2k2)θ1(iπ2(2k1),eπ),

and

π21Γ(34)2=k=θ4(ikπ,eπ)e2πk2,

where Template:Math and Template:Math are two of the Jacobi theta functions.

There also exist a number of Malmsten integrals for certain values of the gamma function:[4]

1lnlnt1+t2=π4(2ln2+3lnπ4Γ(14))
1lnlnt1+t+t2=π63(8ln23ln3+8lnπ12Γ(13))

Products

Some product identities include:

r=12Γ(r3)=2π33.6275987284684357012 Template:OEIS2C
r=13Γ(r4)=2π37.8748049728612098721 Template:OEIS2C
r=14Γ(r5)=4π2517.6552850814935242483
r=15Γ(r6)=4π5340.3993191220037900785
r=16Γ(r7)=8π3793.7541682035825037970
r=17Γ(r8)=4π7219.8287780169572636207

In general:

r=1nΓ(rn+1)=(2π)nn+1

From those products can be deduced other values, for example, from the former equations for r=13Γ(r4), Γ(14) and Γ(24), can be deduced:

Γ(34)=(π2)14AGM(2,1)12

Other rational relations include

Γ(15)Γ(415)Γ(13)Γ(215)=232056575+6654
Γ(120)Γ(920)Γ(320)Γ(720)=54(1+5)2[5]
Γ(15)2Γ(110)Γ(310)=1+5271054

and many more relations for Template:Math where the denominator d divides 24 or 60.[6]

Gamma quotients with algebraic values must be "poised" in the sense that the sum of arguments is the same (modulo 1) for the denominator and the numerator.

A more sophisticated example:

Γ(1142)Γ(27)Γ(121)Γ(12)=8sin(π7)sin(π21)sin(4π21)sin(5π21)21423928713[7]

Imaginary and complex arguments

The gamma function at the imaginary unit Template:Math gives Template:OEIS2C, Template:OEIS2C:

Γ(i)=(1+i)!0.15490.4980i.

It may also be given in terms of the [[Barnes G-function|Barnes Template:Mvar-function]]:

Γ(i)=G(1+i)G(i)=elogG(i)+logG(1+i).

Curiously enough, Γ(i) appears in the below integral evaluation:[8]

0π/2{cot(x)}dx=1π2+i2log(πsinh(π)Γ(i)2).

Here {} denotes the fractional part.

Because of the Euler Reflection Formula, and the fact that Γ(z¯)=Γ¯(z), we have an expression for the modulus squared of the Gamma function evaluated on the imaginary axis:

|Γ(iκ)|2=πκsinh(πκ)

The above integral therefore relates to the phase of Γ(i).

The gamma function with other complex arguments returns

Γ(1+i)=iΓ(i)0.4980.155i
Γ(1i)=iΓ(i)0.498+0.155i
Γ(12+12i)0.81816399950.7633138287i
Γ(1212i)0.8181639995+0.7633138287i
Γ(5+3i)0.01604188279.4332932898i
Γ(53i)0.0160418827+9.4332932898i.

Other constants

The gamma function has a local minimum on the positive real axis

xmin=1.4616321449683623412626595423 Template:OEIS2C

with the value

Γ(xmin)=0.8856031944108887002788159005 Template:OEIS2C.

Integrating the reciprocal gamma function along the positive real axis also gives the Fransén–Robinson constant.

On the negative real axis, the first local maxima and minima (zeros of the digamma function) are:

Approximate local extrema of Template:Math
Template:Mvar Template:Math OEIS
Template:Val Template:Val Template:OEIS2C
Template:Val Template:0Template:Val Template:OEIS2C
Template:Val Template:Val Template:OEIS2C
Template:Val Template:0Template:Val Template:OEIS2C
Template:Val Template:Val Template:OEIS2C
Template:Val Template:0Template:Val Template:OEIS2C
Template:Val Template:Val Template:OEIS2C
Template:Val Template:0Template:Val Template:OEIS2C
Template:Val Template:Val Template:OEIS2C
Template:Val Template:0Template:Val Template:OEIS2C

The only values of Template:Math for which Template:Math are Template:Math and Template:Math... Template:OEIS2C.

See also

References

Further reading