OFFSET
0,3
FORMULA
log(2)=cot(sum(n>=0,n,(-1)^n*acot(a(n))); let b(0)=log(2), 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]=log(2); ?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
AUTHOR
Benoit Cloitre (abmt(AT)wanadoo.fr), Apr 10 2003
STATUS
approved