login
A081789
Continued cotangent for cos(1).
0
0, 1, 3, 31, 1896, 4466398, 51545652359935, 7892797673015743066290382811, 79208317703482281896053478218775642522933780224074809198, 195824124605155340760338534459689354174423866004140103635128465652592661899008938733448659610261796713081609084855
OFFSET
0,3
REFERENCES
D. H. Lehmer, A cotangent analogue of continued fractions, Duke Math. J., 4 (1935), 323-340.
FORMULA
cos(1)=cot(sum(n>=0, n, (-1)^n*acot(a(n))); let b(0)=cos(1), 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]=cos(1); ?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