I ran into the problem discussed in this thread just the other day, and my solution was to use sage.symbolic.ring. How does this solution compare to the others posted in this thread? Here's (basically) what I did:
from sage.symbolic.ring import var as symbvar
def foo(n, k):
t = symbvar('t')
return exp(t^k).series(t, n+1).coefficient(t, n)*factorial(n)
Comments?
Dan
--
--- Dan Drake <dr...@kaist.edu>
----- KAIST Department of Mathematical Sciences
------- http://mathsci.kaist.ac.kr/~drake
signature.asc
Description: Digital signature
