De grandes avancées supplémentaires n'ont été réalisées qu'à partir du XVe siècle (Al-Kashi). Les premiers mathématiciens modernes ont atteint une précision de 35 décimales au début du XVIIe siècle (Ludolph van Ceulen) et 126 chiffres au XIXe siècle (Jurij Vega), dépassant la précision requise pour toute application concevable en dehors des mathématiques pures.
Le record de l'approximation manuelle de π est détenu par William Shanks, qui a calculé 527 décimales correctes vers 1873. Depuis le milieu du XXe siècle, l'approximation de π est effectuée sur ordinateurs par des logiciels spécifiques.
Le , le record est établi avec cent mille milliards de décimales par Emma Haruka Iwao, travaillant sur Google Cloud durant 157 jours[1].
Histoire
Antiquité
Les approximations les plus connues de π datant avant l'ère commune étaient exactes à deux décimales ; cela a été amélioré avec les mathématiques chinoises en particulier au milieu du premier millénaire, à une précision de sept décimales. Après cela, aucun progrès n'a été réalisé jusqu'à la fin de la période médiévale.[pas clair]
Les mathématiques babyloniennes évaluaient généralement π à 3, ce qui suffisait pour les projets architecturaux de l'époque (notamment dans la description du temple de Salomon dans la Bible hébraïque). Les Babyloniens savaient qu'il s'agissait d'une approximation, et une vieille tablette mathématique babylonienne exhumée près de Suse en 1936 (datée entre le XIXe et le XVIIe siècle avant notre ère) apporte une meilleure approximation de π, 25⁄8 = 3,125, soit proche de 0,5% de la valeur exacte[2],[3],[4],[5].
À peu près au même moment, le papyrus Rhind (daté de la Deuxième Période intermédiaire égyptienne, soit 1600 avant notre ère) donne l'approximation 256⁄81 ≈ 3,16 (précis à 0,6 %) en calculant l'aire d'un octogone régulier approximant un cercle[6].
Les calculs astronomiques dans le Shatapatha Brahmana (vers le VIe siècle avant notre ère) utilisent l'approximation fractionnaire 339/108 ≈ 3,139[7].
Au IIIe siècle avant notre ère, Archimède a prouvé les inégalités strictes 223⁄71 < π < 22⁄7, au moyen d'un 96-gone (précision de 2 × 10−4 et 4 × 10−4, respectivement). Il est aussi possible de prouver que 22 / 7 dépasse π grâce à un calcul d'intégrales élémentaire[8], mais cette technique n'existait pas à l'époque d'Archimède.
Au IIe siècle, Ptolémée a utilisé la valeur 377⁄120, la première approximation connue précise à trois décimales[9].
Le mathématicien chinois Liu Hui, en 263, a calculé que π serait compris entre 3,141024 et 3,142708, valeurs correspondant respectivement aux constructions d'un 96-gone et d'un 192-gone ; la moyenne de ces deux valeurs est 3,141866 (précision de 9 × 10−5). Il a également suggéré que 3,14 était une approximation suffisante pour des applications pratiques. Il a souvent été crédité d'un résultat plus précis π ≈ 3927/1250 = 3,1416 (précision de 2 × 10−6), bien que certains chercheurs pensent plutôt que celui-ci est dû au mathématicien chinois Zu Chongzhi (429-500)[10]. Zu Chongzhi est connu pour avoir donné l'encadrement de π entre 3,1415926 et 3,1415927, correct à sept décimales près. Il a donné deux autres approximations de π : π ≈ 22/7 et π ≈ 355/113. La dernière fraction est la meilleure approximation rationnelle possible de π en utilisant moins de cinq chiffres décimaux au numérateur et au dénominateur. Le résultat de Zu Chongzhi restera l'approximation la plus juste pendant près d'un millénaire.
À l'époque de l'empire Gupta (VIe siècle), le mathématicien Aryabhata, dans son traité astronomique Āryabhaṭīya, a calculé la valeur de π à cinq chiffres significatifs (π ≈ 62832/20000 = 3,1416)[11]. Il l'a utilisée pour calculer une approximation de la circonférence de la Terre[12]. Aryabhata a affirmé que son résultat donnait « approximativement » (āsanna) la circonférence d'un cercle. Nilakantha Somayaji (école du Kerala) a fait remarquer, au XVe siècle, que le mot signifie non seulement que c'est une approximation, mais que la valeur est « incommensurable » (c.-à-d. irrationnelle)[13].
Au Ve siècle, π était connu à environ sept décimales dans les mathématiques chinoises, et à environ cinq dans les mathématiques indiennes. Les progrès suivants n'ont été réalisés que plus d'un millénaire après, au XIVe siècle, lorsque le mathématicien et astronome indien Madhava de Sangamagrama, fondateur de l'école du Kerala d'astronomie et de mathématiques, a découvert une série infinie pour π, maintenant connue sous le nom de formule de Madhava-Leibniz[14],[15], et a donné deux méthodes pour calculer la valeur de π. Ce faisant, il a obtenu la série infinie
et a utilisé les 21 premiers termes pour calculer les 11 premières décimales de π.
L'autre méthode qu'il utilisait était d'ajouter un terme restant à la série originale de π. Il a utilisé le terme restant
dans l'extension de série infinie de π⁄4 pour améliorer l'approximation de π à 13 décimales, lorsque n = 75.
Il a atteint ce niveau de précision en calculant le périmètre d'un polygone régulier à 3 × 228 côtés[17].
Du XVIe au XIXe siècle
Dans la seconde moitié du XVIe siècle, le mathématicien français François Viète a découvert un produit infini qui converge vers π, connu sous le nom de formule de Viète.
Le mathématicien germano-néerlandais Ludolph van Ceulen (vers 1600) a calculé les 35 premières décimales de π à l'aide d'un 262-gone. Il était si fier de ce résultat qu'il l'a fait inscrire sur sa pierre tombale.
Dans Cyclometricus (1621), Willebrord Snell affirme que le périmètre du polygone inscrit approche la circonférence deux fois plus vite que le périmètre du polygone circonscrit correspondant. La preuve sera apportée par Christiaan Huygens en 1654. Snell a obtenu 7 décimales de π avec un 96-gone(en)[18].
En 1789, le mathématicien slovèneJurij Vega a calculé les 140 premières décimales de π, cependant les 14 dernières étaient incorrectes[19] ; ce résultat a néanmoins été le record du monde pendant 52 ans, jusqu'en 1841, lorsque William Rutherford(en) a calculé 208 décimales dont les 152 premières étaient correctes. Vega a ensuite amélioré la formule de John Machin à partir de 1706 et sa méthode est encore mentionnée aujourd'hui.
Le mathématicien amateur anglais William Shanks a passé plus de 20 ans à calculer les 707 premières décimales de π. Lorsqu'il eut terminé sa tâche en 1873, les 527 premières décimales étaient correctes. Il calculait de nouvelles décimales toute la matinée et passait l'après-midi à vérifier le travail de la matinée. Il s'agissait de l'approximation la plus précise de π jusqu'à l'avènement de l'ordinateur, un siècle plus tard.
XXe siècle
En 1910, le mathématicien indien Srinivasa Ramanujan a trouvé plusieurs séries infinies de π à convergence rapide, comme
qui calcule à chaque nouvelle itération huit décimales supplémentaires de π. Ses séries sont maintenant la base des algorithmes les plus rapides actuellement utilisés pour calculer π.
À partir du milieu du XXe siècle, tous les calculs de π ont été effectués à l'aide de calculatrices ou d'ordinateurs.
En 1944, D. F. Ferguson, à l'aide d'une calculatrice mécanique, a constaté que William Shanks avait commis une erreur à la 528e décimale et que toutes les décimales suivantes étaient incorrectes.
Dans les premières années de l'ordinateur, un développement de π à 100 000 décimales[20] a été calculé par le mathématicien de l'université du MarylandDaniel Shanks (n'ayant aucun rapport avec William Shanks mentionné ci-dessus) et son équipe au Naval Research Laboratory à Washington. En 1961, Shanks et son équipe ont utilisé deux séries différentes pour calculer les décimales de π. Pour l'une, on savait que toute erreur produirait une valeur légèrement trop élevée, et pour l'autre, on savait que toute erreur produirait une valeur légèrement trop faible par rapport à la valeur réelle ; par conséquent, tant que les deux séries produisaient les mêmes chiffres, il n'y avait pratiquement aucun doute sur la justesse des décimales calculées. Les premiers 100 265 chiffres après la virgule de π ont été publiés en 1962[20].
En 1989, les frères Chudnovsky ont calculé correctement π à 4,8 × 108 décimales grâce au superordinateurIBM 3090 en utilisant la variation suivante de la série infinie de Ramanujan :
.
En 1999, Yasumasa Kanada et son équipe de l'université de Tokyo ont calculé π à 2,061 584 3 × 1011 décimales, grâce au superordinateur Hitachi SR8000/MPP en utilisant une autre variante de la série infinie de Ramanujan.
XXIe siècle
Certains records depuis 2009 ont été accomplis par des calculs sur des ordinateurs personnels en utilisant l'algorithme de Chudnovski
Le 6 décembre 2002, Yasumasa Kanada et une équipe de 9 autres personnes, utilisant le Hitachi SR8000, un supercalculateur avec 1 téraoctet de mémoire, ont calculé π à environ 1,241 1 × 1012 décimales en environ 600 heures[21],[22].
Le 29 avril 2009, Daisuke Takahashi et al., de l'université de Tsukuba, avec un supercalculateur T2K, ont plus que doublé le record précédent, en calculant π à plus de 2,57 × 1012 décimales, en moins de 74 heures au total (calcul et vérification)[23].
Le 31 décembre 2009, Fabrice Bellard a utilisé un ordinateur personnel pour calculer un peu moins de 2,7 × 1012 décimales de π. Les calculs ont été effectués en binaire, puis le résultat a été converti en base 10. Les étapes de calcul, de conversion et de vérification ont pris au total 131 jours[24].
De 2010 à 2020, tous les records ont été obtenus à l'aide du logiciel « y-cruncher » d'Alexander Yee :
En août 2010, Shigeru Kondo a calculé 5 × 1012 décimales de π[25]. Le calcul a été effectué entre le 4 mai et le 3 août, les vérifications primaires et secondaires ayant respectivement duré 64 et 66 heures[26].
Le 16 octobre 2011, Shigeru Kondo a doublé son propre record en calculant 1013 décimales en utilisant la même méthode, mais assisté d'un meilleur matériel[27],[28],[29].
En décembre 2013, Kondo a battu son propre record pour la seconde fois lorsqu'il a calculé 1,21 × 1013 décimales de π. La limitation était principalement due à l'espace de stockage[30].
Le 8 octobre 2014, un individu non identifié, connu sous le pseudonyme de « houkouonchi », a calculé 1,33 × 1013 décimales de π[31].
Le 15 novembre 2016, Peter Trueb et ses sponsors ont calculé, toujours grâce au programme y-cruncher, 2,245 915 771 836 1 × 1013 (πe × 1012) décimales de π. Le calcul a pris (avec 3 interruptions) 105 jours[31],[32].
Le 14 mars 2019, Google, en la personne d'Emma Haruka Iwao, rend public[33] le nouveau record qui s'établit à 31 415 926 535 897 (la partie entière de π e13) décimales de π.
Le 29 janvier 2020 un nouveau record est établi avec 50 billions de décimales (50 000 000 000 000) par Timothy Mullican[34],[35], après 303 jours de calcul[36].
Le , ce record est à nouveau battu par Emma Haruka Iwao, calculant cette fois cent billions de décimales en 157 jours[38].
Développement de formules efficaces
Approximation polygonale d'un cercle
Archimède, dans De la mesure du cercle, a créé le premier algorithme pour le calcul de π, basé sur l'idée que le périmètre de n'importe quel polygone inscrit dans un cercle est inférieur à la circonférence du cercle qui, à son tour, est inférieur au périmètre de tout polygone circonscrit à ce cercle. Il a commencé par des hexagones réguliers inscrits et circonscrits, dont les périmètres sont facilement déterminés. Il montre alors comment calculer les périmètres des polygones réguliers ayant deux fois plus de côtés qui sont inscrits et circonscrits autour du même cercle. Il s'agit d'une procédure récursive qui serait décrite aujourd'hui comme suit :
Soit pk et Pk les périmètres de polygones réguliers à k côtés, inscrits et circonscrits autour du même cercle, respectivement. Alors,
Archimède a utilisé le processus décrit pour calculer successivement P12, p12, P24, p24, P48, p48, P96 et p96[39]. Il a obtenu l'encadrement suivant en partant des trois dernières valeurs
.
On ne sait pas pourquoi Archimède s'est arrêté à un polygone de 96 côtés. Héron rapporte dans son Metrica qu'Archimède a continué le calcul dans un livre maintenant perdu, mais lui attribue une valeur incorrecte[40].
Archimède n'a pas utilisé la trigonométrie dans ce calcul et la difficulté d'appliquer la méthode réside dans l'obtention de bonnes approximations pour les racines carrées des calculs impliqués. La trigonométrie a probablement été utilisée par Claudius Ptolémée pour obtenir la valeur de π donnée dans l'Almageste (environ 150 après J.C.)[41].
La formule de Viète, publiée par François Viète en 1593, a été obtenue en utilisant une méthode polygonale, avec des aires (plutôt que des périmètres) de polygones dont les nombres de côtés sont des puissances de 2[42].
Une amélioration trigonométrique par Willebrord Snell (1621) a permis d'obtenir de meilleures limites à partir de la méthode du polygone. Ainsi, des résultats plus précis ont été obtenus, nécessitant moins de côtés[43]. La dernière tentative majeure de calculer π par cette méthode a été effectuée en 1630 par Christopher Grienberger, qui a calculé 39 décimales en utilisant l'amélioration de Snell.
((x, y) = (239, 132) est une solution de l'équation de Pellx2 − 2y2 = −1.)
Les formules de ce type sont connues sous le nom de formules du type de Machin. Elles ont bien été utilisées dans l'ère informatique pour calculer un nombre record de décimales de π mais plus récemment, d'autres formules similaires ont été utilisées.
Par exemple, Shanks et son équipe ont utilisé la formule suivante de Machin en 1961 pour calculer les 100 000 premiers chiffres de π :
,
ainsi que
,
servant de vérification.
Le record enregistré en décembre 2002 par Yasumasa Kanada (voir supra) utilisait les formules suivantes de Machin :
K. Takano (1982).
F. C. W. Störmer (1896).
Autres formules classiques
D'autres formules qui ont été utilisées pour calculer les estimations de π sont entre autres :
où , la suite converge de manière quartique vers π et donne 171 décimales en trois étapes[45],[46], donc plus de 1012 décimales en 20 étapes. Cependant, l'algorithme de Chudnovsky (qui converge linéairement) est plus rapide[réf. nécessaire] que ces formules itératives.
Le premier million de décimales de π et 1/π sont disponibles grâce au projet Gutenberg (voir les liens externes ci-dessous).
Approximations diverses
Historiquement, la base 60 a été utilisée pour les calculs de π. Dans cette base, π peut être approché à huit chiffres (décimaux) significatifs avec le nombre 3:8:29:4460, c'est-à-dire :
(Le chiffre sexagésimal suivant est 0, ce qui permet de tronquer ici pour obtenir une approximation relativement bonne.)
En outre, les expressions suivantes peuvent être utilisées pour estimer π :
précision de trois décimales :
précision de trois décimales :
Karl Popper a spéculé[47] sur l'hypothèse que cette approximation, mentionnée par Émile Borel[48], aurait pu avoir déjà été remarquée par Platon.
La fraction continue de π peut être utilisée pour générer des approximations rationnelles successives. Ce sont les meilleures approximations rationnelles possibles de π par rapport à la taille de leurs dénominateurs. Voici une liste des treize premières fractions :
.
Approximation de l'aire d'un disque par quadrillage
π peut être obtenu à partir d'un cercle, si son rayon et son aire sont connus, en utilisant la relation :
.
Si on considère un cercle de rayon r avec son centre situé en (0, 0), tout point dont la distance à l'origine est inférieure à r sera situé dans le disque. Le théorème de Pythagore donne la distance de n'importe quel point (x, y) au centre :
Les carrés dont le centre se trouve à l'intérieur ou exactement à la limite du cercle peuvent ensuite être comptés en testant si, pour chaque point (x, y),
Le nombre total de point satisfaisant cette condition se rapproche ainsi de l'aire du cercle, qui peut ensuite être utilisée pour calculer une approximation de π.
Cette formule s'écrit :
Autrement dit, on commence par choisir une valeur pour r. En considérant tous les points (x, y) dans lesquels x et y sont des entiers compris entre −r et r. Il suffit alors de diviser la somme représentant l'aire d'un cercle de rayon r, par r2 pour trouver une approximation de π. Par exemple, si r est égal à 5, alors les points considérés sont :
(−5,5)
(−4,5)
(−3,5)
(−2,5)
(−1,5)
(0,5)
(1,5)
(2,5)
(3,5)
(4,5)
(5,5)
(−5,4)
(−4,4)
(−3,4)
(−2,4)
(−1,4)
(0,4)
(1,4)
(2,4)
(3,4)
(4,4)
(5,4)
(−5,3)
(−4,3)
(−3,3)
(−2,3)
(−1,3)
(0,3)
(1,3)
(2,3)
(3,3)
(4,3)
(5,3)
(−5,2)
(−4,2)
(−3,2)
(−2,2)
(−1,2)
(0,2)
(1,2)
(2,2)
(3,2)
(4,2)
(5,2)
(−5,1)
(−4,1)
(−3,1)
(−2,1)
(−1,1)
(0,1)
(1,1)
(2,1)
(3,1)
(4,1)
(5,1)
(−5,0)
(−4,0)
(−3,0)
(−2,0)
(−1,0)
(0,0)
(1,0)
(2,0)
(3,0)
(4,0)
(5,0)
(−5,−1)
(−4,−1)
(−3,−1)
(−2,−1)
(−1,−1)
(0,−1)
(1,−1)
(2,−1)
(3,−1)
(4,−1)
(5,−1)
(−5,−2)
(−4,−2)
(−3,−2)
(−2,−2)
(−1,−2)
(0,−2)
(1,−2)
(2,−2)
(3,−2)
(4,−2)
(5,−2)
(−5,−3)
(−4,−3)
(−3,−3)
(−2,−3)
(−1,−3)
(0,−3)
(1,−3)
(2,−3)
(3,−3)
(4,−3)
(5,−3)
(−5,−4)
(−4,−4)
(−3,−4)
(−2,−4)
(−1,−4)
(0,−4)
(1,−4)
(2,−4)
(3,−4)
(4,−4)
(5,−4)
(−5,−5)
(−4,−5)
(−3,−5)
(−2,−5)
(−1,−5)
(0,−5)
(1,−5)
(2,−5)
(3,−5)
(4,−5)
(5,−5)
Les 12 points (0, ±5), (±5, 0), (±3, ±4), (±4, ±3) sont exactement sur le cercle, tandis que 69 points sont à l'intérieur. Par conséquent, l'aire est 81, et π est approximé à environ 3,24 car 81/52 = 3,24. Les résultats pour certaines valeurs de r sont présentés dans le tableau ci-dessous :
mais avec une convergence presque aussi lente que la précédente (10 termes supplémentaires donnent trois décimales de plus).
Arcsinus 1/2
En notant que
nous avons
avec une convergence telle que pour cinq termes supplémentaires, on obtient au moins trois décimales de plus.
Calcul de la n-ième décimale de π
La formule Bailey-Borwein-Plouffe (BBP) pour calculer π a été découverte en 1995 par Simon Plouffe. En utilisant le système hexadécimal, la formule peut calculer n'importe quelle décimale de π sans avoir à calculer les décimales situées avant elle.
.
En 1996, Simon Plouffe en a déduit un algorithme pour extraire le n-ième chiffre décimal de π, en base 10, en un temps de calcul amélioré d'un facteur O(n3(log n)3). L'algorithme ne nécessite presque aucune mémoire de stockage[55].
.
Le temps de calcul de la formule de Plouffe a été amélioré de O(n2) par Fabrice Bellard, qui a démontré une formule alternative pour calculer π[56].
.
Logiciels pour calculer π
Au cours des années, plusieurs logiciels ont été écrits pour calculer π sur les ordinateurs personnels. Ils impliquent généralement une vérification et un échange de disque efficace pour faciliter des calculs extrêmement coûteux en ressources.
Super PI par Kanada Laboratory à l'université de Tokyo. Il pouvait calculer un million de chiffres en 40 minutes, deux millions de chiffres en 90 minutes et quatre millions de chiffres en 220 minutes sur un Pentium 90 MHz.
PiFast par Xavier Gourdon était le programme le plus rapide sous Windows en 2003. Selon son auteur, il peut calculer un million de décimales en 3,5 secondes sur un Pentium de 2,4 GHz (ordinateur personnel haut-de-gamme à l'époque). PiFast peut également calculer d'autres nombres irrationnels comme e et √2.
QuickPi par Steve Pagliarulo. Comme PiFast, QuickPi peut également calculer d'autres nombres irrationnels comme √2 ou √3. Le logiciel est disponible sur Stu's Pi.
y-cruncher[57] par Alexander Yee est le programme que Shigeru Kondo a utilisé pour calculer le record mondial de 2020 du nombre de décimales et qui a été utilisé pour la plupart des records subséquents. y-cruncher peut également être utilisé pour calculer d'autres constantes, et atteindre des records mondiaux pour plusieurs d'entre elles.
↑Āryabhaṭīya (gaṇitapāda 10), cité et traduit dans (en) Harold R. Jacobs, Geometry: Seeing, Doing, Understanding, 2003 (3e éd.), New York, W. H. Freeman and Company, p. 70.
↑(en) R. C. Gupta, « On the remainder term in the Madhava–Leibniz's series », Ganita Bharati, vol. 14, nos 1-4, , p. 68-71.
↑Al-Kashi, author: Adolf P. Youschkevitch, chief editor: Boris A. Rosenfeld, p. 256
↑(en) Mohammad K. Azarian, « al-Risāla al-muhītīyya: A Summary », Missouri Journal of Mathematical Sciences, vol. 22, no 2, , p. 64-85 (DOI10.35834/mjms/1312233136).
↑(en) Edward Sandifer, « Why 140 Digits of Pi Matter », dans Jurij baron Vega in njegov čas: Zbornik ob 250-letnici rojstva, Ljubljana, DMFA, (ISBN978-961-6137-98-0, OCLC448882242, LCCN2008467244, lire en ligne), p. 17 : « We should note that Vega's value contains an error in the 127th digit. Vega gives a 4 where there should be an, and all digits after that are incorrect. »
↑(en) J. M. Borwein, P. B. Borwein, D. H. Bailey, « Ramanujan, Modular Equations, and Approximations to Pi or How to Compute One Billion Digits of Pi », Amer. Math. Monthly, vol. 96, no 3, , p. 201-219 (DOI10.1080/00029890.1989.11972169)
↑(en) Ronald Mak, Java Number Cruncher : The Java Programmer's Guide to Numerical Computing, Prentice Hall, (lire en ligne), p. 353.