arxiv: v5 [math.nt] 17 Jul 2017

Transcription

1 π-formulas AND GRAY CODE. PIERLUIGI VELLUCCI AND ALBERTO MARIA BERSANI arxiv: v5 [math.nt] 7 Jul 07 Abstract. In previous papers we introduced a class of polynomials which follow the same recursive formula as the Lucas-Lehmer numbers, studying the distribution of their zeros and remarking that this distribution follows a sequence related to the binary Gray code. It allowed us to give an order for all the zeros of every polynomial L n. In this paper, the zeros, expressed in terms of nested radicals, are used to obtain two formulas for π: the first can be seen as a generalization of the known formula π = lim n n } {{ } n related to the smallest positive zero of L n; the second is an exact formula for π achieved thanks to some identities valid for L n. 00 MSC: 40A05, Y60. Keywords: π formulas; Gray code; Continued roots; Nested square roots; zeros of Chebyshev polynomials.. Introduction. In this paper we obtain π as the limit of a sequence related to the zeros of the class of polynomials L n (x) = L n (x) created by means of the same iterative formula used to build the well-known Lucas-Lehmer sequence, employed in primality tests Bressoud (989); Finch (003); Koshy (00); Lehmer (930); Lucas (878a,c,b); Ribenboim (988). This class of polynomials was introduced in previous papers Vellucci and Bersani (06a), Vellucci and Bersani (06c). The results obtained here are based on the placement of the zeros of the polynomials L n (x), studied in Vellucci and Bersani (06b), and generalize the well-known relation () π = lim n n }{{} n Zeros have a structure of nested radicals, by means of which we can build infinite sequences of nested radicals converging to π. The ordering of the zeros follows the sequence of the binary Gray code, which is very useful in computer science and in telecommunications Vellucci and Bersani (06b). Starting from the seminal papers by Ramanujan (Ramanujan (000), Berndt (989) pp. 08-), there is a vast literature studying the properties of the so-called continued radicals as, for example: Herschfeld (935); Jonathan M. Borwein (99); Sizer (986); Johnson and Richmond (008); Efthimiou (03); Lynd (04). Other authors investigated the properties of more general continued operations and their convergence. For a nice review of these results, see for example Jones (05), which focuses mainly on continued reciprocal roots. The nested square roots of are a special case of the wide class of continued radicals. They have been studied by several authors. In particular, Cipolla Cipolla (908) obtained,

2 PIERLUIGI VELLUCCI, ALBERTO MARIA BERSANI a very elegant formula for π + i n + i n + + i, in terms of cos ( ) k n, where n+ ik {, } and k n is a constant depending on i,..., i n. Servi Servi (003) rediscovered and extended Cipolla s formula, tying the evaluation of nested square roots of the form () R(b k,..., b ) = b ( ) k + b k b π + bk b + sin 4 where b i {, 0, } for i, to expression ( (3) b k 4 b kb k... b ) kb k...b 8 k+ π to obtain, amongst other results, some nested square roots representations of π: (4) π = lim R,,,,...,,, b b }{{} k terms k+ k where b. Nyblom Nyblom (005), citing Servi s work, ) derived a closed-form expression for () with a generic x that replaces sin in (). Efthimiou Efthimiou (0) proved that the radicals + a a 0 ( b π 4 + a + a 3 +, ai {, } have limits two times the fixed points of the Chebyshev polynomials T n(x), unveiling an interesting relation between these topics. The previous formula is equivalent to (5) ± ± ± ± ±... ±. In Moreno and García-Caballero (0, 03), the authors report a relation between the nested square roots of depth n as ± ± ± ± ± + z, z C and the Chebyshev polynomials of degree n in a complex variable, generalizing and unifying Servi and Nyblom s formulas. In Moreno and García-Caballero (0), the authors propose an ordering of the continued roots (6) b k + b k + + b + ξ, where ξ = and each b i is either or, according to formula ( (7) j(b k,..., b ) = k j ) k k j b k i +, j= for each positive integer k. Formula (6) expresses the nested square roots of in (5), and in Vellucci and Bersani (06b) we gave an alternative ordering for them involving the i=0

3 π-formulas AND GRAY CODE. 3 so-called Gray code which, to the best of our knowledge, is applied for the first time to these topics. Actually, there is a strong connection between Moreno and García-Caballero (0) and Vellucci and Bersani (06b). From (7), we have, for example: j(,, ) =, j(,, ) =, j(,, ) = 3, j(,, ) = 4, j(,, ) = 5, j(-,-,-) =6, j(,, ) = 7 and j(,, ) = 8. If we associate bit 0 to number b i = and bit to number b i =, in the expression of index j, we obtain (,, ) (0, 0, 0) (,, ) (0, 0, ) (,, ) (0,, ) (,, ) (0,, 0) (,, ) (,, 0) (,, ) (,, ) (,, ) (, 0, ) (,, ) (, 0, 0), which are just an example of Gray code. In this paper, after having recalled in Section the most important definitions and properties of the Lucas-Lehmer polynomials, in Section 3 we extend formula (4), using the properties of the zeros of these polynomials (shown in Vellucci and Bersani (06b)), stating and proving Theorem 3.3, which produces infinite numerical sequences converging to π. Besides, under suitable assumptions, Proposition 3. simplifies the expression of (3) listed in Servi s Theorem (Servi (003), formula (8)). We also show that the generalizations of the Lucas-Lehmer map, M a n for a > 0 introduced in Vellucci and Bersani (06a), have the same properties of L n, for what concerns the distribution of the zeros and the approximations of π. We also obtain π not as the limit of a sequence, but equal to an expression involving the zeros of the polynomials L n and M a n for a > 0. Some perspectives of future applications of our results are reported in Section 4.. Preliminaries. In this section we recall properties and useful results from our previous papers (Vellucci and Bersani (06c), Vellucci and Bersani (06a), Vellucci and Bersani (06b)), and therefore we will list them without proofs... The class of Lucas-Lehmer polynomials. We recall below some basic facts about Lucas-Lehmer polynomials (8) L 0 (x) = x ; L n (x) = L n n taken from Vellucci and Bersani (06a). respect to the weight function 4 4 x defined on x (, ). Besides, for each n we have ( ) x (9) L n (x) = T n The polynomials L n (x) are orthogonal with

4 4 PIERLUIGI VELLUCCI, ALBERTO MARIA BERSANI where the Tchebycheff polynomials of first kind Bateman et al. (955); Rivlin (990) satisfy the recurrence relation { T n (x) = xt n (x) T n (x) n T 0 (x) =, T (x) = x from which it easily follows that for the n-th term: ( x ) n ( x + x + ) n x (0) T n (x) =. This formula is valid in R for x ; here we assume instead that T n, defined in R, can take complex values, too. Let x = cos θ, then the polynomials L n (x) admit the representation () L n ( cos θ) = cos ( n θ). When x, we can write x = cos(ϑ), thus x = cos(ϑ); hence, for x =, we can also put ( ) x () ϑ(x) = arctan x + bπ where b is a binary digit; thus, using (), we obtain (3) L n (x) = cos ( n ϑ(x)). Moreover, since L (± ) = 0 ; L (± ) = ; L n (± ) = argument of L n (± ) is 0 for every n 3. By setting further ( ) (4) θ(x) = arctan x x n 3, then the we can write: (5) L n (x) = cos ( n θ(x) + n bπ) = cos ( n θ(x)). As those of the first kind, Tchebycheff polynomial of second kind are defined by a recurrence relation Bateman et al. (955); Rivlin (990): { U 0 (x) =, U (x) = x U n (x) = xu n (x) U n (x) n, which is satisfied by (6) U n (x) = n (x + x ) k (x x ) n k x [, ]. k=0 This relation is equivalent to ( x + ) n+ ( x x ) n+ x (7) U n (x) = x for each x (, ). From continuity of function (6), we observe that (7) can be extended by continuity in x = ±, too.

5 π-formulas AND GRAY CODE. 5 It can therefore be put U n (±) = (±) n (n + ) in (7). From Vellucci and Bersani (06a), for each n we have (8) n ( ) x L i (x) = U n i= and (9) n L i ( cos θ) = sin ( n+ θ ) sin θ i=... An ordering for zeros of Lucas-Lehmer polynomials using Gray code. Given a binary code, its order is the number of bits with which the code is built, while its length is the number of strings that compose it. The celebrated Gray code Gardner (986); Nijenhuis and Wilf (978) is a binary code of order n and length n. We recall below how encapsulated sub-code encapsulated sub-code Figure. Sub-codes for m =, m = 3. a Gray Code is generated; if the code for n bit is formed by binary strings g n,. g n, n g n, n,

6 6 PIERLUIGI VELLUCCI, ALBERTO MARIA BERSANI the code for n bits is built from the previous one in the following way: 0g n,. 0g n, n 0g n, n g n, n g n, n. g n,. Just as an example, we have, for n = : g, = 0 ; g, = ; for n = : g, = 00 ; g, = 0 ; g,3 = ; g,4 = 0; for n = 3: (0) g 3, = 000 g 3, = 00 g 3,3 = 0 g 3,4 = 00 g 3,5 = 0 g 3,6 = g 3,7 = 0 g 3,8 = 00, and so on. Following the notation introduced in Vellucci and Bersani (06b), we recall some preliminaries about Gray code. Definition.. Let us consider a Gray code of order n and length n. A sub-code is a Gray code of order m < n and length m. Definition.. Let us consider a Gray code of order n and length n. An encapsulated sub-code is a sub-code built starting from the last string of Gray code of order n that contains it. Figure () contains some examples of encapsulated sub-codes inside a Gray code (with order 4 and length 6). Let us consider the signs plus and minus in the nested form that expresses generic zeros of L n, as follows: () ± ± ± ±... ± ±. }{{} Obviously the underbrace encloses n signs plus or minus, each one placed before each nested radical. Starting from the first nested radical we apply a code (i.e., a system of rules) that associates bits 0 and to plus and minus signs, respectively. Let us define with {ω(g n, j} j=,..., n the set of all the n nested radicals of the form () ± ± ± ±... ± ± }{{} n signs = ω(g n, n ), where each element of the set differs from the others for the sequence of plus and minus signs. Then: (3) ± ± ± ±... ± ± = ω(g n, n ), }{{} n signs

7 π-formulas AND GRAY CODE. 7 Figure. A possible subcode (orange), where the meaning of the limit (37) is highlighted: in this way the number of symbols 0, on the left of the sub-code, increases. where the notation n, n means that it is possible to obtain n strings formed by n bits. 3. Main results: π formulas involving nested radicals. 3.. Infinite sequences tending to π. Let us consider the writing: ω(... }{{} g m,h ) = ω n m (... g m,h ) n m where the asterisks represent n m bits 0 and. Then we give the following results. Lemma 3.. For all m N one has: ( ) h + (4) ω (0g m,h+ ) = sin m+4 π h [0, m ] Proof. We proceed with induction principle for m to prove (4). If m = : ( ) h + (5) ω (0g,h+ ) = sin 5 π h [0, ] i.e. for h = 0, and ( π ) ω (0g, ) = sin 5 ( ) 3π ω (0g, ) = sin 5 for h =, where g, = 0 and g, =. These formulas are easy to check. Now we are going to check (4) for m + : ( ) h + (6) ω (0g m+,h+ ) = sin m+5 π h [0, m+ ] having assumed it true for m. From Gray Code s definition we have that either a) g m+,h+ = (0, g m,h+ ) or b) g m+,h+ = (, g m, m h). In the former case: ω (0g m+,h+ ) = ω (00g m,h+ )

8 8 PIERLUIGI VELLUCCI, ALBERTO MARIA BERSANI where (7) ω (00g m,h+ ) = = ω (00g m,h+ ) + ω (0g m,h+ ). But in fact: ω (0g m,h+ ) = ω (0g m,h+ ), so (7) becomes ω (00g m,h+ ) = + ω (0g m,h+ ) = 4 ω (0g m,h+ ) ( ) = h sin m+4 π (8) ( ) h + = cos m+4 π ( ) h + = sin m+5 π. Therefore (6) is proved for the case a). Now we assume that g m+,h+ = (, g m, m h): ω (0g m+,h+ ) = ω (0g m, m h) thus (9) ω (0g m, m h) = ω (0g m, m h) = + = + ω (g m, m h) ω (g m, m h) Noting that it follows that (30) ω (0g m, m h) = + ω (g m, m h) ω (0g m, m h) = + 4 ω (0g m, m h) From ω (0g m, m h) = ω (0g m, m h), equation (30) becomes (3) ω (0g m, m h) = + 4 [ ω (0g m, m h)] From (4), we have ( m+ ) (h + ) ω (0g m, m h) = sin m+4 π

9 and equation (3) can be rewritten ω (0g m, m h) = + = + 4 (3) Accordingly: π-formulas AND GRAY CODE. 9 = + 4 [ ω (0g m, m h)] [ ( 4 sin m+ )] (h + ) m+4 π 4 4 cos ( m+ (h + ) m+3 π ( = + sin m+ ) (h + ) m+3 π = + cos = cos ( ) π m+ (h + ) m+3 π ( ) π 4 m+ (h + ) m+4 π. ( (h + m ) ) + ω (0g m, m h) = sin m+5 π. Since the term h + m [ m, m+ ] for h [0, m ], then (6) is fully shown and, with it, the whole proposition. Proposition 3.. For each n m +, h N such that h [0, m ]: ( ) h + (33) ω n m (0...0g m,h+ ) = sin n+ π. Proof. Put n m = δ 0, n m = δ, n m = δ,..., n m k = δ k for 0 k n m. Let us proceed by means of induction principle on n. Fixing m, suppose formula (33) to be true for a generic index δ, (34) ω δ (0...0g m,h+ ) = sin ( h + n+ π and proceed to check the case δ 0. We work on both sides of (34): ( ) h + ω δ (0...0g m,h+ ) = 4 sin n+ π ( h + ω δ (0...0g m,h+ ) = 4 4 cos (35) ) ) n+ π ) ( h + ω δ (0...0g m,h+ ) = 4 cos n+ π ( ) h + + ω δ (0...0g m,h+ ) = 4 cos n+ π ( ) h + ω δ (0...0g m,h+ ) = 4 cos n+ π )

10 0 PIERLUIGI VELLUCCI, ALBERTO MARIA BERSANI whence Thus: or ( ) ( ) ω δ (0...0g m,h+ ) = h + h + cos n+ π = cos n+ π ( ( )) h + ω δ (0...0g m,h+ ) = sin n+ π ( h + ω δ (0...0g m,h+ ) = 4 sin n+ π and ( ) h + ω δ0 (0...0g m,h+ ) = 4 sin n+ π, hence ( ) ( ) ω δ0 (0...0g m,h+ ) = h + h + sin n+ π = sin n+ π. The absolute value can be removed by the proposition s assumptions. Therefore, the inductive step is proved. Let us consider the base step: δ 0 =. Indeed: ( ) h + ω (0g m,h+ ) = sin n m+ m π or, (36) ( ) h + ω(0g m,h+ ) = sin m+4 π which is proved, for all m N, in Lemma 3.. Theorem 3.3. (37) lim n n+ h + for every h N such that h [0, m ] and n > m +. Proof. From Proposition 3., we have (38) n+ h + ), h [0, m ] ω n m (0...0g m,h+ ) = π ( ) ω n m (0...0g m,h+ ) = n+ h + h + sin n+ π that, for a well-know limit, tends to π for n. Example 3.4. With the help of computational tools we show below some iterations of a sequence described by n+ ω n m (0...0g m,h+ ). h + Let us consider m = 3; then g 3, = 000 ; g 3, = 00 ; g 3,3 = 0 ; g 3,4 = 00 ; g 3,5 = 0 ; g 3,6 = ; g 3,7 = 0 ; g 3,8 = 00. We choose the binary string g 3,6 = ; in this case, if m = 3, one has h + = 6 and so h = 5. This means that we are iterating Hence, for n = 8: n+ ω n 3 (0...0 }{{} n 3 9 ) = n+ ω(0000) = ω(0...0).

11 9 π-formulas AND GRAY CODE For n = : 3 ω( ) = and so on The generalized map Mn a = a ( Mn ) a a. In Vellucci and Bersani (06a) we introduced an extension of the map L n, obtained through the iterated formula Mn a = a ( Mn ) a a, a > 0, with M 0 a (x) = x. It follows that (39) M0 a (x) = x ; M a (x) = ax ; M a (x) = 8a 3 x 4 8ax a a Note that the map L n is a particular case of Mn, a obtained by setting a = /. We briefly show that the map Mn a leads to the same π formulas stated in the previous sections. Proposition 3.5. For n we have (40) M a n(x) = a cos(a n x) + o(x ). Proof. We must show that: (4) M a n(x) = a an x + o(x ) where we take into account the McLaurin polynomial of cosine. We proceed by induction. For n = : ( (4) M a (x) = a ax ) a a = a 8ax + o(x ) Let us consider the second order McLaurin polynomial of a cos(4ax): it is just a 8ax + o(x ), thus verifying the relation for n =. Let us now assume (40) is true for a generic n, and deduce that it is also true for n + : Mn+ a = a (Mn) a [ ] a = a a an x + o(x ) a = (43) = a an+ x + o(x ) which is in fact the McLaurin polynomial of a cos(a n+ x). Proposition 3.6. At each iteration the zeros of the map Mn(n a ) have the form (44) ± a ± ± ± ±... ±.

12 PIERLUIGI VELLUCCI, ALBERTO MARIA BERSANI Proof. It is obvious that at n = this statement is valid. Now assume that the (44) is valid for n. We have to prove that it is valid for n +. (45) ax a = ± a ± ± ± or (46) x = a ± 4a ± ± ± and placing under the radical sign (47) x = ± a ± 4a ± ± ± the thesis is obtained. ±... ± ±... ± ±... ± It is possible to prove that zeros of the map M a n+ are related to those of M a n, n π-formulas: not only approximations. From (4) and (5) we obtained Vellucci and Bersani (06a) the following formula: ( ) x (48) L n (x) = cos n arctan x, valid for x [, ] and x ±. This expression is equivalent to (49) L n (x) = ( ((x ) )... ) Moreover, we already observed that, for x =, we have (50) L 0 ( ) = ; L ( ) = 0 ; L ( ) = ; L n ( ) = n 3. The right hand side of (48) vanishes when ( ) (5) n x arctan x = ±π (h + ) ; h N ; x ± i.e., ( ) (5) π < arctan x x = ± π n (h + ) < π, x ± whence (53) ( ) x ( ) x = T n,h

13 π-formulas AND GRAY CODE. 3 where T n,h = tan [ ± π (h + ) ], for h = 0,,..., h n max, and h max defined in this way: from (5) and boundedness of inverse tangent function we have π n (h + ) < π from which h < n therefore h max = n, for n. If the factor T n,h is negative, the solutions of (53) belong to the interval (, ); otherwise x [, ) (, ], if T n,h > 0. We have: ( ) x (54) = ( ) x T n,h x = ± + Tn,h Therefore we can write the zeros of L n in the form (55) x n h = ± ± + tan [ π (h + ) ], n ; 0 h n n Moreover, we know that, for every n, the h-th positive zero of L n (x) has the form: (56) ω(g n, n h) where 0 h n. Equating the two expressions, one finds: (57) whence (58) π = n h + arctan [ ] + tan [ π (h + ) ] = ω(g n, n h) n [ ω(g n, n h) ]. In this way we obtain infinite formulas giving π not as the limit of a sequence, but through an equality involving the zeros of the polynomials L n which is true for every choice of n and h as in (55). Similar considerations can be made for the polynomials Mn. a Since, for x a, (59) Mn(x) a = a cos (a n arctan x ) a x vanishes if (a (60) n arctan x ) a x = ± π (h + ) i.e., (6) arctan (a x ) a x = ± π (h + ), n

14 4 PIERLUIGI VELLUCCI, ALBERTO MARIA BERSANI then (6) (a x ) = ( a x ) T n,h where T n,h = tan [ ± π (h + ) ], with h = 0,,..., n. n Furthermore: ( (63) a x ) T n,h > 0 Inequality a x > 0 is verified for ( x < a x > a. If T n,h is negative, the ) [ solutions of (6) belong to the interval a, a, otherwise x ) ( ] a, a a, a, if T n,h > 0. On the other hand: (64) ( a x ) ( = T n,h a x ) a x = ± + Tn,h from which: (65) x n h = ± ± a + tan [ π (h + ) ] n Since, from (44), the zeros of Mn(x) a are proportional to the zeros of L n (x), we can say that also the n positive zeros of Mn, a in decreasing order, follow the order given by the Gray code: (66) ω(g a n, n h) Equating the two expressions we find again the identity: (67) π = n h + arctan [ ω(g n, n h) ] 4. Discussion and perspectives. In previous papers (Vellucci and Bersani (06a) and Vellucci and Bersani (06b)) we introduced a class of polynomials which follow the same recursive formula as the Lucas- Lehmer numbers, studying the distribution of their zeros and remarking that this distributions follows a sequence related to the binary Gray code. It allowed us to give an order for all the zeros of every polynomial L n, Vellucci and Bersani (06b). In this paper, the zeros, expressed in terms of nested radicals, are used to obtain two formulas for π: the first (i.e., formula (37)) can be seen as a generalization of the known formula (), because the latter can be seen as the case related to the smallest positive zero of L n ; the second (i.e., formula (58)) gives infinite formulas reproducing π not as the limit of a sequence, but through an equality involving the zeros of the polynomials L n. The proof of the π-formulas is based on Proposition 3.. Actually, Proposition 3. can be fundamental for further studies, too. In fact, it not only allows to get the main results of this paper, but also allows the evaluation of nested square roots of as: ω n m (0...0g m,h+ ) = ± ± ± for each n m +, h N such that h [0, m ]. This is a result to put in evidence and to generalize in future researches, for example following interesting insights suggested by

15 π-formulas AND GRAY CODE. 5 paper Zimmerman and Ho (008), where the authors defined the set S of all continued radicals of the form + a + a + a a 0 (with a 0 =, a k {, } for k = 0,,..., n ) and investigated some of its properties by assuming that the limit of the sequence of radicals exists. Further studies can involve more general classes of polynomials, as suggested by Pincherle (98), who focused on the realness of the roots in L n = L n a. References Bateman, Harry, Arthur Erdélyi, W Magnus, Fritz Oberhettinger, and Francesco Giacomo Tricomi (955), Higher transcendental functions, volume. McGraw-Hill New York. Berndt, Bruce C. (989), Ramanujan s notebooks. Part II. Springer-Verlag, New York. Bressoud, David M. (989), Factorization and primality testing. Undergraduate Texts in Mathematics, Springer-Verlag, New York. Cipolla, Michele (908), Intorno ad un radicale continuo. Period. mat. l insegnamento second., Ser, 3, Efthimiou, Costas J. (0), A class of periodic continued radicals. Am. Math. Mon., 9, Efthimiou, Costas J. (03), A class of continued radicals. Am. Math. Mon., 0, Finch, Steven R. (003), Mathematical constants, volume 94 of Encyclopedia of Mathematics and its Applications. Cambridge University Press, Cambridge. Gardner, M. (986), The binary Gray code. In Knotted Doughnuts and Other Mathematical Entertainments, chapter, Freeman, New York. Herschfeld, Aaron (935), On infinite radicals. Am. Math. Mon., 4, Johnson, Jamie and Tom Richmond (008), Continued radicals. Ramanujan J., 5, Jonathan M. Borwein, G. de Barra (99), Nested radicals. Am. Math. Mon., 98, Jones, Dixon J. (05), Continued reciprocal roots. Ramanujan J., 38, Koshy, Thomas (00), Fibonacci and Lucas numbers with applications. Pure and Applied Mathematics (New York), Wiley-Interscience, New York. Lehmer, D. H. (930), An extended theory of Lucas functions. Ann. Math., 3, Lucas, Edouard (878a), Theorie des Fonctions Numeriques Simplement Periodiques. Amer. J. Math.,, Lucas, Edouard (878b), Theorie des Fonctions Numeriques Simplement Periodiques. Amer. J. Math.,, Lucas, Edouard (878c), Theorie des Fonctions Numeriques Simplement Periodiques. [Continued]. Amer. J. Math.,, Lynd, Chris D. (04), Using difference equations to generalize results for periodic nested radicals. Amer. Math. Monthly,, Moreno, Samuel G. and Esther M. García-Caballero (0), Chebyshev polynomials and nested square roots. J. Math. Anal. Appl., 394, Moreno, Samuel G. and Esther M. García-Caballero (03), On Viète-like formulas. J. Approx. Theory, 74, 90. Nijenhuis, Albert and Herbert S Wilf (978), Combinatorial algorithms: for computers and calculators. Academic Press, New York. Nyblom, M. A. (005), More nested square roots of. Amer. Math. Monthly,, 8 85.

16 6 PIERLUIGI VELLUCCI, ALBERTO MARIA BERSANI Pincherle, S (98), Sulle radici reali delle equazioni iterate di una equazione quadratica. Rendiconti delle sedute dell Accademia Nazionale dei Lincei, 7, Ramanujan, Srinivasa (000), Collected papers of Srinivasa Ramanujan. AMS Chelsea Publishing, Providence, RI. Ribenboim, Paulo (988), The book of prime number records. Springer-Verlag, New York. Rivlin, Theodore J. (990), Chebyshev polynomials, second edition. Pure and Applied Mathematics (New York), John Wiley & Sons, Inc., New York. Servi, L. D. (003), Nested square roots of. Amer. Math. Monthly, 0, Sizer, Walter S. (986), Continued roots. Math. Mag., 59, 3 7. Vellucci, Pierluigi and Alberto Maria Bersani (06a), The class of Lucas-Lehmer polynomials. Rend. Mat. Appl., 37, Vellucci, Pierluigi and Alberto Maria Bersani (06b), Ordering of nested square roots of according to Gray code. Ramanujan J. Published online, in press. Vellucci, Pierluigi and Alberto Maria Bersani (06c), Orthogonal polynomials and Riesz bases applied to the solution of Love s equation. Math. Mech. Complex Syst., 4, Zimmerman, Seth and Chungwu Ho (008), On infinitely nested radicals. Math. Mag., 8, 3 5. Pierluigi Vellucci: Department of Economics, University of Roma TRE, via Silvio D Amico 77, 0045 Rome, Italy pierluigi.vellucci@uniroma3.it Alberto Maria Bersani: Department of Mechanical and Aerospace Engineering, Sapienza University, Via Eudossiana n. 8, 0084 Rome, Italy alberto.bersani@sbai.uniroma.it