|
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]}
for(n=0, 10, for(k=0, n, print1(T(n, k), ", ")); print(""))
/* Print this sequence, which is column 2 */
for(n=0, 20, print1(T(n+2, 2), ", "))
|