OFFSET
0,2
COMMENTS
PROG
(PARI) {T(n, k)=local(F=x,
LW=serreverse(x*exp(-x+x*O(x^(n+2)))), M, N, P, m=max(n, k));
M=matrix(m+3, m+3, r, c, F=x; for(i=1, r+c-2, F=subst(F, x, LW)); polcoeff(F, c));
N=matrix(m+1, m+1, r, c, M[r, c]);
P=matrix(m+1, m+1, r, c, M[r+1, c]);
(n-k)!*(P~*N~^-1)[n+1, k+1]}
/* Print triangle : */
for(n=0, 10, for(k=0, n, print1(T(n, k), ", ")); print(""))
/* Print this sequence, which is central terms */
for(n=0, 20, print1(T(2*n, n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jun 28 2016
STATUS
approved