Este artículo utiliza notación matemática técnica para logaritmos. Todos los casos de log( x ) sin una base de subíndice deben interpretarse como un logaritmo natural , también comúnmente escrito como ln( x ) o log e ( x ) .
Derivada por Daniel Bernoulli , para números complejos con parte real positiva, la función gamma se define mediante una integral impropia convergente :
Existen otras extensiones de la función factorial, pero la función gamma es la más popular y útil. Es un componente de varias funciones de distribución de probabilidad y, como tal, es aplicable en los campos de la probabilidad y la estadística , así como en la combinatoria .
Motivación
La función gamma puede verse como una solución al problema de interpolación de encontrar una curva suave que conecte los puntos de la secuencia factorial: para todos los valores enteros positivos de . La fórmula simple para el factorial, x ! = 1 × 2 × ⋯ × x solo es válido cuando x es un entero positivo, y ninguna función elemental tiene esta propiedad, pero una buena solución es la función gamma . [1]
La función gamma no sólo es fluida sino también analítica (excepto en los números 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 números enteros positivos, como para un número entero . [1] Esta función se conoce como función pseudogamma , siendo la más famosa la función de Hadamard . [2]
Un requisito más restrictivo es la ecuación funcional que interpola el factorial desplazado : [3] [4]
Pero esto todavía no da una solución única, ya que permite la multiplicación por cualquier función periódica con y , como por ejemplo . Una forma de resolver la ambigüedad es el teorema de Bohr-Mollerup , que muestra que la única función de interpolación sobre los reales positivos es logarítmicamente convexa (superconvexa), [5] lo que significa que es convexa (donde está el logaritmo natural ). [6]
Definición
Definición principal
La notación se debe a Legendre . [1] Si la parte real del número complejo z es estrictamente positiva ( ), entonces la integral converge absolutamente , y se conoce como integral de Euler de segunda especie . (La integral de Euler de primer tipo es la función beta . [1] ) Usando la integración por partes , se ve que:
Reconociendo que como
Podemos calcular :
Así podemos demostrar que para cualquier entero positivo n por inducción . Específicamente, el caso base es ese y el paso de inducción es que
La identidad se puede usar (o, dando el mismo resultado, se puede usar la continuación analítica ) para extender de manera única la formulación integral de a una función meromórfica definida para todos los números complejos z , excepto los enteros menores o iguales a cero. [1] Es esta versión extendida la que comúnmente se conoce como función gamma. [1]
Definiciones alternativas
Hay muchas definiciones equivalentes.
La definición de Euler como producto infinito.
Para un número entero fijo , a medida que el número entero aumenta, tenemos que [7]
Si no es un número entero, entonces no es posible decir si esta ecuación es verdadera porque todavía no hemos definido (en esta sección) la función factorial para números no enteros. Sin embargo, obtenemos una extensión única de la función factorial a los números no enteros al insistir en que esta ecuación continúe siendo válida cuando el número entero arbitrario se reemplaza por un número complejo arbitrario .
Multiplicar ambos lados por da
este producto infinito , que se debe a Euler, [8] converge para todos los números complejos excepto los enteros no positivos, que fallan debido a una división por cero. Intuitivamente, esta fórmula indica que es aproximadamente el resultado de calcular un número entero grande , multiplicar por para aproximar y usar la relación hacia atrás para obtener una aproximación ; y además que esta aproximación se vuelve exacta a medida que aumenta hasta el infinito.
El producto infinito del recíproco
es una función entera , convergente para todo número complejo 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 excepto los enteros no positivos:
donde está la constante de Euler-Mascheroni . [1] Este es el producto de Hadamard en una forma reescrita. No se puede exagerar la utilidad de esta definición, ya que aparece en una determinada identidad que involucra a pi. [ cita necesaria ]
Una propiedad simple pero útil, que puede verse en la definición del límite, es:
En particular, con z = a + bi , este producto es
Si la parte real es un número entero o un semientero, esto se puede expresar de forma finita en forma cerrada :
Quizás el valor más conocido de la función gamma en un argumento no entero es
el que se puede encontrar estableciendo las fórmulas de reflexión o duplicación, usando la relación con la función beta que se proporciona a continuación con , o simplemente haciendo la sustitución en el Definición integral de la función gamma, lo que da como resultado una integral gaussiana . En general, para valores enteros no negativos de tenemos:
donde el doble factorial . Consulte Valores particulares de la función gamma para conocer los valores calculados.
Podría resultar tentador generalizar el resultado de que al buscar una fórmula para otros valores individuales donde sea 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, no se sabe que estos números sean expresables por sí mismos en términos de funciones elementales. Se ha demostrado que es un número trascendental y algebraicamente independiente de cualquier número entero y de cada una de las fracciones . [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 , ψ (0) ( z ) :
Para un entero positivo m, la derivada de la función gamma se puede calcular de la siguiente manera:
Para la derivada enésima de la función gamma es:
(Esto se puede derivar diferenciando la forma integral de la función gamma con respecto a y usando la técnica de diferenciación bajo el signo integral ).
Cuando se restringe a los números reales positivos, la función gamma es una función estrictamente logarítmicamente convexa . Esta propiedad podrá expresarse de cualquiera de las tres formas equivalentes siguientes:
Para dos números reales positivos cualesquiera y , y para cualquiera ,
Para dos números reales positivos cualesquiera y , y >
Para cualquier número real positivo ,
La última de estas afirmaciones es, esencialmente por definición, la misma que la afirmación de que , donde está la función poligamma de orden 1. Para demostrar la convexidad logarítmica de la función gamma, basta observar que tiene una representación en serie que, por positivo real x , consta únicamente de términos positivos.
La convexidad logarítmica y la desigualdad de Jensen juntas implican, para cualquier número real positivo y ,
También existen límites en las proporciones de funciones gamma. La más conocida es la desigualdad de Gautschi , que dice que para cualquier número real positivo x y cualquier s ∈ (0, 1) ,
La fórmula de Stirling.
El comportamiento de para una variable real positiva creciente viene dado por la fórmula de Stirling
donde el símbolo significa convergencia asintótica: la relación de los dos lados converge a 1 en el límite . [1] Este crecimiento es más rápido que exponencial, , para cualquier valor fijo de .
Otro límite útil para aproximaciones asintóticas es:
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]
Residuos
El comportamiento de los no positivos es más complejo. La integral de Euler no converge para , pero la función que define en el semiplano complejo positivo tiene una continuación analítica única en el semiplano negativo. Una forma de encontrar esa continuación analítica es utilizar 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]
eligiendo los que sean positivos. El producto en el denominador es cero cuando es igual a cualquiera de los números enteros . Así, la función gamma debe estar indefinida en esos puntos para evitar la división por cero ; es una función meromórfica con polos simples en los números enteros no positivos. [1]
Para una función de una variable compleja , en un polo simple , el residuo de viene dado por:
Para el polo simple reescribimos la fórmula de recurrencia como:
El numerador en es
y el denominador
Entonces los residuos de la función gamma en esos puntos son: [13] La función gamma es distinta de cero en todas partes a lo largo de la línea real, aunque se acerca arbitrariamente a cero cuando z → −∞ . De hecho, no existe ningún número complejo para el cual y, por tanto, la función gamma recíproca es una función completa , con ceros en . [1]
Mínimos y máximos
En la línea real, la función gamma tiene un mínimo local en z min ≈ +1.46163 21449 68362 34126 [14] donde alcanza el valor Γ( z min ) ≈ +0.88560 31944 10888 70027 . [15] La función gamma se eleva a ambos lados de este mínimo. La solución a Γ( z − 0,5) = Γ( z + 0,5) es z = +1,5 y el valor común es Γ(1) = Γ(2) = +1 . La solución positiva a Γ( 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 . [dieciséis]
La función gamma debe alternar el signo entre sus polos en los números enteros no positivos porque el producto en la recurrencia directa contiene un número impar de factores negativos si el número de polos entre y es impar, y un número par si el número de polos es par. . [13] Los valores en los extremos locales de la función gamma a lo largo del eje real entre los números enteros no positivos son:
Hay muchas fórmulas, además de la integral de Euler de segundo tipo, que expresan la función gamma como una integral. Por ejemplo, cuando la parte real de z es positiva, [22]
y [23]
donde las tres integrales se derivan respectivamente de las sustituciones , [24] y [25] en la segunda integral de Euler. La última integral en particular deja clara la conexión entre la función gamma en argumentos semienteros y la integral gaussiana : si dejamos, obtenemos .
La primera fórmula integral de Binet para la función gamma establece que, cuando la parte real de z es positiva, entonces: [26]
La integral del lado derecho puede interpretarse como una transformada de Laplace . Eso es,
La segunda fórmula integral de Binet establece que, nuevamente, cuando la parte real de z es positiva, entonces: [27]
Sea C un contorno de Hankel , es decir, un camino que comienza y termina en el punto ∞ de la esfera de Riemann , cuyo vector unitario tangente converge a −1 al inicio del camino y a 1 al final, que tiene el devanado número 1 alrededor 0 , y que no cruza [0, ∞) . Fije una rama de tomando una rama cortada a lo largo de [0, ∞) y tomando como real cuando t está en el eje real negativo. Supongamos que z no es un número entero. Entonces la fórmula de Hankel para la función gamma es: [28]
donde se interpreta como . La fórmula de reflexión conduce a la expresión estrechamente relacionada que
vuelve a ser válida siempre que z no sea un número entero.
Representación de fracción continua
La función gamma también se puede representar mediante la suma de dos fracciones continuas : [29] [30]
donde .
Expansión en serie de Fourier
El logaritmo de la función gamma tiene la siguiente expansión en serie de Fourier
, que durante mucho tiempo se atribuyó a Ernst Kummer , quien la derivó en 1847. [31] [32] Sin embargo, Iaroslav Blagouchine descubrió que Carl Johan Malmsten derivó por primera vez esta serie en 1842. [33] [34]
Este último se puede derivar tomando el logaritmo en la fórmula de multiplicación anterior, que da una expresión para la suma de Riemann del integrando. Tomando el límite para da la fórmula.
función pi
Una notación alternativa que fue introducida originalmente por Gauss es la función -, que, en términos de la función gamma, es
tal que para todo número entero no negativo .
Usando la función pi, la fórmula de reflexión toma la forma
donde sinc es la función sinc normalizada , mientras que el teorema de la multiplicación toma la forma
A veces también encontramos
que es una función completa , definida para cada número complejo, al igual que la función gamma recíproca . Eso es completo implica que no tiene polos, por lo que no tiene ceros .
En la primera integral anterior, que define la función gamma, los límites de integración son fijos. Las funciones gamma incompletas superior e inferior son las funciones que se obtienen al permitir que varíe el límite de integración inferior o superior (respectivamente).
La función gamma está relacionada con la función beta mediante la fórmula
La función gamma también aparece en una relación importante con la función zeta de Riemann . También aparece en la siguiente fórmula: que es válida sólo para .El logaritmo de la función gamma satisface la siguiente fórmula debida a Lerch: donde es la función zeta de Hurwitz , es la función zeta de Riemann y el primo ( ′ ) denota 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
Valores particulares
Incluyendo hasta los primeros 20 dígitos después del punto decimal, algunos valores particulares de la función gamma son:
(Estos números se pueden encontrar en la OEIS . [35] [36] [37] [38] [39] [40] Los Los valores presentados aquí están truncados en lugar de redondeados). La función gamma de valores complejos no está definida para números enteros no positivos, pero en estos casos el valor se puede definir en la esfera de Riemann como ∞ . La función gamma recíproca está bien definida y es analítica en estos valores (y en todo el plano complejo ):
Función log-gamma
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 se le da el nombre lgammaen lngammaentornos de programación o gammalnen hojas de cálculo); esto crece mucho más lentamente y, para cálculos combinatorios, permite sumar y restar registros en lugar de multiplicar y dividir valores muy grandes. A menudo se define como [41]
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
se utiliza a menudo ya que permite determinar valores de función en una franja de ancho 1 en z de la franja vecina. En particular, comenzando con una buena aproximación para az con una parte real grande, se puede ir paso a paso hasta la z deseada . Siguiendo una indicación de Carl Friedrich Gauss , Rocktaeschel (1922) propuso para logΓ( z ) una aproximación para Re( z ) grande :
Esto se puede utilizar para aproximar con precisión logΓ( z ) para z con un Re( z ) más pequeño mediante (PEBöhmer, 1939)
Se puede obtener una aproximación más precisa utilizando más términos de las expansiones asintóticas de logΓ( z ) y Γ( z ) , que se basan en la aproximación de Stirling.
como | z | → ∞ en constante | arg( z ) | < π . (Ver secuencias A001163 y A001164 en OEIS ).
En una presentación más "natural":
como | z | → ∞ en constante | arg( z ) | < π . (Ver secuencias A046968 y A046969 en OEIS ).
Los coeficientes de los términos con k > 1 de z 1− k en la última expansión son simplemente
donde los B k son los números de Bernoulli .
La función gamma también tiene una Serie de Stirling (derivada por Charles Hermite en 1900) igual a [42]
La función gamma es la única función que satisface simultáneamente
,
para todos los números complejos excepto los enteros no positivos, y,
para entero n , para todos los números complejos . [1]
En cierto sentido, la función log-gamma es la forma más natural; aclara algunos atributos intrínsecos de la función. Un ejemplo sorprendente es la serie de Taylor de logΓ alrededor de 1:
donde ζ ( k ) denota la función zeta de Riemann en k .
Entonces, usando la siguiente propiedad:
podemos encontrar una representación integral para la función log-gamma:
o, estableciendo z = 1 para obtener una integral para γ , podemos reemplazar el término γ con su integral e incorporarlo a la fórmula anterior, Llegar:
También existen fórmulas especiales para el logaritmo de la función gamma para z racional . Por ejemplo, si y son números enteros con y entonces [43] Esta fórmula se utiliza a veces para cálculos numéricos, ya que el integrando disminuye muy rápidamente.
Cuando se sigue eso y esto también
es una consecuencia de la fórmula de Raabe . O. Espinosa y V. Moll derivaron una fórmula similar para la integral del cuadrado de : [48]
donde es .
DH Bailey y sus coautores [49] evaluaron
cuándo en términos de la función zeta de Tornheim-Witten y sus derivadas.
Además, también se sabe que [50]
Aproximaciones
Los valores complejos de la función gamma se pueden aproximar utilizando la aproximación de Stirling o la aproximación de Lanczos .
Esto es preciso en el sentido de que la relación entre la aproximación y el valor verdadero se acerca a 1 en el límite como | z | va al infinito.
La función gamma se puede calcular con una precisión fija aplicando la integración por partes a la integral de Euler. Para cualquier número positivo x, la función gamma se puede escribir
Cuando Re( z ) ∈ [1,2] y , el valor absoluto de la última integral es menor que . Al elegir un valor lo suficientemente grande , esta última expresión se puede hacer más pequeña que cualquier valor deseado . Por lo tanto, la función gamma se puede evaluar con precisión con la serie anterior.
EA Karatsuba construyó un algoritmo rápido para calcular la función gamma de Euler para cualquier argumento algebraico (incluido el racional). [51] [52] [53]
A diferencia de muchas otras funciones, como la Distribución Normal , no se encuentra fácilmente una implementación obvia, rápida y precisa que sea fácil de implementar para la Función Gamma . 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 se pueden encontrar fácilmente mediante una búsqueda en Internet, como la Biblioteca Wiley en línea. Estas tablas se pueden utilizar con interpolación lineal . Se puede obtener una mayor precisión con el uso de la interpolación cúbica a costa de una mayor sobrecarga computacional. Dado que las tablas generalmente se publican para valores de argumento entre 1 y 2, la propiedad se puede usar para traducir rápida y fácilmente todos los valores reales al rango , de modo que solo sea necesario usar valores tabulados entre 1 y 2. [55]
Si las tablas de interpolación no son deseables, entonces la aproximación de Lanczos mencionada anteriormente funciona bien para obtener una precisión de 1 a 2 dígitos para valores pequeños y comúnmente utilizados de z. Si la aproximación de Lanczos no es lo 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 es posible evitar algunas de ellas al mantenerse alejado de muchas temas matemáticos especializados. Por otro lado, la función gamma Γ( z ) es la más difícil de evitar." [56]
La razón principal de la utilidad de la función gamma en tales contextos es la prevalencia de expresiones del tipo que describen procesos que decaen exponencialmente en el tiempo o el espacio. En ocasiones, las integrales de tales expresiones se pueden resolver en términos de la función gamma cuando no existe una solución elemental. Por ejemplo, si f es una función potencia y g es una función lineal, un simple cambio de variables da la evaluación
El hecho de que la integración se realice a lo largo de toda la línea 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 de 0 y ∞ para describir la acumulación de un proceso finito, en cuyo caso la función gamma ordinaria ya no es una solución; la solución entonces se llama función gamma incompleta . (La función gamma ordinaria, obtenida integrando toda la línea 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
y sus integrales, como la función de error . Existen muchas interrelaciones entre estas funciones y la función gamma; en particular, el factor obtenido al evaluar es el "mismo" que se encuentra en el factor de normalización de la función de error y la distribución normal .
Las integrales que hemos comentado 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 se puede utilizar para calcular el "volumen" y el "área" de hiperesferas de n dimensiones .
Calcular 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 números enteros sucesivos se pueden escribir como alguna combinación de factoriales, siendo quizás el ejemplo más importante el del coeficiente binomial . Por ejemplo, para cualquier número complejo z y n , con | z | < 1 , podemos escribir
cuál se parece mucho al coeficiente binomial cuando n es un número entero no negativo,
El ejemplo de los coeficientes binomiales motiva por qué las propiedades de la función gamma cuando se extiende a números negativos son naturales. Un coeficiente binomial indica el número de formas de elegir k elementos de un conjunto de n elementos; Si k > n , por supuesto no hay formas. Si k > n , ( n − k )! es el factorial de un entero negativo y, por tanto, infinito si utilizamos la definición de factoriales de la función gamma: dividir por infinito da el valor esperado de 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. Si P y Q son polinomios mónicos de grado myn con raíces respectivas p 1 , …, p m y q 1 , …, q n , tenemos
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 sólo del grado de los polinomios, por lo que no importa si b − a es igual a 5 o 10 5 . Al tomar los límites apropiados, también se puede hacer que la ecuación se cumpla incluso cuando el producto de la izquierda contenga 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. Debido al teorema de factorización de Weierstrass , las funciones analíticas se pueden escribir como productos infinitos y, en ocasiones, se pueden representar como productos finitos o cocientes de la función gamma. Ya hemos visto un ejemplo sorprendente: la fórmula de reflexión esencialmente representa la función seno como el producto de dos funciones gamma. A partir de esta fórmula, la función exponencial, así como todas las funciones trigonométricas e hiperbólicas, se pueden expresar en términos de la función gamma.
Esto proporciona, entre otras cosas, 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 llamados "triviales" en la recta real. Borwein et al. llama a esta fórmula "uno de los hallazgos más bellos de las matemáticas". [59] Otro contendiente para ese título podría ser
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 ocurridos en las 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
El problema de extender el factorial a argumentos no enteros aparentemente fue 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 [60]
que está bien definida para valores reales de x distintos de los números enteros negativos.
Leonhard Euler dio más tarde dos definiciones diferentes: la primera no era su integral sino un producto infinito bien definido para todos los números complejos n distintos de los enteros negativos,
de lo que informó a Goldbach en una carta fechada el 13 de octubre de 1729. Le escribió a Goldbach. nuevamente el 8 de enero de 1730, para anunciar su descubrimiento de la representación integral
que es válida cuando la parte real del número complejo n es estrictamente mayor que −1 (es decir, ). Por el cambio de variables t = −ln s , esto se convierte en la familiar integral de Euler. Euler publicó sus resultados en el artículo "De progresivoibus trascendentibus seu quarum termini generales algebraice dari nequeunt" ("Sobre progresiones trascendentales, es decir, aquellas cuyos términos generales no pueden darse algebraicamente"), presentado a la Academia de San Petersburgo el 28 de noviembre de 1729. [ 61] 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 ideó lo que hoy se conoce como fórmula de Stirling . Aunque la fórmula de Stirling da una buena estimación de n ! , además, para números no enteros, no proporciona el valor exacto. Extensiones de su fórmula que corrigen el error fueron dadas por el propio Stirling y por Jacques Philippe Marie Binet .
Siglo XIX: Gauss, Weierstrass y Legendre
Carl Friedrich Gauss reescribió el producto de Euler
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 lo hizo Gauss por primera vez. [62] 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 , a partir de otra representación del producto,
donde γ es la constante de Euler-Mascheroni . Weierstrass originalmente escribió su producto como uno para 1/Γ , en cuyo caso se toman los ceros de la función en lugar de sus polos. Inspirado por este resultado, demostró lo que se conoce como teorema de factorización de Weierstrass : que cualquier función completa puede escribirse como un producto de sus ceros en el plano complejo; una generalización del teorema fundamental del álgebra .
El nombre función gamma y el símbolo Γ fueron introducidos 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 en mayúsculas, no existe un estándar aceptado sobre si el nombre de la función debe escribirse "función gamma" o "función gamma" (algunos autores simplemente escriben " Γ -función"). La notación alternativa de "función pi" Π( z ) = z ! debida a Gauss se encuentra a veces en la literatura más antigua, pero la notación de Legendre es dominante en las obras modernas.
Está justificado preguntar por qué distinguimos entre el "factorial ordinario" y la función gamma usando símbolos distintos, y particularmente por qué la función gamma debería normalizarse a Γ( n + 1) = n . en lugar de simplemente usar " Γ( n ) = n ! ". Considere que la notación para exponentes, x n , se ha generalizado desde números enteros hasta números complejos x z sin ningún cambio. La motivación de Legendre para la normalización no parece ser conocida, y ha sido criticada como engorrosa por algunos (el matemático del siglo XX Cornelius Lanczos , por ejemplo, la llamó "desprovista de cualquier racionalidad" y en su lugar usaría z ! ). [63] 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 aditivo e − x contra el carácter multiplicativo x z con respecto a la medida de Haar en el grupo de Lie R + . Así, esta normalización deja más claro que la función gamma es un análogo continuo de una suma de Gauss . [64]
Siglos XIX y XX: caracterización de la función gamma
Es algo problemático que se hayan dado tantas definiciones para la función gamma. Aunque describen la misma función, no es del todo sencillo demostrar la equivalencia. Stirling nunca demostró que su fórmula extendida corresponda exactamente a la función gamma de Euler; Charles Hermite dio por primera vez una prueba en 1900. [65] En lugar de encontrar una prueba especializada para cada fórmula, sería deseable tener 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 tal ecuación no podría satisfacer la fórmula de recurrencia de la función gamma, convirtiéndola en una función trascendentalmente trascendental . Este resultado se conoce como teorema de Hölder .
Hasta 1922 no se dio una caracterización definitiva y generalmente aplicable de la función gamma. Harald Bohr y Johannes Mollerup demostraron entonces lo que se conoce como teorema de Bohr-Mollerup : que la función gamma es la única solución de la relación de recurrencia factorial que es positiva y logarítmicamente convexa para z positiva y cuyo valor en 1 es 1 (una función es logarítmicamente convexa si su logaritmo es convexo). Otra caracterización viene dada por el teorema de Wielandt .
El teorema de Bohr-Mollerup es útil porque es relativamente fácil demostrar la convexidad logarítmica de cualquiera de las diferentes fórmulas utilizadas para definir la función gamma. Yendo más allá, en lugar de definir la función gamma mediante una fórmula particular, podemos elegir las condiciones del teorema de Bohr-Mollerup como definición y luego elegir cualquier fórmula que nos guste que satisfaga las condiciones como punto de partida para estudiar la función gamma. . Este enfoque fue utilizado por el grupo Bourbaki .
Borwein y Corless [66] revisan tres siglos de trabajo sobre la función gamma.
Tablas de referencia y software.
Aunque la función gamma se puede calcular prácticamente tan fácilmente como cualquier función matemáticamente más simple con una computadora moderna (incluso con una calculadora de bolsillo programable), por supuesto, este no siempre fue así. Hasta mediados del siglo XX, los matemáticos dependían de tablas hechas a mano; en el caso de la función gamma, en particular una tabla calculada por Gauss en 1813 y otra calculada por Legendre en 1825. [67]
Las tablas de valores complejos de la función gamma, así como gráficos dibujados a mano, se dieron en Tablas de funciones con fórmulas y curvas de Jahnke y Emde [de] , publicado por primera vez en Alemania en 1909. Según Michael Berry , "la publicación en J&E de un gráfico tridimensional que muestra los polos de la función gamma en el plano complejo adquirió un estatus casi icónico." [68]
De hecho, había poca necesidad práctica de algo más que valores reales de la función gamma hasta la década de 1930, cuando se descubrieron aplicaciones de la función gamma compleja en la física teórica. A medida que las computadoras electrónicas estuvieron disponibles para la producción de tablas en la década de 1950, se publicaron varias tablas extensas para la función gamma compleja para satisfacer la demanda, incluida una tabla con una precisión de 12 decimales de la Oficina Nacional de Estándares de EE. UU . [1]
Las implementaciones de punto flotante de doble precisión de la función gamma y su logaritmo ahora están disponibles en la mayoría de los software 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 informática , como Mathematica y Maple . PARI/GP , MPFR y MPFUN contienen implementaciones gratuitas de precisión arbitraria. En algunas calculadoras de software , por ejemplo, la Calculadora de Windows y la Calculadora de GNOME , la función factorial devuelve Γ( x + 1) cuando la entrada x es un valor no entero. [69] [70]
^ abcdefghijklmn Davis, PJ (1959). "La integral de Leonhard Euler: un perfil histórico de la función gamma". Mensual Matemático Estadounidense . 66 (10): 849–869. doi :10.2307/2309786. JSTOR 2309786. Archivado desde el original el 7 de noviembre de 2012 . Consultado el 3 de diciembre de 2016 .
^ "¿Está mal definida la función Gamma? O: Hadamard versus Euler: ¿quién encontró la mejor función Gamma?".
^ Beals, Richard; Wong, Roderick (2010). Funciones especiales: un texto de posgrado. Prensa de la Universidad de Cambridge. pag. 28.ISBN978-1-139-49043-6.Extracto de la página 28
^ Ross, arcilla C. (2013). Ecuaciones diferenciales: una introducción a Mathematica (edición ilustrada). Medios de ciencia y negocios de Springer. pag. 293.ISBN978-1-4757-3949-7.Expresión G.2 en la página 293
^ Kingman, JFC (1961). "Una propiedad de convexidad de matrices positivas". La Revista Trimestral de Matemáticas . 12 (1): 283–284. Código Bib : 1961QJMat..12..283K. doi :10.1093/qmath/12.1.283.
^ Waldschmidt, M. (2006). "Trascendencia de los períodos: el estado del arte" (PDF) . Pura aplicación. Matemáticas. Cuarto . 2 (2): 435–463. doi : 10.4310/pamq.2006.v2.n2.a3 . Archivado (PDF) desde el original el 6 de mayo de 2006.
^ "¿Cómo obtener la expansión de Laurent de la función gamma alrededor de $z=0$?". Intercambio de pilas de matemáticas . Consultado el 17 de agosto de 2022 .
^ Artín, Emil (2015). La función gamma . Dover. pag. 24.
^ Gradshteyn, ES; Ryzhik, IM (2007). Tabla de integrales, series y productos (Séptima ed.). Prensa académica. pag. 893.ISBN978-0-12-373637-6.
^ Whittaker y Watson, 12.2 ejemplo 1.
^ Detlef, Gronau. "¿Por qué la función gamma es así?" (PDF) . Imsc.uni-graz.at .
^ Pascal Sebah, Xavier Gourdon. «Introducción a la Función Gamma» (PDF) . Computación de números . Archivado desde el original (PDF) el 30 de enero de 2023 . Consultado el 30 de enero de 2023 .
^ Whittaker y Watson, 12.31.
^ Whittaker y Watson, 12.32.
^ Whittaker y Watson, 22.12.
^ "Integral exponencial E: representaciones de fracciones continuas (Fórmula 06.34.10.0005)".
^ "Integral exponencial E: representaciones de fracciones continuas (Fórmula 06.34.10.0003)".
^ Srivastava, HM; Choi, J. (2001). Series Asociadas al Zeta y Funciones Afines . Países Bajos: Kluwer Academic.
^ Blagouchine, Iaroslav V. (2014). "Redescubrimiento de las integrales de Malmsten, su evaluación mediante métodos de integración de contornos y algunos resultados relacionados". Ramanujan J. 35 (1): 21-110. doi :10.1007/s11139-013-9528-5. S2CID 120943474.
^ Blagouchine, Iaroslav V. (2016). "Erratum y Anexo a "Redescubrimiento de las integrales de Malmsten, su evaluación mediante métodos de integración de contornos y algunos resultados relacionados"". Ramanujan J. 42 ( 3): 777–781. doi :10.1007/s11139-015-9763-z. S2CID 125198685.
^ "Función de registro gamma". Wolfram MathWorld . Consultado el 3 de enero de 2019 .
^ "Integral de Leonhard Euler: un perfil histórico de la función gamma" (PDF) . Archivado (PDF) desde el original el 12 de septiembre de 2014 . Consultado el 11 de abril de 2022 .
^ Blagouchine, Iaroslav V. (2015). "Un teorema para la evaluación en forma cerrada de la primera constante de Stieltjes generalizada en argumentos racionales y algunas sumas relacionadas". Revista de teoría de números . 148 : 537–592. arXiv : 1401.3724 . doi :10.1016/j.jnt.2014.08.009.
^ 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 Weidmanncshe Buchhandluns . 46 : 268–275.
^ Barnes, EW (1899). "La teoría de la función G ". Cuarto de galón. J. Matemáticas . 31 : 264–314.
^ Adamchik, Víctor S. (1998). "Funciones poligamma de orden negativo". J. Computación. Aplica. Matemáticas . 100 (2): 191–199. doi : 10.1016/S0377-0427(98)00192-7 .
^ Gosper, RW (1997). " en funciones especiales, q -series y temas afines". Mermelada. Matemáticas. Soc . 14 .
^ Espinosa, Olivier; Moll, Víctor H. (2002). "Sobre algunas integrales que involucran la función Hurwitz Zeta: Parte 1". El diario Ramanujan . 6 (2): 159–188. doi :10.1023/A:1015706300169. S2CID 128246166.
^ Bailey, David H.; Borwein, David; Borwein, Jonathan M. (2015). "Sobre las integrales log-gamma de Euler y las funciones zeta de Tornheim-Witten". El diario Ramanujan . 36 (1–2): 43–68. doi :10.1007/s11139-012-9427-1. S2CID 7335291.
^ Amdeberhan, T.; Coffey, Mark W.; Espinosa, Olivier; Koutschan, Christoph; Maná, Dante V.; Moll, Víctor H. (2011). "Integrales de potencias de logamma". Proc. América. Matemáticas. Soc . 139 (2): 535–545. doi : 10.1090/S0002-9939-2010-10589-0 .
^ EA Karatsuba, Evaluación rápida de funciones trascendentales. Problema. inf. Trans. Vol.27, No.4, págs. 339–360 (1991).
^ EA Karatsuba, Sobre un nuevo método para la evaluación rápida de funciones trascendentales. Ruso. Matemáticas. Sobrevivir. Vol.46, n.º 2, págs. 246–247 (1991).
^ EA Karatsuba "Algoritmos rápidos y el método FEE".
^ Borwein, JM; Zucker, IJ (1992). "Evaluación rápida de la función gamma para pequeñas fracciones racionales utilizando integrales elípticas completas de primer tipo". Revista IMA de Análisis Numérico . 12 (4): 519–526. doi :10.1093/IMANUM/12.4.519.
^ Werner, Helmut; Collinge, Robert (1961). "Aproximaciones de Chebyshev a la función gamma". Matemáticas. Computación . 15 (74): 195-197. JSTOR 2004230.
^ Michon, GP "Trigonometría y funciones básicas Archivado el 9 de enero de 2010 en Wayback Machine ". Numérica . Consultado el 5 de mayo de 2007.
^ Chaudry, MA y Zubair, SM (2001). Sobre una clase de funciones gamma incompletas con aplicaciones . pag. 37
^ Arroz, JA (1995). Estadística Matemática y Análisis de Datos (Segunda Edición). pag. 52–53
^ Borwein, J.; Bailey, DH y Girgensohn, R. (2003). Experimentación en Matemáticas . AK Peters. pag. 133.ISBN978-1-56881-136-9.
^ "Interpolando el factorial natural n! o El nacimiento de la función factorial real (1729 - 1826)".
^ 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.
^ Remmert, R. (2006). Temas clásicos de la teoría de funciones complejas . Traducido por Kay, LD Springer. ISBN978-0-387-98221-2.
^ Lanczos, C. (1964). "Una aproximación de precisión de la función gamma". J. SIAM Número. Anal. Ser. B . 1 (1): 86. Código bibliográfico : 1964SJNA....1...86L. doi :10.1137/0701008.
^ Ilker Inam; Engin Büyükaşşk (2019). Notas de la Escuela Internacional de Otoño sobre Teoría Computacional de Números. Saltador. pag. 205.ISBN978-3-030-12558-5.Extracto de la página 205
^ Knuth, DE (1997). El arte de la programación informática, volumen 1 (algoritmos fundamentales) . Addison-Wesley.
^ Borwein, Jonathan M .; Corless, Robert M. (2017). "Gamma y Factorial en la Mensual". Mensual Matemático Estadounidense . 125 (5). Asociación Matemática de América: 400–24. arXiv : 1703.05349 . Código Bib : 2017arXiv170305349B. doi :10.1080/00029890.2018.1420983. S2CID 119324101.
^ "¿Cuál es la historia de Gamma_function?". añois.com . Consultado el 5 de noviembre de 2022 .
^ Berry, M. (abril de 2001). "¿Por qué las funciones especiales son especiales?". Física hoy .
^ "microsoft/calculadora". GitHub . Consultado el 25 de diciembre de 2020 .
^ "calculadora de gnomos". GNOME.org . Consultado el 3 de marzo de 2023 .
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.
Artín, 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: Sociedad Matemática Estadounidense.
Birkhoff, George D. (1913). "Nota sobre la función gamma". Toro. América. Matemáticas. Soc . 20 (1): 1–10. doi : 10.1090/s0002-9904-1913-02429-7 . SEÑOR 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". Mensual Matemático Estadounidense . 66 (10): 849–869. doi :10.2307/2309786. JSTOR 2309786.
Correo, Emil (1919). "Las funciones gamma generalizadas". Anales de Matemáticas . Segunda Serie. 20 (3): 202–217. doi :10.2307/1967871. JSTOR 1967871 . Consultado el 5 de marzo de 2021 .
Prensa, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007). "Sección 6.1. Función Gamma". Recetas numéricas: el arte de la informática científica (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 .
Temmé, Nico M. (1996). Funciones especiales: introducción a las funciones clásicas de la física matemática . Nueva York: John Wiley & Sons. ISBN 978-0-471-11313-3.