Approximation de π

Graphique montrant l'évolution historique de la précision record des approximations numériques de Modèle:Math, mesurée en décimales (représentée sur une échelle logarithmique).

Dans l'histoire des mathématiques, les approximations de la constante [[Pi|Modèle:Math]] ont atteint une précision de 0,04 % de la valeur réelle avant le début de notre ère (Archimède). Au Modèle:S-, des mathématiciens chinois les ont améliorées jusqu'à sept décimales.

De grandes avancées supplémentaires n'ont été réalisées qu'à partir du Modèle:S- (Al-Kashi). Les premiers mathématiciens modernes ont atteint une précision de Modèle:Nb au début du Modèle:S- (Ludolph van Ceulen) et Modèle:Nobr au Modèle:S- (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 Modèle:Math est détenu par William Shanks, qui a calculé 527 décimales correctes vers 1873. Depuis le milieu du Modèle:S-, l'approximation de Modèle:Math est effectuée sur ordinateurs par des logiciels spécifiques.

Le Modèle:Date-, le record est établi avec cent mille milliards de décimales par Emma Haruka Iwao, travaillant sur Google Cloud durant 157 jours.

Histoire

Antiquité

Les approximations les plus connues de Modèle:Math 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.Modèle:Pas clair

Les mathématiques babyloniennes évaluaient généralement Modèle:Math à 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 Modèle:S mini- et le Modèle:Sav- avant notre ère) apporte une meilleure approximation de Modèle:Math, Modèle:Formule, soit proche de 0,5% de la valeur exacteModèle:,Modèle:,Modèle:,.

À 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 une procédure algorithmique pour calculer l'aire d'un disque qui correspond, dans le cas du calcul d'une aire, à l'approximation Modèle:Math (précise à 0,6 %). Le papyrus Rhind n'utilise cependant pas de constante multiplicative analogue à Modèle:Math pour passer de l'aire du carré du rayon (ou du diamètre) à celle du disque. Aucune relation n'est faite avec le rapport de la la circonférence du cercle) au diamètre, ni avec le calcul de cette circonférence. En fait, il n'est jamais question de circonférence du cercle dans les textes mathématiques écrits en hiératique qui nous sont parvenus.

Les calculs astronomiques dans le Śatapatha brāhmaṇa (vers le Modèle:Sav- avant notre ère) utilisent l'approximation fractionnaire Modèle:Formule.

Au Modèle:Sav- avant notre ère, Archimède a prouvé les inégalités strictes Modèle:Fraction < Modèle:Math < Modèle:Fraction, au moyen d'un 96-gone (précision de Modèle:Nb et Modèle:Nb, respectivement). Il est aussi possible de prouver que 22 / 7 dépasse π grâce à un calcul d'intégrales élémentaire, mais cette technique n'existait pas à l'époque d'Archimède.

Au Modèle:S-, Ptolémée a utilisé la valeur Modèle:Fraction, la première approximation connue précise à trois décimales.

Modèle:Article détaillé Le mathématicien chinois Liu Hui, en 263, a calculé que Modèle:Math 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 Modèle:Nb). 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 Modèle:Formule (précision de Modèle:Nb), bien que certains chercheurs pensent plutôt que celui-ci est dû au mathématicien chinois Zu Chongzhi (429-500). Zu Chongzhi est connu pour avoir donné l'encadrement de Modèle:Math entre 3,1415926 et 3,1415927, correct à sept décimales près. Il a donné deux autres approximations de Modèle:Math : Modèle:Formule et Modèle:Formule. La dernière fraction est la meilleure approximation rationnelle possible de Modèle:Math 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 (Modèle:S-), le mathématicien Aryabhata, dans son traité astronomique Āryabhaṭīya, a calculé la valeur de Modèle:Math à cinq chiffres significatifs (Modèle:Formule). Il l'a utilisée pour calculer une approximation de la circonférence de la Terre. 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 Modèle:S-, que le mot signifie non seulement que c'est une approximation, mais que la valeur est « incommensurable » (Modèle:C.-à-d. irrationnelle).

Moyen Âge

Au Modèle:S-, Modèle:Math é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 Modèle:S-, 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 Modèle:Math, maintenant connue sous le nom de formule de Madhava-LeibnizModèle:,, et a donné deux méthodes pour calculer la valeur de Modèle:Math. Ce faisant, il a obtenu la série infinie

