OFFSET
1,2
LINKS
Patrick McKinley, Table of n, a(n) for n = 1..12175
MATHEMATICA
ContinuedFraction[x/.First[Solve[x^3-x^2-1==0, x]], 92] (* Stefano Spezia, Jan 21 2024 *)
PROG
(bc)
/* The "test" calculation evaluates the cubic to confirm the calculation of the root. */
define iter(frac)
{j = 0
while(frac > 1){
frac -= 1;
j+=1}
j
return 1/frac}
scale=12578
f=(1+(e(l(((29+3*sqrt(93))/2))/3))+(e(l(((29-3*sqrt(93))/2))/3)))/3
psi=f
test=(psi-1)*psi*psi-1
for(i=0; i<12175; i++)f=iter(f)
(PARI)
\p100 \\ realprecision
contfrac(solve(x = 1, 2, x^3 - x^2 - 1), , 80) \\ Hugo Pfoertner, Jan 21 2024
CROSSREFS
KEYWORD
nonn,cofr
AUTHOR
Patrick McKinley, Jan 20 2024
STATUS
approved