π and e, and the most beautiful theorem in mathematics Robin Wilson Leonhard Euler ( )

Transcription

1 π and e, and the most beautiful theorem in mathematics Robin Wilson Leonhard Euler ( )

2 The circle number π π = (< 22 / 7 ) is the ratio of the circumference C of a circle to its diameter d: π = C/d, so C = π d = 2 π r (r = radius) the ratio of the area A of a circle to the square of its radius r: π = A/r 2, so A = π r 2 = π d 2 /4

3 π to 500 decimal places

4 π to 2000 decimal places

5 The Vienna metro (Karlplatz)

6 From 2πr to πr 2 Sato Moshun s Tengen Shinan, 1698

7 Mesopotamian value of π (c.1800 BC) Ratio of the perimeter of the hexagon to the circumference of the circle is 0; 57, 36 6r / 2πr = 3 / π = 57 / / 3600 So π = 3 1 / 8 = 3.125

8 An Egyptian problem in geometry Problem 50. Example of a round field of diameter 9 khet. What is its area? Answer: Take away 1 / 9 of the diameter, namely 1; the remainder is 8. Multiply 8 times 8; it makes 64. So it contains 64 setat of land. Area = (d 1 / 9 d) 2 = ( 8 / 9 d) 2 = 256 / 81 r 2 which is about r 2

9 Using polygons (Antiphon & Bryson) Antiphon: π > 2 π > Bryson: π < 4 π < 3.32

10 Archimedes value of π perimeter/area of inside 6-gon < circumference/area of circle < perimeter/area of outside 6-gon Now double the number of sides: 6, 12, 24, 48, 96 Result: 3 10 / 71 < π < 3 1 /

11 Chinese values for π Zhang Heng (AD 100) π = 10 ( 3.162) Liu Hui (AD 263) π = (3072 sides) Zu Chongzhi (AD 500) π = (24,576 sides) and π = 355 / 113

12 More polygons Fibonacci (Italy, 1220) π = (96 sides) Al-Kashi (Samarkand, 1424) 16 d. p. ( sides) François Viète (French, 1579) 9 d. p. ( sides) Valentin Otho (Germany, 1573) π = 355 / 113 Adriaan Anthonisz (Dutch, 1585) 333 / 106 < π < 377 / 120 and 355 / 113 (10 9 sides) Adriaan van Roomen (Dutch, 1593) 15 d. p. (2 30 sides)

13 Ludolph van Ceulen (Dutch) ( ) 1596: 20 d. p. 515,396,075,520 = sides 1610: 35 d. p. 4,611,686,018,427,387,904 = 2 62 sides

14 François Viète & John Wallis 1579: Viète 2/π = cos π/4 cos π/8 cos π/16 cos π /32... = ½ 2 ½ (2 + 2) ½ (2 + (2 + 2)) ½ (2 + (2 + (2 + 2))) : Wallis 4 = π

15 The tan 1 (or arctan) function tan θ = a/b θ = tan -1 (a/b) tan π/4 = 1, so tan -1 (1) = π/4 tan π/6 = 1/ 3, so tan -1 (1/ 3) = π/6 tan -1 ( 1 / 2 ) + tan -1 ( 1 / 3 ) = π/4 In general: tan 1 x + tan 1 y = tan 1 x + y 1 xy

16 The series for tan 1 x tan 1 x = x 1 / 3 x / 5 x 5 1 / 7 x So, with x = 1 we have: π/4 = 1 1 / / 5 1 / But this converges extremely slowly: 300 terms give only two decimal places of π. Much better is: π/4 = tan 1 ( 1 / 2 ) + tan 1 ( 1 / 3 ) = { 1 / 2 1 / 3 ( 1 / 2 ) / 5 ( 1 / 2 ) 5 1 / 7 ( 1 / 2 ) } { 1 / 3 1 / 3 ( 1 / 3 ) / 5 ( 1 / 3 ) 5 1 / 7 ( 1 / 3 ) }, which converges much faster.

17 Machin s tan 1 formula (1706) π = 16 tan 1 ( 1 / 5 ) 4 tan 1 ( 1 / 239 ) = 16 { 1 / 5 1 / 3 ( 1 / 5 ) / 5 ( 1 / 5 ) 5 1 / 7 ( 1 / 7 ) } 4 {( 1 / 239 ) 1 / 3 ( 1 / 239 ) 3 + ( 1 / 239 ) 5 1 / 7 ( 1 / 239 ) } Machin used this series to calculate π to 100 decimal places John Machin was Gresham Professor of Astronomy from 1713 to 1751

18 William Jones introduces π (1706)