$\pi =\sqrt{12}\sum _{k=0}^{\mathrm{\infty }}\frac{(-3{)}^{-k}}{2k+1}=\sqrt{12}\sum _{k=0}^{\mathrm{\infty }}\frac{(-\frac{1}{3}{)}^k}{2k+1}=\sqrt{12}(1-\frac{1}{3\cdot 3}+\frac{1}{5\cdot 3^2}-\frac{1}{7\cdot 3^3}+\cdots )$

et a utilisé les 21 premiers termes pour calculer les 11 premières décimales de Modèle:Math.

L'autre méthode qu'il utilisait était d'ajouter un terme restant à la série originale de Modèle:Math. Il a utilisé le terme restant

$\frac{n^2+1}{4n^3+5n}$

dans l'extension de série infinie de Modèle:Frac pour améliorer l'approximation de Modèle:Math à 13 décimales, lorsque Modèle:Formule = 75.

En 1424, Al-Kashi, un astronome et mathématicien perse, a calculé correctement Modèle:Math à 9 chiffres sexagésimaux. Convertie en base dix, cette approximation est donnée par Modèle:Nobr :

$2\pi \approx 6,28318530717958648$,

ce qui équivaut à

$\pi \approx 3,14159265358979324$.

Il a atteint ce niveau de précision en calculant le périmètre d'un polygone régulier à 3 × 2Modèle:Exp côtés.

Du Modèle:Sp-

Le mathématicien germano-néerlandais Ludolph van Ceulen (vers 1600) a calculé les 35 premières décimales de Modèle:Math à l'aide d'un 2Modèle:Exp-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 Modèle:Math avec un Modèle:Lien.

En 1789, le mathématicien slovène Jurij Vega a calculé les 140 premières décimales de Modèle:Math, cependant les 14 dernières étaient incorrectes ; ce résultat a néanmoins été le record du monde pendant Modèle:Nobr, jusqu'en 1841, lorsque Modèle:Lien 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 Modèle:Nobr à calculer les 707 premières décimales de Modèle:Math. 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 Modèle:Math jusqu'à l'avènement de l'ordinateur, un siècle plus tard.

Modèle:S-

En 1910, le mathématicien indien Srinivasa Ramanujan a trouvé plusieurs séries infinies de Modèle:Math à convergence rapide, comme

$\frac{1}{\pi }=\frac{2\sqrt{2}}{9801}\sum _{k=0}^{\mathrm{\infty }}\frac{(4k)!(1103+26390k)}{(k!{)}^{4}396^{4k}}$

qui calcule à chaque nouvelle itération huit décimales supplémentaires de Modèle:Math. Ses séries sont maintenant la base des algorithmes les plus rapides actuellement utilisés pour calculer Modèle:Math.

À partir du milieu du Modèle:S-, tous les calculs de Modèle:Math 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 Modèle: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 Modèle:Math à Modèle:Nb décimales a été calculé par le mathématicien de l'université du Maryland Daniel 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 Modèle:Math. 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 Modèle:Nb chiffres après la virgule de Modèle:Math ont été publiés en 1962.

En 1989, les frères Chudnovsky ont calculé correctement Modèle:Math à Modèle:Nb décimales grâce au superordinateur IBM 3090 en utilisant la variation suivante de la série infinie de Ramanujan :

$\frac{1}{\pi }=12\sum _{k=0}^{\mathrm{\infty }}\frac{(-1{)}^k(6k)!(13591409+545140134k)}{(3k)!(k!{)}^{3}640320^{3k+3/2}}$.

En 1999, Yasumasa Kanada et son équipe de l'université de Tokyo ont calculé Modèle:Math à Modèle:Nb décimales, grâce au superordinateur Hitachi SR8000/MPP en utilisant une autre variante de la série infinie de Ramanujan.

Modèle:S-

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é Modèle:Math à environ Modèle:Nb décimales en environ 600 heuresModèle:,.

Le 29 avril 2009, Daisuke Takahashi Modèle:Et al., de l'université de Tsukuba, avec un supercalculateur T2K, ont plus que doublé le record précédent, en calculant Modèle:Math à plus de Modèle:Nb décimales, en moins de 74 heures au total (calcul et vérification).

