// Algorithm wide_floating_point e() { unsigned n = 33; my_uint128 num = 1; my_uint128 den = 1; for (unsigned i=1; i <= n; i++) { num = num*i + 1; den = den*i; } return (wide_floating_point) num/fact; }