원주율 계산하기

from decimal import Decimal, getcontext

def calculate_pi(n):
    getcontext().prec = n + 1  # set precision
    pi = Decimal(0)
    k = 0
    for i in range(n):
        pi += (Decimal(1) / (16**i)) * (
            Decimal(4) / (8*k+1) - Decimal(2) / (8*k+4) - Decimal(1) / (8*k+5) - Decimal(1) / (8*k+6))
        k += 1
    return str(pi)[:n+1]

n = int(input("Enter the number of digits of pi you want to calculate: "))
pi = calculate_pi(n)
print(f"Pi upto {n} digits: {pi}")

• 1/16^1 [4/(8x1+1) - 2/(8x1+4) - 1/(8x1+5) - 1/(8x1+6)]

• 1/16^2 [4/(8x2+1) - 2/(8x2+4) - 1/(8x2+5) - 1/(8x2+6)]

• ...