OFFSET
0,3
MATHEMATICA
NSolve[x^4 - x^3 - 1 == 0, x] b01 = 1.380277569097614 b11 = -0.8191725133961645 b21 = 0.2194474721492751 - 0.9144736629677264*I b31 = 0.2194474721492751 + 0.9144736629677264*I f[n_] := N[(b01^n - b11^n - b21^n - b31^n)/(b01 - b11 - b21 - b31)]; digits = 200 c = Table[Floor[Re[f[n]]], {n, 1, digits}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Roger L. Bagula, Sep 21 2004; corrected Oct 06 2004
STATUS
approved