Articulo de referencia

Función gamma incompleta

La función gamma incompleta superior para algunos valores de s: 0 (azul), 1 (rojo), 2 (verde), 3 (naranja), 4 (púrpura). Representación gráfica de la función gamma incompleta re...

La función gamma incompleta superior para algunos valores de s: 0 (azul), 1 (rojo), 2 (verde), 3 (naranja), 4 (púrpura).
Representación gráfica de la función gamma incompleta regularizada Q(2,z) en el plano complejo desde -2-2i hasta 2+2i con colores creados con la función ComplexPlot3D de Mathematica 13.1.
Representación gráfica de la función gamma incompleta regularizada Q(2,z) en el plano complejo desde -2-2i hasta 2+2i con colores creados con la función ComplexPlot3D de Mathematica 13.1.

En matemáticas , las funciones gamma incompletas superior e inferior son tipos de funciones especiales que surgen como soluciones a diversos problemas matemáticos, como ciertas integrales .

Sus nombres respectivos derivan de sus definiciones integrales, que se definen de forma similar a la función gamma, pero con límites integrales diferentes o "incompletos". La función gamma se define como una integral de cero a infinito. Esto contrasta con la función gamma incompleta inferior, que se define como una integral de cero a un límite superior variable. De manera similar, la función gamma incompleta superior se define como una integral de un límite inferior variable a infinito.

Definición

La función gamma incompleta superior se define como: Γ(s,incógnita)=incógnitats1mitdt,{\displaystyle \Gamma (s,x)=\int _{x}^{\infty }t^{s-1}\,e^{-t}\,dt,} mientras que la función gamma incompleta inferior se define como: γ(s,incógnita)=0incógnitats1mitdt.{\displaystyle \gamma (s,x)=\int _{0}^{x}t^{s-1}\,e^{-t}\,dt.} En ambos casos, s es un parámetro complejo, de modo que la parte real de s es positiva.

Propiedades

Mediante la integración por partes encontramos las relaciones de recurrencia.Γ(s+1,incógnita)=sΓ(s,incógnita)+incógnitasmiincógnita{\displaystyle \Gamma (s+1,x)=s\Gamma (s,x)+x^{s}e^{-x}} y γ(s+1,incógnita)=sγ(s,incógnita)incógnitasmiincógnita.{\displaystyle \gamma (s+1,x)=s\gamma (s,x)-x^{s}e^{-x}.} Dado que la función gamma ordinaria se define como Γ(s)=0ts1mitdt{\displaystyle \Gamma (s)=\int _{0}^{\infty }t^{s-1}\,e^{-t}\,dt} tenemos Γ(s)=Γ(s,0)=límiteincógnitaγ(s,incógnita){\displaystyle \Gamma (s)=\Gamma (s,0)=\lim _{x\to \infty }\gamma (s,x)} y γ(s,incógnita)+Γ(s,incógnita)=Γ(s).{\displaystyle \gamma (s,x)+\Gamma (s,x)=\Gamma (s).}

Continuación de valores complejos

La función gamma incompleta inferior y la función gamma incompleta superior, definidas anteriormente para s y x reales positivos , pueden desarrollarse en funciones holomorfas , con respecto a x y s , definidas para casi todas las combinaciones de x y s complejos . [ 1 ] El análisis complejo muestra cómo las propiedades de las funciones gamma incompletas reales se extienden a sus contrapartes holomorfas.

Función gamma incompleta inferior

Extensión holomorfa

La aplicación repetida de la relación de recurrencia para la función gamma incompleta inferior conduce a la expansión en serie de potencias : [ 2 ]γ(s,incógnita)=k=0incógnitasmiincógnitaincógnitaks(s+1)(s+k)=incógnitasΓ(s)miincógnitak=0incógnitakΓ(s+k+1).{\displaystyle \gamma (s,x)=\sum _{k=0}^{\infty }{\frac {x^{s}e^{-x}x^{k}}{s(s+1)\cdots (s+k)}}=x^{s}\,\Gamma (s)\,e^{-x}\sum _{k=0}^{\infty }{\frac {x^{k}}{\Gamma (s+k+1)}}.} Dado el rápido crecimiento en valor absoluto de Γ( z + k ) cuando k → ∞ , y el hecho de que el recíproco de Γ( z ) es una función entera , los coeficientes en la suma más a la derecha están bien definidos, y localmente la suma converge uniformemente para todos los complejos s y x . Por un teorema de Weierstrass , [ 3 ] la función límite, a veces denotada comoγ{\displaystyle \gamma ^{*}}, [ 4 ]γ(s,z):=mizk=0zkΓ(s+k+1){\displaystyle \gamma ^{*}(s,z):=e^{-z}\sum _{k=0}^{\infty }{\frac {z^{k}}{\Gamma (s+k+1)}}} es entera con respecto a z (para s fijo ) y s (para z fijo ), [ 1 ] y, por lo tanto, holomorfa en C × C por el teorema de Hartogs . [ 5 ] Por lo tanto, la siguiente descomposición [ 1 ]γ(s,z)=zsΓ(s)γ(s,z),{\displaystyle \gamma (s,z)=z^{s}\,\Gamma (s)\,\gamma ^{*}(s,z),} extiende la función gamma incompleta inferior real como una función holomorfa , tanto conjunta como separadamente en z y s . Se deduce de las propiedades dezs{\displaystyle z^{s}}y la función Γ , que los dos primeros factores capturan las singularidades deγ(s,z){\displaystyle \gamma (s,z)}(en z = 0 o s un entero no positivo), mientras que el último factor contribuye a sus ceros.

Multivaloración

El logaritmo complejo log z = log | z | + i arg z se determina salvo un múltiplo de 2 πi , lo que lo convierte en multivaluado . Las funciones que involucran el logaritmo complejo suelen heredar esta propiedad. Entre ellas se encuentran la potencia compleja y, dado que z s aparece en su descomposición, también la función γ .

La indeterminación de las funciones multivaluadas introduce complicaciones, ya que debe especificarse cómo seleccionar un valor. Las estrategias para abordar esto son:

  • (La forma más general) reemplazar el dominio C de las funciones multivaluadas por una variedad adecuada en C × C llamada superficie de Riemann . Si bien esto elimina la multivaluación, es necesario conocer la teoría subyacente; [ 6 ]
  • restringir el dominio de tal manera que una función multivaluada se descomponga en ramas univaluadas separadas , que pueden manejarse individualmente.

Las siguientes reglas pueden utilizarse para interpretar correctamente las fórmulas de esta sección. Salvo que se indique lo contrario, se asume lo siguiente:

Sectores

Los sectores en C con vértice en z = 0 suelen ser dominios apropiados para expresiones complejas. Un sector D consta de todos los z complejos que satisfacen z ≠ 0 y αδ < arg z < α + δ con algún α y 0 < δπ . A menudo, α puede elegirse arbitrariamente y no se especifica. Si δ no se da, se asume que es π , y el sector es de hecho todo el plano C , con la excepción de una semirrecta que parte de z = 0 y apunta en la dirección de α , que suele servir como corte de rama . Nota: En muchas aplicaciones y textos, α se toma implícitamente como 0, lo que centra el sector alrededor del eje real positivo.

