OFFSET
0,1
FORMULA
sqrt(5) = cot(Sum_{n>=0} (-1)^n*acot(a(n))).
Let b(0) = sqrt(5), b(n) = (b(n-1)*floor(b(n-1))+1)/(b(n-1)-floor(b(n-1))) then a(n) = floor(b(n)).
MATHEMATICA
Floor[NestList[(#*Floor[#]+1)/(#-Floor[#]) &, Sqrt[5], 7]] (* Stefano Spezia, Apr 24 2025 *)
PROG
(PARI) bn=vector(100);
b(n)=if(n<0, 0, bn[n]);
bn[1]=sqrt(5);
for(n=2, 10, bn[n]=(b(n-1)*floor(b(n-1))+1)/(b(n-1)-floor(b(n-1))));
a(n) = floor(b(n+1));
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Apr 10 2003
STATUS
approved
