pi_errors.py

import math,sys

try:
    N = int(sys.argv[1]) + 1
except:
    print('enter the last value')
    sys.exit()

R = range(2,N)

rL = list()

for d in R:
    m = int(math.pi * d)
    n = m + 1
    p = m/d
    q = n/d
    assert p < math.pi < q
    
    if math.gcd(n,d) == 1:
        a = q - math.pi
        rL.append((d,math.log(a,10)))
    if math.gcd(m,d) == 1:
        b = math.pi - p
        rL.append((d,math.log(b,10)))

for x,y in rL:
    print(str(x) + ',' + str(y)[:10])