Reciprocal gamma function

Last updated January 03, 2024

Reciprocal gamma function
1/G(z) in the complex plane, plotted using domain coloring. Reciprocal-gamma-domain-coloring.png — Reciprocal gamma function $1 / Γ(z)$ in the complex plane, plotted using domain coloring.

In mathematics, the reciprocal gamma function is the function

where $Γ(z)$ denotes the gamma function. Since the gamma function is meromorphic and nonzero everywhere in the complex plane, its reciprocal is an entire function. As an entire function, it is of order 1 (meaning that $log log | 1 / Γ(z) |$ grows no faster than $log | z |$ ), but of infinite type (meaning that $log | 1 / Γ(z) |$ grows faster than any multiple of $| z |$ , since its growth is approximately proportional to $| z | log | z |$ in the left-half plane).

The reciprocal is sometimes used as a starting point for numerical computation of the gamma function, and a few software libraries provide it separately from the regular gamma function.

Karl Weierstrass called the reciprocal gamma function the "factorielle" and used it in his development of the Weierstrass factorization theorem.

Infinite product expansion

Following from the infinite product definitions for the gamma function, due to Euler and Weierstrass respectively, we get the following infinite product expansion for the reciprocal gamma function:

{\begin{aligned}{\frac {1}{\Gamma (z)}}&=z\prod _{n=1}^{\infty }{\frac {1+{\frac {z}{n}}}{\left(1+{\frac {1}{n}}\right)^{z}}}\\{\frac {1}{\Gamma (z)}}&=ze^{\gamma z}\prod _{n=1}^{\infty }\left(1+{\frac {z}{n}}\right)e^{-{\frac {z}{n}}}\end{aligned}}

where $γ = 0.577216...$ is the Euler–Mascheroni constant. These expansions are valid for all complex numbers $z$ .

Taylor series

Taylor series expansion around 0 gives:^[1]

{\frac {1}{\ \Gamma (z)\ }}=z+\gamma \ z^{2}+\left({\frac {\gamma ^{2}}{2}}-{\frac {\pi ^{2}}{12}}\right)\ z^{3}+\left({\frac {\gamma ^{3}}{6}}-{\frac {\gamma \pi ^{2}}{12}}+{\frac {\zeta (3)}{3}}\ \right)z^{4}+\cdots \

where $γ$ is the Euler–Mascheroni constant. For $n > 2$ , the coefficient $a n$ for the $z n$ term can be computed recursively as^[2]^[3]

a_{n}={\frac {\ {a_{2}\ a_{n-1}+\sum _{j=2}^{n-1}(-1)^{j+1}\ \zeta (j)\ a_{n-j}}\ }{n-1}}={\frac {\ \gamma \ a_{n-1}-\zeta (2)\ a_{n-2}+\zeta (3)\ a_{n-3}-\cdots \ }{n-1}}

where $ζ$ is the Riemann zeta function. An integral representation for these coefficients was recently found by Fekih-Ahmed (2014):^[3]

a_{n}={\frac {(-1)^{n}}{\pi n!}}\int _{0}^{\infty }e^{-t}\ \operatorname {Im} {\Bigl [}\ {\bigl (}\log(t)-i\pi {\bigr )}^{n}\ {\Bigr ]}\ \mathrm {d} t~.

For small values, these give the following values:

$n$	$a n$
1	+1.0000000000000000000000000000000000000000
2	+0.5772156649015328606065120900824024310422
3	−0.6558780715202538810770195151453904812798
4	−0.0420026350340952355290039348754298187114
5	+0.1665386113822914895017007951021052357178
6	−0.0421977345555443367482083012891873913017
7	−0.0096219715278769735621149216723481989754
8	+0.0072189432466630995423950103404465727099
9	−0.0011651675918590651121139710840183886668
10	−0.0002152416741149509728157299630536478065
11	+0.0001280502823881161861531986263281643234
12	−0.0000201348547807882386556893914210218184
13	−0.0000012504934821426706573453594738330922
14	+0.0000011330272319816958823741296203307449
15	−0.0000002056338416977607103450154130020573
16	+0.0000000061160951044814158178624986828553
17	+0.0000000050020076444692229300556650480600
18	−0.0000000011812745704870201445881265654365
19	+0.0000000001043426711691100510491540332312
20	+0.0000000000077822634399050712540499373114
21	−0.0000000000036968056186422057081878158781
22	+0.0000000000005100370287454475979015481323
23	−0.0000000000000205832605356650678322242954
24	−0.0000000000000053481225394230179823700173
25	+0.0000000000000012267786282382607901588938
26	−0.0000000000000001181259301697458769513765
27	+0.0000000000000000011866922547516003325798
28	+0.0000000000000000014123806553180317815558
29	−0.0000000000000000002298745684435370206592
30	+0.0000000000000000000171440632192733743338

