OFFSET
1,3
COMMENTS
Recipe: find the roots of z^3-z^2-z-1=0. The real one is the tribonacci constant (A058265) and is of no interest here. The other two are complex conjugates; denote their shared magnitude b. Now this constant is the solution of 2*b^(3x)+b^(4x) = 1.
LINKS
Stanislav Sykora, Table of n, a(n) for n = 1..2000
Wikipedia, Rauzy fractal
EXAMPLE
1.09336416428230663992244770117307346816995623374111426411497299420...
PROG
(PARI) \\ Using 2010 digits realprecision:
b=abs(polroots(Pol([1, -1, -1, -1]))[2]);
a=solve(x=1, 2, 2*b^(3*x)+b^(4*x)-1)
CROSSREFS
KEYWORD
AUTHOR
Stanislav Sykora, Apr 29 2016
STATUS
approved