19 A new value for π Other formulas included one of Leonhard Euler (1755): π = 20 tan 1 ( 1 / 7 ) + 8 tan 1 ( 3 / 79 ). This was later used by Jurij Vega to find π to 136 decimal places, for many years the best value known. But persistent rumours pointed to an earlier value correct to 152 decimal places, hidden in Oxford s Bodleian library. This remained undiscovered until 2014 when Benjamin Wardhaugh found it an an obscure book of It used 314 terms of the tan 1 series in the form π = 6 tan 1 ( 1 / 3 ) = 12 {1 1 / 3 ( 1 / 3 ) + 1 / 5 ( 1 / 3 ) 2 1 / 7 ( 1 / 3 ) }

20 William Shanks s value of π In 1873, using Machin s formula, William Shanks calculated π to 707 decimal places, which can be viewed in the Palais de la Découverte in Paris: Unfortunately,...

21 Some weird results Ramanujan (1914): The Chudnovsky brothers (1989): Bailey, Borwein and Plouffe (1995):

22 Enter the computer 1949 ENIAC (Maryland) hours 1955 NORC minutes 1957 PEGASUS (England) 10, hours 1958 IBM704 (Paris) 10, minutes 1961 IBM7090 (London) 20, minutes 1961 Shanks/Wrench (New York) 100, hours 1973 CDC million 23.3 hours 1989 Chudnovskys (New York) 1 billion 2002 Kanada 1 trillion 2011 Shigeni Kondo etc. 10 trillion...

23 Circling the earth The earth s circumference is about 25,000 miles (132 million feet) Tie a string tightly around the earth. Then extend the string by just 2π ( 6.3) feet, and prop it up equally all around the earth. How high above the ground is the string?

24 Circling the earth The earth s circumference is about 25,000 miles. Tie a string tightly around the earth, extend the string by 2π ( 6.3) feet, and prop it up equally all around the earth. How high above the ground is the string? If the earth s radius is r, then the original string has length 2πr. When we extend the string by 2π, the new circumference is 2πr + 2π = 2π (r + 1), so the new radius is r + 1 so everywhere the string is now one foot off the ground.

25 The exponential number e e = The letter e was first used in an unpublished paper of Euler around 1727, and in a letter of It first appeared in print in 1736 in Volume 1 of his Mechanica: where e denotes the number whose hyperbolic logarithm is 1.

26 A chessboard problem The wealthy king of a certain country was so impressed by the new game of chess that he offered the wise man who invented it any reward he wished. The wise man replied: My prize is for you to give me 1 grain of wheat for the first square of the chessboard, 2 grains for the second square, 4 grains for the third square, and so on, doubling the number of grains on each successive square until the chessboard is filled.

27 A chessboard problem Total number of grains: = = 18,446,744,073,709,551,615 enough wheat to form a pile the size of Mount Everest.

28 Exponential growth polynomial growth n : 1, 2, 3, 4, 5,... n 2 : 1, 4, 9, 16, 25,... n 3 : 1, 8, 27, 64, 125,... exponential growth 2 n : 1, 2, 4, 8, 16, 32,... 3 n : 1, 3, 9, 27, 81, 243,... n = 10 n = 30 n = 50 polynomial n seconds seconds seconds n seconds seconds seconds n seconds seconds seconds n seconds 24.3 seconds 5.2 minutes exponential 2 n seconds 17.9 minutes 35.7 years 3 n seconds 6.5 years years In the long run, exponential growth usually exceeds polynomial growth.

29 An interest-ing problem Jakob Bernoulli asked: How much is earned when compound interest is calculated yearly? twice a year?..., every day? n times a year?..., continuously? Invest 1 at 100% interest per year: twice a year: 1 (1 + 1 / 2 ) = 1.50 (1 + 1 / 2 ) 2 = 2.25 n times a year: 1 (1 + 1 / n ) (1 + 1 / n ) 2... (1 + 1 / n ) n As n increases, these numbers tend to the limiting value of e.

30 Leonhard Euler ( )

31 Some properties of e e = lim n (1 + 1/n) n e x = lim n (1 + x/n) n e = / 1! + 1 / 2! + 1 / 3! + 1 / 4! +... e x = 1 + x/1! + x 2 /2! + x 3 /3! +... The slope of the curve y = e x at the value x is y = e x : dy/dx = y

32 John Napier s logarithms (1614) Aim: to turn multiplication and division problems into addition and subtraction ones that is, to turn geometric progressions into arithmetic ones. log (a b) = log a + log b log 1 log 1 = 161,180,956

33 Henry Briggs s logarithms (1617) [Napier has] set my Head and hands a Work with his new and remarkable logarithms... I never saw a book which pleased me better or made me more wonder. I myself, when expounding this doctrine to my auditors at Gresham College, remarked that it would be much more convenient that 0 should be kept for the logarithm of the whole sine [log 1 = 0] Briggs introduced logs to base 10, which became widely used (for example, by mariners) and for which: log (a b) = log a + log b

