OFFSET
0,2
FORMULA
sqrt(3) = cot(Sum_{n>=0} (-1)^n*acot(a(n))).
Let b(0) = sqrt(3), 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[3], 8]] (* Stefano Spezia, Apr 24 2025 *)
PROG
(PARI) \p500
bn=vector(100); bn[1]=sqrt(3); b(n)=if(n<0, 0, bn[n]);
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