Sucursales

En particular, existe un logaritmo holomorfo unívoco en cualquier sector D cuya parte imaginaria está acotada al intervalo ( αδ , α + δ ) . Basándose en este logaritmo restringido, z s y las funciones gamma incompletas se reducen a funciones holomorfas unívocas en D (o C × D ), llamadas ramas de sus contrapartes multivaluadas en D. Si se añade un múltiplo de 2 π a α, se obtiene un conjunto diferente de ramas correlacionadas en el mismo conjunto D . Sin embargo, en cualquier contexto dado, se supone que α es fijo y todas las ramas involucradas están asociadas a él. Si | α | < δ , las ramas se denominan principales , porque son iguales a sus análogos reales en el eje real positivo. Nota: En muchas aplicaciones y textos, las fórmulas solo son válidas para las ramas principales.

Relación entre ramas

Los valores de las diferentes ramas tanto de la función de potencia compleja como de la función gamma incompleta inferior se pueden derivar entre sí mediante la multiplicación demi2πiks{\displaystyle e^{2\pi iks}}, [ 1 ] para k un entero adecuado.

Comportamiento cerca del punto de ramificación

La descomposición anterior muestra además que γ se comporta cerca de z = 0 asintóticamente como: γ(s,z)zsΓ(s)γ(s,0)=zsΓ(s)/Γ(s+1)=zs/s.{\displaystyle \gamma (s,z)\asymp z^{s}\,\Gamma (s)\,\gamma ^{*}(s,0)=z^{s}\,\Gamma (s)/\Gamma (s+1)=z^{s}/s.}

Para valores reales positivos de x , y y s , x y /y → 0 , cuando ( x , y ) → (0, s ) . Esto parece justificar establecer γ ( s , 0) = 0 para s real > 0 . Sin embargo, las cosas son algo diferentes en el ámbito complejo. Solo si (a) la parte real de s es positiva, y (b) los valores u v se toman de un conjunto finito de ramas, se garantiza que convergen a cero cuando ( u , v ) → (0, s ) , y también γ ( u , v ) . En una sola rama de γ ( b ) se cumple naturalmente, por lo que γ ( s , 0) = 0 para s con parte real positiva es un límite continuo . También hay que tener en cuenta que tal continuación no es en absoluto analítica .

Relaciones algebraicas

Todas las relaciones algebraicas y ecuaciones diferenciales observadas por la función real γ ( s , z ) también se cumplen para su contraparte holomorfa. Esto es consecuencia del teorema de identidad , que establece que las ecuaciones entre funciones holomorfas válidas en un intervalo real se cumplen en todas partes. En particular, la relación de recurrencia [ 2 ] y ∂γ ( s , z )/ ∂z = z s⁻¹ e⁻ z [ 2 ] se conservan en las ramas correspondientes .

Representación integral

La última relación nos dice que, para un s fijo , γ es una primitiva o antiderivada de la función holomorfa z s −1 e z . En consecuencia, para cualquier complejo u , v ≠ 0 , vts1mitdt=γ(s,v)γ(s,){\displaystyle \int _{u}^{v}t^{s-1}\,e^{-t}\,dt=\gamma (s,v)-\gamma (s,u)} se cumple siempre que la trayectoria de integración esté completamente contenida en el dominio de una rama del integrando. Si, además, la parte real de s es positiva, entonces se aplica el límite γ ( s , u ) → 0 para u → 0 , llegando finalmente a la definición de integral compleja de γ [ 1 ].γ(s,z)=0zts1mitdt,(s)>0.{\displaystyle \gamma (s,z)=\int _{0}^{z}t^{s-1}\,e^{-t}\,dt,\,\Re (s)>0.}

Cualquier camino de integración que contenga 0 solo al principio, y que de otro modo esté restringido al dominio de una rama del integrando, es válido aquí, por ejemplo, la línea recta que conecta 0 y z .

Límite para z → +∞
Valores reales

Dada la representación integral de una rama principal de γ , la siguiente ecuación se cumple para todos los s y x reales positivos : [ 7 ]Γ(s)=0ts1mitdt=límiteincógnitaγ(s,incógnita){\displaystyle \Gamma (s)=\int _{0}^{\infty }t^{s-1}\,e^{-t}\,dt=\lim _{x\to \infty }\gamma (s,x)}

complejo s

Este resultado se extiende a s complejos . Supongamos primero que 1 ≤ Re( s ) ≤ 2 y 1 < a < b . Entonces |γ(s,b)γ(s,a)|ab|ts1|mitdt=abts1mitdtabtmitdt{\displaystyle \left|\gamma (s,b)-\gamma (s,a)\right|\leq \int _{a}^{b}\left|t^{s-1}\right|e^{-t}\,dt=\int _{a}^{b}t^{\Re s-1}e^{-t}\,dt\leq \int _{a}^{b}te^{-t}\,dt} donde [ 8 ]|zs|=|z|smisargz{\displaystyle \left|z^{s}\right|=\left|z\right|^{\Re s}\,e^{-\Im s\arg z}} se ha utilizado en el medio. Dado que la integral final se vuelve arbitrariamente pequeña si solo a es suficientemente grande, γ ( s , x ) converge uniformemente para x → ∞ en la franja 1 ≤ Re(s) ≤ 2 hacia una función holomorfa, [ 3 ] que debe ser Γ(s) debido al teorema de identidad. Tomando el límite en la relación de recurrencia γ ( s , x ) = ( s − 1) γ ( s − 1, x ) − x s − 1 e x y observando que lim x n e x = 0 para x → ∞ y todo n , se muestra que γ ( s , x ) también converge fuera de la franja hacia una función que obedece la relación de recurrencia de la función Γ. Se sigue Γ(s)=límiteincógnitaγ(s,incógnita){\displaystyle \Gamma (s)=\lim _{x\to \infty }\gamma (s,x)} para todo s complejo que no sea un entero no positivo, x es real y γ es principal.

Convergencia sectorial

Ahora sea u del sector | arg z | < δ < π /2 con algún δ fijo ( α = 0 ), sea γ la rama principal en este sector, y observemos Γ(s)γ(s,)=Γ(s)γ(s,||)+γ(s,||)γ(s,).{\displaystyle \Gamma (s)-\gamma (s,u)=\Gamma (s)-\gamma (s,|u|)+\gamma (s,|u|)-\gamma (s,u).}

