Inverse gamma function: Difference between revisions

Revision as of 07:01, 3 May 2023

Latest revision as of 19:18, 20 August 2025

(44 intermediate revisions by 21 users not shown)

Line 1:

{{Short description|Inverse of the gamma function}}

{{Draft topics|mathematics}}

{{AfC topic|stem}}

{{AfC submission|||ts=20230503070148|u=Onlineuser577215|ns=118}}

{{AFC submission|d|reason|Please discuss at the talk page for [[Gamma function]] whether this should be a section on that page or split off.|u=Onlineuser577215|ns=118|decliner=AngusWOOF|declinets=20230427192227|ts=20230427182236}} <!-- Do not remove this line! -->

{{AFC submission|d|v|u=Onlineuser577215|ns=118|decliner=Majash2020|declinets=20230412023636|small=yes|ts=20221206210833}} <!-- Do not remove this line! -->

{{AFC comment|1=Users on the Gamma function page suggested that this should have its own page instead of being added as a new section on that page. See the talk page on [[Gamma function]]. [[User:Onlineuser577215|Onlineuser577215]] ([[User talk:Onlineuser577215|talk]]) 8:58, 3 May 2023 (UTC)}}

----

{{Distinguish|Inverse-gamma distribution|Reciprocal gamma function}}

{{multiple image

| total_width = 500

| image1 = Inverse Gamma Function.png

| caption1 = Graph of an inverse gamma function

| image2 = Inverse gamma function in complex plane.png

| caption2 = Plot of inverse gamma function in the complex plane

}}

