Approximations de Pi: Formule de Lagny

Soyez le premier à donner votre avis sur cette source.

Vue 2 422 fois - Téléchargée 122 fois

Description

Bonjour,

Thomas Fantet de Lagny a trouvé cette formule en 1682:
Pi = 6 * [1/(1*3°) - 1/(3*3¹) + 1/(5*3²) - 1/(7*3³) + ...] / racine(3)
Il l'utilisa en 1719 pour calculer 127 décimales de Pi: les 112 premiers étaient corrects !
C'est un exemple de série alternée basée sur le développement de arctangente avec l'argument 1/racine(3).

Cet article fait partie de la série CodeS-SourceS: Approximations de Pi.

Voici deux codes "équivalents", dont le second est plus condensé et évite de calculer explicitement (2*i + 1):
double Lagny_A(int n) {
  double pi=1, a=1;
  for (int i=1; i<n; ++i) {
    a /= -3;
    pi += a/(2*i+1);
  }
  return 6*pi/1.73205080756887729;
}

double Lagny_B(int n) { // condensé
  double pi=1, a=1;
  for (int i=3, N=2*n+1; i<N; i+=2) pi += (a /= -3) / i;
  return 6*pi/1.73205080756887729;
}
La convergence est assez rapide puisqu'il ne faut 30 itérations pour calculer pi avec la précision double.

Le Zip contient le seul fichier source Lagny.cpp dont voici l'Output:
Lagny_A:
Pi = 6 * [1/(1*3^0) - 1/(3*3^1) + 1/(5*3^2) - 1/(7*3^3) + ...] / racine(3)
  n =  3: Pi=3.15618147156995
  n =  6: Pi=3.14130878546288
  n =  9: Pi=3.14159977381151
  n = 12: Pi=3.14159245428765
  n = 15: Pi=3.14159265952171
  n = 18: Pi=3.14159265340617
  n = 21: Pi=3.14159265359564
  n = 24: Pi=3.14159265358960
  n = 27: Pi=3.14159265358980
  n = 30: Pi=3.14159265358979
  precis: pi=3.14159265358979

Lagny_B:
Pi = 6 * [1/(1*3^0) - 1/(3*3^1) + 1/(5*3^2) - 1/(7*3^3) + ...] / racine(3)
  n =  3: Pi=3.15618147156995
  n =  6: Pi=3.14130878546288
  n =  9: Pi=3.14159977381151
  n = 12: Pi=3.14159245428765
  n = 15: Pi=3.14159265952171
  n = 18: Pi=3.14159265340617
  n = 21: Pi=3.14159265359564
  n = 24: Pi=3.14159265358960
  n = 27: Pi=3.14159265358980
  n = 30: Pi=3.14159265358979
  precis: pi=3.14159265358979


Voir aussi: Thomas Fantet de Lagny.
 
 
Bonne lecture ...

Codes Sources

A voir également

Ajouter un commentaire

Commentaire

bg62
Messages postés
23253
Date d'inscription
samedi 22 octobre 2005
Statut
Modérateur
Dernière intervention
27 août 2019
-
juste par curiosité, j'ai trouvé ceci :
Les 704 premières décimales de ? affichées dans la salle du Palais de la decouverte sont :

3.14159265358979323846264338327950288419716939937510582097
4944592307816406286208998628034825342117067982148086513282
3066470938446095505822317253594081284811174502841027019385
2110555964462294895493038196442881097566593344612847564823
3786783165271201909145648566923460348610454326648213393607
2602491412737245870066063155881748815209209628292540917153
6436789259036001133053054882046652138414695194151160943305
7270365759591953092186117381932611793105118548074462379962
7495673518857527248912279381830119491298336733624406566430
8602139494639522473719070217986094370277053921717629317675
2384674818467669405132000568127145263560827785771342757789
6091736371787214684409012249534301465495853710507922796892
5892354201...

Les 527 décimales soulignées sont celles que William Shanks calcula « à la main » (soit sans calculatrice ni ordinateur) en 1873. En fait, il en calcula 707 mais seules les 527 premières étaient exactes. Joli travail tout de même !


alors d'où et comment, et par qui proviennent ces ... 704 décimales ?
;)

Vous n'êtes pas encore membre ?

inscrivez-vous, c'est gratuit et ça prend moins d'une minute !

Les membres obtiennent plus de réponses que les utilisateurs anonymes.

Le fait d'être membre vous permet d'avoir un suivi détaillé de vos demandes et codes sources.

Le fait d'être membre vous permet d'avoir des options supplémentaires.