E calculator

e.py

from gmpy2 import mpz, get_context

# 19 loops maxes out 64-bit floating point number
# ~~ 2.718281828459045 == math.e ~~
def E(prec=19):
    num, den = 0, 1
    for i in range(1, prec):
        # add fraction keeping GCD
        # allows for greater precision
        num = (num + 1) * i
        den *= i
    return num / den

def gmp_E(prec=19):
    num, den = mpz(0), mpz(1)
    for i in range(1, prec):
        num = (num + 1) * i
        den *= i
    return num / den