Le 31 décembre 2009, Fabrice Bellard a utilisé un ordinateur personnel pour calculer un peu moins de Modèle:Nb décimales de Modèle:Math. 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 Modèle:Nobr.

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é Modèle:Nb décimales de Modèle:Math. 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.
• Le 16 octobre 2011, Shigeru Kondo a doublé son propre record en calculant Modèle:Nb décimales en utilisant la même méthode, mais assisté d'un meilleur matérielModèle:,Modèle:,.
• En décembre 2013, Kondo a battu son propre record pour la seconde fois lorsqu'il a calculé Modèle:Nb décimales de Modèle:MathPi. La limitation était principalement due à l'espace de stockage.
• Le 8 octobre 2014, un individu non identifié, connu sous le pseudonyme de « houkouonchi », a calculé Modèle:Nb décimales de π.
• Le 15 novembre 2016, Peter Trueb et ses sponsors ont calculé, toujours grâce au programme y-cruncher, Modèle:Nb (Modèle:MathPiModèle:Exp × Modèle:Nb) décimales de Modèle:MathPi. Le calcul a pris (avec 3 interruptions) Modèle:NobrModèle:,.
• Le 14 mars 2019, Google, en la personne d'Emma Haruka Iwao, rend public le nouveau record qui s'établit à Modèle:Nb (la partie entière de Modèle:Nb) décimales de Modèle:Math.
• Le 29 janvier 2020 un nouveau record est établi avec 50 billions de décimales (50 000 000 000 000) par Timothy MullicanModèle:,, après Modèle:Nobr de calcul.
• Le Modèle:Date- un nouveau record est établi avec 62,8 billions de décimales par Thomas Keller et Heiko Rölke de la Modèle:Lien des Grisons en Suisse. Le temps de calcul a été d'un peu plus de Modèle:Nobr.
• Le Modèle:Date-, ce record est à nouveau battu par Emma Haruka Iwao, calculant cette fois cent billions de décimales en 157 jours.

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 Modèle:Math, 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érieure 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 Modèle:Formule et Modèle:Formule les périmètres de polygones réguliers à Modèle:Formule côtés, inscrits et circonscrits autour du même cercle, respectivement. Alors,

$P_{2n}=\frac{2p_n{P}_n}{p_n+P_n},p_{2n}=\sqrt{p_n{P}_{2n}}.$

Archimède a utilisé le processus décrit pour calculer successivement Modèle:Formule et Modèle:Formule. Il a obtenu l'encadrement suivant en partant des trois dernières valeurs

$3+\frac{10}{71}<\pi <3+\frac{1}{7}$.

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.

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 Modèle:Math donnée dans l'Almageste (environ 150 après J.C.).

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.

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. La dernière tentative majeure de calculer Modèle:Math 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.

Formules du type de Machin

Modèle:Article détaillé Pour des calculs rapides, on peut utiliser des formules similaires à celle de John Machin :

$\frac{\pi }{4}=4\arctan \frac{1}{5}-\arctan \frac{1}{239}$

avec le développement en série de Taylor de la fonction arc tangente. Cette formule se vérifie plus facilement à l'aide de coordonnées polaires de nombres complexes, produisant :

$(5+\mathrm{i}{)}^4\cdot (239-\mathrm{i})=2^2\cdot 13^4(1+\mathrm{i})$.

