Métodos para calcular - mat.usal.esmat.usal.es/pi/pdf/espiral-pi.pdf · p 2 2 p 2+ p 2 2 q 2+ p 2+...

1
Métodos para calcular π Series infinitas MADHAVA ( C . 1400), L EIBNIZ (1674): π 4 = 1 - 1 3 + 1 5 - 1 7 + ··· Palabras clave: arco tangente. NEWTON (1665): π = 3 3 4 + 2 - 3 4 X n=0 ( 2n n ) (n + 1)(2n + 5)16 n Palabras clave: áreas de sectores circulares. E ULER (1736) : π 2 6 = X n=1 1 n 2 , π 4 90 = X n=1 1 n 4 , π 6 945 = X n=1 1 n 6 , π 8 9450 = X n=1 1 n 8 ,... Palabras clave: función zeta de Riemann. R AMANUJAN (1914) : 1 π = X n=0 2n n 3 42n + 5 2 12n+4 = 2 2 9801 X n=0 (4n)! n! 4 26390n + 1103 396 4n , etc. Palabras clave: formas modulares. D. Y G. C HUDNOVSKY (1988) : 1 π = 12 640320 3 X n=0 (6n)! (3n)! n! 3 (-1) n (545140134n + 13591409) 640320 3n Palabras clave: Ramanujan. B AILEY ,B ORWEIN Y P LOUFFE (1997) : π = X n=0 1 16 n 4 8n + 1 - 2 8n + 4 - 1 8n + 5 - 1 8n + 6 Palabras clave: bloques independientes de dígitos. GUILLERA (2002,2003) : 1 π 2 = X n=0 2n n 5 (-1) n (820n 2 + 180n + 13) 2 20n+7 ? = 1 128 5 X n=0 (-1) n 2880 3n (6n)! n! 6 (5418n 2 + 693n + 29), etc. Palabras clave: Ramanujan, método WZ, algoritmo PSLQ, relaciones conjeturadas. Productos infinitos VIÈTE (1593): 2 π = 2 2 · p 2 + 2 2 · q 2 + p 2 + 2 2 ··· WALLIS (1655): 2 π = 1 · 3 2 · 2 · 3 · 5 4 · 4 · 5 · 7 6 · 6 · 7 · 9 8 · 8 ··· E ULER (1737): 1 π = v u u t 1 6 Y p primo 1 - 1 p 2 = 4 v u u t 1 90 Y p primo 1 - 1 p 4 = 6 v u u t 1 945 Y p primo 1 - 1 p 6 = ··· Fracciones continuas B ROUNCKER (1665): 4 π = 1 + 1 2 2 + 3 2 2 + 5 2 2 + 7 2 2 + . . . L ANGE (1999): π = 3 + 1 2 6 + 3 2 6 + 5 2 6 + 7 2 6 + . . . P ICKETT Y C OLEMAN (2008): π 2 = 1 + 1 1 + 1 1 2 + 1 1 3 + 1 1 4 + . . . Para más información ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ ππ Integrales ÁREA DE UN DISCO: π = 2 Z 1 -1 p 1 - x 2 dx DISTRIBUCIÓN DE C AUCHY: π = Z -dx x 2 + 1 S ENO INTEGRAL : π = Z -sen x x dx DISTRIBUCIÓN NORMAL : π = Z -e -x 2 dx «Un matemático es alguien para quien esto es tan obvio como que dos veces dos son cuatro» — Lord Kelvin. Métodos empíricos L A AGUJA DE B UFFON: La aguja de Buffon Una aguja de longitud L se lanza sobre una serie de líneas paralelas espaciadas por una distancia λ, con L<λ. La proporción p de veces que la aguja cruza una de las líneas satisface π 2L Por ejemplo, si L = λ 2 se obtiene π 1/p. E L PÉNDULO SIMPLE : La relación entre el periodo de oscilación de un péndu- lo y su longitud es, para amplitudes pequeñas, π T 2 r g L (T = periodo, L = longitud, g = gravedad) Puntos enteros en un disco P UNTOS DE UN RETÍCULO DENTRO DE UN DISCO: Sea D R el disco cerrado con centro el origen y radio R y N(R) el número de puntos con coordenadas enteras dentro de D R . En- tonces π = ım RN(R) R 2 Por ejemplo, N(10)= 317 y N(100000)= 31415925457. S UMAS DE DOS CUADRADOS : Sea r(n) el número de representaciones de un natural n como suma de dos cuadrados enteros, contando signos y permuta- ciones como representaciones distintas. Por ejemplo, r(5)= 8 pues 5 =(±1) 2 +(±2) 2 =(±2) 2 +(±1) 2 . Entonces ım nr(1)+ ··· + r(n) n = π Es decir, un número tiene de media π representaciones como suma de dos cuadrados. Este resultado es la versión aritmética del problema de contar puntos. La distribución de los dígitos de π Se espera que un número irracional e incluso trascendente como π sea normal, lo cual, para sus dígitos decimales, significa que cada uno de ellos, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, se halla estadísticamente en la misma proporción: 10 %. En general, dada una base b, entre los primeros k dígitos del desarrollo de π en esa base, cada una de las b cifras c = 0, 1, . . . , b - 1 aparece un número N b (c, k) de veces. Se conjetura que para cualquier base debería ser ım kN b (c, k) k = 1 b . De hecho, cualquier sucesión finita σ de dígitos, con longitud fija , debería aparecer con frecuencia 1/b , correspon- diente a la equidistribución de dichas sucesiones. Borel demostró que casi todo número real, en el sentido de la medida de Lebesgue, cumple la propiedad de equi- distribución de sucesiones finitas respecto a toda base, pero, a día de hoy, ninguno de los enunciados anteriores ha sido demostrado para una constante concreta «natural» como π, e2 (para esta última hubo una supuesta demos- tración en 2005, pero no avalada por expertos), aunque sí existen ejemplos un tanto «artificiales» de tales números, como la constante de Champernowne 0,123456789101112 . . . . La normalidad de π a menudo se intenta divulgar de manera filosófica o poética: cualquier mensaje escrito, pasa- do o futuro, la letra de cualquier canción, internet al completo, la demostración de la conjetura de Riemann —o su refutación— etc., codificados como sucesiones binarias, por ejemplo, se hallarán (de ser cierta la normalidad) en algún lugar entre los dígitos binarios de π, incluso infinitas veces. Debe tenerse en cuenta, sin embargo, que esto incluye, como en la Biblioteca de Babel de Jorge Luis Borges, sucesiones sin sentido, libros que constan sólo de la letra «a», todo enunciado y su contrario, la paradoja del mentiroso, etc. Tampoco puede predecirse dónde entre los dígitos hay que buscar una sucesión dada y, aunque hallásemos alguna información, nada nos libraría de tener que decidir acerca de su valor, con lo cual es un pobre método de hallar verdades. Más prosaicamente, podemos afirmar que la sucesión 999999 ocurre por primera vez en la posición 762 (aparece en el mural en la vuelta número 9, abajo a la izquierda), pero hay que esperar a la posición 193034 para verla otra vez. Añadiendo otro 9, la sucesión 9999999 aparece por primera vez en la posición 1722776. Existen diversas páginas web dedicadas a hallar sucesiones finitas dadas entre los dígitos de π calculados hasta hoy en día.

