OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..1000
FORMULA
a(n) = [x^n] x/((1-x-x^4)*(1-x)^(n-1)).
MAPLE
a:= n-> coeftayl(x/(1-x-x^4)/(1-x)^(n-1), x=0, n):
seq(a(n), n=0..30);
# second Maple program:
a:= proc(n) option remember; `if`(n<5, n*(n+1)*(n^2-4*n+6)/6,
a(n-4)+(2*(35*n^3-207*n^2+310*n-78)*a(n-1)-(203*n^3
-1244*n^2+1891*n-130)*a(n-2)+(2*n-7)*(7*n-19)*n*
(10*a(n-3)-2*a(n-5)))/((7*n-26)*(n-1)^2))
end:
seq(a(n), n=0..30); # Alois P. Heinz, Aug 18 2019
MATHEMATICA
a[n_] := SeriesCoefficient[x/(1-x-x^4)/(1-x)^(n-1), {x, 0, n}];
Table[a[n], {n, 0, 30}] (* Jean-François Alcover, May 10 2022 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Oct 03 2008
STATUS
approved