login
A134874
Expansion of 1/(1 - 137 x - 133 x^4).
0
1, 137, 18769, 2571353, 352275494, 48261760899, 6611863739440, 905825674293229, 124098164230813075, 17001454918435590842, 2329200203203553290874, 319100548313701481849195
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...
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
Sequence in context: A357131 A103878 A068154 * A064104 A300407 A122703
KEYWORD
nonn
AUTHOR
Roger L. Bagula, Jan 29 2008
STATUS
approved