OFFSET
0,1
FORMULA
sqrt(5)=cot(sum(n>=0, n, (-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))
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