gauss-legendre pi calculation algorithm

gauss-legendre-pi.py

def pi_digits():
    q, r, t, k, n, l = 1, 0, 1, 1, 3, 3
    decimal = []

    for _ in range(100):
        if 4 * q + r - t < n * t:
            decimal.append(n)
            print (n)
            time.sleep(1)
            q, r, t, k, n, l = (
                10 * q, 10 * (r - n * t), t, k, (10 * (3 * q + r)) // t - 10 * n, l
            )
        else:
            q, r, t, k, n, l = (
                q * k, (2 * q + r) * l, t * l, k + 1, (q * (7 * k + 2) + r * l) // (t * l), l + 2
            )

    return decimal

pi_digits()