Re: [sage-support] Fast Ramanujan's tau function

2010-01-17 Thread William Stein
On Sun, Jan 17, 2010 at 10:11 PM, Rolandb wrote: > Hi, > > Consider: > > def tau(m): >    q = PowerSeriesRing(QQ, 'q', default_prec=12).gen() >    pq=prod([(1-q^k)^24 for k in range(1,m)]) >    return pq.coefficients()[:m] > > tau(20) > [1, -24, 252, -1472, 4830, -6048, -16744, 84480, -113643, -11

[sage-support] Fast Ramanujan's tau function

2010-01-17 Thread Rolandb
Hi, Consider: def tau(m): q = PowerSeriesRing(QQ, 'q', default_prec=12).gen() pq=prod([(1-q^k)^24 for k in range(1,m)]) return pq.coefficients()[:m] tau(20) [1, -24, 252, -1472, 4830, -6048, -16744, 84480, -113643, -115920, 534612, -370944, -577738, 401856, 1217160, 987136, -6905934,