OFFSET
0,15
MATHEMATICA
p[x_] = (x^(10) + x^9 - x^7 - x^6 - x^5 - x^4 - x^3 + x + 1);
q[x_] = Expand[x^10*p[1/x]];
a = Table[SeriesCoefficient[Series[1/q[x], {x, 0, 50}], n], {n, 0, 20}];
b[0] := 0;
b[n_] := b[n] = If[n <= 20, a[[n]], b[n - 1 - Mod[n, 20]] + a[[1 + Mod[n, 20]]]];
Table[b[n], {n, 0, 100}]
CROSSREFS
KEYWORD
sign,base,less
AUTHOR
Roger L. Bagula, Feb 06 2010
STATUS
approved