0% found this document useful (0 votes)
16 views16 pages

Eisenstein Series and Approximations To: Bruce C. Berndt and Heng Huat Chan

About Pi

Uploaded by

ambuj.kumar
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
16 views16 pages

Eisenstein Series and Approximations To: Bruce C. Berndt and Heng Huat Chan

About Pi

Uploaded by

ambuj.kumar
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 16

Illinois Journal of Mathematics

Volume 45, Number 1, Spring 2001, Pages 75–90


S 0019-2082

EISENSTEIN SERIES AND APPROXIMATIONS TO π

BRUCE C. BERNDT AND HENG HUAT CHAN

Dedicated to K. Venkatachaliengar

1. Introduction
On page 211 in his lost notebook, in the pagination of [19], Ramanujan
listed eight integers, 11, 19, 27, 43, 67, 163, 35, and 51 at the left margin. To
the right of each integer, Ramanujan recorded a linear equation in Q3 and
R2 . Although Ramanujan did not indicate the definitions of Q and R, we can
easily (and correctly) ascertain that Q and R are the Eisenstein series

X n3 q n
(1.1) Q(q) :=1 + 240
n=1
1 − qn

and

X n5 q n
(1.2) R(q) :=1 − 504 ,
n=1
1 − qn
where |q| < 1. To the right of each equation in Q3 and R2 , Ramanujan
entered an equality involving π and square roots. (For the integer 51, the
linear equation and the equality involving π, in fact, are not recorded by
Ramanujan.)
The equations in Q3 and R2 cannot possibly hold for all values of q with
|q| < 1. Thus, the first task was to find the correct value
√ of q for each equation.
After trial and error we found that q = − exp(−π n), where n is the integer
at the left margin. (We later read that√ K. Venkatachaliengar [20, p. 135] had
also discovered that q = − exp(−π n).) The equalities in the third column
lead to approximations to π that are remindful of approximations given by
Ramanujan in his famous paper on modular equations and approximations
to π [15], [18, p. 33] and studied extensively by J. M. and P. B. Borwein [6,
Chap. 5]. As will be seen, this page in the lost notebook is closely connected
with theorems connected with the modular j-invariant stated by Ramanujan

Received June 18, 1999; received in final form February 16, 2000.
2000 Mathematics Subject Classification. 33C05, 33E05, 11F11, 11R29.

c 2001 University of Illinois

75
76 BRUCE C. BERNDT AND HENG HUAT CHAN

on the last two pages of his third notebook [10] and proved by the authors
[4], [3, pp. 309–322].
In Section 2, we prove a general theorem from which the linear equations
in Q3 and R2 in the second column follow as corollaries. In Sections 3 and 4,
we offer two methods for proving the equalities in the third column and show
how they lead to approximations to π. Ramanujan’s equations for π lead in
Section 4 to certain numbers tn (defined by (4.28)). The numbers tn guide us
in Section 5 to a proof of a general series formula for 1/π, which is equivalent
to a formula found by D. V. and G. V. Chudnovsky [13] and the Borweins [9].
The first series representations for 1/π of this type were found by Ramanujan
[15], [18, pp. 23–39], but first proved by the Borweins [8]. The work of the
Borweins [6], [7], [8] and Chudnovskys [13] significantly extends Ramanujan’s
work. One of Ramanujan’s series for 1/π yields 8 digits of π per term, while
one of the Borweins’ [7] gives 50 digits of π per term. Our simpler method
enables us in Section 5 to determine a series for 1/π which yields about 73 or
74 digits of π per term.

2. Eisenstein series and the modular j–invariant


Recall the definition of the modular j–invariant j(τ ),
Q3 (q)
(2.1) j(τ ) = 1728 , q = e2πiτ , Im τ > 0.
Q3 (q)
− R2 (q)
In particular, if n is a positive integer,
√ 
Q3