((x, y) = (239, 13Modèle:2) est une solution de l'équation de Pell xModèle:2 − 2yModèle:2 = −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 Modèle:Math 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 Modèle:Nb premiers chiffres de Modèle:Math :

$\frac{\pi }{4}=6\arctan \frac{1}{8}+2\arctan \frac{1}{57}+\arctan \frac{1}{239}$,

ainsi que

$\frac{\pi }{4}=12\arctan \frac{1}{18}+8\arctan \frac{1}{57}-5\arctan \frac{1}{239}$,

servant de vérification.

Le record enregistré en décembre 2002 par Yasumasa Kanada Modèle:Supra utilisait les formules suivantes de Machin :

$\frac{\pi }{4}=12\arctan \frac{1}{49}+32\arctan \frac{1}{57}-5\arctan \frac{1}{239}+12\arctan \frac{1}{110443}$

K. Takano (1982).

$\frac{\pi }{4}=44\arctan \frac{1}{57}+7\arctan \frac{1}{239}-12\arctan \frac{1}{682}+24\arctan \frac{1}{12943}$

F. C. W. Störmer (1896).

Autres formules classiques

D'autres formules qui ont été utilisées pour calculer les estimations de Modèle:Math sont entre autres :

Liu Hui (voir aussi Formule de Viète) :

$\begin{array}{cc}\pi & \approx 768\sqrt{2-\sqrt{2+\sqrt{2+\sqrt{2+\sqrt{2+\sqrt{2+\sqrt{2+\sqrt{2+\sqrt{2+1}}}}}}}}}\\ & \approx 3,141590463236763.\end{array}$

Madhava :

$\pi =\sqrt{12}\sum _{k=0}^{\mathrm{\infty }}\frac{(-3{)}^{-k}}{2k+1}=\sqrt{12}\sum _{k=0}^{\mathrm{\infty }}\frac{(-\frac{1}{3}{)}^k}{2k+1}=\sqrt{12}(\frac{1}{1\cdot 3^0}-\frac{1}{3\cdot 3^1}+\frac{1}{5\cdot 3^2}-\frac{1}{7\cdot 3^3}+\cdots )$

Euler :

$\pi =20\arctan \frac{1}{7}+8\arctan \frac{3}{79}$

Transformation de Newton/Euler :

$\frac{\pi }{2}=\sum _{k=0}^{\mathrm{\infty }}\frac{k!}{(2k+1)!!}=\sum _{k=0}^{\mathrm{\infty }}\frac{2^{k}k{!}^2}{(2k+1)!}=1+\frac{1}{3}(1+\frac{2}{5}(1+\frac{3}{7}(1+\cdots )))$

où (2k + 1)!! désigne le produit des entiers impairs jusqu'à 2k + 1.

Ramanujan :

$\frac{1}{\pi }=\frac{2\sqrt{2}}{9801}\sum _{k=0}^{\mathrm{\infty }}\frac{(4k)!(1103+26390k)}{(k!{)}^{4}396^{4k}}$

David et Gregory Chudnovsky :

$\frac{1}{\pi }=12\sum _{k=0}^{\mathrm{\infty }}\frac{(-1{)}^k(6k)!(13591409+545140134k)}{(3k)!(k!{)}^{3}640320^{3k+3/2}}$

Le travail de Ramanujan est à la base de l'algorithme de Chudnovsky, les algorithmes les plus rapides utilisés pour calculer Modèle:Math.

Algorithmes modernes

Les développements décimaux extrêmement longs de Modèle:Math sont généralement calculés avec des algorithmes tels que l'algorithme de Gauss-Legendre et l'Modèle:Lien. Ce dernier, trouvé en 1985 par Jonathan et Peter Borwein, converge très rapidement :

Pour $y_0=\sqrt{2}-1,a_0=6-4\sqrt{2}$ et

$y_{k+1}=(1-f(y_k))/(1+f(y_k)),a_{k+1}=a_k(1+y_{k+1}{)}^4-2^{2k+3}{y}_{k+1}(1+y_{k+1}+y_{k+1}^2)$

où $f(y)=(1-y^4{)}^{1/4}$, la suite $1/a_k$ converge de manière quartique vers Modèle:Math et donne 171 décimales en trois étapesModèle:,, donc plus de Modèle:Nb décimales en 20 étapes. Cependant, l'algorithme de Chudnovsky (qui converge linéairement) est Modèle:Refnec que ces formules itératives.

Le premier million de décimales de Modèle:Math et Modèle:Math sont disponibles grâce au projet Gutenberg (voir les liens externes ci-dessous).

Approximation de l'aire d'un disque par quadrillage

Approximation numérique de Modèle:Math : les points sont dispersés au hasard dans le carré unité, certains tombent dans le cercle unité. La fraction de points dans le cercle approche Modèle:Math.

Modèle:Math peut être obtenu à partir d'un cercle, si son rayon et son aire sont connus, en utilisant la relation :

$A=\pi r^2$.

Si on considère un cercle de rayon Modèle:Formule avec son centre situé en (0, 0), tout point dont la distance à l'origine est inférieure à Modèle:Formule sera situé dans le disque. Le théorème de Pythagore donne la distance de n'importe quel point (Modèle:Formule, Modèle:Formule) au centre :

$d=\sqrt{x^2+y^2}.$

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 (Modèle:Formule, Modèle:Formule),

$\sqrt{x^2+y^2}\le r.$

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 Modèle:Math.

Cette formule s'écrit :

$\pi =\underset{r\to \mathrm{\infty }}{\lim }\frac{1}{r^2}\sum _{x=-r}^r\sum _{y=-r}^r\{\begin{array}{cc}1& \text{si }\sqrt{x^2+y^2}\le r\\ 0& \text{si }\sqrt{x^2+y^2}>r.\end{array}$

Autrement dit, on commence par choisir une valeur pour Modèle:Formule. En considérant tous les points (Modèle:Formule, Modèle:Formule) dans lesquels Modèle:Formule et Modèle:Formule sont des entiers compris entre −Modèle:Formule et Modèle:Formule. Il suffit alors de diviser la somme représentant l'aire d'un cercle de rayon Modèle:Formule, par Modèle:Formule pour trouver une approximation de Modèle:Math. Par exemple, si Modèle:Formule 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 Modèle:Math est approximé à environ 3,24 car Modèle:Nobr Les résultats pour certaines valeurs de Modèle:Formule sont présentés dans le tableau ci-dessous :

Ce cercle tel qu'il serait tiré sur un graphique de coordonnées cartésiennes. Les points de coordonnées (±3, ±4) et (±4, ±3) sont marqués.

Modèle:Math	Aire	Approximation de Modèle:Math
2	13	3,25
3	29	3,22222
4	49	3,0625
5	81	3,24
10	317	3,17
20	1257	3,1425
100	31417	3,1417
1000	3141549	3,141549

Fractions continues

Modèle:Math a de nombreuses représentations en fraction continue généralisée générées par une règle simple:

$\pi =3+\frac{1^2}{6+\frac{3^2}{6+\frac{5^2}{6+\ddots }}}.$

$\pi =\frac{4}{1+\frac{1^2}{3+\frac{2^2}{5+\frac{3^2}{7+\ddots }}}}.$

Trigonométrie

Série de Gregory-Leibniz

La série de Gregory-Leibniz

$\pi =4\sum _{n=0}^{\mathrm{\infty }}\frac{(-1{)}^n}{2n+1}=4(\frac{1}{1}-\frac{1}{3}+\frac{1}{5}-\frac{1}{7}+-\cdots )=\frac{4}{1+\frac{1^2}{2+\frac{3^2}{2+\frac{5^2}{2+\ddots }}}}$

est la série entière pour arctan(x) en posant Modèle:Formule = 1. Elle converge trop lentement pour avoir un intérêt pratique.

Arcsinus 1

Sachant que Modèle:Math, le développement en série entière de arc sin peut être simplifié pour obtenir :

$\begin{array}{cc}\pi & =2(1+\frac{1}{3}+\frac{1\cdot 2}{3\cdot 5}+\frac{1\cdot 2\cdot 3}{3\cdot 5\cdot 7}+\frac{1\cdot 2\cdot 3\cdot 4}{3\cdot 5\cdot 7\cdot 9}+\frac{1\cdot 2\cdot 3\cdot 4\cdot 5}{3\cdot 5\cdot 7\cdot 9\cdot 11}+\cdots )\\ & =2\sum _{n=0}^{\mathrm{\infty }}\frac{n!}{(2n+1)!!}=\sum _{n=0}^{\mathrm{\infty }}\frac{2^{n+1}n{!}^2}{(2n+1)!}=\sum _{n=0}^{\mathrm{\infty }}\frac{2^{n+1}}{(\genfrac{}{}{0ex}{}{2n}{n})(2n+1)}\\ & =2+\frac{2}{3}+\frac{4}{15}+\frac{4}{35}+\frac{16}{315}+\frac{16}{693}+\frac{32}{3003}+\frac{32}{6435}+\frac{256}{109395}+\frac{256}{230945}+\cdots \end{array}$,

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

$\sin \frac{\pi }{6}=\frac{1}{2}$

nous avons

$\begin{array}{cc}\pi & =6\arcsin \left(\frac{1}{2}\right)=6(\frac{1}{2}+\frac{1}{2\cdot 3\cdot 2^3}+\frac{1\cdot 3}{2\cdot 4\cdot 5\cdot 2^5}+\frac{1\cdot 3\cdot 5}{2\cdot 4\cdot 6\cdot 7\cdot 2^7}+\cdots )\\ & =\frac{3}{16^0\cdot 1}+\frac{6}{16^1\cdot 3}+\frac{18}{16^2\cdot 5}+\frac{60}{16^3\cdot 7}+\cdots =\sum _{n=0}^{\mathrm{\infty }}\frac{3\cdot (\genfrac{}{}{0ex}{}{2n}{n})}{16^n(2n+1)}\\ & =3+\frac{1}{8}+\frac{9}{640}+\frac{15}{7168}+\frac{35}{98304}+\frac{189}{2883584}+\frac{693}{54525952}+\frac{429}{167772160}+\cdots \end{array}$

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 Modèle:Math

La formule Bailey-Borwein-Plouffe (BBP) pour calculer Modèle:Math 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 Modèle:Math sans avoir à calculer les décimales situées avant elle.

$\pi =\sum _{n=0}^{\mathrm{\infty }}(\frac{4}{8n+1}-\frac{2}{8n+4}-\frac{1}{8n+5}-\frac{1}{8n+6}){\left(\frac{1}{16}\right)}^n$.

En 1996, Simon Plouffe en a déduit un algorithme pour extraire le n-ième chiffre décimal de Modèle:Math, en base 10, en un temps de calcul amélioré d'un facteur Modèle:Formule. L'algorithme ne nécessite presque aucune mémoire de stockage.

$\pi +3=\sum _{n=1}^{\mathrm{\infty }}\frac{n{2}^{n}n{!}^2}{(2n)!}$.

Le temps de calcul de la formule de Plouffe a été amélioré de Modèle:Formule par Fabrice Bellard, qui a démontré une formule alternative pour calculer Modèle:Math.

$\pi =\frac{1}{2^6}\sum _{n=0}^{\mathrm{\infty }}\frac{(-1{)}^n}{2^{10n}}(-\frac{2^5}{4n+1}-\frac{1}{4n+3}+\frac{2^8}{10n+1}-\frac{2^6}{10n+3}-\frac{2^2}{10n+5}-\frac{2^2}{10n+7}+\frac{1}{10n+9})$.

Approximations diverses

Un certain nombre de formules ne visent pas à donner une approximation de plus en plus précise de Pi, mais expriment une approximation sous forme compacte, à partir de nombres entiers.

Historiquement, la base 60 a été utilisée pour les calculs de Modèle:Math. Dans cette base, Modèle:Math peut être approché à huit chiffres (décimaux) significatifs avec le nombre 3:8:29:44₆₀, c'est-à-dire :

$3+\frac{8}{60}+\frac{29}{60^2}+\frac{44}{60^3}=3,14159259^{+}$

(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 Modèle:Math :

• précision de trois décimales :

$\frac{22}{7}=3,143^{+}$

• précision de trois décimales :

$\sqrt{2}+\sqrt{3}=3,146^{+}$

Karl Popper a spéculé sur l'hypothèse que cette approximation, mentionnée par Émile Borel, aurait pu avoir déjà été remarquée par Platon.

$\sqrt{15}-\sqrt{3}+1=3,140^{+}$

• précision de quatre décimales:

$\sqrt[3]{31}=3,1413^{+}$

• précision de quatre décimales :

$\sqrt{7+\sqrt{6+\sqrt{5}}}=3,1416^{+}$

• une approximation par Ramanujan, précise à 4 chiffres :

$\frac{9}{5}+\sqrt{\frac{9}{5}}=3,1416^{+}$

• précision de cinq décimales :

$\frac{7^7}{4^9}=3,14156^{+}$

• précision de six décimales :

$\frac{355}{113}=3,141592^{+}$

Découverte de Zu Chongzhi.

• précision de neuf décimales :

$\sqrt[4]{3^4+2^4+\frac{1}{2+(\frac{2}{3}{)}^2}}=\sqrt[4]{\frac{2143}{22}}=3,141592652^{+}$

Provenant de Ramanujan, celui-ci a affirmé que la déesse de Namagiri lui était apparue dans un rêve et lui avait dit la vraie valeur de Modèle:Math.

• précision de dix décimales :

$\frac{63}{25}\times \frac{17+15\sqrt{5}}{7+15\sqrt{5}}=3,1415926538^{+}$

• précision de dix décimales :

$\sqrt[193]{\frac{10^{100}}{11222.11122}}=3,1415926536^{+}$

Cette approximation curieuse résulte de l'observation que la Modèle:193e puissance de Modèle:Math, produisant la suite 1122211125...Modèle:Précision nécessaire

• précision de 18 décimales :

$\frac{80\sqrt{15}(5^4+53\sqrt{89}{)}^{\frac{3}{2}}}{3308(5^4+53\sqrt{89})-3\sqrt{89}}$

• précision à 30 chiffres :

$\frac{\ln (640320^3+744)}{\sqrt{163}}=3,141592653589793238462643383279^{+}$

Dérivé du nombre de Heegner 640320Modèle:3+744. Cette méthode n'admet pas de généralisations évidentes, car il n'y a qu'un nombre fini de nombres de Heegner.

• précision de 52 décimales :

$\frac{\ln (5280^3(236674+30303\sqrt{61}{)}^3+744)}{\sqrt{427}}$

Comme celui ci-dessus, une conséquence du j-invariant.

• précision de 161 décimales :

$\frac{\ln ((2u{)}^6+24)}{\sqrt{3502}}$

où u est un produit de quatre unités quartiques,

$u=(a+\sqrt{a^2-1}{)}^2(b+\sqrt{b^2-1}{)}^2(c+\sqrt{c^2-1})(d+\sqrt{d^2-1})$

et

$\begin{array}{cc}a& =\frac{1}{2}(23+4\sqrt{34})\\ b& =\frac{1}{2}(19\sqrt{2}+7\sqrt{17})\\ c& =(429+304\sqrt{2})\\ d& =\frac{1}{2}(627+442\sqrt{2})\end{array}$.

Méthode basée sur le travail de Daniel Shanks. Semblable aux deux précédentes, mais qui est cette fois un quotient de forme modulaire, à savoir la fonction êta de Dedekind, et où l'argument implique $\tau =\sqrt{-3502}$.

• La fraction continue de Modèle:Math peut être utilisée pour générer des approximations rationnelles successives. Ce sont les meilleures approximations rationnelles possibles de Modèle:Math par rapport à la taille de leurs dénominateurs. Voici une liste des treize premières fractions :

$\frac{3}{1},\frac{22}{7},\frac{333}{106},\frac{355}{113},\frac{103993}{33102},\frac{104348}{33215},\frac{208341}{66317},\frac{312689}{99532},\frac{833719}{265381},\frac{1146408}{364913},\frac{4272943}{1360120},\frac{5419351}{1725033}$.

Logiciels pour calculer Modèle:Math

Au cours des années, plusieurs logiciels ont été écrits pour calculer Modèle:Math 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 [[Yasumasa Kanada|Modèle:Langue]] à l'université de Tokyo. Il pouvait calculer un million de chiffres en Modèle:Nombre, deux millions de chiffres en Modèle:Nombre et quatre millions de chiffres en Modèle:Nombre sur un Pentium Modèle:Nombre.
• 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 Modèle:Nombre sur un Pentium de Modèle:Nombre (ordinateur personnel haut-de-gamme à l'époque). PiFast peut également calculer d'autres nombres irrationnels comme [[E (nombre)|Modèle:Math]] et Modèle:Racine.
• QuickPi par Steve Pagliarulo. Comme PiFast, QuickPi peut également calculer d'autres nombres irrationnels comme Modèle:Racine ou [[Racine carrée de trois|Modèle:Racine]]. Le logiciel est disponible sur Stu's Pi.
• Modèle:Formule-cruncher 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. Modèle:Formule-cruncher peut également être utilisé pour calculer d'autres constantes, et atteindre des records mondiaux pour plusieurs d'entre elles.

Voir aussi

Bibliographie

• Modèle:Article.
• Modèle:Ouvrage
• Modèle:Ouvrage
• Modèle:Ouvrage.
• Modèle:Ouvrage.

Articles connexes

• Journée de pi
• Point de Feynman
• Mémorisation des décimales de pi
• Modèle:Lien

Lien externe

Modèle:Lien web (listes d'approximations de Modèle:Math et liens vers d'autres telles listes)