OFFSET
1,13
FORMULA
See program.
EXAMPLE
MAPLE
row:= proc(n) option remember; local f, i, x; f:= unapply(simplify(sum('cat(a||i) *x^i', 'i'=0..n-1) ), x); unapply(subs(solve({seq(f(i+1)= coeftayl(x/ (1-x-x^4)/ (1-x)^i, x=0, n), i=0..n-1)}, {seq(cat(a||i), i=0..n-1)}), sum('cat(a||i) *x^i', 'i'=0..n-1) ), x); end: T:= (n, k)-> coeff(row(n)(x), x, k): seq(seq(numer(T(n, k)), k=0..n-1), n=1..14);
MATHEMATICA
row[n_] := Module[{f, eq}, f = Function[x, Sum[a[k]*x^k, {k, 0, n-1}]]; eq = Table[f[k+1] == SeriesCoefficient[x/(1-x-x^4)/(1-x)^k, {x, 0, n}], {k, 0, n-1}]; Table[a[k], {k, 0, n-1}] /. Solve[eq] // First]; Table[row[n] // Numerator, {n, 1, 14}] // Flatten (* Jean-François Alcover, Feb 04 2014, after Alois P. Heinz *)
CROSSREFS
KEYWORD
AUTHOR
Alois P. Heinz, Oct 03 2008
STATUS
approved