OFFSET
0,2
PROG
(PARI) {a(n)=if(n<0, 0, sum(k=1, n+1, (matrix(n+k, n+k, m, j, if(m==j, 1, if(m==j+1, -m+1, -(m-j-1)*polcoeff(log(sum(i=0, m, (i+3)!/3!*x^i)), m-j-1))))^-1)[n+1, k]))}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Aug 07 2005
STATUS
approved