La formule de Chudnovsky
Découverte en 1988 par les frères David et Gregory Chudnovsky, elle
fournit 14 décimales exactes par terme. C'est encore
aujourd'hui l'une des formules les plus utilisées pour calculer π sur
un ordinateur personnel.
1/π = 12 · Σk≥0 (−1)k (6k)! (13591409 + 545140134·k)
/ [ (3k)! (k!)³ · 6403203k+3/2 ]
Le script utilise la variante dite binary splitting : les
termes sont regroupés par paires, ce qui transforme une addition
naïve en une multiplication d'entiers efficaces. Résultat : environ
45 secondes en Python pur pour un million de
décimales, sans aucune bibliothèque externe.
La vérification
Un calcul de π n'a de valeur que s'il est contrôlé. Le script
procède à un test systématique avant d'écrire quoi que ce soit :
- il compare les 100 premières décimales à une référence connue et
s'arrête avec une erreur en cas d'écart ;
- le module JavaScript expose sa taille réelle
(PI_DATA.count), ce qui permet de
vérifier que l'intégralité du million est bien chargée ;
- les positions annoncées dans la page
Anecdotes ont toutes été
mesurées sur ce même million de décimales.
Limite connue : cette méthode de contrôle est empirique, pas une
preuve formelle. Chaque terme de la série majore l'erreur, et la
précision est fixée à plusieurs centaines de chiffres de marge —
largement suffisant pour un usage pédagogique.