Como se muestra arriba, la primera diferencia puede hacerse arbitrariamente pequeña, si | u | es suficientemente grande. La segunda diferencia permite la siguiente estimación: |γ(s,||)γ(s,)||||zs1miz|dz=|||z|s1misargzmizdz,{\displaystyle \left|\gamma (s,|u|)-\gamma (s,u)\right|\leq \int _{u}^{|u|}\left|z^{s-1}e^{-z}\right|dz=\int _{u}^{|u|}\left|z\right|^{\Re s-1}\,e^{-\Im s\,\arg z}\,e^{-\Re z}\,dz,} donde utilizamos la representación integral de γ y la fórmula sobre | z s | anterior. Si integramos a lo largo del arco con radio R = | u | alrededor de 0 que conecta u y | u | , entonces la última integral es R|arg|Rs1mis|arg|miRporqueargδRsmisδmiRporqueδ=METRO(Rporqueδ)smiRporqueδ{\displaystyle \leq R\left|\arg u\right|R^{\Re s-1}\,e^{\Im s\,|\arg u|}\,e^{-R\cos \arg u}\leq \delta \,R^{\Re s}\,e^{\Im s\,\delta }\,e^{-R\cos \delta }=M\,(R\,\cos \delta )^{\Re s}\,e^{-R\cos \delta }} donde M = δ (cos δ ) −Re s e Im es una constante independiente de u o R . Refiriéndose nuevamente al comportamiento de x n e x para valores grandes de x , vemos que la última expresión se aproxima a 0 a medida que R aumenta hacia . En total, ahora tenemos: Γ(s)=límite|z|γ(s,z),|argz|<π/2ϵ,{\displaystyle \Gamma (s)=\lim _{|z|\to \infty }\gamma (s,z),\quad \left|\arg z\right|<\pi /2-\epsilon ,} Si s no es un entero no negativo, 0 < ε < π /2 es arbitrariamente pequeño, pero fijo, y γ denota la rama principal en este dominio.

Descripción general

γ(s,z){\displaystyle \gamma (s,z)}es:

  • entero en z para un entero positivo fijo s ;
  • holomorfo multivaluado en z para s fijo que no es un entero, con un punto de ramificación en z = 0 ;
  • en cada rama meromorfa en s para z ≠ 0 fijo , con polos simples en enteros no positivos s.

Función gamma incompleta superior

En cuanto a la función gamma incompleta superior , una extensión holomorfa , con respecto a z o s , viene dada por [ 1 ].Γ(s,z)=Γ(s)γ(s,z){\displaystyle \Gamma (s,z)=\Gamma (s)-\gamma (s,z)} en los puntos ( s , z ) , donde existe el lado derecho. Dado queγ{\displaystyle \gamma }es multivaluado, lo mismo ocurre conΓ{\displaystyle \Gamma }, pero una restricción a valores principales solamente produce la rama principal de valor único deΓ{\displaystyle \Gamma }.

Cuando s es un entero no positivo en la ecuación anterior, ninguna parte de la diferencia está definida, y un proceso límite , desarrollado aquí para s → 0 , completa los valores faltantes. El análisis complejo garantiza la holomorfía , porqueΓ(s,z){\displaystyle \Gamma (s,z)}Se demuestra que está acotado en un entorno de ese límite para un z fijo .

Para determinar el límite, la serie de potencias deγ{\displaystyle \gamma ^{*}}en z = 0 es útil. Al reemplazarmiincógnita{\displaystyle e^{-x}}por su serie de potencias en la definición integral deγ{\displaystyle \gamma }, se obtiene (supongamos que x y s son números reales positivos por ahora): γ(s,incógnita)=0incógnitats1mitdt=0incógnitak=0(1)kts+k1k¡dt=k=0(1)kincógnitas+kk¡(s+k)=incógnitask=0(incógnita)kk¡(s+k){\displaystyle {\begin{aligned}\gamma (s,x)&=\int _{0}^{x}t^{s-1}e^{-t}\,dt=\int _{0}^{x}\sum _{k=0}^{\infty }\left(-1\right)^{k}\,{\frac {t^{s+k-1}}{k!}}\,dt\\[1ex]&=\sum _{k=0}^{\infty }\left(-1\right)^{k}\,{\frac {x^{s+k}}{k!(s+k)}}=x^{s}\,\sum _{k=0}^{\infty }{\frac {(-x)^{k}}{k!(s+k)}}\end{aligned}}} o [ 4 ]γ(s,incógnita)=k=0(incógnita)kk¡Γ(s)(s+k),{\displaystyle \gamma ^{*}(s,x)=\sum _{k=0}^{\infty }{\frac {(-x)^{k}}{k!\,\Gamma (s)(s+k)}},} que, como representación en serie de todoγ{\displaystyle \gamma ^{*}}función, converge para todo x complejo (y todo s complejo que no sea un entero no positivo).

Al eliminarse su restricción a valores reales, la serie permite la expansión: γ(s,z)1s=1s+zsk=0(z)kk¡(s+k)=zs1s+zsk=1(z)kk¡(s+k),(s)>1,s0.{\displaystyle {\begin{aligned}\gamma (s,z)-{\frac {1}{s}}&=-{\frac {1}{s}}+z^{s}\,\sum _{k=0}^{\infty }{\frac {(-z)^{k}}{k!(s+k)}}\\[1ex]&={\frac {z^{s}-1}{s}}+z^{s}\,\sum _{k=1}^{\infty }{\frac {\left(-z\right)^{k}}{k!(s+k)}},&\Re (s)>-1,\,s\neq 0.\end{aligned}}}

Cuando s → 0 : [ 9 ]zs1sln(z),Γ(s)1s=1sγ+O(s)1sγ,{\displaystyle {\frac {z^{s}-1}{s}}\to \ln(z),\quad \Gamma (s)-{\frac {1}{s}}={\frac {1}{s}}-\gamma +O(s)-{\frac {1}{s}}\to -\gamma ,} (γ{\displaystyle \gamma }es la constante de Euler-Mascheroni aquí), por lo tanto, Γ(0,z)=límites0(Γ(s)1s(γ(s,z)1s))=γln(z)k=1(z)kk(k¡){\displaystyle {\begin{aligned}\Gamma (0,z)&=\lim _{s\to 0}\left(\Gamma (s)-{\tfrac {1}{s}}-\left(\gamma (s,z)-{\tfrac {1}{s}}\right)\right)\\&=-\gamma -\ln(z)-\sum _{k=1}^{\infty }{\frac {\left(-z\right)^{k}}{k\,(k!)}}\end{aligned}}} es la función límite de la función gamma incompleta superior cuando s → 0 , también conocida como la integral exponencial.mi1(z){\displaystyle E_{1}(z)}. [ 10 ]

Mediante la relación de recurrencia, los valores deΓ(norte,z){\displaystyle \Gamma (-n,z)}para enteros positivos n se puede derivar de este resultado, [ 11 ]Γ(norte,z)=1norte¡(mizznortek=0norte1(1)k(nortek1)¡zk+(1)norteΓ(0,z)){\displaystyle \Gamma (-n,z)={\frac {1}{n!}}\left({\frac {e^{-z}}{z^{n}}}\sum _{k=0}^{n-1}(-1)^{k}(n-k-1)!\,z^{k}+\left(-1\right)^{n}\Gamma (0,z)\right)} Por lo tanto, se demuestra que la función gamma incompleta superior existe y es holomorfa, con respecto a z y s , para todo s y z ≠ 0 .