Transcript of Métodos para calcular - mat.usal.esmat.usal.es/pi/pdf/espiral-pi.pdf · p 2 2 p 2+ p 2 2 q 2+ p 2+...

Page 1: Métodos para calcular - mat.usal.esmat.usal.es/pi/pdf/espiral-pi.pdf · p 2 2 p 2+ p 2 2 q 2+ p 2+ 2 2 WALLIS (1655): 2 ... EL PÉNDULO SIMPLE: La relación entre el periodo de oscilación

Métodos para calcular πSeries infinitas

MADHAVA (C. 1400), LEIBNIZ (1674):π

4= 1−

13+

15−

17+ · · ·

Palabras clave: arco tangente.

NEWTON (1665):

π =3√3

4+ 2−

34

∞∑n=0

(2nn

)(n+ 1)(2n+ 5)16n

Palabras clave: áreas de sectores circulares.

EULER (1736):π2

6=

∞∑n=1

1n2 ,

π4

90=

∞∑n=1

1n4 ,

π6

945=

∞∑n=1

1n6 ,

π8

9450=

∞∑n=1

1n8 , . . .

Palabras clave: función zeta de Riemann.

RAMANUJAN (1914):

1π=

∞∑n=0

(2nn

)3 42n+ 5212n+4 =

2√2

9801

∞∑n=0

(4n)!n!4

26390n+ 11033964n

, etc.

Palabras clave: formas modulares.

D. Y G. CHUDNOVSKY (1988):

1π=

12√6403203

∞∑n=0

(6n)!(3n)!n!3

(−1)n(545140134n+ 13591409)6403203n

Palabras clave: Ramanujan.

BAILEY, BORWEIN Y PLOUFFE (1997):

π =

∞∑n=0

116n

(4

8n+ 1−

28n+ 4

−1

8n+ 5−

18n+ 6

)Palabras clave: bloques independientes de dígitos.

GUILLERA (2002,2003):

1π2

=

∞∑n=0

(2nn

)5(−1)n(820n2 + 180n+ 13)

220n+7?=

1128√5

∞∑n=0

(−1)n

28803n(6n)!n!6

(5418n2 + 693n+ 29), etc.

Palabras clave: Ramanujan, método WZ, algoritmo PSLQ, relaciones conjeturadas.

Productos infinitos

VIÈTE (1593):

2π=

√22·√

2+√2

√2+

√2+√2

2· · ·

WALLIS (1655):2π=

1 · 32 · 2

· 3 · 54 · 4

· 5 · 76 · 6

· 7 · 98 · 8· · ·

EULER (1737):

1π=

√√√√16

∏p primo

(1−

1p2

)= 4

√√√√ 190

∏p primo

