Pi_Bailey_Borwein_Plouffe.jl

# Pi_Bailey_Borwein_Plouffe.jl
# Juliaohjelma laskee Piin likiarvon halutulla tarkkuudella
# Perustuu Bailey-Borwein-Plouffe kaavaan vuodelta 1995
# Juhani Kaukoranta 12.9.2019

function Pii(tarkkuus,bitprecision)
	setprecision(BigFloat,bitprecision+40)

    s = BigFloat(0.0)
    for k = 0:tarkkuus+40
        s = s + BigFloat(1.0)/big(16)^k*(BigFloat(4.0)/(8*k+1)-BigFloat(2.0)/(8*k+4)-BigFloat(1.0)/(8*k+5)-BigFloat(1.0)/(8*k+6))
    end
    return s
end

println("Kuinka monen numeron tarkkuudella haluat Piin likiarvon? ")
tarkkuustxt = readline()
tarkkuus = parse(Int,tarkkuustxt)
bitprecision = round(Int,log2(BigFloat(10.0)^tarkkuus))
println("tarkkuus = ",tarkkuus)
println("bitprecision =", bitprecision)
piinlikiarvo = Pii(tarkkuus,bitprecision)
# asetetaan tulostustarkkuus halutuksi
setprecision(BigFloat,bitprecision)
println("Piin likiarvo ",tarkkuus," tarkkuudella = ",BigFloat(piinlikiarvo))