Γ(s,z){\displaystyle \Gamma (s,z)}es:

  • entero en z para una integral positiva fija s ;
  • holomorfo multivaluado en z para s fijo distinto de cero y no un entero positivo, con un punto de ramificación en z = 0 ;
  • igual aΓ(s){\displaystyle \Gamma (s)}para s con parte real positiva y z = 0 (el límite cuando(si,zi)(s,0){\displaystyle (s_{i},z_{i})\to (s,0)}), pero esta es una extensión continua, no analíticano se cumple para s real < 0 !);
  • en cada rama entera en s para z ≠ 0 fijo .

Valores especiales

  • Γ(s+1,1)=mis¡mi{\displaystyle \Gamma (s+1,1)={\frac {\lfloor es!\rfloor }{e}}}si s es un entero positivo ,
  • Γ(s,incógnita)=(s1)¡miincógnitak=0s1incógnitakk¡{\displaystyle \Gamma (s,x)=(s-1)!\,e^{-x}\sum _{k=0}^{s-1}{\frac {x^{k}}{k!}}}si s es un entero positivo , [ 12 ]
  • Γ(s,0)=Γ(s),(s)>0{\displaystyle \Gamma (s,0)=\Gamma (s),\Re (s)>0},
  • Γ(1,incógnita)=miincógnita{\displaystyle \Gamma (1,x)=e^{-x}},
  • γ(1,incógnita)=1miincógnita{\displaystyle \gamma (1,x)=1-e^{-x}},
  • Γ(0,incógnita)=Ei(incógnita){\displaystyle \Gamma (0,x)=-\operatorname {Ei} (-x)}paraincógnita>0{\displaystyle x>0},
  • Γ(s,incógnita)=incógnitasmi1s(incógnita){\displaystyle \Gamma (s,x)=x^{s}\operatorname {E} _{1-s}(x)},
  • Γ(12,incógnita)=πerfc(incógnita){\displaystyle \Gamma \left({\tfrac {1}{2}},x\right)={\sqrt {\pi }}\operatorname {erfc} \left({\sqrt {x}}\right)},
  • γ(12,incógnita)=πterreno(incógnita){\displaystyle \gamma \left({\tfrac {1}{2}},x\right)={\sqrt {\pi }}\operatorname {erf} \left({\sqrt {x}}\right)}.

Aquí,Ei{\displaystyle \operatorname {Ei} }es la integral exponencial ,minorte{\displaystyle \operatorname {E} _{n}}es la integral exponencial generalizada ,terreno{\displaystyle \operatorname {erf} }es la función de error yerfc{\displaystyle \operatorname {erfc} }es la función de error complementaria ,erfc(incógnita)=1terreno(incógnita){\displaystyle \operatorname {erfc} (x)=1-\operatorname {erf} (x)}.

Comportamiento asintótico

  • γ(s,incógnita)incógnitas1s{\displaystyle {\frac {\gamma (s,x)}{x^{s}}}\to {\frac {1}{s}}}comoincógnita0{\displaystyle x\to 0},
  • Γ(s,incógnita)incógnitas1s{\displaystyle {\frac {\Gamma (s,x)}{x^{s}}}\to -{\frac {1}{s}}}comoincógnita0{\displaystyle x\to 0}y(s)<0{\displaystyle \Re (s)<0}(para s real , el error de Γ( s , x ) ~ − x s / s es del orden de O ( x min{ s + 1, 0} ) si s ≠ −1 y O (ln( x )) si s = −1 ),
  • Γ(s,incógnita)Γ(s)norte=0(1)norteincógnitas+nortenorte¡(s+norte){\displaystyle \Gamma (s,x)\sim \Gamma (s)-\sum _{n=0}^{\infty }(-1)^{n}{\frac {x^{s+n}}{n!(s+n)}}}como una serie asintótica dondeincógnita0+{\displaystyle x\to 0^{+}}ys0,1,2,{\displaystyle s\neq 0,-1,-2,\dots }. [ 13 ]
  • Γ(norte,incógnita)donorte+(1)norte+1norte¡lnincógnitanorte=0,nortenorte(1)norteincógnitanortenortenorte¡(nortenorte){\displaystyle \Gamma (-N,x)\sim C_{N}+{\frac {(-1)^{N+1}}{N!}}\ln x-\sum _{n=0,n\neq N}^{\infty }(-1)^{n}{\frac {x^{n-N}}{n!(n-N)}}}como una serie asintótica dondeincógnita0+{\displaystyle x\to 0^{+}}ynorte=1,2,{\displaystyle N=1,2,\dots }, dóndedonorte=(1)norte+1norte¡(γnorte=1norte1norte){\textstyle C_{N}={\frac {(-1)^{N+1}}{N!}}\left(\gamma -\displaystyle \sum _{n=1}^{N}{\frac {1}{n}}\right)}, dóndeγ{\displaystyle \gamma }es la constante de Euler-Mascheroni . [ 13 ]
  • γ(s,incógnita)Γ(s){\displaystyle \gamma (s,x)\to \Gamma (s)}comoincógnita{\displaystyle x\to \infty }. [ 14 ]
  • Γ(s,incógnita)incógnitas1miincógnita1{\displaystyle {\frac {\Gamma (s,x)}{x^{s-1}e^{-x}}}\to 1}comoincógnita{\displaystyle x\to \infty }[ 15 ] ,
  • sis{\displaystyle s}≥ 0 es de valor real,Γ(s,z)zs1mizk=0Γ(s)Γ(sk)zk{\displaystyle \Gamma (s,z)\sim z^{s-1}e^{-z}\sum _{k=0}{\frac {\Gamma (s)}{\Gamma (s-k)}}z^{-k}}como una serie asintótica donde|z|{\displaystyle |z|\to \infty }y|argz|<32π{\displaystyle \left|\arg z\right|<{\tfrac {3}{2}}\pi }. [ 16 ]

Fórmulas de evaluación

La función gamma inferior se puede evaluar utilizando la expansión en serie de potencias: [ 17 ]γ(s,z)=k=0zsmizzks(s+1)(s+k)=zsmizk=0zksk+1¯{\displaystyle \gamma (s,z)=\sum _{k=0}^{\infty }{\frac {z^{s}e^{-z}z^{k}}{s(s+1)\dots (s+k)}}=z^{s}e^{-z}\sum _{k=0}^{\infty }{\dfrac {z^{k}}{s^{\overline {k+1}}}}} dóndesk+1¯{\displaystyle s^{\overline {k+1}}}es el símbolo de Pochhammer .

