OFFSET
1,2
LINKS
Iain Fox, Table of n, a(n) for n = 1..20000 (first 1000 terms from Vincenzo Librandi).
Wikipedia, Pisot number
FORMULA
A root of x^5-x^4-x^3+x^2-1.
EXAMPLE
1.443268791270373107628127607386...
MAPLE
fsolve(x^5-x^4-x^3+x^2-1, x, 1.4..1.5) ;
MATHEMATICA
Root[Function[x, x^5-x^4-x^3+x^2-1], 1] // RealDigits[#, 10, 90]& // First (* Jean-François Alcover, Feb 20 2014 *)
PROG
(PARI) default(realprecision, 20080); x=solve(x=1, 2, x^5 - x^4 - x^3 + x^2 - 1); for (n=1, 20000, d=floor(x); x=(x-d)*10; write("b228777.txt", n, " ", d)); \\ Iain Fox, Oct 23 2017
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
R. J. Mathar, Sep 04 2013
STATUS
approved