OFFSET
0,2
PROG
(PARI) {a(n)=if(n<0, 0, (matrix(n+3, n+3, 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+3, 3])}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Aug 07 2005
STATUS
approved