Una expansión alternativa es γ(s,z)=k=0(1)kk¡zs+ks+k=zssMETRO(s,s+1,z),{\displaystyle \gamma (s,z)=\sum _{k=0}^{\infty }{\frac {(-1)^{k}}{k!}}{\frac {z^{s+k}}{s+k}}={\frac {z^{s}}{s}}M(s,s+1,-z),} donde M es la función hipergeométrica confluente de Kummer .

Conexión con la función hipergeométrica confluente de Kummer

Cuando la parte real de z es positiva, γ(s,z)=s1zsmizMETRO(1,s+1,z){\displaystyle \gamma (s,z)=s^{-1}z^{s}e^{-z}M(1,s+1,z)} dóndeMETRO(1,s+1,z)=1+z(s+1)+z2(s+1)(s+2)+z3(s+1)(s+2)(s+3)+{\displaystyle M(1,s+1,z)=1+{\frac {z}{(s+1)}}+{\frac {z^{2}}{(s+1)(s+2)}}+{\frac {z^{3}}{(s+1)(s+2)(s+3)}}+\cdots }tiene un radio de convergencia infinito.

Nuevamente con funciones hipergeométricas confluentes y empleando la identidad de Kummer, Γ(s,z)=mizU(1s,1s,z)=zsmizΓ(1s)0mis(z+)d=mizzsU(1,1+s,z)=miz0mi(z+)s1d=mizzs0miz(1+)s1d.{\displaystyle {\begin{aligned}\Gamma (s,z)&=e^{-z}U(1{-}s,1{-}s,z)={\frac {z^{s}e^{-z}}{\Gamma (1-s)}}\int _{0}^{\infty }{\frac {e^{-u}}{u^{s}(z+u)}}du\\&=e^{-z}z^{s}U(1,1{+}s,z)=e^{-z}\int _{0}^{\infty }e^{-u}(z+u)^{s-1}du\\&=e^{-z}z^{s}\int _{0}^{\infty }e^{-zu}(1+u)^{s-1}du.\end{aligned}}}

Para el cálculo real de valores numéricos, la fracción continua de Gauss proporciona una expansión útil: γ(s,z)=zsmizsszs+1+zs+2(s+1)zs+3+2zs+4(s+2)zs+5+3zs+6.{\displaystyle \gamma (s,z)={\cfrac {z^{s}e^{-z}}{s-{\cfrac {sz}{s+1+{\cfrac {z}{s+2-{\cfrac {(s+1)z}{s+3+{\cfrac {2z}{s+4-{\cfrac {(s+2)z}{s+5+{\cfrac {3z}{s+6-\ddots }}}}}}}}}}}}}}.}

Esta fracción continua converge para todo z complejo , siempre que s no sea un entero negativo.

La función gamma superior tiene la fracción continua [ 18 ]Γ(s,z)=zsmizz+1s1+1z+2s1+2z+3s1+{\displaystyle \Gamma (s,z)={\cfrac {z^{s}e^{-z}}{z+{\cfrac {1-s}{1+{\cfrac {1}{z+{\cfrac {2-s}{1+{\cfrac {2}{z+{\cfrac {3-s}{1+\ddots }}}}}}}}}}}}} yΓ(s,z)=zsmiz1+zs+s13+zs+2(s2)5+zs+3(s3)7+zs+4(s4)9+zs+{\displaystyle \Gamma (s,z)={\cfrac {z^{s}e^{-z}}{1+z-s+{\cfrac {s-1}{3+z-s+{\cfrac {2(s-2)}{5+z-s+{\cfrac {3(s-3)}{7+z-s+{\cfrac {4(s-4)}{9+z-s+\ddots }}}}}}}}}}}

Teorema de la multiplicación

El siguiente teorema de multiplicación es verdadero : Γ(s,z)=1tsi=0(11t)ii¡Γ(s+i,tz)=Γ(s,tz)(tz)smitzi=1(1t1)iiLi1(si)(tz).{\displaystyle {\begin{aligned}\Gamma (s,z)&={\frac {1}{t^{s}}}\sum _{i=0}^{\infty }{\frac {\left(1-{\frac {1}{t}}\right)^{i}}{i!}}\Gamma (s+i,tz)\\&=\Gamma (s,tz)-(tz)^{s}e^{-tz}\sum _{i=1}^{\infty }{\frac {\left({\frac {1}{t}}-1\right)^{i}}{i}}L_{i-1}^{(s-i)}(tz).\end{aligned}}}

Implementación de software

Las funciones gamma incompletas están disponibles en varios sistemas de álgebra computacional .

Sin embargo, aunque no estén disponibles directamente, los valores incompletos de las funciones pueden calcularse utilizando funciones comunes en las hojas de cálculo (y en los paquetes de álgebra computacional). En Excel , por ejemplo, se pueden calcular utilizando la función gamma combinada con la función de distribución gamma .

  • La función incompleta inferior:γ(s,incógnita){\displaystyle \gamma (s,x)} = EXP(GAMMALN(s))*GAMMA.DIST(x,s,1,TRUE).
  • La función incompleta superior:Γ(s,incógnita){\displaystyle \Gamma (s,x)} = EXP(GAMMALN(s))*(1-GAMMA.DIST(x,s,1,TRUE)).

Estas se derivan de la definición de la función de distribución acumulativa de la distribución gamma .

En Python , la biblioteca Scipy proporciona implementaciones de funciones gamma incompletas scipy.special, sin embargo, no admite valores negativos para el primer argumento. La función gammaincde la biblioteca mpmath admite todos los argumentos complejos.

Funciones gamma regularizadas y variables aleatorias de Poisson

Dos funciones relacionadas son las funciones gamma regularizadas: PAG(s,incógnita)=γ(s,incógnita)Γ(s),Q(s,incógnita)=Γ(s,incógnita)Γ(s)=1PAG(s,incógnita).{\displaystyle {\begin{aligned}P(s,x)&={\frac {\gamma (s,x)}{\Gamma (s)}},\\[1ex]Q(s,x)&={\frac {\Gamma (s,x)}{\Gamma (s)}}=1-P(s,x).\end{aligned}}}PAG(s,incógnita){\displaystyle P(s,x)}es la función de distribución acumulativa para variables aleatorias gamma con parámetro de formas{\displaystyle s}y parámetro de escala 1.

Cuandos{\displaystyle s}es un número entero,Q(s+1,λ){\displaystyle Q(s+1,\lambda )}es la función de distribución acumulativa para variables aleatorias de Poisson : Siincógnita{\displaystyle X}es unPAGoi(λ){\displaystyle \mathrm {Poi} (\lambda )}variable aleatoria entonces Pr(incógnitas)=ismiλλii¡=Γ(s+1,λ)Γ(s+1)=Q(s+1,λ).{\displaystyle \Pr(X\leq s)=\sum _{i\leq s}e^{-\lambda }{\frac {\lambda ^{i}}{i!}}={\frac {\Gamma (s+1,\lambda )}{\Gamma (s+1)}}=Q(s+1,\lambda ).}

