import sys
import math
a = 1.0
b = 1.0/math.sqrt(2)
t = 1.0/4.0
p = 1.0
while True:
at = (a+b)/2
bt = math.sqrt(a*b)
tt = t - p*(a-at)**2
pt = 2*p
a = at;b = bt;t = tt;p = pt
my_pi = (a+b)**2/(4*t)
print("pi to the last digit = " + str(my_pi))
aW1wb3J0IHN5cwppbXBvcnQgbWF0aAoKYSA9IDEuMApiID0gMS4wL21hdGguc3FydCgyKQp0ID0gMS4wLzQuMApwID0gMS4wCgp3aGlsZSBUcnVlOgogICAgYXQgPSAoYStiKS8yCiAgICBidCA9IG1hdGguc3FydChhKmIpCiAgICB0dCA9IHQgLSBwKihhLWF0KSoqMgogICAgcHQgPSAyKnAKICAgIGEgPSBhdDtiID0gYnQ7dCA9IHR0O3AgPSBwdAogICAgbXlfcGkgPSAoYStiKSoqMi8oNCp0KQoKcHJpbnQoInBpIHRvIHRoZSBsYXN0IGRpZ2l0ID0gIiArIHN0cihteV9waSkp