3 + −n
(2.2) j = 1728 3 n 2 ,
2 Qn − Rn
where, for brevity, we set
√ √
(2.3) Qn := Q(−e−π n
) and Rn := R(−e−π n
).
In his third notebook, at the top of page 392 in the pagination of [17], Ramanu-
jan defined a certain function Jn of singular moduli, which, as the authors [4]
easily showed, has the representation
s  √ 
1 3 3 + −n
(2.4) Jn = − j .
32 2

Hence, from (2.2) and (2.4),


Q3n
(2.5) (−32Jn )3 = 1728 .
Q3n− Rn2
After a simple manipulation of (2.5), we deduce the following theorem.
EISENSTEIN SERIES AND APPROXIMATIONS TO π 77

Theorem 2.1. For each positive integer n,


 3 !  3
8 8
(2.6) Jn + 1 Q3n − Jn Rn2 = 0,
3 3

where Jn is defined by (2.4), and Qn and Rn are defined by (2.3).

Examples 2.2 [19, p. 211]. We have


539Q311 − 512R11
2
=0,
(83 + 1)Q319 − 83 R19
2
=0,
(403 + 9)Q327 − 403 R27
2
=0,
(803 + 1)Q343 − 803 R43
2
=0,
(4403 + 1)Q367 − 4403 R67
2
=0,
(533603 + 1)Q3163 − 533603 R163
2
=0,
√ 3 3
√ 3 2
((60 + 28 5) + 27)Q35 − (60 + 28 5) R35 =0,
and
√ √ √ √
((4(4 + 17)2/3 (5 + 17))3 + 1)Q351 − (4(4 + 17)2/3 (5 + 17))3 R51
2
= 0.

Proof. In [4], [3, pp. 310, 311], we showed that


J11 =1,


 J19 =3,
 1/3
J =5 · 3 , J43 =30,

 27



(2.7) J67 =165, J163 =20, 010,
 √ !4 √ !


 √ 1+ 5 √ 2/3 5+ 17
 J35 = 5 , J51 =3(4 + 17) .



2 2

Using (2.7) in (2.6), we readily deduce all eight equations in Qn and Rn . 

3. Eisenstein series and equations in π – first method


Define, after Ramanujan,

X nq n
(3.1) P (q) := 1 − 24 , |q| < 1,
n=1
1 − qn
and put

(3.2) Pn := P (−e−π n
).
Next, set
1/2
(3.3) bn = {n(1728 − jn )}
78 BRUCE C. BERNDT AND HENG HUAT CHAN

and
  
1 Qn 6
(3.4) an = bn 1 − Pn − √ .
6 Rn π n
The numbers an and bn arise in series representations for 1/π proved by the
Chudnovskys [13] and the Borweins [8], namely,

1 1 X (6k)! an + kbn
(3.5) =√ ,
π −jn k=0 (3k)!(k!)3 jnk
where (c)0 = 1, (c)k = c(c + 1) · · · (c + k − 1), for k ≥ 1, and
 √ 
3 + −n
jn = j .
2
These authors have calculated an and bn for several values of n, but we are
uncertain if these calculations are theoretically grounded. We show how (3.3)
and (3.4) lead to a formula from which Ramanujan’s equalities in the third
column on page 211 of [19] follow.
From (2.6), we easily see that
3 !−1/2
8
Qn 1 3 Jn +1
(3.6) =√ 3 ,
Rn Qn 8
Jn
3
and from (3.4), we find that
Qn 6 √ √ an √
 
(3.7) − nPn = 6 n − n.
Rn π bn
The substitution of (3.6) into (3.7) leads to the following theorem.

Theorem 3.1. If Pn , bn , an , and Jn are defined by (3.2)–(3.4) and (2.4),


respectively, then
3 !1/2
8
√ √
   
1 6 an 3 Jn +1
(3.8) √ nPn − = n 1−6 3 .
Qn π bn 8
Jn
3

Examples 3.2 [19, p. 211]. We have


√ √
 
1 6
√ 11P11 − = 2,
Q11 π
√ √
 
1 6
√ 19P19 − = 6,
Q19 π

  r
1 6 6
√ 27P27 − =3 ,
Q27 π 5

  r
