OFFSET
0,2
COMMENTS
(x*exp(x)/(exp(x)-1))^m = 1+sum(n>0, ((-1)^n*sum(k=1..n, (stirling1(m+k,m) *stirling2(n,k))/binomial(m+k,k)))*x^n/n!).
FORMULA
a(n) = 2*(-1)^n*(n+1)!*sum(k=1..n, (stirling1(k+2,2) *stirling2(n,k))/((k+1)*(k+2))), a(0)=1.
PROG
(Maxima) a(n):=if(n=0 then 1 else 2*(-1)^n*(n+1)!* sum((stirling1(k+2, 2) *stirling2(n, k))/((k+1)*(k+2)), k, 1, n);
CROSSREFS
KEYWORD
sign
AUTHOR
Vladimir Kruchinin, Jun 07 2011
STATUS
approved