In [[mathematics]], the '''inverse gamma function''' <math>\Gamma^{-1}(x)</math> is the [[inverse function]] of the [[gamma function]]. In other words, it<math>y is= the function\Gamma^{-1}(x)</math> satisfyingwhenever <math display="inline">\Gamma(y)=x</math>. For example, <math>\Gamma^{-1}(24)=5</math> .<ref>{{Citecite journal |lastlast1=Borwein, |first1= Jonathan M. |last2=Corless |first2=Robert M.|title=Gamma and Factorial in the Monthly |journal=The American Mathematical Monthly |year=2017 |volume=125 |issue=5 |pages= 400–424 |doi= 10.1080/00029890.2018.1420983 |arxiv=1703.05349 |jstor=48663320 |s2cid=119324101}}</ref>. Usually, the inverse gamma function refers to the principal branch with ___domain on the real interval <math>\left([\Gamma(beta, +\alphainfty\right)=</math> 0.8856031...and image on the real interval <math>\left[\alpha, +\infty\right)</math>, where <math>\alphabeta =1.4616321.. 0.8856031\ldots</math><ref>{{oeis|A030171}}</ref> is the uniqueminimum positivevalue numberof suchthe thatgamma function on the positive real axis and <math>\Psialpha = \Gamma^{-1}(x\beta) =0 1.4616321\ldots</math><ref>{{oeis|A030169}}</ref> is the ___location of that minimum.<ref>{{cite journal |last1=Uchiyama |first1=MITSURUMitsuru |title=The principal inverse of the gamma function |date=April 2012 |url= https://www.jstor.org/stable/41505586 |journal=Proceedings of the American Mathematical Society|volume=140 |issue=4 |pages=1347 |doi= 10.1090/S0002-9939-2011-11023-2

|jstor=41505586 |s2cid=85549521 |doi-access=free }}</ref>

== Definition ==

In [[mathematics]], the inverse gamma function <math>\Gamma^{-1}(x)</math> is the [[inverse function]] of the [[gamma function]]. In other words, it is the function satisfying <math display="inline">\Gamma(y)=x</math>. For example, <math>\Gamma^{-1}(24)=5</math> <ref>{{Cite journal |last=Borwein, Corless |title=Gamma and Factorial in the Monthly |journal=|year=2017 |arxiv=1703.05349 }}</ref>. Usually, the inverse gamma function refers to the principal branch on the interval <math>\left(\Gamma(\alpha)= 0.8856031..., \infty\right)</math> where <math>\alpha=1.4616321...</math> is the unique positive number such that <math>\Psi(x)=0</math> <ref>{{cite journal |last1=Uchiyama |first1=MITSURU |title=The principal inverse of the gamma function |date=April 2012 |url= https://www.jstor.org/stable/41505586 |journal=Proceedings of the American Mathematical Society|volume=140 |issue=4 |pages=1347 |doi= 10.1090/S0002-9939-2011-11023-2

The inverse gamma function may be defined by the following integral representation<ref>{{cite journal |last1=PedersePedersen |first1=Henrik |title="Inverses of gamma functions" |journal=Constructive Approximation |date=9 SepSeptember 2013 |pagesvolume=7 |issue=2 |pages=251–267 |doi=10.1007/s00365-014-9239-1 }}<|arxiv=1309.2167 |s2cid=253898042 |url=https:/ref><math>\Gamma^{/link.springer.com/article/10.1007/s00365-1}(x)=a+bx+\int_{014-\infty}^{\Gamma(\alpha)}\left(\frac{1}{x-t}-\frac{t}{t^{2}9239-1}\right)d\mu(t) }</mathref>

|jstor=41505586 |s2cid=85549521 |access-date=20 March 2023}}</ref> where <math>\Psi(x)</math> is the [[digamma function]].

<math display="block">\Gamma^{-1}(x)=a+bx+\int_{-\infty}^{\Gamma(\alpha)}\left(\frac{1}{x-t}-\frac{t}{t^{2}-1}\right)d\mu(t)\,,</math>

where <math>\mu (t)</math> is a [[Borel measure]] such that <math display="block">\int_{-\infty}^{\Gamma\left(\alpha\right)}\left(\frac{1}{t^{2}+1}\right)d\mu(t)<\infty \,,</math> and <math>a</math> and <math>b</math> are real numbers with <math>b \geqq 0</math>.

==== DefinitionApproximation ====

To compute the branches of the inverse gamma function one can first compute the [[Taylor series]] of <math>\Gamma(x)</math> near <math>\alpha</math>. The series can then be truncated and inverted, which yields successively better approximations to <math>\Gamma^{-1}(x)</math>. For instance, we have the quadratic approximation:<ref>{{cite journalconference |first1=Robert M.|last1=Corless |last2first2=Folitse Komla|last2=Amenyou |last3=Jeffrey |first3=David |book-title=2017 19th International Symposium on Symbolic and Numeric Algorithms for Scientific Computing (SYNASC) |title=Properties and Computation of the Functional Inverse of Gamma |journalconference=International Symposium on Symbolic and Numeric Algorithms for Scientific Computing (SYNASC) |date=2017 |pages=65 |doi=10.1109/SYNASC.2017.00020|isbn=978-1-5386-2626-9 |s2cid=53287687 }}</ref>

The inverse gamma function may be defined by the following integral representation<ref>{{cite journal |last1=Pederse |first1=Henrik |title=Inverses of gamma functions |journal=Constructive Approximation |date=9 Sep 2013 |pages=7 |doi=10.1007/s00365-014-9239-1 }}</ref><math>\Gamma^{-1}(x)=a+bx+\int_{-\infty}^{\Gamma(\alpha)}\left(\frac{1}{x-t}-\frac{t}{t^{2}-1}\right)d\mu(t) </math>

Where <math display="block"> \int_Gamma^{-1}\infty}^left(x\right)\approx\alpha+\sqrt{\frac{2\left(x-\Gamma\left(\alpha\right)}\left(\frac{1right)}{t\psi^{2}+\left(1} \right)d}\muleft(t)<\infty</math>,alpha and a and b are real numbers with <math>b\geqq0</math>, and <math>right)\mu Gamma\left(t\alpha\right)}}.</math> is the [[Borel measure|Borel Meausure]].

