|
Search: id:A147652
|
|
|
| A147652 |
|
G.f.: 1/(1 - x^4 - x^5 - x^6 + x^10). |
|
+0
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; listen)
|
|
|
OFFSET
|
0,10
|
|
|
LINKS
|
Index entries for sequences related to linear recurrences with constant coefficients
|
|
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: A029115 A029101 A029080 this_sequence A058360 A098527 A035635
Adjacent sequences: A147649 A147650 A147651 this_sequence A147653 A147654 A147655
|
|
KEYWORD
|
nonn
|
|
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
|
|
|
Search completed in 0.002 seconds
|