34 exp and log are inverse to each other so If y = e x, then x = log y log e x = x and e log y = y So the exponential law e x e y = e x+y and the logarithm law log (x y) = log x + log y are essentially the same result.

35 Derangements Given any n letters a, b, c, d, e,..., in how many ways can we rearrange them so that none is in its original position? n = 4: badc, bcda, bdac, cadb, cdab, cdba, dabc, dcab, dcba n D n ,833 De Moivre/Euler: D n = n! {1 1/1! + 1/2! 1/3! +... ± 1/n!}, which is very close to n! e 1 = n!/e. For example, when n = 8, n!/e = 14,832.9 and D n = 14,833.

36 Exponential growth as predicted by Thomas Malthus (1796) If N(t), the size of a population at time t, grows at a fixed rate k proportional to its size, then dn/dt = kn or dn/n = k dt. This can be integrated to give log N = kt + const. or N = N 0 e kt, where N 0 is the initial population. Similarly for the decay of radium, the cooling of a cup of tea, etc.

37 Euler s equation the most beautiful theorem in mathematics Five important constants: 1 the counting number 0 the nothingness number π the circle number e the exponential number i the imaginary number e iπ + 1 = 0 (or e iπ = 1) Feynman (aged 14): the most remarkable formula in math M. Atiyah: the mathematical equivalent of Hamlet s To be or not to be: very succinct, but at the same time very deep Featured in The Simpsons and in a criminal court case

38 A near-miss: Johann Bernoulli What is log ( 1)? 2 log ( 1) = log ( 1) + log ( 1) = log ( 1 1) = log (1) = 0, so log ( 1) = 0 ; similarly, log ( x) = log x, for all x. Bernoulli also found the area A of the sector on the right to be A = (a 2 /4i) log (x + iy)/(x iy). Euler then put x = 0: A = (a 2 /4i) log ( 1) = πa 2 /4. So log ( 1) 0, and we deduce log ( 1) = iπ Taking exponentials then gives Euler s equation in the form e iπ = 1

39 Euler s identity There s no real reason why there should be any connection between e x and cos x and sin x. e x = 1 + x/1! + x 2 /2! + x 3 /3! + x 4 /4! + x 5 /5! +... cos x = 1 x 2 /2! + x 4 /4!... sin x = x x 3 /3! + x 5 /5!... so e ix = 1 + ix/1! + (ix) 2 /2! + (ix) 3 /3! + (ix) 4 /4! + (ix) 5 /5! +... = (1 x 2 /2! + x 4 /4!...) + i (x x 3 /3! + x 5 /5!... ) = cos x + i sin x e ix = cos x + i sin x

40 Euler s identity: e iv = cos v + i sin v

41 Two postage stamps

42 Another near-miss: Roger Cotes Roger Cotes introduced radian measure for angles, and worked with Newton on the Principia Mathematica. He also investigated the surface area of an ellipsoid, and found two expressions for it involving logarithms and trigonometry, and both involving an angle φ. Equating these, he found that log (cos φ + i sin φ) = iφ. Taking exponentials then gives Euler s identity

43 Consequences of Euler s identity e ix = cos x + i sin x 1. Put x = π (= 180 ): then e iπ = cos π + i sin π = i = 1, which is Euler s equation. 2. For any number n, (cos x + i sin x) n = (e ix ) n = e i (nx) = cos nx + i sin nx : this is De Moivre s theorem. 3. Since e ix = cos x + i sin x and e ix = cos x i sin x, cos x = (e ix + e ix )/2 and sin x = (e ix e ix )/2i.

44 Who discovered Euler s equation? Bernoulli/Euler: iπ = log ( 1) Roger Cotes: log (cos φ + i sin φ) = iφ Euler seems never to have written down e iπ + 1 = 0 explicitly though he surely realised that it follows from his identity, e ix = cos x + i sin x. We don t know who first stated it explicitly though there s an early appearance in Most people attribute it to Euler, to honour this great mathematical pioneer.

45 Who discovered Euler s equation? Bernoulli/Euler: iπ = log ( 1) Roger Cotes: log (cos φ + i sin φ) = iφ Euler seems never to have written down e iπ + 1 = 0 explicitly though he surely realised that it follows from his identity, e ix = cos x + i sin x. We don t know who first stated it explicitly though there s an early appearance in Most people attribute it to Euler, to honour this great mathematical pi-0ne-e-r.

46 Euler s pioneering equation (Oxford University Press) Gresham College book launch at Barnard s Inn Hall on 15 February 2018 at 6 pm