Pi_Bailey_Borwein_Plouffe.py

# Ohjelma laskee Piin likiarvon halutulla tarkkuudella
# Perustuu Bailey-Borwein-Plouffe kaavaan vuodelta 1995
# Kaava antaa k:nnen numeron 16-järjestelmän avulla ilman että tarvitsee laskea
# aikaisempia numeroita
# Tämä ohjelma laskee samaan kaavaan nojautuen halutun määrän Pi:n desimaaleja.
# Juhani Kaukoranta 14.3.2022

from decimal import *
import time

def pii(precision):
    getcontext().prec=precision
    s = Decimal(0)
    for k in range(0,precision):
        s = s + 1/Decimal(16)**k*(
        Decimal(4)/(8*k+1) -
        Decimal(2)/(8*k+4) -
        Decimal(1)/(8*k+5) -
        Decimal(1)/(8*k+6))
    return s

tarkkuus = int(input("Kuinka monen numeron tarkkuudella haluat Piin likiarvon? "))
time0 = time.perf_counter()
piinlikiarvo = pii(tarkkuus)
time1 = time.perf_counter()
print ("Aikaa kului ",time1-time0," sekuntia")
print ("Piin likiarvo ",tarkkuus," tarkkuudella = ",piinlikiarvo)