OFFSET
1,3
COMMENTS
Also decimal expansion of (1+(1+(1+ ... )^(1/k))^(1/k))^(1/k), with k integer and k<0. Case k=9.
LINKS
Paolo P. Lava, Table of n, a(n) for n = 1..1000
EXAMPLE
1.0850702454914508283368958640973142340506536310308965814...
MAPLE
with(numtheory); P:=proc(q, h) local a, n; a:=(q+1)^(1/h);
for n from q by -1 to 1 do a:=(1+a)^(1/h); od;
print(evalf(a, 1000)); end: P(1000, 9);
MATHEMATICA
Root[(#^9-#-1)&, 1] // RealDigits[#, 10, 100]& // First (* Jean-François Alcover, Feb 18 2014 *)
RealDigits[x/.FindRoot[x^9-x-1==0, {x, 1}, WorkingPrecision->100]][[1]] (* Harvey P. Dale, Jul 31 2017 *)
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Paolo P. Lava, Oct 11 2013
STATUS
approved