OFFSET
1,8
MATHEMATICA
(* Adamson's matrix functions alternating golden mean and minimal Pisot*)
digits = 20
NSolve[x^3 - x - 1 == 0, x]
k = -1.32471795724474605
q = k - 1/k
m0 = {{0, 1}, {1, q}}
NSolve[x^2 - x - 1 == 0, x]
k1 = 1.6180339887498949
q1 = k1 - 1/k1
m1 = {{0, 1}, {1, q1}}
m[n_Integer?Positive] := If[Mod[n, 2] == 0, m[n - 1].m0, m[n - 1].m1]
m[0] = {{1, 0}, {0, 1}}
a = Table[Floor[m[n][[2, 2]]], {n, 1, digits}]
CROSSREFS
KEYWORD
sign,uned,obsc,less
AUTHOR
Roger L. Bagula, Dec 02 2003
STATUS
approved