OFFSET
1,2
COMMENTS
From characteristic polynomial of the matrix: {{0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 1}, {133, 0, 0, 137}}.
Polynomial base is Pisot and the large positiuve root is near 1/alpha: 137.00005172368373...
LINKS
Index entries for linear recurrences with constant coefficients, signature (137, 0, 0, 133).
FORMULA
a(n) =Expansion(1/(1 - 137 x - 133 x^4))
MATHEMATICA
mc = {{0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 1}, {133, 0, 0, 137}}; f[x_] = CharacteristicPolynomial[mc, x]; g[x_] = Expand[x^4*f[1/x]]; Table[ SeriesCoefficient[Series[1/g[x], {x, 0, 30}], n], {n, 0, 30}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Roger L. Bagula, Jan 29 2008
STATUS
approved