A veces se la llama serie Madhava-Leibniz, ya que fue descubierta por primera vez por el matemático indio Madhava de Sangamagrama o sus seguidores en los siglos XIV-XV (ver serie Madhava ), [1] y luego fue redescubierta de forma independiente por James Gregory en 1671 y Leibniz en 1673. [2] La serie de Taylor para la función tangente inversa , a menudo llamada serie de Gregory , es
Considerando sólo la integral en el último término, tenemos:
Por lo tanto, según el teorema de compresión , como n → ∞ , nos queda la serie de Leibniz:
Prueba 2
Sea , cuando , la serie converge uniformemente, entonces
Por lo tanto, si se aproxima de manera que sea continua y converja uniformemente, la prueba está completa, donde la serie converge según la prueba de Leibniz y también se aproxima desde dentro del ángulo de Stolz, por lo que, según el teorema de Abel, esto es correcto.
Convergencia
La fórmula de Leibniz converge extremadamente lentamente: presenta una convergencia sublineal . Calcular π con 10 decimales correctos usando la suma directa de la serie requiere precisamente cinco mil millones de términos porque 4/2k + 1 < 10 −10 para k > 2 × 10 10 − 1/2 (es necesario aplicar el límite de error calabrese ). Para obtener 4 decimales correctos (error de 0,00005) se necesitan 5000 términos. [4] Se encuentran disponibles límites de error incluso mejores que los de Calabrese o Johnsonbaugh. [5]
Sin embargo, la fórmula de Leibniz se puede utilizar para calcular π con alta precisión (cientos de dígitos o más) utilizando varias técnicas de aceleración de convergencia . Por ejemplo, la transformación de Shanks , la transformada de Euler o la transformación de Van Wijngaarden , que son métodos generales para series alternas, se pueden aplicar eficazmente a las sumas parciales de la serie de Leibniz. Además, la combinación de términos por pares da la serie no alterna
Si la serie se trunca en el momento adecuado, la expansión decimal de la aproximación coincidirá con la de π para muchos más dígitos, excepto para dígitos aislados o grupos de dígitos. Por ejemplo, tomando cinco millones de términos se obtiene
donde los dígitos subrayados son incorrectos. De hecho, los errores se pueden predecir; son generados por los números de Euler E n según la fórmula asintótica
donde N es un número entero divisible por 4. Si se elige que N sea una potencia de diez, cada término de la suma correcta se convierte en una fracción decimal finita. La fórmula es un caso especial de la fórmula de suma de Euler-Boole para series alternas, y proporciona otro ejemplo más de una técnica de aceleración de convergencia que se puede aplicar a la serie de Leibniz. En 1992, Jonathan Borwein y Mark Limber utilizaron los primeros mil números de Euler para calcular π con 5263 decimales con la fórmula de Leibniz. [6]
Producto Euler
La fórmula de Leibniz se puede interpretar como una serie de Dirichlet utilizando el carácter único no principal de Dirichlet módulo 4. Al igual que con otras series de Dirichlet, esto permite convertir la suma infinita en un producto infinito con un término para cada número primo . Un producto de este tipo se denomina producto de Euler . Es:
En este producto, cada término es una razón superparticular , cada numerador es un número primo impar y cada denominador es el múltiplo de 4 más cercano al numerador. [7]
^ Plofker, Kim (noviembre de 2012), " Tantrasaṅgraha de Nīlakaṇṭha Somayājī por K. Ramasubramanian y MS Sriram", The Mathematical Intelligencer , 35 (1): 86–88, doi :10.1007/s00283-012-9344-6, S2CID 124507583
^ Roy, Ranjan (1990). "El descubrimiento de la fórmula en serie para π por Leibniz, Gregory y Nilakantha" (PDF) . Revista Matemáticas . 63 (5): 291–306. doi :10.1080/0025570X.1990.11977541.Horvath, Miklos (1983). «Sobre la cuadratura leibniziana del círculo» (PDF) . Annales Universitatis Scientiarum Budapestiensis (Sectio Computatorica) . 4 : 75–83.
^ Villarino, Mark B. (21 de abril de 2018). "El error en una serie alterna". El Mensual Matemático Estadounidense . 125 (4): 360–364. doi :10.1080/00029890.2017.1416875. hdl : 10669/75532 . ISSN 0002-9890. S2CID 56124579.
^ Rattaggi, Diego (30 de agosto de 2018). "Estimaciones de error para la serie de Gregory-Leibniz y la serie armónica alterna utilizando integrales de Dalzell". arXiv : 1809.00998 [matemáticas.CA].
^ Borwein, Jonathan ; Bailey, David ; Girgensohn, Roland (2004), "1.8.1: Serie de Gregory reexaminada", Experimentación en matemáticas: caminos computacionales hacia el descubrimiento , AK Peters, págs. 28–30, ISBN1-56881-136-5, señor 2051473
^ Debnath, Lokenath (2010), El legado de Leonhard Euler: un tributo al tricentenario, World Scientific, p. 214, ISBN9781848165267.
Enlaces externos
Fórmula de Leibniz en C, ensamblaje x86 FPU, ensamblaje x86-64 SSE3 y ensamblaje DEC Alpha