The idea: trap the circle between two polygons

Nobody can lay a ruler along a circle, but a polygon has straight sides. Around 250 BC Archimedes put a regular polygon inside a circle, with its corners on the circle, and another one outside, with its sides touching the circle. The inside polygon is shorter than the circle and the outside one is longer, so the circle's length is trapped between two numbers he could work out. Take a circle of radius 1: half of it is π long, so half of each perimeter gives a bound on π.

A circle of radius 1 with a regular hexagon inside it and another outside it. The inside hexagon has half-perimeter 3, the outside one 2 times the square root of 3, about 3.4641, so 3 is less than pi, which is less than 3.4641
With six sides: 3 < π < 2√3 ≈ 3.4641. The circle always runs between the two polygons.

One side: sin against tan

Cut a regular n-gon into n slices from the centre. Each slice has the angle 360°/n at the centre; halve it to get θ = 180°/n. For the inside polygon, half a side is sin θ. For the outside polygon, half a side is tan θ, because it stands on the radius at a right angle. Adding up the 2n half-sides and halving:

a = n·sin(180°/n) < π < n·tan(180°/n) = b.

This is the same as sin θ < θ < tan θ, the inequality that also gives sin x / x → 1 on the L'Hôpital's Rule page and appears as areas on Proofs Without Words. On the canvas one side of each polygon is red when it is drawn. (Demo: hexagon gives 3 < π < 2√3)

One slice of the polygons with centre O, half-angle theta, radius 1. The inside half-side is sin theta and the outside half-side, which touches the circle at a right angle to the radius, is tan theta. The arc between them has length theta
Half a side inside is sin θ, the arc is θ, half a side outside is tan θ, and sin θ < θ < tan θ.

Doubling the sides

Press Double n and every arc is cut in half. The new points do two jobs. The inside polygon bends out to them, so each side becomes two shorter sides that together are longer: the inside perimeter grows. The tangent line at each new point cuts a corner off the outside polygon: the outside perimeter shrinks. Both stay on their side of the circle, so the bracket only gets tighter.

Archimedes had no sine table. He worked from the two numbers he already had. With a and b for n sides, the values for 2n sides are

New valueFormulaName6 → 12
outside b′2ab / (a + b)harmonic mean of a and b3.21539
inside a′√(a · b′)geometric mean of a and b′3.10583

Only adding, multiplying, dividing and square roots. Archimedes did the square roots by hand with fractions, always rounding the inside value down and the outside value up so the bounds stayed true. On the canvas each formula appears under the table and flies into the new row. (Demo: Archimedes, hexagon → 96 sides)

How fast the bracket closes

For a small angle, sin θ ≈ θ − θ³/6 and tan θ ≈ θ + θ³/3. So the inside value is short of π by about π³/(6n²) and the outside value is over by about twice that. Both errors fall like 1/n²: each doubling makes the gap about 4 times smaller, which the gap ÷ column shows.

ninside aoutside bgap b − agap ÷
63.000003.464100.46410—
123.105833.215390.109564.24
243.132633.159660.027034.05
483.139353.146090.006744.01
963.141033.142710.001684.00
Bars for the gap between the outside and the inside half-perimeter for 6, 12, 24, 48 and 96 sides: 0.46410, 0.10956, 0.02703, 0.00674 and 0.00168, each about a quarter of the one before
Each doubling cuts the gap by about 4, so every two doublings give roughly one more correct digit.

Two number lines on the canvas show the bracket [a, b] around π: the whole view, and a zoom ×20 that only catches the bracket from 24 sides on. A start from a square or a pentagon is further away at first, but the same factor 4 takes over after a doubling or two. (Demo: the gap shrinks 4× per doubling; Demo: start from a square, 4 → 128; Demo: pentagon, 5 → 80)

Archimedes' 96-gon and later records

Archimedes stopped at 96 sides and wrote the bounds as fractions: 3 10/71 < π < 3 1/7, that is 3.14085 < π < 3.14286. His fractions are a little wider than the exact 96-gon values because he rounded safely at every square root. 22/7 is still a common school value for π. Liu Hui in China ran the same doubling to 3072 sides, Zu Chongzhi found 3.1415926 < π < 3.1415927 in the 5th century, and Ludolph van Ceulen spent years on polygons with 2⁶² sides to get 35 digits. After that, infinite series were far faster.

For a way to estimate π with random dots instead of polygons, see Monte Carlo π: there the error shrinks only like 1/√n. Measuring a curve with chords in general is on Arc Length, and sin and tan as lengths on the Unit Circle page.

What the page leaves out

Archimedes' own fractions and his hand-rounded square roots (such as 265/153 < √3 < 1351/780); the areas of the polygons, which trap π the same way; the proof that the inside perimeter always grows and the outside one always shrinks; Richardson extrapolation, which combines two rows to cancel the 1/n² error; and the infinite series and algorithms that later computed π to trillions of digits.

On the mind map

This page on the Trigonometry & Geometry mind map · 34 pages in 7 branches · open it to see the whole map and test yourself in recall mode.