Esta fórmula se puede obtener mediante integración sucesiva por partes.

PAG(s,incógnita){\displaystyle P(s,x)}yQ(s,incógnita){\displaystyle Q(s,x)}se implementan como gammainc[ 19 ] y gammaincc[ 20 ] en scipy .

Derivadas parciales

Utilizando la representación integral anterior, la derivada de la función gamma incompleta superiorΓ(s,incógnita){\displaystyle \Gamma (s,x)}con respecto a x es Γ(s,incógnita)incógnita=incógnitas1miincógnita{\displaystyle {\frac {\partial \Gamma (s,x)}{\partial x}}=-x^{s-1}e^{-x}} La derivada con respecto a su primer argumentos{\displaystyle s}está dado por [ 21 ]Γ(s,incógnita)s=lnincógnitaΓ(s,incógnita)+incógnitaT(3,s,incógnita){\displaystyle {\frac {\partial \Gamma (s,x)}{\partial s}}=\ln x\Gamma (s,x)+x\,T(3,s,x)} y la segunda derivada por 2Γ(s,incógnita)s2=ln2incógnitaΓ(s,incógnita)+2incógnita[lnincógnitaT(3,s,incógnita)+T(4,s,incógnita)]{\displaystyle {\frac {\partial ^{2}\Gamma (s,x)}{\partial s^{2}}}=\ln ^{2}x\Gamma (s,x)+2x\left[\ln x\,T(3,s,x)+T(4,s,x)\right]} donde la funciónT(metro,s,incógnita){\displaystyle T(m,s,x)}es un caso especial de la función G de MeijerT(metro,s,incógnita)=GRAMOmetro1,metrometro,0(0,0,,0s1,1,,1|incógnita).{\displaystyle T(m,s,x)=G_{m-1,\,m}^{\,m,\,0}\!\left(\left.{\begin{matrix}0,0,\dots ,0\\s-1,-1,\dots ,-1\end{matrix}}\;\right|\,x\right).} Este caso especial en particular tiene propiedades de cierre interno propias porque puede usarse para expresar todas las derivadas sucesivas. En general, metroΓ(s,incógnita)smetro=lnmetroincógnitaΓ(s,incógnita)+metroincógnitanorte=0metro1PAGnortemetro1lnmetronorte1incógnitaT(3+norte,s,incógnita){\displaystyle {\frac {\partial ^{m}\Gamma (s,x)}{\partial s^{m}}}=\ln ^{m}x\Gamma (s,x)+mx\,\sum _{n=0}^{m-1}P_{n}^{m-1}\ln ^{m-n-1}x\,T(3+n,s,x)} dóndePAGjnorte{\displaystyle P_{j}^{n}}es la permutación definida por el símbolo de Pochhammer : PAGjnorte=(nortej)j¡=norte¡(nortej)¡.{\displaystyle P_{j}^{n}={\binom {n}{j}}j!={\frac {n!}{(n-j)!}}.} Todos estos derivados pueden generarse sucesivamente a partir de: T(metro,s,incógnita)s=lnincógnita T(metro,s,incógnita)+(metro1)T(metro+1,s,incógnita){\displaystyle {\frac {\partial T(m,s,x)}{\partial s}}=\ln x~T(m,s,x)+(m-1)T(m+1,s,x)} y T(metro,s,incógnita)incógnita=T(metro1,s,incógnita)+T(metro,s,incógnita)incógnita{\displaystyle {\frac {\partial T(m,s,x)}{\partial x}}=-{\frac {T(m-1,s,x)+T(m,s,x)}{x}}} Esta funciónT(metro,s,incógnita){\displaystyle T(m,s,x)}se puede calcular a partir de su representación en serie válida para|z|<1{\displaystyle |z|<1}, T(metro,s,z)=(1)metro1(metro2)¡dmetro2dtmetro2[Γ(st)zt1]|t=0+norte=0(1)nortezs1+nortenorte¡(snorte)metro1{\displaystyle T(m,s,z)=-{\frac {\left(-1\right)^{m-1}}{(m-2)!}}\left.{\frac {d^{m-2}}{dt^{m-2}}}\left[\Gamma (s-t)z^{t-1}\right]\right|_{t=0}+\sum _{n=0}^{\infty }{\frac {\left(-1\right)^{n}z^{s-1+n}}{n!\left(-s-n\right)^{m-1}}}} con el entendimiento de que s no es un entero negativo ni cero. En tal caso, se debe usar un límite. Resultados para|z|1{\displaystyle |z|\geq 1}se puede obtener mediante continuación analítica . Algunos casos especiales de esta función se pueden simplificar. Por ejemplo,T(2,s,incógnita)=Γ(s,incógnita)/incógnita{\displaystyle T(2,s,x)=\Gamma (s,x)/x},incógnitaT(3,1,incógnita)=mi1(incógnita){\displaystyle x\,T(3,1,x)=\mathrm {E} _{1}(x)}, dóndemi1(incógnita){\displaystyle \mathrm {E} _{1}(x)}es la integral exponencial . Estas derivadas y la funciónT(metro,s,incógnita){\displaystyle T(m,s,x)}proporcionar soluciones exactas a varias integrales mediante la diferenciación repetida de la definición integral de la función gamma incompleta superior. [ 22 ] [ 23 ] Por ejemplo, incógnitats1lnmetrotmitdt=metrosmetroincógnitats1mitdt=metrosmetroΓ(s,incógnita){\displaystyle \int _{x}^{\infty }{\frac {t^{s-1}\ln ^{m}t}{e^{t}}}dt={\frac {\partial ^{m}}{\partial s^{m}}}\int _{x}^{\infty }{\frac {t^{s-1}}{e^{t}}}dt={\frac {\partial ^{m}}{\partial s^{m}}}\Gamma (s,x)} Esta fórmula puede ampliarse o generalizarse a una amplia gama de transformadas de Laplace y Mellin . Al combinarse con un sistema de álgebra computacional , el uso de funciones especiales proporciona un método eficaz para resolver integrales definidas, en particular las que se presentan en aplicaciones prácticas de ingeniería (véase Integración simbólica para más detalles).

Integrales indefinidas y definidas

