Y llego el día. En la fecha de hoy, 14 de marzo, se ha celebrado desde hace ya unos años el Día de Pi, y desde 2019 se celebra el Día Internacional de las Matemáticas.
Llamadme nostálgico si queréis, pero, para mí, este día será siempre el Día de Pi, la constante matemática que más me ha fascinado desde siempre. Y, como suelo hacer en este día tan especial, os traigo una curiosidad sobre ella que me pasaron hace ya tiempo por Twitter (por desgracia, no recuerdo quién fue).
Si preguntas a alguien sobre si recuerda el valor de Pi, muchos te responderan , y otros muchos te dirán
. Teniendo en cuenta que Pi es irracional (por lo que tiene infinitos decimales que no siguen ningún patrón periódico de repetición), es evidente que ambas son aproximaciones («redondeos») del valor de Pi (por cierto, aquí otra demostración de la irracionalidad de Pi). Ahora, ¿son buenas o malas aproximaciones? Preguntado de otra forma: ¿cuántos decimales de Pi necesitamos en nuestra aproximación para que podamos considerar que nuestros cálculos son «prácticamente correctos» (en el sentido de que el error sea «despreciable»)?
Se sabe que solamente hacen falta 10 decimales de Pi para calcular la circunferencia de la Tierra con un error de un milímetro. Ahora, hay situaciones en las que hay que realizar cálculos con distancias mucho mayores que las que determinan el tamaño de nuestro planeta: los cálculos en Astronomía.
Por ello, creo que tiene interés saber cúantos decimales utilizan en la NASA para sus cálculos. Y va a ser Marc Rayman, director e ingeniero jefe de la misión Dawn de la propia NASA, quien nos lo cuente.
Digo que va a ser él porque es él mismo quien responde a una pregunta sobre ello en la web del Laboratorio de Propulsión a Chorro (JPL) de la NASA. Hace unos años, les llego la siguiente cuestión a través de Facebook:
¿JPL utiliza solamente 3,14 para sus cálculos con Pi? ¿O se usan más decimales, digamos algo como 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117067982148086513282306647093844609550582231725359408128481117450284102701938521105559644622948954930381964428810975665933446128475648233786783165271201909145648566923460348610454326648213393607260249141273724587006606315588174881520920962829254091715364367892590360?
La respuesta de Marc Rayman fue la siguiente: en la NASA usamos 15 decimales de Pi. Es decir, usan la siguiente aproximación de Pi:
El propio Marc justifica esta decisión de la siguiente forma:
Para apoyar su respuesta, Marc aporta los siguientes ejemplos:
Vamos a hacer nosotros los cálculos con las unidades que más familiares nos son. Se tiene que 12500 millones de millas es, aproximadamente, 20116,8 millones de kilómetros. Si queremos calcular la longitud de una circunferencia que tenga ese radio (recuerdo que la fórmula sería ), obtenemos los siguientes resultados:
- Usando el valor real de Pi:
- Usando la aproximación de Pi con 15 decimales:
Como podéis ver, la primera diferencia se produce en el cuarto decimal (de un 3 a un 2). Mirando el siguiente decimal, vemos que la diferencia sería de menos de un centímetro. Considerando que el valor inicial del radio de nuestra circunferencia es una aproximación, podríamos decir que si tenemos un círculo cuyo radio es de unos 20116,8 millones de kilómetros, usando una aproximación de Pi con 15 decimales podemos calcular la longitud de su circunferencia con un error de unos pocos centímetros.
Segundo ejemplo de Marc:
Hagamos de nuevo los cálculos con kilómetros. Tenemos que 7926 millas son 12755,66 kilómetros. Como ése sería el diámetro de la circunferencia, se tiene que la longitud de la misma se calcula con la fórmula . Vamos con ello:
- Usando el valor real de Pi:
- Usando la aproximación de Pi de 15 decimales:
Como se ve, la primera diferencia viene en el undécimo decimal. Esto es: la diferencia es bastante menor que un nanómetro de unos 10 nanómetros (gracias, Ángel). Para intentar hacernos una idea de cuánto es esto, os dejo este dato: el diámetro de un átomo de helio mide 0,1 nanómetros el diámetro de un glóbulo rojo está, aproximadamente, entre 6 y 8 nanómetros. Vamos, nah.
Marc nos deja un comentario más al respecto:
Todo esto ejemplifica bastante bien el porqué de que no necesitemos más de esos 15 decimales para obtener buenos resultados en los cálculos con nuestra aproximación de Pi.
Curioso este asunto de la aproximación de Pi que usa la NASA. Tengo que reconocer que yo pensaba que usaban una aproximación con mayor cantidad de decimales. ¿Y vosotros? Contádnoslo en los comentarios.
Por cierto, según cuenta Marc Rayman, no es la primera vez que llega esa pregunta a la NASA. Hace unos cuantos años, la propuso un estudiante de ciencias de sexto grado (no sé a qué nivel corresponde en España) y entusiasta del espacio. Más adelante, ese estudiante tuvo la oportunidad de hacer un doctorado en Física y de involucrarse en la exploración espacial. ¿Su nombre? Marc Rayman.
Fuente: How many decimals of Pi do we really need?.
Y para finalizar, os dejo un buen puñado de enlaces a artículos de Gaussianos relacionados con Pi que, creo, merecen mucho la pena:
- Construyendo Pi con regla y compás de forma exacta (bueno, casi…)
- Calculando Pi con volúmenes de hiperfsferas
- Cinco cosas que posiblemente no sepas sobre Pi (y varias que seguro que ya sabías)
- Un límite inocente y una cuestión poco conocida (y sin demostrar) del número Pi
- Pi y la fórmula de sumación de Poisson
- Echegaray y la trascendencia de Pi: no lo cuento, lo hago
- Pi no siempre vale 3,14159…
- ¿Quién fue el primero que probó que “la constante del círculo” (Pi) es constante?
- El primer producto infinito con Pi como protagonista
- Un día de Pi muy especial
- (Vídeo) Singing Pi-gram
- La cuestión más importante que aún no se ha respondido sobre el número Pi
- El algoritmo de Chudnovsky, o cómo se calculan los decimales de Pi en el siglo XXi
- El desarrollo más bello de Pi como suma infinita
- Celebrando el día de Pi con matemáticos nacidos ese mismo día
- Pi y el conjunto de Mandelbrot
- Celebrando infinitamente el día de Pi
- Probabilidad de escoger dos números coprimos
- Celebrando el día de Pi con una aguja y una medusa
- Cómo demostrar que Pi es irracional (II)
- Cadaeic Cadenza
- El problema de Basilea
- La identidad de Euler
¿Te ha gustado la entrada? Puedes invitarme a un café, Gauss te lo agradecerá 😉
El problema real es: cuántas operaciones encadenadas se realizan? El error aumenta con la acumulación de operaciones posteriores al redondeo.
Interesante cuestión. Habría que preguntarle a Marc sobre ello, a ver qué nos responde.
Pido disculpas por escribir para corregir un dato, pero creo que hay que hacerlo.
Ha calculado el valor de la longitud de la circunferencia de la Tierra con dos valores de pi, y concluye que la diferencia es del orden de 0,1 nm.
Yo no he hecho el calculo, por lo que me fío de sus resultados, y veo que la diferencia es de 10 nm.
Ángel, tienes toda la razón, me confundí en mi observación sobre los resultados.
Y nada de pedir disculpas por ello, era totalmente necesario hacer esa corrección. Te lo agradezco, y te animo a volver a hacerlo cuando lo creas oportuno.
Lo de usar 15 dígitos para PI es porque las variables utilizadas por los softwares son de doble precisión (15 cifras significativas), una longitud más que suficiente para la mayoría de los cálculos físicos.
Esto dicho de manera muy general porque el error final depende del volumen de cálculos y los algoritmos.