La fórmula de Stirling es una buena aproximación del factorial bien conocida por, al menos, los lectores más antiguos de este blog (seguro que por mucha más gente). Dicha fórmula dice que
entendiendo como «equivalente». De hecho, cuando mayor es
más equivalentes son los valores de esas dos expresiones. Esto significa, grosso modo, que el límite del cociente de esas dos expresiones cuando
tiende a infinito es 1. Vamos, que en términos de límites las dos son «iguales» (no es exactamente así, pero nos vale para que quede más o menos clara la relación entre ellas).
Como hemos dicho, esta fórmula de Stirling da aproximaciones buenas del factorial en el sentido descrito anteriormente, pero según qué precisión necesitemos las aproximaciones obtenidas quizás no sean tan buenas. Por ejemplo, para :
Con esta aproximación estaríamos cometiendo un error relativo del 2%:
que no está mal, pero quizás en ocasiones sea muy grande. Por ejemplo, para que el error relativo baje del 1% debemos irnos hasta , y el primer número entero positivo para el cual la aproximación por la fórmula de Stirling de su factorial nos da un error relativo menor del 0.1% es
.
La cuestión es que hay varias expresiones que dan mejores aproximaciones que la de Stirling. Una de estas mejoras, debida a Srinivasa Ramanujan, es la que nos ocupa hoy. Aquí la tenéis:
Nadie duda ahora que es de Ramanujan, ¿verdad?
Con esta fórmula, para obtenemos el valor
, mucho más aproximado al valor real, que es 24. El error relativo en este caso es 0.000282%, muchísimo más pequeño que el que da la fórmula de Stirling. Y por poner otro ejemplo, para el caso
el error relativo de la fórmula de Stirling era un pelín mayor del 0.1%. Con la fórmula de Ramanujan es 0.000000001987%. La diferencia es sustancial, como se puede ver.
La fórmula de Stirling tiene la característica de que el error relativo va disminuyendo conforme aumenta el valor de . Pues el caso es que la fórmula de Ramanujan también posee esa interesante propiedad: cuanto mayor es
menor es el error relativo de la aproximación de Ramanujan. Al ser éste mucho más bajo que el que da la fórmula de Stirling, es bastante evidente que la de Ramanujan es mucho mejor aproximación.
¿Hay demostración de todo esto? Pues sí, pero no por el propio Ramanujan, que lo dejó sin demostrar, sino por Ekatherina Karatsuba, que lo probó en el año 2000 en su trabajo On the asymptotic representation of the Euler gamma function by Ramanujan. Lo que demostró Karatsuba fue que la función
con es monótona creciente y está acotada «en infinito» por
. Como ya sabemos que
para valores enteros positivos de
(vamos, que la función gamma es una generalización del factorial), este resultado demuestra que la aproximación de Ramanujan comentada anteriormente es correcta.
Y para terminar comentar que ni mucho menos estas dos aproximaciones para el factorial, la de Stirling y la de Ramanujan, son las únicas interesantes. Os dejo por aquí tres más que he encontrado junto la aproximación que dan para , para que podáis comparar con ellas:
- Fórmula de Burnside:
Para
su valor es 24.2226179.
- Fórmula de Gosper:
Para
su valor es 23.9908895.
- Y una tremenda de Necdet Batir que aparece en el último enlace al final de este post:
Su valor para
es 24.00000001332. Buenísima, ¿verdad? Por ser ésta tan buena vamos a dar un par de valores más. Para
tenemos que
, y nuestra fórmula da como resultado
. Y para
tenemos que
, y nuestra fórmula nos da
. Podéis comparar estos valores con los que dan las otras fórmulas para convenceros de que ésta es realmente impresionante.
Y ahora os toca a vosotros. ¿Conocéis más aproximaciones buenas del factorial tipo las que hemos descrito por aquí? Los comentarios son vuestros.
Fuentes y enlaces relacionados:
- On the asymptotic representation of the Euler gamma function by Ramanujan, de Ekatherina Karatsuba.
- The Ramanujan asymptotic formula for the Euler Gamma function, de Ekatherina Karatsuba (básicamente igual que el anterior, con algún detalle más).
- Ramanujan’s Factorial Approximation, en The Endeavour.
- On Gosper’s formula for the Gamma function (pdf).
- A Refinement of Ramanujan’s Factorial Approximation.
- Very accurate approximation for the factorial function (pdf).
¿Te ha gustado la entrada? Puedes invitarme a un café, Gauss te lo agradecerá 😉
se supone qeu la formula de EULER MACLAURIN aplicada a la suma de logaritmos predice una formula asintotica exacta para el factorial n!
Información Bitacoras.com…
Valora en Bitacoras.com: La fórmula de Stirling es una buena aproximación del factorial bien conocida por, al menos, los lectores más antiguos de este blog (seguro que por mucha más gente). Dicha fórmula dice que entendiendo como âequivalenteâ…..
Aplicando la fórmula de Euler McLaurin se obtiene una expresión asintótica, pero no exacta, pues la serie resulta divergente.
Creo que es mejor decir que la función Gamma interpola al factorial (es decir su gráfica pasa por los puntos de la gráfica de n!) en lugar de generaliza al factorial. ¿Qué sería, por ejemplo, el factorial de cualquier real x con 3< x < 4? La definición de factorial necesita de “sucesores” que se dan en la relación de orden usual de los reales pero sólo en el subconjunto Z y no en los racionales y a fortiori no en los reales (lo que se puede hacer es definir para cada real x una función F “pseudo- factorial” por… Lee más »
Entiendo que es una generalización del factorial en el sentido de que restringida a los naturales coincide con el factorial.
Sobre otras generalizaciones del factorial quizás hable en los próximos tiempos.
¡Dios!, ¿pero de dónde sacó Ramanujan esa equivalencia?…y su manía de dejar las cosas sin demostrar…¿¡por qué!?, jajaja…
Disculpen, me parecio tan maravilloso que lo introduje en mi calculadora, y me dio algo mas de error.
Por las dudas lo meti en el wolfram alpha, y coincidio con la calculadora.
http://www.wolframalpha.com/input/?i=n^n+e^%28-n%29+%282*pi*%28n%2B1%2F4%2B1%2F%2872n%29-31%2F%286480n^2%29-139%2F%28155520n^3%29%2B9871%2F%286531840n^4%29%29%29^%281%2F2%29%2Cn%3D5
¿Podrian decirme donde copie mal? (Si es que paso… sino corrijan el post por favor xD)
Edit: Revise el enlace del fin del post, y en vez de 1/4 dice 1/6, ahora la calculadora me arroja el resultado adecuado 😀
Como dice Luis Cardenas, en la fórmula de Necdet Batir el término 1/4 debería ser 1/6 según este enlace: http://files.ele-math.com/articles/jmi-04-30.pdf (que pones al final del artículo)
Saludos
Cierto, es 1/6. En cuanto pueda lo cambio. Gracias por el aviso :-).
Este es mi aporte al calculo de factoriales, en especial a un caso específico de coeficiente binomial. Ver en https://docs.google.com/file/d/0B0aVaPjfCtc3ZktZNUQ0TkVXSnM/edit