OFFSET
0,3
LINKS
Eric Weisstein's World of Mathematics, Lehmer Cotangent Expansion
FORMULA
Gamma=cot(sum(n>=0, n, (-1)^n*acot(a(n))); let b(0)=Gamma, 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]=Euler; ?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