|
| |
|
|
A147652
|
|
G.f.: 1/(1 - x^4 - x^5 - x^6 + x^10).
|
|
1
| |
|
|
1, 0, 0, 0, 1, 1, 1, 0, 1, 2, 2, 2, 2, 3, 4, 5, 5, 7, 8, 10, 12, 15, 18, 22, 26, 32, 40, 48, 58, 70, 86, 105, 128, 154, 188, 229, 279, 339, 412, 501, 610, 742, 902, 1098, 1335, 1624, 1975, 2403, 2923, 3556, 4324
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 0,10
|
|
|
LINKS
| Index entries for sequences related to linear recurrences with constant coefficients, signature (0,0,0,1,1,1,0,0,0,-1).
|
|
|
MATHEMATICA
| f[x_] = 1 - x^4 - x^5 - x^6 + x^10; g[x] = ExpandAll[x^10*f[1/x]]; a = Table[SeriesCoefficient[Series[1/g[x], {x, 0, 50}], n], {n, 0, 50}] (*vector Matrix Markov*); CompanionMatrix[p_, x_] := Module[{cl = CoefficientList[p, x], deg, m}, cl = Drop[cl/Last[cl], -1]; deg = Length[cl]; If[deg == 1, {-cl}, m = RotateLeft[ IdentityMatrix[deg]]; m[[ -1]] = -cl; Transpose[m]]]; M = Transpose[CompanionMatrix[g[x], x]]; v[0] = Table[a[[n]], {n, 1, 10}]; v[n_] := v[n] = M.v[n - 1]; Table[v[n][[1]], {n, 0, 50}]
|
|
|
CROSSREFS
| Sequence in context: A029101 A194815 A029080 * A058360 A193942 A098527
Adjacent sequences: A147649 A147650 A147651 * A147653 A147654 A147655
|
|
|
KEYWORD
| nonn,easy
|
|
|
AUTHOR
| Roger L. Bagula (rlbagulatftn(AT)yahoo.com), Nov 09 2008
|
|
|
EXTENSIONS
| Definition corrected by N. J. A. Sloane (njas(AT)research.att.com), Nov 10 2008
|
| |
|
|