Fekih-Ahmed (2014)^[3] also gives an approximation for $a_{n}$ :

a_{n}\approx {\frac {(-1)^{n}}{\ (n-1)!\ }}\ {\sqrt {{\frac {2}{\ \pi n\ }}\ }}\ \operatorname {Im} \left({\frac {\ z_{0}^{\left(1/2-n\right)}\ e^{-nz_{0}}\ }{\sqrt {1+z_{0}\ }}}\right)\ ,

where $z_{0}=-{\frac {1}{n}}\exp \!{\Bigl (}W_{-1}(-n){\Bigr )}\ ,$ and $W_{-1}$ is the minus-first branch of the Lambert W function.

The Taylor expansion around $1$ has the same (but shifted) coefficients, i.e.:

{\frac {1}{\Gamma (1+z)}}={\frac {1}{z\Gamma (z)}}=1+\gamma \ z+\left({\frac {\gamma ^{2}}{2}}-{\frac {\pi ^{2}}{12}}\right)\ z^{2}+\left({\frac {\gamma ^{3}}{6}}-{\frac {\gamma \pi ^{2}}{12}}+{\frac {\zeta (3)}{3}}\ \right)z^{3}+\cdots \

(the reciprocal of Gauss' pi-function).

Asymptotic expansion

As $| z |$ goes to infinity at a constant $arg(z)$ we have:

\ln(1/\Gamma (z))\sim -z\ln(z)+z+{\tfrac {1}{2}}\ln \left({\frac {z}{2\pi }}\right)-{\frac {1}{12z}}+{\frac {1}{360z^{3}}}-{\frac {1}{1260z^{5}}}\qquad {\text{for}}~\left|\arg(z)\right|<\pi

Contour integral representation

An integral representation due to Hermann Hankel is

{\frac {1}{\Gamma (z)}}={\frac {i}{2\pi }}\oint _{H}(-t)^{-z}e^{-t}\,dt,

where $H$ is the Hankel contour, that is, the path encircling 0 in the positive direction, beginning at and returning to positive infinity with respect for the branch cut along the positive real axis. According to Schmelzer & Trefethen,^[4] numerical evaluation of Hankel's integral is the basis of some of the best methods for computing the gamma function.

Integral representations at the positive integers

For positive integers $n\geq 1$ , there is an integral for the reciprocal factorial function given by^[5]

{\frac {1}{n!}}={\frac {1}{2\pi }}\int _{-\pi }^{\pi }e^{-nit}e^{e^{it}}\ dt.

Similarly, for any real $c>0$ and $z\in \mathbb {C}$ we have the next integral for the reciprocal gamma function along the real axis in the form of:^[6]

{\frac {1}{\Gamma (z)}}={\frac {1}{2\pi }}\int _{-\infty }^{\infty }(c+it)^{-z}e^{c+it}dt,

where the particular case when $z=n+1/2$ provides a corresponding relation for the reciprocal double factorial function, ${\frac {1}{(2n-1)!!}}={\frac {\sqrt {\pi }}{2^{n}\cdot \Gamma \left(n+{\frac {1}{2}}\right)}}.$

Integral along the real axis

Integration of the reciprocal gamma function along the positive real axis gives the value

\int _{0}^{\infty }{\frac {1}{\Gamma (x)}}\,dx\approx 2.80777024,

which is known as the Fransén–Robinson constant.

We have the following formula (^[7] chapter 9, exercise 100)

\int _{0}^{\infty }{\dfrac {a^{x}}{\Gamma (x)}}\,dx=ae^{a}+a\int _{0}^{\infty }{\dfrac {e^{-ax}}{\log ^{2}(x)+\pi ^{2}}}\,dx

Related Research Articles

In mathematics, the gamma function 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. For every positive integer $n$ ,

The Riemann zeta function or Euler–Riemann zeta function, denoted by the Greek letter $ζ$ (zeta), is a mathematical function of a complex variable defined as

<span class="mw-page-title-main">Euler's constant</span> Constant value used in mathematics

Euler's constant is a mathematical constant, usually denoted by the lowercase Greek letter gamma, defined as the limiting difference between the harmonic series and the natural logarithm, denoted here by $log$ :

<span class="mw-page-title-main">Harmonic number</span> Sum of the first n whole number reciprocals; 1/1 + 1/2 + 1/3 + ... + 1/n

In mathematics, the $n$ -th harmonic number is the sum of the reciprocals of the first $n$ natural numbers:

In mathematics, the Clausen function, introduced by Thomas Clausen (1832), is a transcendental, special function of a single variable. It can variously be expressed in the form of a definite integral, a trigonometric series, and various other forms. It is intimately connected with the polylogarithm, inverse tangent integral, polygamma function, Riemann zeta function, Dirichlet eta function, and Dirichlet beta function.

In mathematics, the polygamma function of order $m$ is a meromorphic function on the complex numbers $defined as the (m + 1) th derivative of the logarithm of the gamma function:$

<span class="mw-page-title-main">Digamma function</span> Mathematical function

In mathematics, the digamma function is defined as the logarithmic derivative of the gamma function:

<span class="mw-page-title-main">Hurwitz zeta function</span> Special function in mathematics

In mathematics, the Hurwitz zeta function is one of the many zeta functions. It is formally defined for complex variables $s$ with $Re(s) > 1$ and $a \neq 0, -1, -2, \dots$ by

In mathematics, in the area of analytic number theory, the Dirichlet eta function is defined by the following Dirichlet series, which converges for any complex number having real part > 0:

In mathematics, the polylogarithm (also known as Jonquière's function, for Alfred Jonquière) is a special function $Li s (z)$ of order $s$ and argument $z$ . Only for special values of $s$ does the polylogarithm reduce to an elementary function such as the natural logarithm or a rational function. In quantum statistics, the polylogarithm function appears as the closed form of integrals of the Fermi–Dirac distribution and the Bose–Einstein distribution, and is also known as the Fermi–Dirac integral or the Bose–Einstein integral. In quantum electrodynamics, polylogarithms of positive integer order arise in the calculation of processes represented by higher-order Feynman diagrams.

In mathematics, the Lerch zeta function, sometimes called the Hurwitz–Lerch zeta function, is a special function that generalizes the Hurwitz zeta function and the polylogarithm. It is named after Czech mathematician Mathias Lerch, who published a paper about the function in 1887.

<span class="mw-page-title-main">Stieltjes constants</span>

In mathematics, the Stieltjes constants are the numbers $that occur in the Laurent series expansion of the Riemann zeta function:$

In mathematics, the Barnes G-functionG(z) is a function that is an extension of superfactorials to the complex numbers. It is related to the gamma function, the K-function and the Glaisher–Kinkelin constant, and was named after mathematician Ernest William Barnes. It can be written in terms of the double gamma function.

In mathematics, the explicit formulae for L-functions are relations between sums over the complex number zeroes of an L-function and sums over prime powers, introduced by Riemann (1859) for the Riemann zeta function. Such explicit formulae have been applied also to questions on bounding the discriminant of an algebraic number field, and the conductor of a number field.

In mathematics, Apéry's constant is the sum of the reciprocals of the positive cubes. That is, it is defined as the number

In q-analog theory, the $-gamma function$ , or basic gamma function, is a generalization of the ordinary gamma function closely related to the double gamma function. It was introduced by Jackson (1905). It is given by

Carl Johan Malmsten was a Swedish mathematician and politician. He is notable for early research into the theory of functions of a complex variable, for the evaluation of several important logarithmic integrals and series, for his studies in the theory of Zeta-function related series and integrals, as well as for helping Mittag-Leffler start the journal Acta Mathematica. Malmsten became Docent in 1840, and then, Professor of mathematics at the Uppsala University in 1842. He was elected a member of the Royal Swedish Academy of Sciences in 1844. He was also a minister without portfolio in 1859–1866 and Governor of Skaraborg County in 1866–1879.

<span class="mw-page-title-main">Wrapped Cauchy distribution</span>

In probability theory and directional statistics, a wrapped Cauchy distribution is a wrapped probability distribution that results from the "wrapping" of the Cauchy distribution around the unit circle. The Cauchy distribution is sometimes known as a Lorentzian distribution, and the wrapped Cauchy distribution may sometimes be referred to as a wrapped Lorentzian distribution.

In mathematics, Ramanujan's Master Theorem, named after Srinivasa Ramanujan, is a technique that provides an analytic expression for the Mellin transform of an analytic function.

In mathematical physics, the Wu–Sprung potential, named after Hua Wu and Donald Sprung, is a potential function in one dimension inside a Hamiltonian $with the potential defined by solving a non-linear integral equation defined by the Bohr-Sommerfeld quantization conditions involving the spectral staircase, the energies and the potential .$

References

↑ Weisstein, Eric W. "Gamma function". mathworld.wolfram.com. Retrieved 2021-06-15.
↑ Wrench, J.W. (1968). "Concerning two series for the gamma function". Mathematics of Computation. 22 (103): 617–626. doi: 10.1090/S0025-5718-1968-0237078-4 . S2CID 121472614. and
Wrench, J.W. (1973). "Erratum: Concerning two series for the gamma function". Mathematics of Computation. 27 (123): 681–682. doi: 10.1090/S0025-5718-1973-0319344-9 .
1 2 3 Fekih-Ahmed, L. (2014). "On the power series expansion of the reciprocal gamma function". HAL archives.
↑ Schmelzer, Thomas; Trefethen, Lloyd N. (2007). "Computing the Gamma function using contour integrals and rational approximations". SIAM Journal on Numerical Analysis. Society for Industrial and Applied Mathematics. 45 (2): 558–571. doi:10.1137/050646342.; "Copy on Trefethen's academic website" (PDF). Mathematics, Oxford, UK. Retrieved 2020-08-03.; "Link to two other copies". CiteSeerX 10.1.1.210.299 .
↑ Graham, Knuth, and Patashnik (1994). Concrete Mathematics. Addison-Wesley. p. 566.{{cite book}}: CS1 maint: multiple names: authors list (link)
↑ Schmidt, Maxie D. (2019-05-19). "A Short Note on Integral Transformations and Conversion Formulas for Sequence Generating Functions". Axioms. 8 (2): 62. arXiv: 1809.03933 . doi: 10.3390/axioms8020062 .
↑ Henri Cohen (2007). Number Theory Volume II: Analytic and Modern Tools. Graduate Texts in Mathematics. Vol. 240. doi:10.1007/978-0-387-49894-2. ISBN 978-0-387-49893-5. ISSN 0072-5285.

Mette Lund, An integral for the reciprocal Gamma function
Milton Abramowitz & Irene A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables
Eric W. Weisstein, Gamma Function , MathWorld

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Weisstein, Eric W. "Gamma function". mathworld.wolfram.com. Retrieved 2021-06-15.

[2] Wrench, J.W. (1968). "Concerning two series for the gamma function". Mathematics of Computation. 22 (103): 617–626. doi: 10.1090/S0025-5718-1968-0237078-4 . S2CID 121472614. and
Wrench, J.W. (1973). "Erratum: Concerning two series for the gamma function". Mathematics of Computation. 27 (123): 681–682. doi: 10.1090/S0025-5718-1973-0319344-9 .

[Fekih-Ahmed2014-3] 1 2 3 Fekih-Ahmed, L. (2014). "On the power series expansion of the reciprocal gamma function". HAL archives.

[4] Schmelzer, Thomas; Trefethen, Lloyd N. (2007). "Computing the Gamma function using contour integrals and rational approximations". SIAM Journal on Numerical Analysis. Society for Industrial and Applied Mathematics. 45 (2): 558–571. doi:10.1137/050646342.; "Copy on Trefethen's academic website" (PDF). Mathematics, Oxford, UK. Retrieved 2020-08-03.; "Link to two other copies". CiteSeerX 10.1.1.210.299 .

[5] Graham, Knuth, and Patashnik (1994). Concrete Mathematics. Addison-Wesley. p. 566.{{cite book}}: CS1 maint: multiple names: authors list (link)

[6] Schmidt, Maxie D. (2019-05-19). "A Short Note on Integral Transformations and Conversion Formulas for Sequence Generating Functions". Axioms. 8 (2): 62. arXiv: 1809.03933 . doi: 10.3390/axioms8020062 .

[7] Henri Cohen (2007). Number Theory Volume II: Analytic and Modern Tools. Graduate Texts in Mathematics. Vol. 240. doi:10.1007/978-0-387-49894-2. ISBN 978-0-387-49893-5. ISSN 0072-5285.

[1]

[2]

[3]

[4]

[5]

[6]

[7]