|
| |
|
|
A084765
|
|
a(n)=2a(n-1)^2-1, a(0)=1, a(1)=5.
|
|
1
| | |
|
|
|
OFFSET
| 0,2
|
|
|
COMMENTS
| Product((1+1/a(k)), k=1,..,n) converges to sqrt(3/2).
A subsequence of A001079 (cf. formula), which must contain any prime occurring in A001079. The initial term a(0)=1 seems rather unnatural; using the recurrence relation it would yield the constant sequence 1,1,1,... Note that this sequence corresponds to sequence b(n) in Shallit's paper, which starts only at offset n=1. [From M. F. Hasler (www.univ-ag.fr/~mhasler), Sep 27 2009]
|
|
|
REFERENCES
| H. S. Wilf, Limit of a sequence, Elementary Problem E 1093, Amer. Math. Monthly 61 (1954), 424-425
|
|
|
LINKS
| J. O. Shallit, Rational numbers with non-terminating, non-periodic modified Engel-type expansions, Fib. Quart., 31 (1993), 37-40.
|
|
|
FORMULA
| With a=5+2sqrt(6), b=5-2sqrt(6): a(n+1)=(a^(2^n)+b^(2^n))/2.
A084765(n+1)=A001079(2^n). [From M. F. Hasler (www.univ-ag.fr/~mhasler), Sep 27 2009]
|
|
|
MATHEMATICA
| For n>0: b[n_] := b[n] = 2 b[n - 1]^2 - 1; b[1] = 5 Table[b[n], {n, 1, 8}]
|
|
|
CROSSREFS
| Cf. A084764.
Sequence in context: A064618 A193199 A075986 * A203411 A082795 A059008
Adjacent sequences: A084762 A084763 A084764 * A084766 A084767 A084768
|
|
|
KEYWORD
| easy,nonn
|
|
|
AUTHOR
| Mario Catalani (mario.catalani(AT)unito.it), Jun 04 2003
|
| |
|
|