(1−

1p4

)= 6

√√√√ 1945

∏p primo

(1−

1p6

)= · · ·

Fracciones continuas

BROUNCKER (1665):4π= 1+

12

2+32

2+52

2+72

2+ . . .

LANGE (1999):

π = 3+12

6+32

6+52

6+72

6+ . . .

PICKETT Y COLEMAN (2008):π

2= 1+

1

1+1

12+

113+

114+ . . .

Para másinformación

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

π

Integrales

ÁREA DE UN DISCO:

π = 2∫ 1−1

√1− x2 dx

DISTRIBUCIÓN DE CAUCHY:π =

∫∞−∞

dx

x2 + 1

SENO INTEGRAL:π =

∫∞−∞

sen xx

dx

DISTRIBUCIÓN NORMAL:√π =

∫∞−∞ e

−x2dx

«Un matemático es alguien para quien esto es tan obvio como que dos veces dos son cuatro» — Lord Kelvin.

Métodos empíricos

LA AGUJA DE BUFFON:

La aguja de Buffon

Una aguja de longitud L se lanza sobre una serie delíneas paralelas espaciadas por una distancia λ, conL < λ. La proporción p de veces que la aguja cruzauna de las líneas satisface

π ≈ 2Lpλ

Por ejemplo, si L = λ2 se obtiene π ≈ 1/p.

EL PÉNDULO SIMPLE:

La relación entre el periodo de oscilación de un péndu-lo y su longitud es, para amplitudes pequeñas,

π ≈ T2

√g

L

(T = periodo, L = longitud, g = gravedad)

Puntos enteros en un disco

PUNTOS DE UN RETÍCULO DENTRO DE UN DISCO:

SeaDR el disco cerrado con centro el origen y radio R yN(R) elnúmero de puntos con coordenadas enteras dentro de DR. En-tonces

π = lımR→∞

N(R)

R2

Por ejemplo, N(10) = 317 y N(100000) = 31415925457.

SUMAS DE DOS CUADRADOS:

Sea r(n) el número de representaciones de un natural n comosuma de dos cuadrados enteros, contando signos y permuta-ciones como representaciones distintas. Por ejemplo, r(5) = 8pues 5 = (±1)2 + (±2)2 = (±2)2 + (±1)2. Entonces

lımn→∞ r(1) + · · ·+ r(n)n

= π

Es decir, un número tiene de media π representaciones comosuma de dos cuadrados. Este resultado es la versión aritmética

del problema de contar puntos.

La distribución de los dígitos de π

Se espera que un número irracional e incluso trascendente como π sea normal, lo cual, para sus dígitos decimales,significa que cada uno de ellos, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, se halla estadísticamente en la misma proporción: 10 %.

En general, dada una base b, entre los primeros k dígitos del desarrollo de π en esa base, cada una de las b cifrasc = 0, 1, . . . ,b− 1 aparece un número Nb(c,k) de veces. Se conjetura que para cualquier base debería ser

lımk→∞

Nb(c,k)k

=1b.

De hecho, cualquier sucesión finita σ de dígitos, con longitud fija `, debería aparecer con frecuencia 1/b`, correspon-diente a la equidistribución de dichas sucesiones.

Borel demostró que casi todo número real, en el sentido de la medida de Lebesgue, cumple la propiedad de equi-distribución de sucesiones finitas respecto a toda base, pero, a día de hoy, ninguno de los enunciados anteriores hasido demostrado para una constante concreta «natural» como π, e, ó

√2 (para esta última hubo una supuesta demos-

tración en 2005, pero no avalada por expertos), aunque sí existen ejemplos un tanto «artificiales» de tales números,como la constante de Champernowne 0,123456789101112 . . . .

La normalidad de π a menudo se intenta divulgar de manera filosófica o poética: cualquier mensaje escrito, pasa-do o futuro, la letra de cualquier canción, internet al completo, la demostración de la conjetura de Riemann —o surefutación— etc., codificados como sucesiones binarias, por ejemplo, se hallarán (de ser cierta la normalidad) en algúnlugar entre los dígitos binarios de π, incluso infinitas veces. Debe tenerse en cuenta, sin embargo, que esto incluye,como en la Biblioteca de Babel de Jorge Luis Borges, sucesiones sin sentido, libros que constan sólo de la letra «a»,todo enunciado y su contrario, la paradoja del mentiroso, etc. Tampoco puede predecirse dónde entre los dígitos hayque buscar una sucesión dada y, aunque hallásemos alguna información, nada nos libraría de tener que decidir acercade su valor, con lo cual es un pobre método de hallar verdades.

Más prosaicamente, podemos afirmar que la sucesión 999999 ocurre por primera vez en la posición 762 (aparece enel mural en la vuelta número 9, abajo a la izquierda), pero hay que esperar a la posición 193034 para verla otra vez.Añadiendo otro 9, la sucesión 9999999 aparece por primera vez en la posición 1722776. Existen diversas páginas webdedicadas a hallar sucesiones finitas dadas entre los dígitos de π calculados hasta hoy en día.