## 1  The geometric period

Up to the Seventeenth Century, approximations of p were obtained by mean of geometrical considerations. Most of the methods were dealing with regular polygons circumscribed about and inscribed in the circle. The perimeter or the area of those polygons were calculated with elementary geometrical rules.

During this period the notation p was not used and it was not yet a constant but just a geometrical ratio or even just implicit.

### 1.1  Ancient estimations

#### 1.1.1  Egypt

In one of the oldest mathematical text, the Rhind papyrus (from the name of the Egyptologist Henry Rhind who purchased this document in 1858 at Luxor), the scribe Ahmes copied, around 1650 B.C.E., eighty-five mathematical problems. Among those is given a rule, the problem 48, to find the area of a circular field of diameter 9: take away 1/9 of the diameter and take the square of the remainder. In modern notation, it becomes
 A= ćč d- 1 9 d öř 2 = ćč 8 9 öř 2 d2,
(A is the area of the field and d it's diameter): so if we use the formula A=pd2/4, comes the following approximation
 p = 4 ćč 8 9 öř 2 = ćč 4 3 öř 4 » 3.1605.

This accuracy is astonishing for such ancient time. See  for a possible justification of this value.

#### 1.1.2  Babylon

On a Babylonian cuneiform tablet from Susa, about 2000 B.C.E., and discovered in 1936, the ratio of the perimeter of the circle to its diameter was founded to be
 p = 3+ 1 8 =3.125,
and this estimation is one of the oldest we know.

### 1.2  Archimedes' method

The famous treatise On the measurement of the Circle from the Greek mathematician and engineer Archimedes of Syracuse (287-212 B.C.E.) is a major step in the knowledge of the circle properties. Among those are given the following numerical bounds for p
 3+ 10 71
 <
 p < 3+ 1 7 = 22 7
 3.14085
 <
 p < 3.142857,
and, that is a novelty, it's the first algorithm which allows, in theory, to compute as many digits of p as required.

The proof uses a recursive algorithm starting with two regular hexagons (see figure 1), one is circumscribed and the other is inscribed. Archimedes then subdivides both polygons by 2, so the number of sides of the polygons is 6.2n (respectively 6,12,24,48,96 sides in his computations). The perimeters of the last polygons gives the announced bounds. Figure 1: Archimedes's hexagons

In modern notations this gives, for the circumscribed polygons, the following recursive sequence (where un is half the length of the side of the polygon and vn the distance from a vertex of the polygon to the center of the circle)
 un+1
 =
 un 1+vn
 vn+1
 =
 Ö 1+un+12
 pn
 =
 6.2nun
starting with the regular hexagon for which:u0=1/Ö3,v0=2/Ö3,p0=2Ö3.

We may observe that Archimedes needed to compute square roots to achieve his computation, and in particular he used a good approximation of Ö3 » 265/153 (this value is one of the partial quotients of the continued fraction of Ö3). The upper bound founded by Archimedes is given with n=4.

Now we observe the value of the first iterates of this sequence:
 p0
 =
 3.(46410161513...)
 p1
 =
 3.(21539030917...)
 p2
 =
 3.1(5965994209...)
 p3
 =
 3.14(608621513...)
 p4
 =
 3.14(271459964...)
 p10
 =
 3.141592(92738...)
 p16
 =
 3.141592653(65...)

The same kind of iteration may be written for the inscribed polygons. Archimedes method is the first known algorithm, in theory, to compute p at a desired accuracy.

For example, p16 is a nine digits approximation and it is the perimeter of a polygon of 393216 sides. This approximation was given much later than Archimedes by François Vičte (1540-1603) after a long calculation. He gave the bounds :

 3.1415926535 < p < 3.1415926537.

#### 1.2.1  Trigonometric formulation

Using trigonometric functions (unknown to Archimedes), it's easy to show that the perimeter of the circumscribed polygon with n sides on a circle of radius r is given by
 Ln=2nrtan ćč p n öř ,
and the perimeter of the inscribed polygon is
 ln=2nrsin ćč p n öř .

It follows the bounds for p
 nsin ćč p n öř < p < ntan ćč p n öř
and, for example, with n=6 (hexagon)
 3 < p < 2Ö3.

If we observe that, when n is large,
 nsin ćč p n öř
 =
 p- p3 6n2 +O ćč 1 n4 öř
 ntan ćč p n öř
 =
 p+ p3 3n2 +O ćč 1 n4 öř ,
this suggests to form (in order to eliminate the term in 1/n2)
 2 3 nsin ćč p n öř + 1 3 ntan ćč p n öř =p+O ćč 1 n4 öř
whose convergence is much faster than Archimedes' method. This was noted by a geometric mean, in 1621, by Willebrod Snellius (1580-1626) in his work Cyclometricus. For example with the hexagons (n=6), it becomes
 p » 2+ 2 3 Ö3=3.1(547...),
and with the last polygons considered by Archimedes (n=96)
 p » 3.141592(833...).

See  for other accelerations of this method.

### 1.3  Pfaff formulation

In 1800, Johann Friedrich Pfaff (1765-1825) gave a modern and analytical formulation of the previous geometric iteration. Let an and bn be respectively the length of the circumscribed and inscribed regular polygon with 6.2n sides, then starting with a0=2Ö3,b0=3 and
 an+1
 =
 2anbn an+bn
 bn+1
 =
 Ö an+1bn
we have
 bn
 <
 bn+1 < an+1 < an
 p
 =
 lim n® Ą an= lim n® Ą bn
so that p is the common limit of the two sequences (an,bn), and the error is O(4-n).

If we set an=1/un and bn=1/vn the iteration becomes (starting with a0=1/(2Ö3),b0=1/3)
 an+1
 =
 an+bn 2
 bn+1
 =
 Ö an+1bn
and
 1 p = lim n® Ą an= lim n®Ą bn.

It's interesting to note that this formulation is not very far from the AGM iteration which is one of the recent and very efficient methods used to compute many digits of p.

### 1.4  Record of computation during the geometric period

Here are some other approximations computed by various mathematicians and mostly by mean of polygons:

• 150 : Claudius Ptolemy (c. 85-165, Egypt) published 3.141666... in his Almagest (an astronomical treatise), the value was given in sexagesimal fractions 3+8/60+30/602.
• 263 : Liu Hui (China) gave 3927/1250=3.1416 with a polygon of 3072 sides .
• 480? : Zu Chongzhi (429-500, China) determined 355/113 and also the very impressive bounds 3.1415926 < p < 3.1415927 probably by Liu's method .
• 499 : Aryabhata (c. 476-550, India) gave 3.1416 maybe by mean of a polygon of 384 sides. In fact, in rule 10 of the Aryabhatiya we are told: add four to one hundred, multiply by eight and then add sixty-two thousand. The result is approximately the circumference of a circle of diameter twenty thousand.
• 830 : Al'Khwarizmi (c. 780-850, Persia): 22/7,Ö{10} and 62832/20000. The term algorithm came from his name.
• 1220 : Leonardo of Pisa (1180-1240, Italy): 3.141818. He is better known as Fibonacci.
• 1424 : Al-Kashi (c. 1380-1429, Samarkand): 14 digits with a polygon of 6.227sides in his Treatise on the Circumference.
• 1579 : François Vičte (1540-1603, France): 9 digits with a polygon of 393216 sides.
• 1593 : Adrianus Romanus (1561-1615, Netherlands): 15 digits with a polygon of 230 sides.
• 1596 : Ludolph van Ceulen (1540-1610, Germany): 20 digits with a polygon of 60.233 sides.
• 1609 : Ludolph van Ceulen: 35 digits with a polygon of 262 sides. He spent a considerable part of his life with such computations. p was often known as the Ludolphine Number in Germany and his digits are engraved on his tombstone.
• 1621 : Willebrod Snellius (1580-1626, Netherlands): 34 digits with a polygon of only 230 sides thanks to an acceleration of Archimedes' method.

A very complete outline of the history of p is given in  and a more modern one are in  and .

## References


Le Petit Archimčde, no. hors série, Le nombre p, (1980)


J. Arndt and C. Haenel, p- Unleashed, Springer, (2001)


J.P. Delahaye, Le fascinant nombre p, Bibliothčque Pour la Science, Belin, (1997)


H. Engels, Quadrature of the Circle in Ancient Egypt, Historia Mathematica, (1977), vol. 4, p. 137-140


T.L. Heath, The Works of Archimedes, Cambridge University Press, (1897)


L.Y. Lam and T.S. Ang, Circle Measurements in Ancient China, Historia Mathematica, (1986), vol. 13, p. 325-340


G.M. Phillips, Archimedes the Numerical Analyst, The American Mathematical Monthly, (1981), vol. 88, p. 165-169


H.C. Schepler, The Chronology of Pi, Mathematics Magazine, (1950)


F. Vičte, Opera Mathematica (reprinted), Georg Olms Verlag, H