1 6 3
√ 43P43 − =6 ,
Q43 π 5
EISENSTEIN SERIES AND APPROXIMATIONS TO π 79


  r
1 6 6
√ 67P67 − =19 ,
Q67 π 55

  r
1 6 3
√ 163P163 − =362 ,
Q163 π 3335
s
√ √
 
1 6 2
√ 35P35 − =(2 + 5) √ ,
Q35 π 5

 
1 6
√ 51P51 − =...
Q51 π

Ramanujan’s formulation of the first of these examples is apparently given


by
√ 6
11 − + · · · √
π
s  3  = 2.
1
1 − 240 √ ···
eπ 11
(The denominator with √ 13 in the numerator√is unreadable.) Further equalities
are even briefer, with Qn replaced by ·. Note that Pn is replaced by
“1 + · · · ” in Ramanujan’s examples. Also observe that Ramanujan did not
record the right side when n = 51. Because it is unwieldy, we also have not
recorded it. However, readers can readily complete the equality, since J51 is
given in (2.7) and a51 and b51 are given in the next table.

Proof. The first six values of an and bn were calculated by the Borweins [8,
pp. 371, 372]. The values for n = 35 and 51 were calculated by the present
authors. We record all 8 pairs of values for an and bn in the following table.

n an bn
11 60 616
19 300 4104
27 1116 18216
43 9468 195048
67 122124 3140424
163 163096908√ 6541681608√
35 1740 + 768 √5 32200 + 14336√5
51 11820 + 2880 17 265608 + 64512 17

If we substitute these values of an and bn in Theorem 3.1, we obtain, after


some calculation and simplification, Ramanujan’s equalities. 
80 BRUCE C. BERNDT AND HENG HUAT CHAN

Theorem 3.1 and the last set of examples yield approximations to π. Let
rn denote the algebraic expression on the right side in (3.8). If we use the
expansions
√ √
Pn = 1 + 24e−π n − · · ·
p
and Qn = 1 − 120e−π n + · · · ,
we easily find that
 √ 
6 24 n + 120rn −π√n
π=√ 1− √ e + ··· .
n − rn n − rn
We thus have proved the following theorem.

Theorem 3.3. We have


6
π≈√ =: An ,
n − rn
with the error approximately equal to

n + 5rn −π√n
144 √ e ,
( n − rn )2
where rn is the algebraic expression on the right side in (3.8).

See Ramanujan’s paper [15], [18, p. 33] for other approximations to π of


this sort.
In the table below, we record the decimal expansion of each approximation
An and the number Nn of digits of π agreeing with the approximation.

n An Nn
11 3.1538 . . . 1
19 3.1423 . . . 2
27 3.1416621 . . . 3
43 3.141593 . . . 5
67 3.14159266 . . . 7
163 3.14159265358980 . . . 12
35 3.141601 . . . 3
51 3.14159289 . . . 6

4. Eisenstein series and equations in π – second method


Set P := P(q) := P (−q), Q := Q(q) := Q(−q), R := R(q)  := R(−q),
∆ := ∆(q) := Q3 (q) − R2 (q), and J := J(q) := 1728/j 3+τ
2 , where q =
e2πiτ . Set
 1/3

(4.1) z 4 := Q = ,
J
EISENSTEIN SERIES AND APPROXIMATIONS TO π 81

by (2.1). Then, by (4.1) and the definition of ∆,


r
p ∆√ √
(4.2) R = Q3 − ∆ = 1 − J = z 6 1 − J.
J
Recall the differential equations [16], [17, p. 142]
(4.3)
dP P 2 (q) − Q(q) dQ P (q)Q(q) − R(q) dR P (q)R(q) − Q2 (q)
q = , q = , q = ,
dq 12 dq 3 dq 2
which yield the associated differential equations
(4.4)
dP P2 (q) − Q(q) dQ P(q)Q(q) − R(q) dR P(q)R(q) − Q2 (q)
q = ,q = ,q = .
dq 12 dq 3 dq 2
Now, by rearranging the second equation in (4.4), with the help of (4.1) and
(4.2), we find that

