Articulo de referencia

Función gamma

\\Gamma(z) = \\int_0^\\infty t^{z-1} e^{-t}\\,dt \n| fields_of_application = Calculus, mathematical analysis, statistics, physics\n}}"}},"i":0}}]}"> En matemáticas , la función ...

En matemáticas , la función gamma (representada por Γ{\displaystyle \Gamma }La función gamma ( letra griega mayúscula gamma ) es la extensión más común de la función factorial a los números complejos . Estudiada por primera vez por Daniel Bernoulli , la función gammaΓ(z){\displaystyle \Gamma (z)}está definido para todos los números complejosz{\displaystyle z}excepto los enteros no positivos, yΓ(norte)=(norte1)¡{\displaystyle \Gamma (n)=(n-1)!}para cada entero positivonorte{\displaystyle n}La función gamma se puede definir mediante una integral impropia convergente para números complejos con parte real positiva: Γ(z)=0tz1mitdt, (z)>0.{\displaystyle \Gamma (z)=\int _{0}^{\infty }t^{z-1}e^{-t}\,dt,\ \qquad \Re (z)>0.}La función gamma se define entonces en el plano complejo como la continuación analítica de esta función integral: es una función meromorfa que es holomorfa excepto en cero y en los enteros negativos, donde tiene polos simples .

Dado que la función gamma no tiene ceros, su recíproca1Γ{\displaystyle {\frac {1}{\Gamma }}}es una función completa . De hecho, la función gamma corresponde a la transformada de Mellin del decaimiento exponencial : Γ(z)=METRO{miincógnita}(z).{\displaystyle \Gamma (z)={\mathcal {M}}\{e^{-x}\}(z)\,.}

Existen otras extensiones de la función factorial, pero la función gamma es la más popular y útil. Aparece como factor en diversas funciones de distribución de probabilidad y otras fórmulas en los campos de la probabilidad , la estadística , la teoría analítica de números y la combinatoria .

Motivación

Γ(incógnita+1){\displaystyle \Gamma (x+1)}interpola la función factorial a valores no enteros.

La función gamma puede considerarse una solución al problema de interpolación para encontrar una curva suave.y=F(incógnita){\displaystyle y=f(x)}que conecta los puntos de la secuencia factorial:(incógnita,y)=(norte,norte¡){\displaystyle (x,y)=(n,n!)}para todos los valores enteros positivos de norte{\displaystyle n} . La fórmula simple para el factorial,incógnita¡=123incógnita{\displaystyle x!=1\cdot 2\cdot 3\cdots x}es válido solo cuandoincógnita{\displaystyle x}es un entero positivo, y ninguna función elemental tiene esta propiedad, pero una buena solución es la función gamma .F(incógnita)=Γ(incógnita+1){\displaystyle f(x)=\Gamma (x+1)} . [ 1 ]

La función gamma no solo es suave sino analítica (excepto en los enteros no positivos), y puede definirse de varias maneras explícitas. Sin embargo, no es la única función analítica que extiende el factorial, ya que se puede agregar cualquier función analítica que sea cero en los enteros positivos, comokpecado(metroπincógnita){\displaystyle k\sin(m\pi x)}para un número enterometro{\displaystyle m} . [ 1 ] Dicha función se conoce como función pseudogamma , siendo la más famosa la función de Hadamard . [ 2 ]

La función gamma, Γ( z ) en azul, se representa junto con Γ( z ) + sin(πz ) en verde. Nótese la intersección en los enteros positivos. Ambas son extensiones válidas de los factoriales a una función meromorfa en el plano complejo.

Un requisito más restrictivo es la ecuación funcional que interpola el factorial desplazado .F(norte)=(norte1)¡{\displaystyle f(n)=(n-1)!}: [ 3 ] [ 4 ]F(incógnita+1)=incógnitaF(incógnita)  a pesar de incógnita>0,F(1)=1.{\displaystyle f(x+1)=xf(x)\ {\text{ para todo }}x>0,\qquad f(1)=1.}

Pero esto aún no proporciona una solución única, ya que permite la multiplicación por cualquier función periódica.gramo(incógnita){\displaystyle g(x)}congramo(incógnita)=gramo(incógnita+1){\displaystyle g(x)=g(x+1)}ygramo(0)=1{\displaystyle g(0)=1} , como por ejemplogramo(incógnita)=mikpecado(metroπincógnita){\displaystyle g(x)=e^{k\sin(m\pi x)}}para cualquier entero m .

Una forma de resolver la ambigüedad es el teorema de Bohr-Mollerup , que muestra queF(incógnita)=Γ(incógnita){\displaystyle f(x)=\Gamma (x)}es la única función de interpolación para el factorial, definida sobre los reales positivos, que es logarítmicamente convexa , [ 5 ] lo que significa quey=registroF(incógnita){\displaystyle y=\log f(x)}es convexa . [ 6 ]

Definición

Definición principal

La notaciónΓ(z){\displaystyle \Gamma (z)}se debe a Legendre . [ 1 ] Si la parte real del número complejoz{\displaystyle z}es estrictamente positivo ( (z)>0{\displaystyle \Re (z)>0} ), entonces la integralΓ(z)=0tz1mitdt{\displaystyle \Gamma (z)=\int _{0}^{\infty }t^{z-1}e^{-t}\,dt}converge absolutamente y se conoce como la integral de Euler de segunda especie . (La integral de Euler de primera especie es la función beta . [ 1 ] )

Valor absoluto (vertical) y argumento (tono) de la función gamma en el plano complejo.
Valor absoluto (vertical) y argumento (tono) de la función gamma en el plano complejo.

El valorΓ(1){\displaystyle \Gamma (1)}se puede calcular comoΓ(1)=0t11mitdt=0mitdt=1.{\displaystyle {\begin{aligned}\Gamma (1)&=\int _{0}^{\infty }t^{1-1}e^{-t}\,dt\\[6pt]&=\int _{0}^{\infty }e^{-t}\,dt=1.\end{aligned}}}

Integrando por partes , se observa queΓ(z+1)=0tzmitdt=[tzmit]0+0ztz1mitdt.{\displaystyle {\begin{aligned}\Gamma (z+1)&=\int _{0}^{\infty }t^{z}e^{-t}\,dt\\[6pt]&={\Bigl [}-t^{z}e^{-t}{\Bigr ]}_{0}^{\infty }+\int _{0}^{\infty }z\,t^{z-1}e^{-t}\,dt.\end{aligned}}}Reconociendo quetzmit0{\displaystyle -t^{z}e^{-t}\to 0}comot0{\displaystyle t\to 0} (siempre y cuando(z)>0{\displaystyle \Re (z)>0}) y comot{\displaystyle t\to \infty } ,Γ(z+1)=z0tz1mitdt=zΓ(z).{\displaystyle {\begin{aligned}\Gamma (z+1)&=z\int _{0}^{\infty }t^{z-1}e^{-t}\,dt\\[6pt]&=z\,\Gamma (z).\end{aligned}}}

Así hemos demostrado queΓ(norte)=(norte1)¡{\displaystyle \Gamma (n)=(n-1)!}para cualquier entero positivo n por inducción .

La identidadzΓ(z)=Γ(z+1){\textstyle z\Gamma (z)=\Gamma (z+1)}se puede utilizar (o, obteniendo el mismo resultado, se puede utilizar la continuación analítica ) para extender de forma única la formulación integral paraΓ(z){\displaystyle \Gamma (z)}a una función meromorfa definida para todos los números complejosz{\displaystyle z} , excepto los enteros menores o iguales a cero. [ 1 ] Es esta versión extendida la que comúnmente se conoce como la función gamma. [ 1 ]

Definiciones alternativas

Existen muchas definiciones equivalentes.

La definición de Euler como un producto infinito

Para un número entero fijometro{\displaystyle m} , como el enteronorte{\displaystyle n}aumenta, tenemos que [ 7 ]límitenortenorte¡(norte+1)metro(norte+metro)¡=1.{\displaystyle \lim _{n\to \infty }{\frac {n!\,\left(n+1\right)^{m}}{(n+m)!}}=1\,.}

Simetro{\displaystyle m}Si no es un número entero, entonces esta ecuación carece de sentido, ya que en esta sección aún no se ha definido el factorial de un no entero. Sin embargo, para definir la función Gamma para no enteros, supongamos que esta ecuación sigue siendo válida cuandometro{\displaystyle m}se reemplaza por un número complejo arbitrario .z{\displaystyle z}:

límitenortenorte¡(norte+1)z(norte+z)¡=1.{\displaystyle \lim _{n\to \infty }{\frac {n!\,\left(n+1\right)^{z}}{(n+z)!}}=1\,.} Multiplicando ambos lados por(z1)¡{\displaystyle (z-1)!}da (z1)¡=1zlímitenortenorte¡z¡(norte+z)¡(norte+1)z=1zlímitenorte(12norte)1(1+z)(norte+z)(2132norte+1norte)z=1znorte=1[11+znorte(1+1norte)z].{\displaystyle {\begin{aligned}(z-1)!&={\frac {1}{z}}\lim _{n\to \infty }n!{\frac {z!}{(n+z)!}}(n+1)^{z}\\[6pt]&={\frac {1}{z}}\lim _{n\to \infty }(1\cdot 2\cdots n){\frac {1}{(1+z)\cdots (n+z)}}\left({\frac {2}{1}}\cdot {\frac {3}{2}}\cdots {\frac {n+1}{n}}\right)^{z}\\[6pt]&={\frac {1}{z}}\prod _{n=1}^{\infty }\left[{\frac {1}{1+{\frac {z}{n}}}}\left(1+{\frac {1}{n}}\right)^{z}\right].\end{aligned}}}Este producto infinito , que se debe a Euler, [ 8 ] converge para todos los números complejos.z{\displaystyle z}excepto los enteros no positivos, que fallan debido a una división por cero. De hecho, la suposición anterior produce una definición única deΓ(z){\displaystyle \Gamma (z)}como(z1)¡{\displaystyle (z-1)!}.

Intuitivamente, esta fórmula indica queΓ(z){\displaystyle \Gamma (z)}es aproximadamente el resultado de la computaciónΓ(norte+1)=norte¡{\displaystyle \Gamma (n+1)=n!}para algún número entero grandenorte{\displaystyle n} , multiplicando por(norte+1)z{\displaystyle (n+1)^{z}}aproximarΓ(norte+z+1){\displaystyle \Gamma (n+z+1)}y luego utilizando la relaciónΓ(incógnita+1)=incógnitaΓ(incógnita){\displaystyle \Gamma (x+1)=x\Gamma (x)}hacia atrásnorte+1{\displaystyle n+1}tiempos para obtener una aproximación para Γ(z){\displaystyle \Gamma (z)} , y además que esta aproximación se vuelve exacta a medida que n aumenta hasta el infinito.

El producto infinito para el recíproco1Γ(z)=znorte=1[(1+znorte)/(1+1norte)z]{\displaystyle {\frac {1}{\Gamma (z)}}=z\prod _{n=1}^{\infty }\left[\left(1+{\frac {z}{n}}\right)/{\left(1+{\frac {1}{n}}\right)^{z}}\right]} es una función entera que converge para cada número complejo .z{\displaystyle z}.

Definición de Weierstrass

La definición de la función gamma debida a Weierstrass también es válida para todos los números complejos. z{\displaystyle z}excepto los enteros no positivos: Γ(z)=miγzznorte=1(1+znorte)1miz/norte,{\displaystyle \Gamma (z)={\frac {e^{-\gamma z}}{z}}\prod _{n=1}^{\infty }\left(1+{\frac {z}{n}}\right)^{-1}e^{z/n},} dóndeγ0,577216{\displaystyle \gamma \approx 0.577216}es la constante de Euler-Mascheroni . [ 1 ] Este es el producto de Hadamard de1Γ(z){\displaystyle \textstyle {\frac {1}{\Gamma (z)}}}en una forma reescrita.

Propiedades

General

Además de la propiedad fundamental mencionada anteriormente,Γ(z+1)=z Γ(z).{\displaystyle \Gamma (z+1)=z\ \Gamma (z).}Otras ecuaciones funcionales importantes para la función gamma son la fórmula de reflexión de Euler ,Γ(1z)Γ(z)=πpecadoπz,zZ,{\displaystyle \Gamma (1-z)\Gamma (z)={\frac {\pi }{\sin \pi z}},\qquad z\not \in \mathbb {Z} ,}lo cual implicaΓ(znorte)=(1)norte1Γ(z)Γ(1+z)Γ(norte+1z),norteZ{\displaystyle \Gamma (z-n)=(-1)^{n-1}\;{\frac {\Gamma (-z)\Gamma (1+z)}{\Gamma (n+1-z)}},\qquad n\in \mathbb {Z} } y la fórmula de duplicación de LegendreΓ(z)Γ(z+12)=212zπΓ(2z).{\displaystyle \Gamma (z)\Gamma \left(z+{\frac {1}{2}}\right)=2^{1-2z}\,{\sqrt {\pi }}\,\Gamma (2z).}

La fórmula de duplicación es un caso especial del teorema de multiplicación (véase [ 9 ] Eq. 5.5.6):   k=0metro1Γ(z+kmetro)=(2π)metro12metro12metrozΓ(metroz).{\displaystyle \prod _{k=0}^{m-1}\Gamma \left(z+{\frac {k}{m}}\right)=(2\pi )^{\frac {m-1}{2}}\;m^{{\frac {1}{2}}-mz}\;\Gamma (mz).}

Una propiedad simple pero útil, que se puede observar en la definición de límite, es: Γ(z)¯=Γ(z¯)Γ(z)Γ(z¯)R.{\displaystyle {\overline {\Gamma (z)}}=\Gamma ({\overline {z}})\;\Rightarrow \;\Gamma (z)\Gamma ({\overline {z}})\in \mathbb {R} .}

En particular, conz=a+bi{\displaystyle z=a+bi} , este producto es |Γ(a+bi)|2=|Γ(a)|2k=011+b2(a+k)2{\displaystyle |\Gamma (a+bi)|^{2}=|\Gamma (a)|^{2}\prod _{k=0}^{\infty }{\frac {1}{1+{\frac {b^{2}}{(a+k)^{2}}}}}}

Si la parte real es un número entero o un semi-entero , esto se puede expresar de forma finita en forma cerrada : |Γ(bi)|2=πbsinhπb|Γ(12+bi)|2=πaporrearπb|Γ(1+bi)|2=πbsinhπb|Γ(1+norte+bi)|2=πbsinhπbk=1norte(k2+b2),nortenorte|Γ(norte+bi)|2=πbsinhπbk=1norte(k2+b2)1,nortenorte|Γ(12±norte+bi)|2=πaporrearπbk=1norte((k12)2+b2)±1,nortenorte{\displaystyle {\begin{aligned}|\Gamma (bi)|^{2}&={\frac {\pi }{b\sinh \pi b}}\\[6pt]\left|\Gamma \left({\tfrac {1}{2}}+bi\right)\right|^{2}&={\frac {\pi }{\cosh \pi b}}\\[6pt]\left|\Gamma \left(1+bi\right)\right|^{2}&={\frac {\pi b}{\sinh \pi b}}\\[6pt]\left|\Gamma \left(1+n+bi\right)\right|^{2}&={\frac {\pi b}{\sinh \pi b}}\prod _{k=1}^{n}\left(k^{2}+b^{2}\right),\quad n\in \mathbb {N} \\[6pt]\left|\Gamma \left(-n+bi\right)\right|^{2}&={\frac {\pi }{b\sinh \pi b}}\prod _{k=1}^{n}\left(k^{2}+b^{2}\right)^{-1},\quad n\in \mathbb {N} \\[6pt]\left|\Gamma \left({\tfrac {1}{2}}\pm n+bi\right)\right|^{2}&={\frac {\pi }{\cosh \pi b}}\prod _{k=1}^{n}\left(\left(k-{\tfrac {1}{2}}\right)^{2}+b^{2}\right)^{\pm 1},\quad n\in \mathbb {N} \\[-1ex]&\end{aligned}}}

Quizás el valor más conocido de la función gamma en un argumento no entero es Γ(12)=π,{\displaystyle \Gamma \left({\tfrac {1}{2}}\right)={\sqrt {\pi }},} que se puede encontrar configurandoz=12{\textstyle z={\frac {1}{2}}}en la fórmula de reflexión, utilizando la relación con la función beta que se muestra a continuación con z1=z2=12{\displaystyle z_{1}=z_{2}={\frac {1}{2}}}, o simplemente haciendo la sustituciónt=2{\displaystyle t=u^{2}}en la definición integral de la función gamma, resultando en una integral gaussiana . En general, para valores enteros no negativos denorte{\displaystyle n}tenemos: Γ(12+norte)=(2norte)¡4nortenorte¡π=(2norte1)¡¡2norteπ=(norte12norte)norte¡πΓ(12norte)=(4)nortenorte¡(2norte)¡π=(2)norte(2norte1)¡¡π=π(1/2norte)norte¡{\displaystyle {\begin{aligned}\Gamma \left({\frac {1}{2}}+n\right)&={(2n)! \over 4^{n}n!}{\sqrt {\pi }}={\frac {(2n-1)!!}{2^{n}}}{\sqrt {\pi }}={\binom {n-{\frac {1}{2}}}{n}}\,n!\,{\sqrt {\pi }}\\[6pt]\Gamma \left({\frac {1}{2}}-n\right)&={(-4)^{n}n! \over (2n)!}{\sqrt {\pi }}={\frac {(-2)^{n}}{(2n-1)!!}}{\sqrt {\pi }}={\frac {\sqrt {\pi }}{{\binom {-1/2}{n}}\,n!}}\end{aligned}}}donde el factorial doble(2norte1)¡¡=(2norte1)(2norte3)(3)(1){\displaystyle (2n-1)!!=(2n-1)(2n-3)\cdots (3)(1)}Consulte los valores particulares de la función gamma para obtener los valores calculados.

Podría ser tentador generalizar el resultado de queΓ(12)=π{\textstyle \Gamma \left({\frac {1}{2}}\right)={\sqrt {\pi }}}buscando una fórmula para otros valores individualesΓ(r){\displaystyle \Gamma (r)}dónder{\displaystyle r}es racional, especialmente porque según el teorema digamma de Gauss , es posible hacerlo para la función digamma estrechamente relacionada en cada valor racional. Sin embargo, estos númerosΓ(r){\displaystyle \Gamma (r)}No se sabe que puedan expresarse por sí mismas en términos de funciones elementales. Se ha demostrado queΓ(norte+r){\displaystyle \Gamma (n+r)}es un número trascendental y algebraicamente independiente deπ{\displaystyle \pi }para cualquier número enteronorte{\displaystyle n}y cada una de las fraccionesr=16,14,13,23,34,56{\displaystyle \textstyle r={\frac {1}{6}},{\frac {1}{4}},{\frac {1}{3}},{\frac {2}{3}},{\frac {3}{4}},{\frac {5}{6}}} . [ 10 ] En general, al calcular los valores de la función gamma, debemos conformarnos con aproximaciones numéricas.

Las derivadas de la función gamma se describen en términos de la función poligamma , ψ(metro)(z){\displaystyle \psi ^{(m)}(z)}:Γ(z)=Γ(z)ψ(0)(z).{\displaystyle \Gamma '(z)=\Gamma (z)\psi ^{(0)}(z).}Para un entero positivometro{\displaystyle m}La derivada de la función gamma se puede calcular de la siguiente manera:

Función gamma en el plano complejo con el color mostrando su argumento.
Tono que muestra el argumento de la función gamma en el plano complejo desde −2 − 2 i hasta 6 + 2 i

Γ(metro+1)=metro¡(γ+k=1metro1k)=metro¡(γ+H(metro)),{\displaystyle \Gamma '(m+1)=m!\left(-\gamma +\sum _{k=1}^{m}{\frac {1}{k}}\right)=m!\left(-\gamma +H(m)\right)\,,} dóndeH(metro){\displaystyle H(m)}es elmetro{\displaystyle m}número armónico yγ{\displaystyle \gamma }es la constante de Euler-Mascheroni .

Para(z)>0{\displaystyle \Re (z)>0}elnorte{\displaystyle n}La derivada de la función gamma es:dnortedznorteΓ(z)=0tz1mit(registrot)nortedt.{\displaystyle {\frac {d^{n}}{dz^{n}}}\Gamma (z)=\int _{0}^{\infty }t^{z-1}e^{-t}(\log t)^{n}\,dt.}(Esto se puede derivar diferenciando la forma integral de la función gamma con respecto a z{\displaystyle z}. )

Utilizando la identidadΓ(norte)(1)=(1)norteBnorte(γ,1¡ζ(2),,(norte1)¡ζ(norte)),{\displaystyle \Gamma ^{(n)}(1)=(-1)^{n}B_{n}(\gamma ,1!\zeta (2),\ldots ,(n-1)!\,\zeta (n)),}dóndeζ(z){\displaystyle \zeta (z)}es la función zeta de Riemann yBnorte{\displaystyle B_{n}}es elnorte{\displaystyle n}polinomio de Bell , tenemos en particular la expansión en serie de Laurent de la función gamma [ 11 ].Γ(z)=1zγ+12(γ2+π26)z16(γ3+γπ22+2ζ(3))z2+O(z3).{\displaystyle \Gamma (z)={\frac {1}{z}}-\gamma +{\frac {1}{2}}\left(\gamma ^{2}+{\frac {\pi ^{2}}{6}}\right)z-{\frac {1}{6}}\left(\gamma ^{3}+{\frac {\gamma \pi ^{2}}{2}}+2\zeta (3)\right)z^{2}+O(z^{3}).}

Desigualdades

Cuando se restringe a los números reales positivos , la función gamma es una función estrictamente logarítmicamente convexa . Esta propiedad puede enunciarse de cualquiera de las siguientes tres formas equivalentes:

  • Para cualesquiera dos números reales positivosincógnita1{\displaystyle x_{1}}yincógnita2{\displaystyle x_{2}} , y para cualquiert[0,1]{\displaystyle t\in [0,1]} ,Γ(tincógnita1+(1t)incógnita2)Γ(incógnita1)tΓ(incógnita2)1t.{\displaystyle \Gamma (tx_{1}+(1-t)x_{2})\leq \Gamma (x_{1})^{t}\Gamma (x_{2})^{1-t}.}
  • Para cualesquiera dos números reales positivosincógnita1{\displaystyle x_{1}}yincógnita2{\displaystyle x_{2}}yincógnita2{\displaystyle x_{2}}>incógnita1{\displaystyle x_{1}}(Γ(incógnita2)Γ(incógnita1))1incógnita2incógnita1>exp(Γ(incógnita1)Γ(incógnita1)).{\displaystyle \left({\frac {\Gamma (x_{2})}{\Gamma (x_{1})}}\right)^{\frac {1}{x_{2}-x_{1}}}>\exp \left({\frac {\Gamma '(x_{1})}{\Gamma (x_{1})}}\right).}
  • Para cualquier número real positivoincógnita{\displaystyle x} ,Γ(incógnita)Γ(incógnita)>Γ(incógnita)2.{\displaystyle \Gamma ''(x)\Gamma (x)>\Gamma '(x)^{2}.}

La última de estas afirmaciones es, esencialmente por definición, la misma que la afirmación de queψ(1)(incógnita)>0{\displaystyle \psi ^{(1)}(x)>0}, dondeψ(1){\displaystyle \psi ^{(1)}}es la función poligamma de orden 1. Para demostrar la convexidad logarítmica de la función gamma, basta con observar queψ(1){\displaystyle \psi ^{(1)}}tiene una representación en serie que, para x real positivo , consta únicamente de términos positivos.

La convexidad logarítmica y la desigualdad de Jensen implican conjuntamente, para cualesquiera números reales positivos :incógnita1,,incógnitanorte{\displaystyle x_{1},\ldots ,x_{n}}ya1,,anorte{\displaystyle a_{1},\ldots ,a_{n}} , Γ(a1incógnita1++anorteincógnitanortea1++anorte)(Γ(incógnita1)a1Γ(incógnitanorte)anorte)1a1++anorte.{\displaystyle \Gamma \left({\frac {a_{1}x_{1}+\cdots +a_{n}x_{n}}{a_{1}+\cdots +a_{n}}}\right)\leq {\bigl (}\Gamma (x_{1})^{a_{1}}\cdots \Gamma (x_{n})^{a_{n}}{\bigr )}^{\frac {1}{a_{1}+\cdots +a_{n}}}.}

También existen límites para las razones de las funciones gamma. El más conocido es la desigualdad de Gautschi , que dice que para cualquier número real positivo x y cualquier s ∈ (0, 1) , incógnita1s<Γ(incógnita+1)Γ(incógnita+s)<(incógnita+1)1s.{\displaystyle x^{1-s}<{\frac {\Gamma (x+1)}{\Gamma (x+s)}}<\left(x+1\right)^{1-s}.}

Fórmula de Stirling

Representación de la función gamma en el plano complejo. Cada puntoz{\displaystyle z}se colorea según el argumento de Γ(z){\displaystyle \Gamma (z)} . El gráfico de contorno del módulo|Γ(z)|{\displaystyle |\Gamma (z)|}También se muestra.
Representación tridimensional del valor absoluto de la función gamma compleja.

El comportamiento deΓ(incógnita){\displaystyle \Gamma (x)}para una variable real positiva creciente viene dada por la fórmula de StirlingΓ(incógnita+1)2πincógnita(incógnitami)incógnita,{\displaystyle \Gamma (x+1)\sim {\sqrt {2\pi x}}\left({\frac {x}{e}}\right)^{x},}donde el símbolo{\displaystyle \sim }significa convergencia asintótica: la razón de los dos lados converge a 1{\displaystyle 1}en el límiteincógnita+{\displaystyle \textstyle x\to +\infty } . [ 1 ] Este crecimiento es más rápido que el exponencial,exp(βincógnita){\displaystyle \exp(\beta x)}, para cualquier valor fijo deβ{\displaystyle \beta }.

Otro límite útil para aproximaciones asintóticas paraincógnita+{\displaystyle x\to +\infty }es:Γ(incógnita+α)Γ(incógnita)incógnitaα,αdo.{\displaystyle {\Gamma (x+\alpha )}\sim {\Gamma (x)x^{\alpha }},\qquad \alpha \in \mathbb {C} .}

Al escribir el término de error como un producto infinito, se puede utilizar la fórmula de Stirling para definir la función gamma: [ 12 ]Γ(incógnita)=2πincógnita(incógnitami)incógnitanorte=0[1mi(1+1incógnita+norte)incógnita+norte+12].{\displaystyle \Gamma (x)={\sqrt {\frac {2\pi }{x}}}\left({\frac {x}{e}}\right)^{x}\prod _{n=0}^{\infty }\left[{\frac {1}{e}}\left(1+{\frac {1}{x+n}}\right)^{x+n+{\frac {1}{2}}}\right].}

Extensión a valores negativos y no enteros.

Aunque la definición principal de la función gamma —la integral de Euler de segundo tipo— solo es válida (en el eje real) para argumentos positivos, su dominio puede extenderse con continuación analítica [ 13 ] a argumentos negativos desplazando el argumento negativo a valores positivos mediante el uso de la fórmula de reflexión de Euler,Γ(incógnita)=1Γ(incógnita+1)πpecado(π(incógnita+1)),{\displaystyle \Gamma (-x)={\frac {1}{\Gamma (x+1)}}{\frac {\pi }{\sin {\big (}\pi (x+1){\big )}}},}o la propiedad fundamental,Γ(incógnita):=1incógnitaΓ(incógnita+1),{\displaystyle \Gamma (-x):={\frac {\,1}{-x}}\,\Gamma (-x+1),}cuandoincógnitaZ{\displaystyle x\not \in \mathbb {Z} } . Por ejemplo,Γ(12)=2Γ(12).{\displaystyle \Gamma \!\left(\!-{\frac {1}{2}}\right)=-2\,\Gamma \!\left({\frac {1}{2}}\right).}

Residuos

El comportamiento para no positivosz{\displaystyle z}es más intrincado. La integral de Euler no converge para (z)0{\displaystyle \Re (z)\leq 0} , pero la función que define en el semiplano complejo positivo tiene una continuación analítica única al semiplano negativo. Una forma de encontrar esa continuación analítica es usar la integral de Euler para argumentos positivos y extender el dominio a números negativos mediante la aplicación repetida de la fórmula de recurrencia, [ 1 ]Γ(z)=Γ(z+norte+1)z(z+1)(z+norte),{\displaystyle \Gamma (z)={\frac {\Gamma (z+n+1)}{z(z+1)\cdots (z+n)}},}elegirnorte{\displaystyle n}de tal manera quez+norte{\displaystyle z+n}es positivo. El producto en el denominador es cero cuandoz{\displaystyle z}es igual a cualquiera de los números enteros0,1,2,{\displaystyle 0,-1,-2,\ldots }Por lo tanto, la función gamma debe ser indefinida en esos puntos para evitar la división por cero ; es una función meromorfa con polos simples en los enteros no positivos. [ 1 ]

Para una funciónF{\displaystyle f}de una variable complejaz{\displaystyle z} , en un poste simple do{\displaystyle c} , el residuo deF{\displaystyle f}está dado por:Res(F,do)=límitezdo(zdo)F(z).{\displaystyle \operatorname {Res} (f,c)=\lim _{z\to c}(z-c)f(z).}

Para el poste simplez=norte{\displaystyle z=-n} , la fórmula de recurrencia se puede reescribir como:(z+norte)Γ(z)=Γ(z+norte+1)z(z+1)(z+norte1).{\displaystyle (z+n)\Gamma (z)={\frac {\Gamma (z+n+1)}{z(z+1)\cdots (z+n-1)}}.}El numerador en z=norte{\displaystyle z=-n} , esΓ(z+norte+1)=Γ(1)=1{\displaystyle \Gamma (z+n+1)=\Gamma (1)=1}y el denominadorz(z+1)(z+norte1)=norte(1norte)(norte1norte)=(1)nortenorte¡.{\displaystyle z(z+1)\cdots (z+n-1)=-n(1-n)\cdots (n-1-n)=(-1)^{n}n!.}Entonces, los residuos de la función gamma en esos puntos son: [ 14 ]Res(Γ,norte)=(1)nortenorte¡.{\displaystyle \operatorname {Res} (\Gamma ,-n)={\frac {(-1)^{n}}{n!}}.}La función gamma es distinta de cero en todos los puntos de la recta real, aunque se aproxima arbitrariamente a cero cuando z{\displaystyle z\rightarrow -\infty } . De hecho, no existe ningún número complejo.z{\displaystyle z}para lo cualΓ(z)=0{\displaystyle \Gamma (z)=0}y , por lo tanto, la función gamma recíproca1Γ(z){\textstyle {\dfrac {1}{\Gamma (z)}}}es una función completa , con ceros en z=0,1,2,{\displaystyle z=0,-1,-2,\ldots } . [ 1 ]

Mínimos y máximos

En la recta real, la función gamma tiene un mínimo local en z min = +1,46163 21449 68362 34126... [ 15 ] donde alcanza el valor Γ( z min ) = +0,88560 31944 10888 70027... . [ 16 ] La función gamma aumenta a ambos lados de este mínimo. La solución de Γ( z − 0,5) = Γ( z + 0,5) es z = +1,5 y el valor común es Γ(1) = Γ(2) = +1 . La solución positiva de Γ( z − 1) = Γ( z + 1) es z = φ ≈ +1,618 , la proporción áurea , y el valor común es Γ( φ − 1) = Γ( φ + 1) = φ ! = +1,44922 96022 69896 60037... . [ 17 ]

La función gamma debe alternar el signo entre sus polos en los enteros no positivos porque el producto en la recurrencia hacia adelante contiene un número impar de factores negativos si el número de polos entrez{\displaystyle z}yz+norte{\displaystyle z+n}es impar, y un número par si el número de polos es par. [ 14 ] Los valores en los extremos locales de la función gamma a lo largo del eje real entre los enteros no positivos son:

Γ( −0,50408 30082 64455 40925... [ 18 ] ) = −3,54464 36111 55005 08912... ,
Γ( −1,57349 84731 62390 45877... [ 19 ] ) = 2,30240 72583 39680 13582... ,
Γ( −2,61072 08684 44144 65000... [ 20 ] ) = −0,88813 63584 01241 92009... ,
Γ( −3,63529 33664 36901 09783... [ 21 ] ) = 0,24512 75398 34366 25043... ,
Γ( −4,65323 77617 43142 44171... [ 22 ] ) = −0,05277 96395 87319 40076... , etc.

Representaciones integrales

Hay muchas fórmulas, además de la integral de Euler de segunda especie, que expresan la función gamma como una integral. Por ejemplo, cuando la parte real dez{\displaystyle z}es positivo, [ 23 ]Γ(z)=miztmitdt{\displaystyle \Gamma (z)=\int _{-\infty }^{\infty }e^{zt-e^{t}}\,dt}y [ 24 ]Γ(z)=01(registro1t)z1dt,{\displaystyle \Gamma (z)=\int _{0}^{1}\left(\log {\frac {1}{t}}\right)^{z-1}\,dt,}Γ(z)=2doz0t2z1midot2dt,do>0{\displaystyle \Gamma (z)=2c^{z}\int _{0}^{\infty }t^{2z-1}e^{-ct^{2}}\,dt\,,\;c>0} donde las tres integrales se obtienen respectivamente de las sustitucionest=miincógnita{\displaystyle t=e^{-x}} ,t=registroincógnita{\displaystyle t=-\log x}[ 25 ] yt=doincógnita2{\displaystyle t=cx^{2}}[ 26 ] en la segunda integral de Euler. La última integral en particular deja clara la conexión entre la función gamma en argumentos semi-enteros y laintegral gaussiana: siz=1/2,do=1{\displaystyle z=1/2,\;c=1}obtenemosΓ(1/2)=20mit2dt=π.{\displaystyle \Gamma (1/2)=2\int _{0}^{\infty }e^{-t^{2}}\,dt={\sqrt {\pi }}\;.}

La primera fórmula integral de Binet para la función gamma establece que, cuando la parte real de z es positiva, entonces: [ 27 ]logramoΓ(z)=(z12)registrozz+12registro(2π)+0(121t+1mit1)mitztdt.{\displaystyle \operatorname {log\Gamma } (z)=\left(z-{\frac {1}{2}}\right)\log z-z+{\frac {1}{2}}\log(2\pi )+\int _{0}^{\infty }\left({\frac {1}{2}}-{\frac {1}{t}}+{\frac {1}{e^{t}-1}}\right){\frac {e^{-tz}}{t}}\,dt.}La integral del lado derecho puede interpretarse como una transformada de Laplace . Es decir,registro(Γ(z)(miz)zz2π)=L(12t1t2+1t(mit1))(z).{\displaystyle \log \left(\Gamma (z)\left({\frac {e}{z}}\right)^{z}{\sqrt {\frac {z}{2\pi }}}\right)={\mathcal {L}}\left({\frac {1}{2t}}-{\frac {1}{t^{2}}}+{\frac {1}{t(e^{t}-1)}}\right)(z).}

La segunda fórmula integral de Binet establece que, nuevamente cuando la parte real de z es positiva, entonces: [ 28 ]logramoΓ(z)=(z12)registrozz+12registro(2π)+20arctan(t/z)mi2πt1dt.{\displaystyle \operatorname {log\Gamma } (z)=\left(z-{\frac {1}{2}}\right)\log z-z+{\frac {1}{2}}\log(2\pi )+2\int _{0}^{\infty }{\frac {\arctan(t/z)}{e^{2\pi t}-1}}\,dt.}

Sea C un contorno de Hankel , es decir, un camino que comienza y termina en el punto en la esfera de Riemann , cuyo vector tangente unitario converge a −1 al inicio del camino y a 1 al final, que tiene número de vueltas 1 alrededor de 0 , y que no cruza [0,){\displaystyle [0,\infty )} . Arreglar una rama deregistro(t){\displaystyle \log(-t)}tomando una rama cortada a lo largo[0,){\displaystyle [0,\infty )}y tomandoregistro(t){\displaystyle \log(-t)}ser real cuandot{\displaystyle t}está en el eje real negativo. Siz{\displaystyle z}Si no es un número entero, entonces la fórmula de Hankel para la función gamma es: [ 29 ]Γ(z)=12ipecadoπzdo(t)z1mitdt,{\displaystyle \Gamma (z)=-{\frac {1}{2i\sin \pi z}}\int _{C}(-t)^{z-1}e^{-t}\,dt,}dónde(t)z1{\displaystyle (-t)^{z-1}}se interpreta comoexp((z1)registro(t)){\displaystyle \exp((z-1)\log(-t))} . La fórmula de reflexión conduce a la expresión estrechamente relacionada1Γ(z)=i2πdo(t)zmitdt,{\displaystyle {\frac {1}{\Gamma (z)}}={\frac {i}{2\pi }}\int _{C}(-t)^{-z}e^{-t}\,dt,}que es válido siempre quezZ{\displaystyle z\not \in \mathbb {Z} }.

Representación de fracciones continuas

La función gamma también puede representarse mediante una suma de dos fracciones continuas : [ 30 ] [ 31 ]Γ(z)=mi12+0z+1z12+2z+2z22+4z+3z32+6z+4z42+8z+5z52+10z++ mi1z+0z+0z+1+1z+2z+1z+3+2z+4z+2z+5+3z+6{\displaystyle {\begin{aligned}\Gamma (z)&={\cfrac {e^{-1}}{2+0-z+1{\cfrac {z-1}{2+2-z+2{\cfrac {z-2}{2+4-z+3{\cfrac {z-3}{2+6-z+4{\cfrac {z-4}{2+8-z+5{\cfrac {z-5}{2+10-z+\ddots }}}}}}}}}}}}\\&+\ {\cfrac {e^{-1}}{z+0-{\cfrac {z+0}{z+1+{\cfrac {1}{z+2-{\cfrac {z+1}{z+3+{\cfrac {2}{z+4-{\cfrac {z+2}{z+5+{\cfrac {3}{z+6-\ddots }}}}}}}}}}}}}}\end{aligned}}}dondezdo{\displaystyle z\in \mathbb {C} }.

expansión de la serie de Fourier

El logaritmo de la función gamma tiene la siguiente expansión en serie de Fourier para0<z<1:{\displaystyle 0<z<1:}logramoΓ(z)=(12z)(γ+registro2)+(1z)registroπ12registropecado(πz)+1πnorte=1registronortenortepecado(2πnortez),{\displaystyle \operatorname {log\Gamma } (z)=\left({\frac {1}{2}}-z\right)(\gamma +\log 2)+(1-z)\log \pi -{\frac {1}{2}}\log \sin(\pi z)+{\frac {1}{\pi }}\sum _{n=1}^{\infty }{\frac {\log n}{n}}\sin(2\pi nz),} que durante mucho tiempo se atribuyó a Ernst Kummer , quien la derivó en 1847. [ 32 ] [ 33 ] Sin embargo, Iaroslav Blagouchine descubrió que Carl Johan Malmsten derivó esta serie por primera vez en 1842. [ 34 ] [ 35 ]

Fórmula de Raabe

En 1840 Joseph Ludwig Raabe demostró queaa+1registroΓ(z)dz=12registro2π+aregistroaa,a>0.{\displaystyle \int _{a}^{a+1}\log \Gamma (z)\,dz={\tfrac {1}{2}}\log 2\pi +a\log a-a,\quad a>0.}En particular, sia=0{\displaystyle a=0}entonces01registroΓ(z)dz=12registro(2π).{\displaystyle \int _{0}^{1}\log \Gamma (z)\,dz={\frac {1}{2}}\log(2\pi ).}

Este último se puede derivar tomando el logaritmo en la fórmula de multiplicación anterior, lo que da una expresión para la suma de Riemann del integrando. Tomando el límite paraa{\displaystyle a\to \infty }da la fórmula.

Función Pi

Una notación alternativa introducida por Gauss es laΠ{\displaystyle \Pi }-función, una versión desplazada de la función gamma: Π(z)=Γ(z+1)=zΓ(z)=0mittzdt,{\displaystyle \Pi (z)=\Gamma (z+1)=z\Gamma (z)=\int _{0}^{\infty }e^{-t}t^{z}\,dt,} de modo queΠ(norte)=norte¡{\displaystyle \Pi (n)=n!}para cada entero no negativonorte{\displaystyle n}.

Utilizando la función pi, la fórmula de reflexión es: Π(z)Π(z)=πzpecado(πz)=1desde(z){\displaystyle \Pi (z)\Pi (-z)={\frac {\pi z}{\sin(\pi z)}}={\frac {1}{\operatorname {sinc} (z)}}}utilizando la función sinc normalizada ; mientras que el teorema de la multiplicación se convierte en: Π(zmetro)Π(z1metro)Π(zmetro+1metro)=(2π)metro12metroz12Π(z) .{\displaystyle \Pi \left({\frac {z}{m}}\right)\,\Pi \left({\frac {z-1}{m}}\right)\cdots \Pi \left({\frac {z-m+1}{m}}\right)=(2\pi )^{\frac {m-1}{2}}m^{-z-{\frac {1}{2}}}\Pi (z)\ .}

La función gamma recíproca desplazada a veces se denota comoπ(z)=1Π(z){\displaystyle \pi (z)={\frac {1}{\Pi (z)}}} , una función completa .

El volumen de un n -elipsoide con radios r 1 , ..., r n se puede expresar comoVnorte(r1,,rnorte)=πnorte2Π(norte2)k=1norterk.{\displaystyle V_{n}(r_{1},\dotsc ,r_{n})={\frac {\pi ^{\frac {n}{2}}}{\Pi \left({\frac {n}{2}}\right)}}\prod _{k=1}^{n}r_{k}.}

Relación con otras funciones

  • En la primera integral que define la función gamma, los límites de integración están fijos. La función gamma incompleta superior se obtiene permitiendo que varíe el límite inferior de integración:Γ(z,incógnita)=incógnitatz1mitdt.{\displaystyle \Gamma (z,x)=\int _{x}^{\infty }t^{z-1}e^{-t}dt.}Existe una función gamma incompleta inferior similar.
  • La función gamma está relacionada con la función beta de Euler mediante la fórmulaB(z1,z2)=01tz11(1t)z21dt=Γ(z1)Γ(z2)Γ(z1+z2).{\displaystyle \mathrm {B} (z_{1},z_{2})=\int _{0}^{1}t^{z_{1}-1}(1-t)^{z_{2}-1}\,dt={\frac {\Gamma (z_{1})\,\Gamma (z_{2})}{\Gamma (z_{1}+z_{2})}}.}
  • La derivada logarítmica de la función gamma se llama función digamma ; las derivadas de orden superior son las funciones poligamma .
  • El análogo de la función gamma sobre un campo finito o un anillo finito son las sumas gaussianas , un tipo de suma exponencial .
  • La función gamma recíproca es una función completa y se ha estudiado como un tema específico.
  • La función gamma también aparece en una relación importante con la función zeta de Riemann ,ζ(z){\displaystyle \zeta (z)}.πz2Γ(z2)ζ(z)=π1z2Γ(1z2)ζ(1z).{\displaystyle \pi ^{-{\frac {z}{2}}}\;\Gamma \left({\frac {z}{2}}\right)\zeta (z)=\pi ^{-{\frac {1-z}{2}}}\;\Gamma \left({\frac {1-z}{2}}\right)\;\zeta (1-z).}También aparece en la siguiente fórmula:ζ(z)Γ(z)=0zmi1d,{\displaystyle \zeta (z)\Gamma (z)=\int _{0}^{\infty }{\frac {u^{z}}{e^{u}-1}}\,{\frac {du}{u}},}que es válido solo para(z)>1{\displaystyle \Re (z)>1}.
    El logaritmo de la función gamma satisface la siguiente fórmula debida a Lerch:logramoΓ(z)=ζH(0,z)ζ(0),{\displaystyle \operatorname {log\Gamma } (z)=\zeta _{H}'(0,z)-\zeta '(0),}dóndeζH{\displaystyle \zeta _{H}}es la función zeta de Hurwitz ,ζ{\displaystyle \zeta }es la función zeta de Riemann y la prima ( ) denota la diferenciación en la primera variable.
  • La función gamma está relacionada con la función exponencial estirada . Por ejemplo, los momentos de esa función son:τnorte0tnorte1mi(tτ)βdt=τnorteβΓ(norteβ).{\displaystyle \langle \tau ^{n}\rangle \equiv \int _{0}^{\infty }t^{n-1}\,e^{-\left({\frac {t}{\tau }}\right)^{\beta }}\,\mathrm {d} t={\frac {\tau ^{n}}{\beta }}\Gamma \left({n \over \beta }\right).}

Valores particulares

Incluyendo hasta los primeros 20 dígitos después del punto decimal, algunos valores particulares de la función gamma son: Γ(32)=4π3+2.36327180120735470306Γ(12)=2π3.54490770181103205459Γ(12)=π+1,77245385090551602729Γ(1)=0¡=+1Γ(32)=π2+0,88622692545275801364Γ(2)=1¡=+1Γ(52)=3π4+1.32934038817913702047Γ(3)=2¡=+2Γ(72)=15π8+3.32335097044784255118Γ(4)=3¡=+6{\displaystyle {\begin{array}{rcccl}\Gamma \left(-{\frac {3}{2}}\right)&=&{\frac {4{\sqrt {\pi }}}{3}}&\approx &+2.36327\,18012\,07354\,70306\\[6pt]\Gamma \left(-{\frac {1}{2}}\right)&=&-2{\sqrt {\pi }}&\approx &-3.54490\,77018\,11032\,05459\\[6pt]\Gamma \left({\frac {1}{2}}\right)&=&{\sqrt {\pi }}&\approx &+1.77245\,38509\,05516\,02729\\[6pt]\Gamma (1)&=&0!&=&+1\\[6pt]\Gamma \left({\frac {3}{2}}\right)&=&{\frac {\sqrt {\pi }}{2}}&\approx &+0.88622\,69254\,52758\,01364\\[6pt]\Gamma (2)&=&1!&=&+1\\[6pt]\Gamma \left({\frac {5}{2}}\right)&=&{\frac {3{\sqrt {\pi }}}{4}}&\approx &+1.32934\,03881\,79137\,02047\\[6pt]\Gamma (3)&=&2!&=&+2\\[6pt]\Gamma \left({\frac {7}{2}}\right)&=&{\tfrac {15{\sqrt {\pi }}}{8}}&\approx &+3.32335\,09704\,47842\,55118\\[6pt]\Gamma (4)&=&3!&=&+6\end{array}}}(Estos números se pueden encontrar en el OEIS . [ 36 ] [ 37 ] [ 38 ] [ 39 ] [ 40 ] [ 41 ] Los valores presentados aquí están truncados en lugar de redondeados.) La función gamma de valor complejo no está definida para enteros no positivos, pero en estos casos el valor se puede definir en la esfera de Riemann como {\displaystyle \infty } . La función gamma recíproca está bien definida y es analítica en estos valores (y en todo el plano complejo ):1Γ(3)=1Γ(2)=1Γ(1)=1Γ(0)=0.{\displaystyle {\frac {1}{\Gamma (-3)}}={\frac {1}{\Gamma (-2)}}={\frac {1}{\Gamma (-1)}}={\frac {1}{\Gamma (0)}}=0.}

función logaritmo-gamma

La función analítica logΓ( z )

Debido a que las funciones gamma y factorial crecen tan rápidamente para argumentos moderadamente grandes, muchos entornos informáticos incluyen una función que devuelve el logaritmo natural de la función gamma, a menudo denominada lgammao lngammaen entornos de programación o gammalnen hojas de cálculo. Esta crece mucho más lentamente y, para cálculos combinatorios, permite sumar y restar valores logarítmicos en lugar de multiplicar y dividir valores muy grandes. A menudo se define como [ 42 ] .logramoΓ(z)=γzregistroz+k=1[zkregistro(1+zk)].{\displaystyle \operatorname {log\Gamma } (z)=-\gamma z-\log z+\sum _{k=1}^{\infty }\left[{\frac {z}{k}}-\log \left(1+{\frac {z}{k}}\right)\right].}

La función digamma , que es la derivada de esta función, también se ve comúnmente. En el contexto de aplicaciones técnicas y físicas, por ejemplo con la propagación de ondas, la ecuación funcional logramoΓ(z)=logramoΓ(z+1)registroz{\displaystyle \operatorname {log\Gamma } (z)=\operatorname {log\Gamma } (z+1)-\log z}

Función gamma logarítmica en el plano complejo desde −2 − 2i hasta 2 + 2i con hue que da el argumento complejo
Función gamma logarítmica en el plano complejo desde −2 − 2i hasta 2 + 2i, donde el tono proporciona el argumento complejo.

se usa a menudo ya que permite determinar valores de función en una franja de ancho 1 enz{\displaystyle z}de la franja vecina. En particular, comenzando con una buena aproximación para unaz{\displaystyle z}Con una parte real grande, uno puede ir paso a paso hasta llegar al valor deseado .z{\displaystyle z} . Siguiendo una indicación de Carl Friedrich Gauss , Rocktaeschel (1922) propuso pararegistroΓ(z){\displaystyle \log \Gamma (z)}una aproximación para grandes (z){\displaystyle \Re (z)}:logramoΓ(z)(z12)registrozz+12registro(2π).{\displaystyle \operatorname {log\Gamma } (z)\approx (z-{\tfrac {1}{2}})\log z-z+{\tfrac {1}{2}}\log(2\pi ).}

Esto se puede utilizar para aproximar con precisión.registroΓ(z){\displaystyle \log \Gamma (z)}paraz{\displaystyle z}con uno más pequeño(z){\displaystyle \Re (z)}vía (PEBöhmer, 1939) logramoΓ(zmetro)=logramoΓ(z)k=1metroregistro(zk).{\displaystyle \operatorname {log\Gamma } (z-m)=\operatorname {log\Gamma } (z)-\sum _{k=1}^{m}\log(z-k).}

Se puede obtener una aproximación más precisa utilizando más términos de las expansiones asintóticas deregistroΓ(z){\displaystyle \log \Gamma (z)}yΓ(z){\displaystyle \Gamma (z)} , que se basan en la aproximación de Stirling.Γ(z)zz12miz2π(1+112z+1288z213951840z35712488320z4){\displaystyle \Gamma (z)\sim z^{z-{\frac {1}{2}}}e^{-z}{\sqrt {2\pi }}\left(1+{\frac {1}{12z}}+{\frac {1}{288z^{2}}}-{\frac {139}{51\,840z^{3}}}-{\frac {571}{2\,488\,320z^{4}}}\right)}

como |z|{\displaystyle |z|\rightarrow \infty }a constante|arg(z)|<π{\displaystyle \vert \arg(z)\vert <\pi }( Véanse las secuencias A001163 y A001164 en el OEIS ).

En una presentación más "natural",registroΓ(z)=zregistrozz12registroz+12registro2π+112z1360z3+11260z5+O(1z5){\displaystyle \log \Gamma (z)=z\log z-z-{\frac {1}{2}}\log z+{\frac {1}{2}}\log 2\pi +{\frac {1}{12z}}-{\frac {1}{360z^{3}}}+{\frac {1}{1260z^{5}}}+O\left({\frac {1}{z^{5}}}\right)}

como|z|{\displaystyle |z|\rightarrow \infty } a constante|arg(z)|<π{\displaystyle \vert \arg(z)\vert <\pi }( Véanse las secuencias A046968 y A046969 en el OEIS ).

Los coeficientes de los términos conk>1{\displaystyle k>1} dez1k{\displaystyle z^{1-k}} en la última expansión son simplementeBkk(k1),{\displaystyle {\frac {B_{k}}{k(k-1)}},}donde elBk{\displaystyle B_{k}}son los números de Bernoulli .

La función gamma también tiene la serie de Stirling (derivada por Charles Hermite en 1900) igual a [ 43 ].logramoΓ(1+incógnita)=incógnita(incógnita1)2¡registro(2)+incógnita(incógnita1)(incógnita2)3¡(registro(3)2registro(2))+,(incógnita)>0.{\displaystyle \operatorname {log\Gamma } (1+x)={\frac {x(x-1)}{2!}}\log(2)+{\frac {x(x-1)(x-2)}{3!}}(\log(3)-2\log(2))+\cdots ,\quad \Re (x)>0.}

Propiedades

El teorema de Bohr-Mollerup establece que, entre todas las funciones que extienden las funciones factoriales a los números reales positivos, solo la función gamma es log-convexa , es decir, su logaritmo natural es convexo en el eje real positivo. Otra caracterización la proporciona el teorema de Wielandt .

La función gamma es la única función que satisface simultáneamente

  1. Γ(1)=1{\displaystyle \Gamma (1)=1} ,
  2. Γ(z+1)=zΓ(z){\displaystyle \Gamma (z+1)=z\Gamma (z)}para todos los números complejosz{\displaystyle z}excepto los enteros no positivos, y,
  3. para entero n ,límitenorteΓ(norte+z)Γ(norte)nortez=1{\textstyle \lim _{n\to \infty }{\frac {\Gamma (n+z)}{\Gamma (n)\;n^{z}}}=1}para todos los números complejosz{\displaystyle z} . [ 1 ]

En cierto sentido, la función log-gamma es la forma más natural; hace más evidentes algunos atributos intrínsecos de la función. Un ejemplo notable es la serie de Taylor de logΓ alrededor de 1: logramoΓ(z+1)=γz+k=2ζ(k)k(z)k|z|<1{\displaystyle \operatorname {log\Gamma } (z+1)=-\gamma z+\sum _{k=2}^{\infty }{\frac {\zeta (k)}{k}}\,(-z)^{k}\qquad \forall \;|z|<1} conζ(k){\displaystyle \zeta (k)} que denota la función zeta de Riemann enk{\displaystyle k}.

Entonces, utilizando la siguiente propiedad: ζ(s)Γ(s)=0tsmit1dtt{\displaystyle \zeta (s)\Gamma (s)=\int _{0}^{\infty }{\frac {t^{s}}{e^{t}-1}}\,{\frac {dt}{t}}} Una representación integral para la función log-gamma es: logramoΓ(z+1)=γz+0mizt1+ztt(mit1)dt{\displaystyle \operatorname {log\Gamma } (z+1)=-\gamma z+\int _{0}^{\infty }{\frac {e^{-zt}-1+zt}{t\left(e^{t}-1\right)}}\,dt} o bien, estableciendoz=1{\displaystyle z=1}para obtener una integral paraγ{\displaystyle \gamma } , podemos reemplazar elγ{\displaystyle \gamma } término con su integral e incorporar eso en la fórmula anterior, para obtener: logramoΓ(z+1)=0miztzmit1+zt(mit1)dt.{\displaystyle \operatorname {log\Gamma } (z+1)=\int _{0}^{\infty }{\frac {e^{-zt}-ze^{-t}-1+z}{t\left(e^{t}-1\right)}}\,dt\,.}

También existen fórmulas especiales para el logaritmo de la función gamma para números racionales .z{\displaystyle z} . Por ejemplo, sik{\displaystyle k}ynorte{\displaystyle n}son números enteros conk<norte{\displaystyle k<n}yknorte/2{\displaystyle k\neq n/2} , entonces [ 44 ]logramoΓ(knorte)=(norte2k)registro2π2norte+12{registroπregistropecadoπknorte}+1πr=1norte1γ+registrorrpecado2πrknorte12πpecado2πknorte0minorteincógnitaregistroincógnitaaporrearincógnitaporque(2πk/norte)dincógnita.{\displaystyle {\begin{aligned}\operatorname {log\Gamma } \left({\frac {k}{n}}\right)={}&{\frac {\,(n-2k)\log 2\pi \,}{2n}}+{\frac {1}{2}}\left\{\,\log \pi -\log \sin {\frac {\pi k}{n}}\,\right\}+{\frac {1}{\pi }}\!\sum _{r=1}^{n-1}{\frac {\,\gamma +\log r\,}{r}}\cdot \sin {\frac {\,2\pi rk\,}{n}}\\&{}-{\frac {1}{2\pi }}\sin {\frac {2\pi k}{n}}\cdot \!\int _{0}^{\infty }\!\!{\frac {\,e^{-nx}\!\cdot \log x\,}{\,\cosh x-\cos(2\pi k/n)\,}}\,{\mathrm {d} }x.\end{aligned}}}Esta fórmula se utiliza a veces para cálculos numéricos, ya que el integrando disminuye muy rápidamente.

Integración sobre log-gamma

El integral 0zlogramoΓ(incógnita)dincógnita{\displaystyle \int _{0}^{z}\operatorname {log\Gamma } (x)\,dx} puede expresarse en términos de la función G de Barnes [ 45 ] [ 46 ] (véase la función G de Barnes para una demostración): 0zlogramoΓ(incógnita)dincógnita=z2registro(2π)+z(1z)2+zlogramoΓ(z)registroGRAMO(z+1){\displaystyle \int _{0}^{z}\operatorname {log\Gamma } (x)\,dx={\frac {z}{2}}\log(2\pi )+{\frac {z(1-z)}{2}}+z\operatorname {log\Gamma } (z)-\log G(z+1)} donde(z)>1{\displaystyle \Re (z)>-1}.

También se puede escribir en términos de la función zeta de Hurwitz : [ 47 ] [ 48 ]0zlogramoΓ(incógnita)dincógnita=z2registro(2π)+z(1z)2ζ(1)+ζ(1,z).{\displaystyle \int _{0}^{z}\operatorname {log\Gamma } (x)\,dx={\frac {z}{2}}\log(2\pi )+{\frac {z(1-z)}{2}}-\zeta '(-1)+\zeta '(-1,z).}

Cuandoz=1{\displaystyle z=1}resulta que 01logramoΓ(incógnita)dincógnita=12registro(2π),{\displaystyle \int _{0}^{1}\operatorname {log\Gamma } (x)\,dx={\frac {1}{2}}\log(2\pi ),} y esto también es consecuencia de la fórmula de Raabe . Espinosa y Moll derivaron una fórmula similar para la integral del cuadrado de logramoΓ{\displaystyle \operatorname {log\Gamma } }: [ 49 ]01registro2Γ(incógnita)dincógnita=γ212+π248+13γL1+43L12(γ+2L1)ζ(2)π2+ζ(2)2π2,{\displaystyle \int _{0}^{1}\log ^{2}\Gamma (x)dx={\frac {\gamma ^{2}}{12}}+{\frac {\pi ^{2}}{48}}+{\frac {1}{3}}\gamma L_{1}+{\frac {4}{3}}L_{1}^{2}-\left(\gamma +2L_{1}\right){\frac {\zeta ^{\prime }(2)}{\pi ^{2}}}+{\frac {\zeta ^{\prime \prime }(2)}{2\pi ^{2}}},} dóndeL1{\displaystyle L_{1}}es12registro(2π){\displaystyle {\frac {1}{2}}\log(2\pi )}.

DH Bailey y sus coautores [ 50 ] dieron una evaluación para Lnorte:=01registronorteΓ(incógnita)dincógnita{\displaystyle L_{n}:=\int _{0}^{1}\log ^{n}\Gamma (x)\,dx} cuandonorte=1,2{\displaystyle n=1,2}en términos de la función zeta de Tornheim-Witten y sus derivadas.

Además, también se sabe que [ 51 ]límitenorteLnortenorte¡=1.{\displaystyle \lim _{n\to \infty }{\frac {L_{n}}{n!}}=1.}

Aproximaciones

Comparación de la función gamma (línea azul) con el factorial (puntos azules) y la aproximación de Stirling (línea magenta).

Los valores complejos de la función gamma pueden aproximarse utilizando la aproximación de Stirling o la aproximación de Lanczos .Γ(z)2πzz1/2mizcomo z en |arg(z)|<π.{\displaystyle \Gamma (z)\sim {\sqrt {2\pi }}z^{z-1/2}e^{-z}\quad {\hbox{as }}z\to \infty {\hbox{ in }}\left|\arg(z)\right|<\pi .}Esto es preciso en el sentido de que la relación entre la aproximación y el valor verdadero se aproxima a1{\displaystyle 1}en el límite como|z|{\displaystyle \vert z\vert \rightarrow \infty }.

La función gamma se puede calcular con precisión fija para(z)[1,2]{\displaystyle \Re (z)\in [1,2]}aplicando la integración por partes a la integral de Euler. Para cualquier número positivo incógnita{\displaystyle x}La función gamma se puede escribir Γ(z)=0incógnitamittzdtt+incógnitamittzdtt=incógnitazmiincógnitanorte=0incógnitanortez(z+1)(z+norte)+incógnitamittzdtt.{\displaystyle {\begin{aligned}\Gamma (z)&=\int _{0}^{x}e^{-t}t^{z}\,{\frac {dt}{t}}+\int _{x}^{\infty }e^{-t}t^{z}\,{\frac {dt}{t}}\\&=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}{t}}.\end{aligned}}}

Cuando(z)[1,2]{\displaystyle \Re (z)\in [1,2]}yincógnita1{\displaystyle x\geq 1} , el valor absoluto de la última integral es menor que(incógnita+1)miincógnita{\displaystyle (x+1)e^{-x}} . Al elegir un tamaño suficientemente grandeincógnita{\displaystyle x} , esta última expresión se puede hacer más pequeña que2norte{\displaystyle 2^{-N}}para cualquier valor deseadonorte{\displaystyle N} . Por lo tanto, la función gamma puede evaluarse comonorte{\displaystyle N}detalles de precisión con la serie anterior.

EA Karatsuba construyó un algoritmo rápido para el cálculo de la función gamma de Euler para cualquier argumento algebraico (incluidos los racionales). [ 52 ] [ 53 ] [ 54 ]

Para argumentos que son múltiplos enteros de 124{\displaystyle {\tfrac {1}{24}}}La función gamma también puede evaluarse rápidamente mediante iteraciones de media aritmético-geométrica (véanse los valores particulares de la función gamma ). [ 55 ]

Implementaciones prácticas

A diferencia de muchas otras funciones, como la distribución normal , no existe una implementación obvia, rápida y precisa que sea fácil de implementar para la función gamma.Γ(z){\displaystyle \Gamma (z)}es fácil de encontrar. Por lo tanto, vale la pena investigar posibles soluciones. En el caso de que la velocidad sea más importante que la precisión, las tablas publicadas paraΓ(z){\displaystyle \Gamma (z)}Se pueden encontrar fácilmente en una búsqueda en Internet, como la Biblioteca Wiley en línea. Dichas tablas se pueden usar con interpolación lineal . Se puede obtener mayor precisión con el uso de interpolación cúbica a costa de una mayor sobrecarga computacional. Dado queΓ(z){\displaystyle \Gamma (z)}Las tablas se publican normalmente para valores de argumentos entre 1 y 2, la propiedadΓ(z+1)=z Γ(z){\displaystyle \Gamma (z+1)=z\ \Gamma (z)}Puede utilizarse para traducir rápida y fácilmente todos los valores reales.z<1{\displaystyle z<1}yz>2{\displaystyle z>2}dentro del rango1z2{\displaystyle 1\leq z\leq 2} , de modo que solo los valores tabulados dez{\displaystyle z}Se debe utilizar entre 1 y 2. [ 56 ]

Si no se desean tablas de interpolación, entonces la aproximación de Lanczos mencionada anteriormente funciona bien para una precisión de 1 a 2 dígitos para valores pequeños y comúnmente utilizados de z{\displaystyle z}Si la aproximación de Lanczos no es suficientemente precisa,se puede utilizar la fórmula de Stirling para la función gamma .

Aplicaciones

Un autor describe la función gamma como "Posiblemente, la función especial más común , o la menos 'especial' de ellas. Las otras funciones trascendentales [...] se llaman 'especiales' porque teóricamente se podrían evitar algunas de ellas manteniéndose alejado de muchos temas matemáticos especializados. Por otro lado, la función gammaΓ(z){\displaystyle \Gamma (z)}es muy difícil de evitar." [ 57 ]

Problemas de integración

La función gamma encuentra aplicación en áreas tan diversas como la física cuántica , la astrofísica y la dinámica de fluidos . [ 58 ] La distribución gamma , que se formula en términos de la función gamma, se utiliza en estadística para modelar una amplia gama de procesos; por ejemplo, el tiempo entre la ocurrencia de terremotos. [ 59 ]

La razón principal de la utilidad de la función gamma en tales contextos es la prevalencia de expresiones del tipo F(t)migramo(t){\displaystyle f(t)e^{-g(t)}} , que describen procesos que decaen exponencialmente en el tiempo o el espacio. Las integrales de tales expresiones pueden resolverse ocasionalmente en términos de la función gamma cuando no existe una solución elemental. Por ejemplo, siF{\displaystyle f}es una función potencial ygramo{\displaystyle g}es una función lineal, un simple cambio de variables:=at{\displaystyle u:=a\cdot t}da la evaluación 0tbmiatdt=1ab0bmid(a)=Γ(b+1)ab+1.{\displaystyle \int _{0}^{\infty }t^{b}\,e^{-at}\,dt={\frac {1}{a^{b}}}\int _{0}^{\infty }u^{b}\,e^{-u}\,d\left({\frac {u}{a}}\right)={\frac {\Gamma (b+1)}{a^{b+1}}}.}

El hecho de que la integración se realice a lo largo de toda la recta real positiva podría significar que la función gamma describe la acumulación de un proceso dependiente del tiempo que continúa indefinidamente, o que el valor podría ser el total de una distribución en un espacio infinito.

Por supuesto, con frecuencia es útil tomar límites de integración distintos de0{\displaystyle 0}y{\displaystyle \infty }Para describir la acumulación de un proceso finito, la función gamma ordinaria deja de ser una solución; en ese caso, se la denomina función gamma incompleta . (La función gamma ordinaria, obtenida mediante la integración a lo largo de toda la recta real positiva, a veces se denomina función gamma completa para contrastar).

Una categoría importante de funciones que decaen exponencialmente es la de las funciones gaussianas.ami(incógnitab)2do2{\displaystyle ae^{-{\frac {(x-b)^{2}}{c^{2}}}}} y sus integrales, como la función de error . Existen muchas interrelaciones entre estas funciones y la función gamma; en particular, el factorπ{\displaystyle {\sqrt {\pi }}}obtenido mediante evaluaciónΓ(12){\textstyle \Gamma \left({\frac {1}{2}}\right)}es "igual" que el que se encuentra en el factor de normalización de la función de error y la distribución normal .

Las integrales analizadas hasta ahora involucran funciones trascendentales , pero la función gamma también surge de integrales de funciones puramente algebraicas. En particular, las longitudes de arco de las elipses y de la lemniscata , que son curvas definidas por ecuaciones algebraicas, vienen dadas por integrales elípticas que, en casos especiales, pueden evaluarse en términos de la función gamma. La función gamma también puede utilizarse para calcular el "volumen" y el "área" denorte{\displaystyle n}Hiperesferas -dimensionales .

Cálculo de productos

La capacidad de la función gamma para generalizar productos factoriales conduce inmediatamente a aplicaciones en muchas áreas de las matemáticas; en combinatoria y, por extensión, en áreas como la teoría de la probabilidad y el cálculo de series de potencias . Muchas expresiones que involucran productos de enteros sucesivos pueden escribirse como alguna combinación de factoriales, siendo quizás el ejemplo más importante el del coeficiente binomial . Por ejemplo, para cualquier número complejoz{\displaystyle z}ynorte{\displaystyle n}, con|z|<1{\displaystyle \vert z\vert <1} , podemos escribir(1+z)norte=k=0Γ(norte+1)k¡Γ(nortek+1)zk,{\displaystyle (1+z)^{n}=\sum _{k=0}^{\infty }{\frac {\Gamma (n+1)}{k!\Gamma (n-k+1)}}z^{k},}que se asemeja mucho al coeficiente binomial cuandonorte{\displaystyle n}es un número entero no negativo, (1+z)norte=k=0nortenorte¡k¡(nortek)¡zk=k=0norte(nortek)zk.{\displaystyle (1+z)^{n}=\sum _{k=0}^{n}{\frac {n!}{k!(n-k)!}}z^{k}=\sum _{k=0}^{n}{\binom {n}{k}}z^{k}.}

El ejemplo de los coeficientes binomiales motiva por qué las propiedades de la función gamma, cuando se extienden a números negativos, son naturales. Un coeficiente binomial da el número de maneras de elegirk{\displaystyle k}elementos de un conjunto denorte{\displaystyle n}elementos; sik>norte{\displaystyle k>n} , por supuesto que no hay maneras. Sik>norte{\displaystyle k>n}, entonces (nortek)¡{\displaystyle (n-k)!}es el factorial de un entero negativo y, por lo tanto, infinito si usamos la definición de factorial de la función gamma; al dividir por infinito se obtiene el valor esperado de 0{\displaystyle 0}.

Podemos reemplazar el factorial por una función gamma para extender cualquier fórmula de este tipo a los números complejos. Generalmente, esto funciona para cualquier producto en el que cada factor sea una función racional de la variable índice, factorizando la función racional en expresiones lineales. SiPAG{\displaystyle P}yQ{\displaystyle Q}son polinomios mónicos de gradometro{\displaystyle m}ynorte{\displaystyle n}con raíces respectivaspag1,,pagmetro{\displaystyle p_{1},\cdots ,p_{m}}yq1,,qmetro{\displaystyle q_{1},\cdots ,q_{m}} , tenemosi=abPAG(i)Q(i)=(j=1metroΓ(bpagj+1)Γ(apagj))(k=1norteΓ(aqk)Γ(bqk+1)).{\displaystyle \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).}

Si tenemos una forma de calcular numéricamente la función gamma, es muy sencillo calcular los valores numéricos de dichos productos. El número de funciones gamma en el lado derecho depende únicamente del grado de los polinomios, por lo que no importa si ba{\displaystyle b-a}es igual a 5 o 10 5 . Tomando los límites apropiados, la ecuación también puede hacerse válida incluso cuando el producto de la izquierda contiene ceros o polos.

Al tomar límites, ciertos productos racionales con infinitos factores también pueden evaluarse en términos de la función gamma. Gracias al teorema de factorización de Weierstrass , las funciones analíticas pueden escribirse como productos infinitos, y estos a veces pueden representarse como productos finitos o cocientes de la función gamma. Ya hemos visto un ejemplo notable: la fórmula de reflexión representa esencialmente la función seno como el producto de dos funciones gamma. Partiendo de esta fórmula, la función exponencial, así como todas las funciones trigonométricas e hiperbólicas, pueden expresarse en términos de la función gamma.

Otras funciones, incluyendo la función hipergeométrica y casos especiales de la misma, pueden representarse mediante integrales de contorno complejas de productos y cocientes de la función gamma, denominadas integrales de Mellin-Barnes .

Teoría analítica de números

Una aplicación de la función gamma es el estudio de la función zeta de Riemann . Una propiedad fundamental de la función zeta de Riemann es su ecuación funcional :Γ(s2)ζ(s)πs2=Γ(1s2)ζ(1s)π1s2.{\displaystyle \Gamma \left({\frac {s}{2}}\right)\,\zeta (s)\,\pi ^{-{\frac {s}{2}}}=\Gamma \left({\frac {1-s}{2}}\right)\,\zeta (1-s)\,\pi ^{-{\frac {1-s}{2}}}.}

Entre otras cosas, esto proporciona una forma explícita para la continuación analítica de la función zeta a una función meromorfa en el plano complejo y conduce a una prueba inmediata de que la función zeta tiene infinitos ceros denominados "triviales" en la recta real. [ 60 ] Otra fórmula poderosa esζ(s)Γ(s)=0tsmit1dtt.{\displaystyle \zeta (s)\;\Gamma (s)=\int _{0}^{\infty }{\frac {t^{s}}{e^{t}-1}}\,{\frac {dt}{t}}.}

Ambas fórmulas fueron derivadas por Bernhard Riemann en su influyente artículo de 1859 " Ueber die Anzahl der Primzahlen unter einer gegebenen Größe " ("Sobre el número de primos menores que una magnitud dada"), uno de los hitos en el desarrollo de la teoría analítica de números , la rama de las matemáticas que estudia los números primos utilizando las herramientas del análisis matemático.

Historia

La función gamma ha despertado el interés de algunos de los matemáticos más destacados de todos los tiempos. Su historia, notablemente documentada por Philip J. Davis en un artículo que le valió el Premio Chauvenet en 1963 , refleja muchos de los principales avances en matemáticas desde el siglo XVIII. En palabras de Davis: «Cada generación ha encontrado algo interesante que decir sobre la función gamma. Quizás la próxima generación también lo haga». [ 1 ]

Siglo XVIII: Euler y Stirling

Carta de Daniel Bernoulli a Christian Goldbach (6 de octubre de 1729)

El problema de extender el factorial a argumentos no enteros fue aparentemente considerado por primera vez por Daniel Bernoulli y Christian Goldbach en la década de 1720. En particular, en una carta de Bernoulli a Goldbach fechada el 6 de octubre de 1729, Bernoulli introdujo la representación del producto [ 61 ].incógnita¡=límitenorte(norte+1+incógnita2)incógnita1k=1nortek+1k+incógnita,{\displaystyle x!=\lim _{n\to \infty }\left(n+1+{\frac {x}{2}}\right)^{x-1}\prod _{k=1}^{n}{\frac {k+1}{k+x}},}lo cual está bien definido para valores reales de x distintos de los enteros negativos.

Leonhard Euler posteriormente dio dos definiciones diferentes: la primera no era su integral, sino un producto infinito que está bien definido para todos los números complejos n distintos de los enteros negativos,norte¡=k=1(1+1k)norte1+nortek,{\displaystyle n!=\prod _{k=1}^{\infty }{\frac {\left(1+{\frac {1}{k}}\right)^{n}}{1+{\frac {n}{k}}}}\,,}De lo cual informó a Goldbach en una carta fechada el 13 de octubre de 1729. Le escribió nuevamente a Goldbach el 8 de enero de 1730 para anunciar su descubrimiento de la representación integral.norte¡=01(registros)norteds,{\displaystyle n!=\int _{0}^{1}(-\log s)^{n}\,ds,} lo cual es válido cuando la parte real del número complejonorte{\displaystyle n}es estrictamente mayor que1{\displaystyle -1}(es decir ,(norte)>1{\displaystyle \Re (n)>-1}) . Mediante el cambio de variablest=lns{\displaystyle t=-\ln s} , esto se convierte en la familiar integral de Euler. Euler publicó sus resultados en el artículo "De progressibus transcendentibus seu quarum termini generales algebraice dari nequeunt" ("Sobre progresiones trascendentales, es decir, aquellas cuyos términos generales no se pueden dar algebraicamente"), presentado a la Academia de San Petersburgo el 28 de noviembre de 1729. [ 62 ] Euler descubrió además algunas de las propiedades funcionales importantes de la función gamma, incluida la fórmula de reflexión.

James Stirling , contemporáneo de Euler, también intentó encontrar una expresión continua para el factorial y llegó a lo que ahora se conoce como la fórmula de Stirling . Aunque la fórmula de Stirling da una buena estimación denorte¡{\displaystyle n!} , incluso para números no enteros, no proporciona el valor exacto. Stirling y Jacques Philippe Marie Binet aportaron extensiones de su fórmula que corrigen el error .

Siglo XIX: Gauss, Weierstrass y Legendre

De progresivoibus trascendentibus, seu quarum termini generales algebraicae dari nequeunt
La primera página del artículo de Euler

Carl Friedrich Gauss reescribió el producto de Euler comoΓ(z)=límitemetrometrozmetro¡z(z+1)(z+2)(z+metro){\displaystyle \Gamma (z)=\lim _{m\to \infty }{\frac {m^{z}m!}{z(z+1)(z+2)\cdots (z+m)}}}y utilizó esta fórmula para descubrir nuevas propiedades de la función gamma. Aunque Euler fue un pionero en la teoría de variables complejas, no parece haber considerado el factorial de un número complejo, como sí lo hizo Gauss. [ 63 ] Gauss también demostró el teorema de multiplicación de la función gamma e investigó la conexión entre la función gamma y las integrales elípticas .

Karl Weierstrass estableció además el papel de la función gamma en el análisis complejo , partiendo de otra representación de producto,Γ(z)=miγzzk=1(1+zk)1mizk,{\displaystyle \Gamma (z)={\frac {e^{-\gamma z}}{z}}\prod _{k=1}^{\infty }\left(1+{\frac {z}{k}}\right)^{-1}e^{\frac {z}{k}},}dóndeγ{\displaystyle \gamma }es la constante de Euler-Mascheroni . Weierstrass originalmente escribió su producto como uno para 1Γ{\displaystyle \textstyle {\frac {1}{\Gamma }}} , en cuyo caso se toma sobre los ceros de la función en lugar de sobre sus polos. Inspirado por este resultado, demostró lo que se conoce como el teorema de factorización de Weierstrass : que cualquier función entera puede escribirse como un producto sobre sus ceros en el plano complejo; una generalización del teorema fundamental del álgebra .

El nombre función gamma y el símboloΓ{\displaystyle \Gamma }fueron introducidas por Adrien-Marie Legendre alrededor de 1811; Legendre también reescribió la definición integral de Euler en su forma moderna. Aunque el símbolo es una gamma griega mayúscula, no existe un estándar aceptado sobre si el nombre de la función debe escribirse "gamma function" o "Gamma function" (algunos autores simplemente escriben " Γ{\displaystyle \Gamma }-función "). La notación alternativa de "función pi"Π(z)=z¡{\displaystyle \Pi (z)=z!}La notación de Gauss se encuentra a veces en la literatura antigua, pero la notación de Legendre predomina en las obras modernas.

Es justificable preguntar por qué distinguimos entre el "factorial ordinario" y la función gamma mediante el uso de símbolos distintos, y en particular por qué la función gamma debe normalizarse aΓ(norte+1)=norte¡{\displaystyle \Gamma (n+1)=n!}en lugar de simplemente usarΓ(norte)=norte¡{\displaystyle \Gamma (n)=n!} . Consideremos que la notación para exponentes,incógnitanorte{\displaystyle x^{n}} , se ha generalizado de números enteros a números complejosincógnitaz{\displaystyle x^{z}}sin ningún cambio. Se desconoce la motivación de Legendre para la normalización, y algunos la han criticado por engorrosa (el matemático del siglo XX Cornelius Lanczos , por ejemplo, la llamó "carente de toda racionalidad" y en su lugar usaría z¡{\displaystyle z!} ). [ 64 ] La normalización de Legendre simplifica algunas fórmulas, pero complica otras. Desde un punto de vista moderno, la normalización de Legendre de la función gamma es la integral del carácter aditivomiincógnita{\displaystyle e^{-x}}contra el carácter multiplicativoincógnitaz{\displaystyle x^{z}} con respecto a la medida Haardincógnitaincógnita{\textstyle {\frac {dx}{x}}}en el grupo de LieR+{\displaystyle \mathbb {R} ^{+}} . Por lo tanto, esta normalización deja más claro que la función gamma es un análogo continuo de una suma de Gauss . [ 65 ]

Siglos XIX-XX: caracterización de la función gamma

Resulta algo problemático que se hayan dado numerosas definiciones para la función gamma. Si bien describen la misma función, no es del todo sencillo demostrar su equivalencia. Stirling nunca demostró que su fórmula extendida se corresponde exactamente con la función gamma de Euler; la primera demostración la dio Charles Hermite en 1900. [ 66 ] En lugar de buscar una demostración específica para cada fórmula, sería deseable contar con un método general para identificar la función gamma.

Una forma de demostrar la equivalencia sería encontrar una ecuación diferencial que caracterice la función gamma. La mayoría de las funciones especiales en matemáticas aplicadas surgen como soluciones de ecuaciones diferenciales, cuyas soluciones son únicas. Sin embargo, la función gamma no parece satisfacer ninguna ecuación diferencial simple. Otto Hölder demostró en 1887 que la función gamma, al menos, no satisface ninguna ecuación diferencial algebraica , al mostrar que una solución a dicha ecuación no podría satisfacer la fórmula de recurrencia de la función gamma, lo que la convierte en una función trascendentalmente trascendental . Este resultado se conoce como el teorema de Hölder .

Una caracterización definitiva y de aplicación general de la función gamma no se dio hasta 1922. Harald Bohr y Johannes Mollerup demostraron entonces lo que se conoce como el teorema de Bohr-Mollerup : que la función gamma es la solución única de la relación de recurrencia factorial que es positiva y logarítmicamente convexa para γ positivo .z{\displaystyle z}y cuyo valor en1{\displaystyle 1}es1{\displaystyle 1}( Una función es logarítmicamente convexa si su logaritmo es convexo). Otra caracterización la proporciona el teorema de Wielandt .

El teorema de Bohr-Mollerup es útil porque permite demostrar con relativa facilidad la convexidad logarítmica de cualquiera de las fórmulas utilizadas para definir la función gamma. Además, en lugar de definir la función gamma mediante una fórmula específica, podemos elegir las condiciones del teorema de Bohr-Mollerup como definición y, a partir de ahí, seleccionar cualquier fórmula que satisfaga dichas condiciones como punto de partida para el estudio de la función gamma. Este enfoque fue utilizado por el grupo de Bourbaki .

Borwein y Corless repasan tres siglos de trabajo sobre la función gamma. [ 67 ]

Tablas de referencia y software

Un gráfico dibujado a mano del valor absoluto de la función gamma compleja, de las Tablas de funciones superiores de Jahnke y Emde.

Aunque la función gamma se puede calcular prácticamente con la misma facilidad que cualquier función matemáticamente más sencilla con una computadora moderna —incluso con una calculadora de bolsillo programable—, por supuesto, no siempre fue así. Hasta mediados del siglo XX, los matemáticos se basaban en tablas hechas a mano; en el caso de la función gamma, en particular una tabla calculada por Gauss en 1813 y otra por Legendre en 1825. [ 68 ]

Las tablas de valores complejos de la función gamma, así como los gráficos dibujados a mano, se presentaron en Tablas de funciones con fórmulas y curvas de Jahnke y Emde , publicadas por primera vez en Alemania en 1909. Según Michael Berry , «la publicación en J&E de un gráfico tridimensional que mostraba los polos de la función gamma en el plano complejo adquirió un estatus casi icónico». [ 69 ]

De hecho, hasta la década de 1930, cuando se descubrieron aplicaciones para la función gamma compleja en la física teórica, existía poca necesidad práctica de utilizar otros valores que no fueran reales de la función gamma. Con la llegada de las computadoras electrónicas para la elaboración de tablas en la década de 1950, se publicaron varias tablas extensas para la función gamma compleja con el fin de satisfacer la demanda, incluyendo una tabla con una precisión de 12 decimales de la Oficina Nacional de Estándares de EE . UU. [ 1 ]

Reproducción de una famosa gráfica compleja Γ(z) de Janhke y Emde de la función gamma para -4,5 < Re(z) < 4,5 y -2,5 < Im(z) < 2,5. (Tablas de funciones con fórmulas y curvas, 4.ª ed., Dover, 1945)
Reproducción de una famosa trama compleja de Γ(z){\displaystyle \Gamma (z)} por Janhke y Emde de la función gamma para4.5<(z)<4.5{\displaystyle -4.5<\Re (z)<4.5}y2.5<(z)<2.5{\displaystyle -2.5<\Im (z)<2.5}( Tablas de funciones con fórmulas y curvas , 4.ª ed., Dover, 1945)

Las implementaciones de punto flotante de doble precisión de la función gamma y su logaritmo están disponibles en la mayoría de los programas de computación científica y bibliotecas de funciones especiales, por ejemplo, TK Solver , Matlab , GNU Octave y la Biblioteca Científica GNU . La función gamma también se agregó a la biblioteca estándar de C ( math.h ). Las implementaciones de precisión arbitraria están disponibles en la mayoría de los sistemas de álgebra computacional , como Mathematica y Maple . PARI/GP , MPFR y MPFUN contienen implementaciones gratuitas de precisión arbitraria. En algunas calculadoras de software , como la Calculadora de Windows y la Calculadora GNOME , la función factorial devuelveΓ(incógnita+1){\displaystyle \Gamma (x+1)}cuando la entradaincógnita{\displaystyle x}es un valor no entero. [ 70 ] [ 71 ]

Véase también

Notas

  1. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 Davis, PJ (1959). " La integral de Leonhard Euler: un perfil histórico de la función gamma" . American Mathematical Monthly . 66 (10): 849– 869. doi : 10.2307/2309786 . JSTOR 2309786. Archivado del original el 7 de noviembre de 2012. Recuperado el 3 de diciembre de 2016 . 
  2. "¿Está mal definida la función Gamma? O: Hadamard contra Euler: ¿Quién encontró la mejor función Gamma?" .
  3. Beals, Richard; Wong, Roderick (2010). Funciones especiales: Un texto para graduados . Cambridge University Press. pág. 28. ISBN  978-1-139-49043-6.Extracto de la página 28
  4. Ross, Clay C. (2013). Ecuaciones diferenciales: Una introducción con Mathematica ( edición ilustrada). Springer Science & Business Media. pág. 293. ISBN   978-1-4757-3949-7.Expresión G.2 en la página 293
  5. Kingman, JFC (1961). "Una propiedad de convexidad de matrices positivas". The Quarterly Journal of Mathematics . 12 (1): 283– 284. Bibcode : 1961QJMat..12..283K . doi : 10.1093/qmath/12.1.283 .
  6. Weisstein, Eric W. "Teorema de Bohr-Mollerup" . MathWorld .
  7. Davis, Philip. "La integral de Leonhard Euler: un perfil histórico de la función gamma" (PDF) . maa.org .
  8. ^ Bonvini, Marco (9 de octubre de 2010). "La función Gamma" (PDF) . Roma1.infn.it .
  9. Askey, RA ; Roy, R. (2010), "Expansiones en serie" , en Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.), NIST Handbook of Mathematical Functions , Cambridge University Press, ISBN 978-0-521-19225-5, MR 2723248 .
  10. Waldschmidt, M. (2006). "Trascendencia de los períodos: estado del arte" (PDF) . Pure Appl. Math. Quart . 2 (2): 435– 463. doi : 10.4310/pamq.2006.v2.n2.a3 . Archivado (PDF) del original el 6 de mayo de 2006.
  11. "¿Cómo obtener la expansión de Laurent de la función gamma alrededor de $z=0$?" . Mathematics Stack Exchange . Consultado el 17 de agosto de 2022 .
  12. Artin, Emil (2015). La función gamma . Dover. pág. 24. 
  13. Oldham, Keith; Myland, Jan; Spanier, Jerome (2010). "Capítulo 43 - La función gammaΓ(ν){\displaystyle \Gamma (\nu )}Atlas de funciones (2.ª  ed.). Nueva York, NY: Springer Science & Business Media . ISBN 9780387488073.
  14. 1 2 Weisstein, Eric W. "Función Gamma" . MathWorld .
  15. Sloane, N. J. A. (ed.). "Secuencia A030169 (Expansión decimal de un número real x tal que y = Gamma(x) es un mínimo)" . La enciclopedia en línea de secuencias de enteros . Fundación OEIS.  
  16. Sloane, N. J. A. (ed.). "Secuencia A030171 (Expansión decimal del número real y tal que y = Gamma(x) es un mínimo)" . La enciclopedia en línea de secuencias de enteros . Fundación OEIS.  
  17. Sloane, N. J. A. (ed.). "Secuencia A178840 (Expansión decimal del factorial de la proporción áurea)" . La enciclopedia en línea de secuencias de enteros . Fundación OEIS.  
  18. Sloane, N. J. A. (ed.). "Secuencia A175472 (Expansión decimal del valor absoluto de la abscisa del máximo local de la función Gamma en el intervalo [ -1,0 ] )" . La enciclopedia en línea de secuencias enteras . Fundación OEIS.  
  19. Sloane, N. J. A. (ed.). "Secuencia A175473 (Expansión decimal del valor absoluto de la abscisa del mínimo local de la función Gamma en el intervalo [ -2,-1 ] )" . La enciclopedia en línea de secuencias de enteros . Fundación OEIS.  
  20. Sloane, N. J. A. (ed.). "Secuencia A175474 (Desarrollo decimal del valor absoluto de la abscisa del máximo local de la función Gamma en el intervalo [ -3,-2 ] )" . La enciclopedia en línea de secuencias enteras . Fundación OEIS.  
  21. Sloane, N. J. A. (ed.). "Secuencia A256681 (Expansión decimal de la abscisa [ negada ] del mínimo local de la función Gamma en el intervalo [ -4,-3 ] )" . La enciclopedia en línea de secuencias de enteros . Fundación OEIS.  
  22. Sloane, N. J. A. (ed.). "Secuencia A256682 (Expansión decimal de la abscisa [ negada ] del máximo local de la función Gamma en el intervalo [ -5,-4 ] )" . La enciclopedia en línea de secuencias de enteros . Fundación OEIS.  
  23. Gradshteyn, IS; Ryzhik, IM (2007). Tabla de integrales, series y productos (séptima ed.). Academic Press. pág. 893. ISBN   978-0-12-373637-6.
  24. Whittaker y Watson, 12.2 ejemplo 1.
  25. Detlef, Gronau. "¿Por qué la función gamma es como es?" (PDF) . Imsc.uni-graz.at .
  26. Pascal Sebah, Xavier Gourdon. "Introducción a la función gamma" (PDF) . Numbers Computation . Archivado del original (PDF) el 30 de enero de 2023. Consultado el 30 de enero de 2023 .
  27. Whittaker y Watson, 12.31.
  28. Whittaker y Watson, 12.32.
  29. Whittaker y Watson, 12.22.
  30. "Integral exponencial E: Representaciones de fracciones continuas (Fórmula 06.34.10.0005)" .
  31. "Integral exponencial E: Representaciones de fracciones continuas (Fórmula 06.34.10.0003)" .
  32. Bateman, Harry; Erdélyi, Arthur (1955). Funciones trascendentales superiores . McGraw-Hill. OCLC 627135 . 
  33. Srivastava, HM; Choi, J. (2001). Series asociadas con la función zeta y funciones relacionadas . Países Bajos: Kluwer Academic. ISBN 0-7923-7054-6.
  34. Blagouchine, Iaroslav V. (2014). "Redescubrimiento de las integrales de Malmsten, su evaluación mediante métodos de integración de contorno y algunos resultados relacionados" . Ramanujan J. 35 ( 1): 21–110 . doi : 10.1007/s11139-013-9528-5 . S2CID 120943474 . 
  35. Blagouchine, Iaroslav V. (2016). "Fe de erratas y adenda a "Redescubrimiento de las integrales de Malmsten, su evaluación mediante métodos de integración de contorno y algunos resultados relacionados"". Ramanujan J . 42 (3): 777– 781. doi : 10.1007/s11139-015-9763-z . S2CID 125198685 . 
  36. Sloane, N. J. A. (ed.). "Secuencia A245886 (Expansión decimal de Gamma(-3/2), donde Gamma es la función gamma de Euler)" . La enciclopedia en línea de secuencias enteras . Fundación OEIS.  
  37. Sloane, N. J. A. (ed.). "Secuencia A019707 (Expansión decimal de sqrt(Pi)/5)" . La enciclopedia en línea de secuencias de enteros . Fundación OEIS.  
  38. Sloane, N. J. A. (ed.). "Secuencia A002161 (Expansión decimal de la raíz cuadrada de Pi)" . La enciclopedia en línea de secuencias de enteros . Fundación OEIS.  
  39. Sloane, N. J. A. (ed.). "Secuencia A019704 (Expansión decimal de sqrt(Pi)/2)" . La enciclopedia en línea de secuencias de enteros . Fundación OEIS.  
  40. Sloane, N. J. A. (ed.). "Secuencia A245884 (Expansión decimal de Gamma(5/2), donde Gamma es la función gamma de Euler)" . La enciclopedia en línea de secuencias enteras . Fundación OEIS.  
  41. Sloane, N. J. A. (ed.). "Secuencia A245885 (Expansión decimal de Gamma(7/2), donde Gamma es la función gamma de Euler)" . La enciclopedia en línea de secuencias enteras . Fundación OEIS.  
  42. "Función logarítmica gamma" . Wolfram MathWorld . Consultado el 3 de enero de 2019 .
  43. "La integral de Leonhard Euler: un perfil histórico de la función gamma" (PDF) . Archivado (PDF) del original el 12 de septiembre de 2014. Consultado el 11 de abril de 2022 .
  44. Blagouchine, Iaroslav V. (2015). "Un teorema para la evaluación en forma cerrada de la primera constante generalizada de Stieltjes en argumentos racionales y algunas sumas relacionadas". Journal of Number Theory . 148 : 537–592 . arXiv : 1401.3724 . doi : 10.1016/j.jnt.2014.08.009 .
  45. ^ Alexejewsky, WP (1894). "Über eine Classe von Funktionen, die der Gammafunktion analog sind" [ Sobre una clase de funciones análogas a la función gamma ] . Leipzig Weidmannsche Buchhandlung . 46 : 268-275 .
  46. Barnes, EW (1899). "La teoría de la función G ". Quart. J. Math . 31 : 264–314 .
  47. Adamchik, Victor S. (1998). "Funciones poligamma de orden negativo" . J. Comput. Appl. Math . 100 (2): 191– 199. doi : 10.1016/S0377-0427(98)00192-7 .
  48. Gosper, RW (1997). "norte/4metro/6registroF(z)dz{\displaystyle \textstyle \int _{n/4}^{m/6}\log F(z)\,dz}en funciones especiales, series q y temas relacionados". J. Am. Math. Soc . 14 .
  49. Espinosa, Olivier; Moll, Victor H. (2002). "Sobre algunas integrales que involucran la función zeta de Hurwitz: Parte 1". The Ramanujan Journal . 6 (2): 159– 188. doi : 10.1023/A:1015706300169 . S2CID 128246166 . 
  50. Bailey, David H.; Borwein, David; Borwein, Jonathan M. (2015). "Sobre las integrales log-gamma eulerianas y las funciones zeta de Tornheim-Witten". The Ramanujan Journal . 36 ( 1–2 ): 43–68 . doi : 10.1007/s11139-012-9427-1 . S2CID 7335291 . 
  51. Amdeberhan, T.; Coffey, Mark W.; Espinosa, Olivier; Koutschan, Christoph; Manna, Dante V.; Moll, Victor H. (2011). "Integrales de potencias de loggamma" . Proc. Amer. Math. Soc . 139 (2): 535– 545. doi : 10.1090/S0002-9939-2010-10589-0 .
  52. EA Karatsuba, Evaluación rápida de funciones trascendentales. Probl. Inf. Transm. Vol.27, No. 4, pp. 339–360 (1991).
  53. EA Karatsuba, Sobre un nuevo método para la evaluación rápida de funciones trascendentales. Russ. Math. Surv. Vol.46, No. 2, pp. 246–247 (1991).
  54. EA Karatsuba " Algoritmos rápidos y el método FEE ".
  55. Borwein, JM; Zucker, IJ (1992). "Evaluación rápida de la función gamma para fracciones racionales pequeñas utilizando integrales elípticas completas de primer tipo". IMA Journal of Numerical Analysis . 12 (4): 519– 526. doi : 10.1093/IMANUM/12.4.519 .
  56. Werner, Helmut; Collinge, Robert (1961). "Aproximaciones de Chebyshev a la función gamma". Math. Comput . 15 (74): 195– 197. doi : 10.1090/S0025-5718-61-99220-1 . JSTOR 2004230 . 
  57. Michon, GP " Trigonometría y funciones básicas Archivado el 9 de enero de 2010 en Wayback Machine ". Numericana . Consultado el 5 de mayo de 2007.
  58. Chaudry, MA; Zubair, SM (2001). Sobre una clase de funciones gamma incompletas con aplicaciones . Boca Raton: CRC Press. pág. 37. ISBN  1-58488-143-7.
  59. Rice, JA (1995). Estadística matemática y análisis de datos (Segunda edición). Belmont: Duxbury Press. págs. 52–53 . ISBN   0-534-20934-3.
  60. Borwein, J.; Bailey, DH y Girgensohn, R. (2003). Experimentación en matemáticas . AK Peters. pág. 133. ISBN  978-1-56881-136-9.
  61. "Interpolación del factorial natural n ! o El nacimiento de la función factorial real (1729–1826)" .
  62. ^ El artículo de Euler fue publicado en Commentarii academiae scientiarum Petropolitanae 5, 1738, 36–57. Véase E19 – De progresivoibus trascendentibus seu quarum termini generales algebraice dari nequeunt , de The Euler Archive, que incluye una copia escaneada del artículo original.
  63. Remmert, R. (2006). Temas clásicos en la teoría de funciones complejas . Traducido por Kay, LD Springer. ISBN 978-0-387-98221-2.
  64. Lanczos, C. (1964). "Una aproximación precisa de la función gamma". Journal of the Society for Industrial and Applied Mathematics, Series B: Numerical Analysis . 1 (1): 86. Bibcode : 1964SJNA....1...86L . doi : 10.1137/0701008 .
  65. Ilker Inam; Engin Büyükaşşk (2019). Notas de la Escuela Internacional de Otoño sobre Teoría Computacional de Números . Springer. pág. 205. ISBN  978-3-030-12558-5.Extracto de la página 205
  66. Knuth, DE (1997). El arte de la programación informática . Vol. 1 (Algoritmos fundamentales). Addison-Wesley. ISBN  0-201-89683-4.
  67. Borwein, Jonathan M. ; Corless, Robert M. (2017). "Gamma y factorial en el mensual". American Mathematical Monthly . 125 (5). Mathematical Association of America: 400– 24. arXiv : 1703.05349 . Bibcode : 2017arXiv170305349B . doi : 10.1080/00029890.2018.1420983 . S2CID 119324101 . 
  68. "¿Cuál es la historia de la función Gamma?" . yearis.com . Consultado el 5 de noviembre de 2022 .
  69. Berry, M. (abril de 2001). "¿Por qué son especiales las funciones especiales?" . Physics Today .
  70. "microsoft/calculator" . GitHub . Consultado el 25 de diciembre de 2020 .
  71. "gnome-calculator" . GNOME.org . Consultado el 3 de marzo de 2023 .
  • Este artículo incorpora material del artículo de Citizendium " Función Gamma ", que está bajo la licencia Creative Commons Attribution-ShareAlike 3.0 Unported , pero no bajo la GFDL .

Lecturas adicionales

  • Abramowitz, Milton; Stegun, Irene A., eds. (1972). "Capítulo 6" . Manual de funciones matemáticas con fórmulas, gráficas y tablas matemáticas . Nueva York: Dover.
  • Andrews, GE ; Askey, R.; Roy, R. (1999). «Capítulo 1 (Funciones gamma y beta)». Funciones especiales . Nueva York: Cambridge University Press. ISBN 978-0-521-78988-2.
  • Artin, Emil (2006). "La función gamma". En Rosen, Michael (ed.). Exposición de Emil Artin: una selección . Historia de las matemáticas. Vol.  30. Providence, RI: American Mathematical Society.
  • Askey, R.; Roy, R. (2010), "Función gamma" , en Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.), NIST Handbook of Mathematical Functions , Cambridge University Press, ISBN 978-0-521-19225-5, MR 2723248 .
  • Birkhoff, George D. (1913). "Nota sobre la función gamma" . Bull. Amer. Math. Soc . 20 (1): 1– 10. doi : 10.1090/s0002-9904-1913-02429-7 . MR 1559418 . 
  • Böhmer, PE (1939). Differenzengleichungen und bestimmte Integrale [ Ecuaciones diferenciales e integrales definidas ] . Leipzig: Köhler Verlag.
  • Davis, Philip J. (1959). "La integral de Leonhard Euler: un perfil histórico de la función gamma". American Mathematical Monthly . 66 (10): 849– 869. doi : 10.2307/2309786 . JSTOR 2309786 . 
  • Post, Emil (1919). "Las funciones gamma generalizadas". Anales de Matemáticas . Segunda serie. 20 (3): 202– 217. doi : 10.2307/1967871 . JSTOR 1967871 . 
  • Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007). «Sección 6.1. Función Gamma» . Numerical Recipes: The Art of Scientific Computing (3.ª  ed.). Nueva York: Cambridge University Press. ISBN 978-0-521-88068-8.
  • Rocktäschel, Oregón (1922). Methoden zur Berechnung der Gammafunktion für komplexes Argument [ Métodos para calcular la función gamma para argumentos complejos ] . Dresde: Universidad Técnica de Dresde .
  • Temme, Nico M. (1996). Funciones especiales: Una introducción a las funciones clásicas de la física matemática . Nueva York: John Wiley & Sons. ISBN 978-0-471-11313-3.
  • Whittaker, ET ; Watson, GN (1927). Un curso de análisis moderno . Cambridge University Press.ISBN 978-0-521-58807-2
  • Li, Xin; Chen, Chao-Ping (2017). " Aproximante de Padé relacionado con la asintótica de la función gamma" . J. Inequal. Applic . 2017 (1): 53. doi : 10.1186/s13660-017-1315-1 . PMC 5331117. PMID 28303079 .  
  • Biblioteca Digital de Funciones Matemáticas del NIST: Función gamma
  • Pascal Sebah y Xavier Gourdon. Introducción a la función gamma . En formatos PostScript y HTML .
  • Referencia de C++ parastd::tgamma
  • gamma()y lgamma()expuesto desde C99 en Postgres versión 18
  • En Exampleproblems.com se pueden encontrar ejemplos de problemas que involucran la función gamma .
  • "Función gamma" , Enciclopedia de Matemáticas , EMS Press , 2001 [1994]
  • Evaluador de la función gamma de Wolfram (precisión arbitraria). Archivado el 28 de octubre de 2019 en Wayback Machine.
  • "Gamma" . Sitio web de Wolfram Functions.
  • Volumen de n-esferas y la función gamma en MathPages