def estimate_euler_constant_factorial(number_of_elements=1000):
    """Estimation of euler's constant using its taylor series."""
    euler_constant = 1
    factorial = 1.
    for number in range(1, number_of_elements):
        factorial *= number
        euler_constant += 1. / factorial
    return euler_constant