R(q) 12q dz
(4.5) P(q) = + .
Q(q) z dq
From the chain rule and (4.5), it follows that, for any positive integer n,
R(q n ) 12q dz(q n )
P(q n ) = + .
Q(q n ) nz(q n ) dq
Subtracting (4.5) from the last equality and setting
z(q)
m := ,
z(q n )
we find that
R(q n ) R(q) q dz(q n ) q dz(q)
nP(q n ) − P(q) =n n
− + 12 n
− 12
Q(q ) Q(q) z(q ) dq z(q) dq
R(q n ) R(q) q dm
(4.6) =n − − 12 .
Q(q n ) Q(q) m dq
dm dm
Our next aim is to replace in (4.6) by (J(q), J(q n )). From (2.1),
dq dJ
the definition of J, (4.4), (4.1), and (4.2), upon differentiation, we find that
dJ (3Q2 Q0 − 2RR0 )Q3 − 3Q2 Q0 (Q3 − R2 )
q =
dq Q6
{Q (PQ − R) − R(PR − Q2 )}Q3 − Q2 (PQ − R)(Q3 − R2 )
2
=
Q6
3
RQ − R 3
R∆ 6
√ J 2

(4.7) = = = z 1 − J = z J 1 − J,
Q4 Q4 Q
82 BRUCE C. BERNDT AND HENG HUAT CHAN

which implies that


1 dJ(q)
(4.8) z 2 (q) = p q .
J(q) 1 − J(q) dq
Replacing q by q n in (4.8) and simplifying, we deduce that
1 dJ(q n )
(4.9) z 2 (q n ) = p q .
nJ(q n ) 1 − J(q n ) dq

Using (4.8) and (4.9), we conclude that


p
2 J(q n ) 1 − J(q n ) dJ(q)
(4.10) m =n .
J(q) 1 − J(q) dJ(q n )
p

It is well known that there is a relation (known as the class equation) be-
tween j(τ ) and j(nτ ) for any integer n [14, p. 231, Theorem 11.18(i)]. With
the definition of J given at the beginning of this section, the class equation
translates to a relation between J(q) and J(q n ). It follows that,
dJ(q)
(4.11) = F (J(q), J(q n )),
dJ(q n )
for some rational function F (x, y). Thus, by (4.10) and (4.11), we may differ-
entiate m with respect to J, and so, by (4.7) and the definition of m(q),
q dm q dm
2 = 2z(q)z(q n ) 2
m(q) dq z (q) dq
dm dJ √ dm2 (q)
= 2z 2 (q n )m(q) q = z 2 (q n )J 1 − J .
dJ dq dJ
Using this in (4.6), we deduce that
nP(q n ) − P(q) R(q n ) R(q) 2 n
p dm2
(4.12) = n − − 6z (q )J(q) 1 − J(q) .
z(q)z(q n ) Q(q n ) Q(q) dJ

If we put q = e−π/ n
, n > 0, (4.12) takes the shape
√ √

−π n

−π/ n R(e−π n ) R(e−π/ n )
nP(e ) − P(e ) =n √ − √
Q(e−π n ) Q(e−π/ n )
√ √ q √
− 6z 2 (e−π n )J(e−π/ n ) 1 − J(e−π/ n )
dm2  −π√n √ 
(4.13) × J(e ), J(e−π/ n ) .
dJ
It is well known that [12, p. 84]
√ √
(4.14) J(e−π/ n
) = J(e−π n
).
EISENSTEIN SERIES AND APPROXIMATIONS TO π 83

Furthermore, if
∞ ∞
X 2 X
(4.15) ϕ(q) = qn and ψ(q) = q n(n+1)/2 , |q| < 1,
n=−∞ n=0

then [2, p. 127, Entries 13(iii), (iv)],


(4.16) Q(q) = z24 (1 + 14x2 + x22 ),

and

(4.17) R(q) = z26 (1 + x2 )(1 − 34x2 + x22 ),


