Example of Bailey–Borwein–Plouffe digit extraction formula

BBP.py

# Example of Bailey–Borwein–Plouffe digit extraction formula
extra = 2

def BBP_8k(k,j):
	return (8*k+j)

def BBP_mod(n, d):
	return (16**n) % d

def BBP_inner(n, j):
	total = 0
	for k in range(n+extra):
		tmp = BBP_8k(k, j)
		total += BBP_mod(n-k, tmp) / tmp
	return total

def BBP_one(n):
	T1 = 4.0 * BBP_inner(n, 1)
	T2 = -2.0 * BBP_inner(n, 4)
	T3 = -1.0 * BBP_inner(n, 5)
	T4 = -1.0 * BBP_inner(n, 6)
	return T1 + T2 + T3 + T4

for n in range(200):
	p = 4
	s = BBP_one(n) % 1
	s *= 16
	print(hex(int(s))[2:], end=' ')
print()