Las siguientes integrales indefinidas se obtienen fácilmente mediante la integración por partes (omitiendo la constante de integración en ambos casos): incógnitab1γ(s,incógnita)dincógnita=1b(incógnitabγ(s,incógnita)γ(s+b,incógnita)),incógnitab1Γ(s,incógnita)dincógnita=1b(incógnitabΓ(s,incógnita)Γ(s+b,incógnita)).{\displaystyle {\begin{aligned}\int x^{b-1}\gamma (s,x)\,dx&={\frac {1}{b}}\left(x^{b}\gamma (s,x)-\gamma (s+b,x)\right),\\[1ex]\int x^{b-1}\Gamma (s,x)\,dx&={\frac {1}{b}}\left(x^{b}\Gamma (s,x)-\Gamma (s+b,x)\right).\end{aligned}}} Las funciones gamma incompletas inferior y superior están conectadas mediante la transformada de Fourier : γ(s2,z2π)(z2π)s2mi2πikzdz=Γ(1s2,k2π)(k2π)1s2.{\displaystyle \int _{-\infty }^{\infty }{\frac {\gamma \left({\frac {s}{2}},z^{2}\pi \right)}{(z^{2}\pi )^{\frac {s}{2}}}}e^{-2\pi ikz}dz={\frac {\Gamma \left({\frac {1-s}{2}},k^{2}\pi \right)}{(k^{2}\pi )^{\frac {1-s}{2}}}}.} Esto se deduce, por ejemplo, de la especialización adecuada de ( Gradshteyn et al. 2015 , §7.642) .

Notas

  1. 1 2 3 4 5 6 "DLMF: §8.2 Definiciones y propiedades básicas ‣ Funciones gamma incompletas ‣ Capítulo 8 Gamma incompleta y funciones relacionadas" . dlmf.nist.gov .
  2. 1 2 3 "DLMF: §8.8 Relaciones de recurrencia y derivadas ‣ Funciones gamma incompletas ‣ Capítulo 8 Gamma incompleta y funciones relacionadas" . dlmf.nist.gov .
  3. 1 2 Donald E. Marshall (otoño de 2009). "Análisis complejo" (PDF) . Matemáticas 534 (material para estudiantes). Universidad de Washington. Teorema 3.9 en la pág. 56. Archivado del original (PDF) el 16 de mayo de 2011. Recuperado el 23 de abril de 2011 .
  4. 1 2 "DLMF: §8.7 Expansiones en serie ‣ Funciones gamma incompletas ‣ Capítulo 8 Gamma incompleta y funciones relacionadas" . dlmf.nist.gov .
  5. Paul Garrett. "Teorema de Hartogs: la analiticidad separada implica la conjunta" (PDF) . cse.umn.edu . Consultado el 21 de diciembre de 2023 .
  6. C. Teleman. "Superficies de Riemann" (PDF) . berkeley.edu . Consultado el 21 de diciembre de 2023 .
  7. "DLMF: §5.2 Definiciones ‣ Propiedades ‣ Capítulo 5 Función Gamma" . dlmf.nist.gov .
  8. "DLMF: §4.4 Valores especiales y límites ‣ Logaritmo, exponencial, potencias ‣ Capítulo 4 Funciones elementales" . dlmf.nist.gov .
  9. ver última ecuación.
  10. "DLMF: §8.4 Valores especiales ‣ Funciones gamma incompletas ‣ Capítulo 8 Gamma incompleta y funciones relacionadas" . dlmf.nist.gov .
  11. "DLMF: 8.4 Valores especiales" .
  12. Weisstein, Eric W. "Función gamma incompleta" . MathWorld .(ecuación 2)
  13. 1 2 Bender & Orszag (1978). Métodos matemáticos avanzados para científicos e ingenieros . Springer. Bibcode : 1978amms.book.....B .
  14. "DLMF: §8.11 Aproximaciones y expansiones asintóticas ‣ Funciones gamma incompletas ‣ Capítulo 8 Gamma incompleta y funciones relacionadas" . dlmf.nist.gov .
  15. Temme, NM (1979). "La expansión asintótica de las funciones gamma incompletas". SIAM Journal on Mathematical Analysis . 10 (4): 757– 766. doi : 10.1137/0510071 .
  16. "DLMF: §8.11 Aproximaciones y expansiones asintóticas ‣ Funciones gamma incompletas ‣ Capítulo 8 Gamma incompleta y funciones relacionadas" . dlmf.nist.gov .
  17. "DLMF: §8.11 Aproximaciones y expansiones asintóticas ‣ Funciones gamma incompletas ‣ Capítulo 8 Gamma incompleta y funciones relacionadas" . dlmf.nist.gov .
  18. Abramowitz y Stegun p.  263, 6.5.31
  19. "scipy.special.gammainc — Manual de SciPy v1.11.4" . docs.scipy.org .
  20. "scipy.special.gammaincc — Manual de SciPy v1.11.4" . docs.scipy.org .
  21. KO Geddes , ML Glasser, RA Moore y TC Scott, Evaluación de clases de integrales definidas que involucran funciones elementales mediante la diferenciación de funciones especiales , AAECC (Álgebra Aplicable en Ingeniería, Comunicación e Informática), vol. 1, (1990), págs. 149–165,
  22. Milgram, MS (1985). "La función integro-exponencial generalizada" . Math. Comp . 44 (170): 443– 458. doi : 10.1090/S0025-5718-1985-0777276-4 . MR 0777276 . 
  23. Mathar (2009). "Evaluación numérica de la integral oscilatoria sobre exp(i*pi*x)*x^(1/x) entre 1 e infinito". arXiv : 0912.3844 [ math.CA ]., Aplicación B