where [2, p. 122–123, Entries 10(i), 11(iii)]
ψ 4 (q 2 )
(4.18) z2 := ϕ2 (q) and x2 := 16q .
ϕ4 (q)
Replacing q by −q in (4.16) and (4.17), and using (4.18), we find that
(4.19) Q(q) = ϕ8 (−q) − 224qϕ4 (−q)ψ 4 (q 2 ) + 162 q 2 ψ 8 (q 2 )

and

R(q) = (ϕ4 (−q) − 16qψ 4 (q 2 ))


(4.20) × (ϕ8 (−q) + 544qϕ4 (−q)ψ 4 (q 2 ) + 162 q 2 ψ 8 (q 2 )).
Using the transformation formula [2, p. 43, Entry 27(ii)],

ϕ(e−π/t ) = 2e−πt/4 tψ(e−2πt )
in (4.19) and (4.20), we deduce that
√ √
(4.21) R(e−π/ n
) = −n3 R(e−π n
)

and
√ √
(4.22) Q(e−π/ n
) = n2 Q(e−π n
).
Using (4.14), (4.21), and (4.22), we may rewrite (4.13) as
√ √
nP(e−π n
) − P(e−π/ n
)

R(e−π n ) √
2 −π n
q dm2
= 2n −π

n
− 6z (e )J n 1 − J1/n (Jn , Jn ) ,
Q(e ) dJ
dm2
 p q  √
(4.23) = 2n 1 − Jn − 6Jn 1 − J1/n (Jn , Jn ) z 2 (e−π n ),
dJ
where

(4.24) Jk = J(e−π k
), k > 0.
√ √
This gives the first relation between P(e−π n
) and P(e−π/ n
).
84 BRUCE C. BERNDT AND HENG HUAT CHAN

Recall the definitions of Ramanujan’s function f (−q) and the Dedekind


eta-function η(τ ), namely,

Y
f (−q) := (1 − q k ) =: e−2πiτ /24 η(τ ), q = e2πiτ , Im τ > 0.
k=1

The function f satisfies the well-known transformation formula [2, p. 43]


√ √ √ √
(4.25) n1/4 e−π n/24
f (e−π n
) = e−π/(24 n)
f (e−π/ n
), n > 0.
Logarithmically differentiating (4.25) with respect to n, multiplying both sides
by 48n3/2 /π, rearranging terms, and employing the definition of P (q) given
in (3.1), we find that
√ √ √
12 n
(4.26) = nP(e−π n ) + P(e−π/ n ).
π
√ √
This gives a second relation between P(e−π n ) and P(e−π/ n ).
Now adding (4.23) and (4.26) and dividing by 2, we arrive at

dm2

 p  √
−π n 6 n q
nP(e )= + n 1 − Jn − 3Jn 1 − J1/n (Jn , Jn ) z 2 (e−π n ),
π dJ
or, by (4.1),
p !
1 − J1/n dm2
  p
1 6
√ Pn − √ = 1 − Jn 1 − 3Jn √ (Jn , Jn ) ,
Qn nπ n 1 − Jn dJ
where Qn is defined by (2.3), and Pn is defined by (3.2). (Be careful; Jn 6= Jn ,
where Jn is defined by (2.4).)
We record the last result in the following theorem, which should be com-
pared with Theorem 3.1.

Theorem 4.1. If Pn , Qn , and Jn are defined by (3.2), (2.3), and (4.24),


respectively, then
 
1 6 p
(4.27) √ Pn − √ = 1 − J n tn ,
Qn nπ
where
p !
1 − J1/n dm2
(4.28) tn := 1 − 3Jn √ (Jn , Jn ) .
n 1 − Jn dJ

Observe that, by (4.1), (4.2), and Theorem 2.1,


3 !1/2
8
p
3 J n + 1
(4.29) 1 − Jn = 3 .
8
3 Jn
EISENSTEIN SERIES AND APPROXIMATIONS TO π 85


