8.6 Une approximation de π

On peut obtenir une approximation du nombre π avec la formule :

      ┌│  ---∘-----------------------
      │∘           ∘ ------∘-----√---
π ≈ 2k   2-   2 +   2 + ...  2 +   2
k est le nombre de racines carrées. Plus k devient grand et plus cette expression se rapproche du nombre π.

La formule est constituée de l’expression 2 + ∘ ------∘-----√---
  2+ ...  2 +   2 qui est clairement récursive d’où le programme :
# k désigne le nombre de racines  
pour approxpi :k  
tape "Approximation:\  ec (puissance 2 :k) * racine (2- racine (calc :k-2))  
ec "-------------------------  
tape "Pi:\  ec pi  
fin  
 
pour calc :p  
si :p=0 [ret 2][ret 2+racine calc :p-1]  
fin  
 
approxpi 10  
Approximation: 3.141591421568446  
-------------------------  
Pi: 3.141592653589793

On a obtenu les 5 premières décimales ! Si l’on souhaite davantage, il faudra éliminer certaines erreurs de calculs dues aux racines carrées imbriquées. Pour cela nous allons augmenter le nombre de décimales avec la primitive fixedecimales.

fixedecimales 100  
approxpi 100  
Approximation: 3.1415926535897932384626433832795028841973393069670160975807684313880468...  
-------------------------  
Pi: 3.141592653589793238462643383279502884197169399375105820974944592307816406....

Et on obtient à présent 39 décimales...