Referencias

  • Abramowitz, Milton ; Stegun, Irene Ann , eds. (1983) [junio de 1964]. «Capítulo 6.5». Manual de funciones matemáticas con fórmulas, gráficas y tablas matemáticas . Serie de Matemáticas Aplicadas. Vol.  55 (novena reimpresión con correcciones adicionales de la décima edición original con correcciones (diciembre de 1972); primera  ed.). Washington D. C.; Nueva York: Departamento de Comercio de los Estados Unidos, Oficina Nacional de Normas; Dover Publications. ISBN 978-0-486-61272-0. LCCN 64-60036 . MR 0167642 . LCCN 65-12253 .   "Función Gamma Incompleta" .§6.5.
  • Allasia, Giampietro; Besenghi, Renata (1986). "Cálculo numérico de funciones gamma incompletas mediante la regla del trapecio". Número. Matemáticas . 50 (4): 419– 428. doi : 10.1007/BF01396662 . S2CID 121964300 . 
  • Amore, Paolo (2005). "Representaciones asintóticas y exactas en serie para la función Gamma incompleta". Europhys. Lett . 71 (1): 1– 7. arXiv : math-ph/0501019 . Bibcode : 2005EL.....71....1A . doi : 10.1209 / epl/i2005-10066-6 . MR 2170316. S2CID 1921569 .  
  • G. Arfken y H. Weber. Métodos matemáticos para físicos . Harcourt/Academic Press, 2000. (Véase el capítulo 10).
  • DiDonato, Armido R.; Morris, Jr., Alfred H. (diciembre de 1986). "Cálculo de las razones de la función gamma incompleta y su inversa". ACM Transactions on Mathematical Software . 12 (4): 377– 393. doi : 10.1145/22721.23109 . S2CID 14351930 . 
  • Barakat, Richard (1961). "Evaluación de la función gamma incompleta de argumento imaginario mediante polinomios de Chebyshev" . Math. Comp . 15 (73): 7–11 . doi : 10.1090/s0025-5718-1961-0128058-1 . MR 0128058 . 
  • Carsky, Petr; Polasek, Martin (1998). "Funciones Gamma F_m(x) incompletas para argumentos reales y complejos". J. Comput. Phys . 143 (1): 259– 265. Bibcode : 1998JCoPh.143..259C . doi : 10.1006/jcph.1998.5975 . MR 1624704 . 
  • Chaudhry, M. Aslam; Zubair, SM (1995). "Sobre la descomposición de funciones Gamma incompletas generalizadas con aplicaciones a transformadas de Fourier" . J. Comput. Appl. Math . 59 (101): 253–284 . doi : 10.1016/0377-0427(94)00026-w . MR 1346414 . 
  • DiDonato, Armido R.; Morris, Jr., Alfred H. (septiembre de 1987). "ALGORITMO 654: Subrutinas FORTRAN para calcular las razones de la función gamma incompleta y su inversa" . ACM Transactions on Mathematical Software . 13 (3): 318–319 . doi : 10.1145/29380.214348 . S2CID 19902932 . (Véase también www.netlib.org/toms/654 ).
  • Früchtl, H.; Otto, P. (1994). "Un nuevo algoritmo para la evaluación de la función Gamma incompleta en computadoras vectoriales" . ACM Trans. Math. Softw . 20 (4): 436– 446. doi : 10.1145/198429.198432 . S2CID 16737306 . 
  • Gautschi, Walter (1998). "La función gamma incompleta desde Tricomi". Atti Convegni Lincei . 147 : 203– 237. SEÑOR 1737497 . 
  • Gautschi, Walter (1999). " Una nota sobre el cálculo recursivo de funciones gamma incompletas" . ACM Trans. Math. Softw . 25 (1): 101– 107. doi : 10.1145/305658.305717 . MR 1697463. S2CID 36469885 .  
  • Gradshteyn, Izrail Solomonovich ; Ryzhik, Iosif Moiseevich ; Geronimus, Yuri Veniaminovich ; Tseytlin, Michail Yulyevich ; Jeffrey, Alan (2015) [octubre de 2014]. "8.35". En Zwillinger, Daniel; Moll, Víctor Hugo (eds.). Tabla de Integrales, Series y Productos . Traducido por Scripta Technica, Inc. (8  ed.). Academic Press, Inc. págs. 908–911 . ISBN  978-0-12-384933-5. LCCN 2014010276 . 
  • Jones, William B.; Thron, WJ (1985). "Sobre el cálculo de funciones gamma incompletas en el dominio complejo" . J. Comput. Appl. Math . 12–13 : 401–417 . doi : 10.1016/0377-0427(85)90034-2 . MR 0793971 . 
  • "Función gamma incompleta" , Enciclopedia de Matemáticas , EMS Press , 2001 [1994]
  • Mathar, Richard J. (2004). "Representación numérica de la función gamma incompleta de argumento de valor complejo". Numerical Algorithms . 36 (3): 247– 264. arXiv : math/0306184 . Bibcode : 2004NuAlg..36..247M . doi : 10.1023/B:NUMA.0000040063.91709.58 . MR 2091195 . S2CID 30860614 .  
  • Miller, Allen R.; Moskowitz, Ira S. (1998). "Sobre ciertas funciones Gamma incompletas generalizadas" . J. Comput. Appl. Math . 91 (2): 179– 190. doi : 10.1016/s0377-0427(98)00031-4 .
  • Paris, RB (2010), "Función gamma incompleta" , 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 .
  • Paris, RB (2002). "Una expansión asintótica uniforme para la función gamma incompleta" . J. Comput. Appl. Math . 148 (2): 323– 339. Bibcode : 2002JCoAM.148..323P . doi : 10.1016/S0377-0427(02)00553-8 . MR 1936142 . 
  • Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007). «Sección 6.2. Función gamma incompleta y función de error» . Numerical Recipes: The Art of Scientific Computing (3.ª  ed.). Nueva York: Cambridge University Press. ISBN 978-0-521-88068-8Archivado del original el 15 de abril de 2021. Consultado el 9 de agosto de 2011 .
  • Takenaga, Roy (1966). "Sobre la evaluación de la función gamma incompleta" . Math. Comp . 20 (96): 606– 610. doi : 10.1090/S0025-5718-1966-0203911-3 . MR 0203911 . 
  • Temme, Nico (1975). " Expansiones asintóticas uniformes de las funciones gamma incompletas y la función beta incompleta" . Math . Comp . 29 (132): 1109–1114 . doi : 10.1090/S0025-5718-1975-0387674-2.MR 0387674 . 
  • Terras, Riho (1979). "La determinación de funciones gamma incompletas mediante integración analítica". J. Comput. Phys . 31 (1): 146– 151. Bibcode : 1979JCoPh..31..146T . doi : 10.1016/0021-9991(79)90066-4 . MR 0531128 . 
  • Tricomi, Francesco G. (1950). "Sulla funzione gamma incompleta". Ana. Estera. Pura Appl . 31 : 263–279.doi : 10.1007 / BF02428264 . SEÑOR 0047834 . S2CID 120404791 .  
  • Tricomi, FG (1950). "Asymptotische Eigenschaften der unvollst. Gammafunktion". Matemáticas. Z.53 (2): 136– 148. doi : 10.1007/bf01162409 . SEÑOR 0045253 . S2CID 121234109 .  
  • van Deun, Joris; Cools, Ronald (2006). " Una recurrencia estable para la función gamma incompleta con segundo argumento imaginario". Numer. Math . 104 (4): 445– 456. doi : 10.1007/s00211-006-0026-1 . MR 2249673. S2CID 43780150 .  
  • Winitzki, Serge (2003). "Cálculo de la función gamma incompleta con precisión arbitraria". En Vipin Kumar; Marina L. Gavrilova ; Chih Jeng Kenneth Tan; Pierre L'Ecuyer (eds.). Ciencia Computacional y sus Aplicaciones — ICSSA 2003. Conferencia Internacional sobre Ciencia Computacional y sus Aplicaciones, Montreal, Canadá, 18-21 de mayo de 2003, Actas, Parte I. Lecture Notes in Computer Science. Vol.  2667. pp. 790-798 . doi : 10.1007/3-540-44839-x_83 . ISBN  978-3-540-40155-1MR 2110953 . 
  • Weisstein, Eric W. "Función gamma incompleta" . MathWorld .
  • PAG(a,incógnita){\displaystyle P(a,x)}— Calculadora de la función gamma incompleta inferior regularizada
  • Q(a,incógnita){\displaystyle Q(a,x)}— Calculadora de la función gamma incompleta superior regularizada
  • γ(a,incógnita){\displaystyle \gamma (a,x)}— Calculadora de la función gamma incompleta inferior
  • Γ(a,incógnita){\displaystyle \Gamma (a,x)}— Calculadora de la función gamma incompleta superior
  • Fórmulas e identidades de la función gamma incompleta functions.wolfram.com