Hence, the values of 1 − Jn for those n given on page 211 of the Lost Note-
book follow immediately from (2.7). In order to rederive Examples 3.2, it
suffices to compute tn .

Theorem 4.2. √ If n > 1 is an odd positive integer, then tn lies in the ring
class field of Z[ −n].

Proof. By differentiating (4.10) with respect to J(q), we conclude that


p
dm2 1 − J(q n )
(4.30) = p G(J(q), J(q n )),
dJ(q) 1 − J(q)
for some rational function G(x, y). Using (4.28) and (4.30), we find that √ tn
can be expressed in terms of Jn . Since Jn is in the ring class field of Z[ −n]
when n > 1 is odd and squarefree [14, p. 220, Theorem 11.1], we complete
our proof. 

An equivalent form of Theorem 4.2 is first mentioned without proof and


conditions on n by the Chudnovskys on page 391 of [13].
Theorem 4.2 allows us √to devise an empirical process for deriving tn when-
ever the class group of Q( −n) is of the type Zr2 , r ∈ N, where r = 0 refers to
imaginary quadratic fields with class number 1. By (3.6), (4.27), and (4.29),
we find that
 
Qn 6
(4.31) tn = Pn − √ .
Rn nπ
When r = 0, we numerically compute the right hand side of (4.31) using the
definitions of Pn , Qn , and Rn , and then using the command “minpoly(tn ,1)”
on the computer algebra system MAPLE V, we derive the values of tn for
n = 3, 7, 11, 19, 27, 43, 67 and 163. We summarize our findings in the following
table.

n tn
3 0
5
7 21
32
11 77
32
19 57
160
27 253
640
43 903
33440
67 43617
77265280
163 90856689
86 BRUCE C. BERNDT AND HENG HUAT CHAN

When r = 1, we are in a situation where tn satisfies a polynomial of degree


2. We use the command “minpoly(tn ,2)” and derive the resulting polynomial
satisfied by tn . Using this method, we could, in particular, determine tn for
n = 35 and 51, namely,
√ √
1504 + 576 5 144 400 17
t35 = and t51 = + .
4123 329 5593
Using the values of tn given in the tables above and Theorem 4.1, we rederive
Examples 3.2.
So far, our method is still similar to that of the Borweins and Chudnovskys.
Of course, we could continue to use “minpoly” and derive minimal polynomials
satisfied by tn for r = 2, 3, and 4. However, this means that we have to solve
polynomials of degress 4, 8, and 16, respectively, which is cumbersome when
the degree of the polynomial exceeds 4. We now describe a new process for
deriving tn , for r > 1. We illustrate our method using n = 21.
When n = 21, we numerically compute t21 and t7/3 and use “minpoly” to
deduce that
30493 3820 √
t21 + t7/3 = + 7
34279 34279
and
√ 28800 1935 √
3(t21 − t7/3 ) = + 7.
34279 34279
Solving for t21 yields
30493 645 √ 4800 √ 1910 √
t21 = + 21 + 3− 7.
68558 68558 34279 34279
The largest n for which tn can be computed using this new empirical
method is n = 3315, namely,
1095255033002752301233099478037584
t3315 :=
2050242335692983321671746996556833

1006588064225996719872149534306400 √ √
+ 17 5
34854119706780716468419698941466161

692779168175128551453280427070000 √
+ 17
34854119706780716468419698941466161

136434536163779492503565618457696 √
− 5
2050242335692983321671746996556833

400179322879781860521299209248000 √
+ 13
26653150364008783181732710955238829
EISENSTEIN SERIES AND APPROXIMATIONS TO π 87

1077564413015882021519209726762688 √ √ √
+ 13 17 5
453103556188149314089456086239060093

120226784218523863048087030809600 √ √
+ 17 13
64729079455449902012779440891294299

239369594240980944219359445009600 √ √
+ 13 5.
26653150364008783181732710955238829
This value is deduced from determining the quadratic numbers
t3315 + t1105/3 + t663/5 + t221/15 ,

221(t3315 − t1105/3 + t663/5 − t221/15 ),

