OFFSET
1,4
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..1000
MATHEMATICA
(* k is a root of x^4 - x^3 - 1: *)
k=1.38027756909761411567330169182273187781662670155876302541177133121124957411;
q=k^2-k-1/k-1/k^2;
m0={{0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 1}, {1, 1, 1, q}};
m[n_]=MatrixPower[m0, n];
Table[Floor[Re[m[n][[4, 4]]]], {n, 1, 1000}]
CROSSREFS
KEYWORD
sign,easy,less
AUTHOR
Roger L. Bagula, Dec 04 2003
EXTENSIONS
Edited by G. C. Greubel, Feb 19 2021
STATUS
approved