from fractions import Fraction as R
def Gregory():
"""
Gregory Coefficients (
http://en.wikipedia.org/wiki/Euler-Mascheroni_constant)
1/2, 1/12, 1/24, 19/720, 3/160, 863/60480...
using fraction module to keep it rational.
"""
G = [R(-1,1)]
def Sum(n):
total = R(0,1)
for k in range(0,n):
total += R(G[k],n+1-k)
return -total
n = 1
while True:
G.append(Sum(n))
yield G[-1]
n += 1
====
>>> from anyproject import Gregory
>>> thegen = Gregory()
>>> next(thegen)
Fraction(1, 2)
>>> next(thegen)
Fraction(1, 12)
>>> next(thegen)
Fraction(1, 24)
>>> next(thegen)
Fraction(19, 720)
>>> next(thegen)
Fraction(3, 160)
>>> next(thegen)
Fraction(863, 60480)
>>> next(thegen)
Fraction(275, 24192)
>>> next(thegen)
Fraction(33953, 3628800)
>>> next(thegen)
Fraction(8183, 1036800)
...
Kirby
_______________________________________________
Edu-sig mailing list
Edu-sig@python.org
http://mail.python.org/mailman/listinfo/edu-sig