13(t3315 + t1105/3 − t663/5 − t221/15 ),
and √
5(t3315 − t1105/3 − t663/5 + t221/15 )

in Q( 85). We show in the next section that knowledge of t3315 and J3315
leads to a new series for 1/π which converges at a rate of 73/74 decimal places
per term. This series appears to be the fastest known convergent series for
1/π which involves quadratic radicals. The previous fastest convergent series,
giving 50 decimal places per term, was obtained by the Borweins√ [8]. Their
series involves quartic radicals which arise from the field Q( −1555), which
has class number 4.

5. tn , Jn and the Ramanujan–Borweins–Chudnovskys series for 1/π


Let ϕ(q) and ψ(q) be given as in (4.15), and recall the definition of the
hypergeometric function 2 F1 ,

X (a)k (b)k uk
2 F1 (a, b; c; u) := , |u| < 1.
(c)k k!
k=0

In [20], using Ramanujan’s differential equations (4.3) and the relations


Q(q 2 ) = z24 (1 − x2 + x22 ) and R(q 2 ) = z26 (1 + x2 )(1 − 2x2 )(1 − x2 /2),
where z2 and x2 are given by (4.18), Venkatachaliengar proved the inversion
formula
z2 = 2 F1 ( 21 , 12 ; 1; x2 ).
His method can be applied whenever we have relations of the form
Q(q) = z 4 F (x) and R(q) = z 6 G(x),
for some function z and x. The relations (4.1) and (4.2) indicate that we are
in such a situation with z := Q and x := J. Invoking Venkatachaliengar’s
88 BRUCE C. BERNDT AND HENG HUAT CHAN

method (see [11], [10] for examples of such calculations), we readily deduce
that
1 5
(5.1) z = 2 F1 ( 12 , 12 ; 1; J).
We are now ready to give a proof of an equivalent form of the general series
for 1/π given by (3.5). As far as we know, a proof of this series has never
been written down in the literature. First, by Clausen’s formula [1, p. 116],
we find that
2
(5.2) z 2 = 2 F1 ( 12
1 5
, 12 ; 1; J) = 3 F2 ( 16 , 56 , 12 ; 1, 1; J),
where

X (a)k (b)k (c)k uk
3 F2 (a, b, c; d, e; u) = , |u| < 1.
(d)k (e)k k!
k=0
This implies, by (5.1), that

dz X
(5.3) 2z = Ak kJk−1 ,
dJ
k=1

where
1 5 1
  
6 k 6 k 2 k
Ak := .
(k!)3
By (4.5) and (4.7), we deduce that
 
dz 1 P
(5.4) 2z = √ − z2 .
dJ 6J 1−J
Substituting (5.4) into (5.3) and using (5.2), we deduce that

P X
(5.5) √ = Ak (6k + 1)Jk .
1 − J k=0

Next, set q = e−π n
and deduce from (5.5) that

Pn X
(5.6) √ = Ak (6k + 1)Jkn .
1 − Jn k=0

On the other hand, by Theorem 4.1 and (4.1),


6 p
(5.7) Pn = √ + zn2 1 − Jn tn .
π n
Using (5.7), (5.2) (with z replaced by zn and J replaced by Jn ) and (5.6), we
conclude that
∞ 1
 5 1
6 1 X
6 k 6 k 2 k k
(5.8) √ √ = Jn (6k + 1 − tn ).
n 1 − Jn π (k!)3
k=0
EISENSTEIN SERIES AND APPROXIMATIONS TO π 89

It can be shown that (5.8) is equivalent to the Ramanujan–Borweins–Chudnov-


skys series. The series (5.8) enables us to write down for each n a series for
1/π if we know the values of tn and Jn .
We have already determined the value of t3315 in Section 4. It suffices
to compute J3315 in order to write down a series for 1/π associated with
n = 3315. We first quote the identity [4] [5],

(λ2 − 1)(9λ2n − 1)3
 