where <math> \psi^{\left(1 \right)} \left(x \right)</math> is the [[trigamma function]]. The inverse gamma function also has the following [[asymptotic formula]]<ref>{{Citecite journalthesis |lasttype=MS |last1=Amenyou |firstfirst1=Folitse Komla |last2=Jeffrey |first2=David |title="Properties and Computation of the Inverseinverse of the Gamma functionFunction" |date=2018 |pages=28 |url=https://ir.lib.uwo.ca/cgi/viewcontent.cgi?article=7340&context=etd |journal=Western:Graduate & Postdoctoral Studies}}</ref>

==== Approximation ====

<math display="block"> \Gamma^{-1}(x)\sim\frac{1}{2}+\frac{\ln\left(\frac{x}{\sqrt{2\pi}}\right)}{W_{0}\left(e^{-1}\ln\left(\frac{x}{\sqrt{2\pi}}\right)\right)}\,,</math>

To compute the branches of the inverse gamma function one can first compute the Taylor series of <math>\Gamma(x)</math> near <math>\alpha</math>. The series can then be truncated and inverted, which yields successively better approximations to <math>\Gamma^{-1}(x)</math>. For instance, we have the quadratic approximation<ref>{{cite journal |last1=Corless |last2=Folitse |last3=Jeffrey |title=Properties and Computation of the Functional Inverse of Gamma |journal=SYNASC |date=2017 |pages=65 |doi=10.1109/SYNASC.2017.00020}}</ref>

Wherewhere <math>W_0(x)</math> is the [[Lambert W function]]. The formula is found by inverting the [[Stirling's approximation|Stirling approximation]], and so can also be expanded into an asymptotic series.

=== Series expansion ===

<math>

To obtain a series expansion of the inverse gamma function one can first compute the series expansion of the [[reciprocal gamma function]] <math>\frac{1}{\Gamma(x)}</math> near the poles at the negative integers, and then invert the series.

\Gamma^{-1}\left(x\right)\approx\alpha+\sqrt{\frac{2\left(x-\Gamma\left(\alpha\right)\right)}{\Psi\left(1,\ \alpha\right)\Gamma\left(\alpha\right)}}.</math>

Setting <math>z=\frac{1}{x}</math> then yields, for the ''n'' th branch <math>\Gamma_{n}^{-1}(z)</math> of the inverse gamma function (<math>n\ge 0</math>)<ref>{{Cite journal |last1=Couto |first1=Ana Carolina Camargos |last2=Jeffrey |first2=David |last3=Corless |first3=Robert |date=November 2020 |title=The Inverse Gamma Function and its Numerical Evaluation |url=https://www.maplesoft.com/mapleconference/2020/highlights.aspx |at=Section 8 |journal=Maple Conference Proceedings}}</ref>

The inverse gamma function also has the following [[asymptotic formula]]<ref>{{Cite journal |last=Amenyou |first=Komla |title=Properties and Computation of the Inverse of the Gamma function |url=https://ir.lib.uwo.ca/cgi/viewcontent.cgi?article=7340&context=etd |journal=Western:Graduate & Postdoctoral Studies}}</ref>

<math display="block"> \Gamma_{n}^{-1}(z)=-n+\frac{\left(-1\right)^{n}}{n!z}+\frac{\psi^{(0)}\left(n+1\right)}{\left(n!z\right)^2}+\frac{\left(-1\right)^{n}\left(\pi^{2}+9\psi^{(0)}\left(n+1\right)^{2}-3\psi^{(1)}\left(n+1\right)\right)}{6\left(n!z\right)^3}+O\left(\frac{1}{z^{4}}\right)\,,</math>

where <math>\psi^{(n)}(x)</math> is the [[polygamma function]].

<math>\Gamma^{-1}(x)\sim\frac{1}{2}+\frac{\ln\left(\frac{x}{\sqrt{2\pi}}\right)}{W_{0}\left(e^{-1}\ln\left(\frac{x}{\sqrt{2\pi}}\right)\right)}</math>

Where <math>W_0(x)</math> is the [[Lambert W function]]. The formula is found by inverting the [[Stirling's approximation|Stirling approximation]], and so can also be expanded into an asymptotic series.

== References ==

<!-- Inline citations added to your article will automatically display here. See en.wikipedia.org/wiki/WP:REFB for instructions on how to add citations. -->

{{reflist}}

[[Category:Gamma and related functions]]

{{Draft topics|mathematics-stub}}