3 + −3n
(5.9) j = −27 n ,
2 λ2n
where √ √
eπ n/2 f 6 (e−π n/3 )
λn = √ √ ,
3 3 f 6 (e−π 3n )
where f (−q) is defined prior to (4.25). Since [5]
√ !24 √ !6
2 5+1 √ 6 15 + 221 √
λ1105 = (4 + 17) (8 + 65)6 ,
2 2
the value of J3315 follows immediately from (5.9). The values J3315 and t3315 ,
when substituted into (5.8), yield the series which we mentioned at the end
of Section 4.
Acknowledgement. We are grateful to Seung Hwan Son for helping us
find the correct value of q for the linear equations in Q3 and R2 . The first
author is grateful to the National University of Singapore for financially sup-
porting his visit there when this research was conducted.

References
[1] G. E. Andrews, R. A. Askey, and R. Roy, Special functions, Cambridge University
Press, Cambridge, 1999.
[2] B. C. Berndt, Ramanujan’s notebooks, part III, Springer–Verlag, New York, 1991.
[3] , Ramanujan’s notebooks, part V, Springer–Verlag, New York, 1998.
[4] B. C. Berndt and H. H. Chan, Ramanujan and the modular j–invariant, Canad. Math.
Bull. 42 (1999), 427–440.
[5] B. C. Berndt, H. H. Chan, S.-Y. Kang, and L.-C. Zhang, A certain quotient of eta-
functions found in Ramanujan’s lost notebook, Pacific J. Math., to appear.
[6] J. M. Borwein and P. M. Borwein, Pi and the AGM, Wiley, New York, 1987.
[7] , Ramanujan’s rational and algebraic series for 1/π, Indian J. Math. 51 (1987),
147–160.
[8] , More Ramanujan-type series for 1/π, in: Ramanujan Revisited (G.E. An-
drews et al., eds.), Academic Press, Boston, 1988, pp. 359–374.
[9] , Some observations on computer aided analysis, Notices Amer. Math. Soc. 39
(1992), 825–829.
[10] H. H. Chan, On Ramanujan’s cubic transformation formula for 2 F1 (1/3, 2/3; 1; z),
Math. Proc. Cambridge Philos. Soc. 124 (1998), 193–204.
[11] H. H. Chan and Y. L. Ong, On Eisenstein series and ∞ m2 +mn+2n2 , Proc.
P
m,n=−∞ q
Amer. Math. Soc. 127 (1999), 1735–1744.
[12] K. Chandrasekharan, Elliptic functions, Springer–Verlag, Berlin, 1985.
90 BRUCE C. BERNDT AND HENG HUAT CHAN

[13] D. V. Chudnovsky and G. V. Chudnovsky, Approximation and complex multiplica-


tion according to Ramanujan, in: Ramanujan Revisited (G.E. Andrews et al., eds.),
Academic Press, Boston, pp. 375–472.
[14] D. A. Cox, Primes of the form x2 + ny 2 , Wiley, New York, 1989.
[15] S. Ramanujan, Modular equations and approximations to π, Quart. J. Math. 45 (1914),
350–372.
[16] , On certain arithmetical functions, Trans. Cambridge Philos. Soc. 22 (1916),
159–184.
[17] , Collected papers, Cambridge University Press, Cambridge, 1927; reprinted by
Chelsea, New York, 1962; reprinted by the American Math. Soc., Providence, 2000.
[18] , Notebooks (2 volumes), Tata Institute of Fundamental Research, Bombay,
1957.
[19] , The lost notebook and other unpublished papers, Narosa, New Delhi, 1988.
[20] K. Venkatachaliengar, Development of elliptic functions according to Ramanujan,
Madurai Kamaraj University, Madurai, 1988, Tech. Rep. 2.

Bruce C. Berndt, Department of Mathematics, University of Illinois, 1409


West Green Street, Urbana, IL 61801
E-mail address: [email protected]

Heng Huat Chan, Department of Mathematics, National University of Singa-


pore, Kent Ridge, Singapore 119260, Republic of Singapore
E-mail address: [email protected]

You might also like