Certain quadratics: The golden ratio (sqrt(5)+1)/2 , and sqrt(2)

Dan Brown's bestseller The Da Vinci Code features one of the most famous algebraic numbers (without, of course, being referred to as such!); unfortunately his presentation contains a number of serious errors. I quote from his Chapter 20:

" [Langdon] felt himself suddenly reeling back to Harvard, standing in front of his 'Symbolism in Art' class, writing his favourite number on the chalkboard: 1.618
Langdon turned to face his sea of eager students. 'Who can tell me what this number is?
A long-legged maths major at the back of the class raised his hand. 'That's the number
PHI.' He pronounced it fee .
'Nice job, Stettner,' Langdon said. 'Everyone, meet PHI.'

... ...
'This number PHI,' Langdon continued, 'one-point-six-one-eight, is a very important number in art. Who can tell me why?'
... ...
'Actually,' Langdon said, ... 'PHI is generally considered the most beautiful number in the universe.'
... ...
As Langdon loaded his slide projector, he explained that the number PHI was derived from the Fibonacci sequence - a progression famous not only because the sum of adjacent terms equalled the next term, but because the
quotients [Brown's italics] of adjacent terms possessed the astonishing property of approaching the number 1.618 - PHI! " [end of quote]

' phi ' (PHI), the so-called golden ratio , is the number (sqrt(5)+1)/2 , and is one of the two solutions of the quadratic equation x^2-x-1 = 0 . This number is so well known that it would be foolish for me to write about it at length; however I will make a few comments.

First, phi is NOT "1.618". It is entirely elementary to argue that (sqrt(5)+1)/2 cannot have a terminating decimal expansion, and not only that, but being irrational it cannot have a periodic decimal expansion. (Of course I can see that Brown could have put some readers off their coffee had he used (sqrt(5)+1)/2 ... , but surely he could have been advised to use '1.618... ' What view could readers form of mathematicians who would 'consider' 1.618 to be 'the most beautiful number in the universe'?)

> solve(x^2 - x - 1 = 0);
alpha := (sqrt(5) + 1)/2;
evalf(alpha);

1/2*sqrt(5)+1/2, 1/2-1/2*sqrt(5)

alpha := 1/2*sqrt(5)+1/2

1.618033989

>

It is well-known that phi is the limit of the infinite sequence { F[n+1]/F[n] }, where F[n] is the n -th Fibonacci number (they, too, make an appearance in the Da Vinci Code).

The Fibonacci numbers are the terms of the sequence 1, 1, 2, 3, 5, 8, 13, ... , F[n] , ... , where F[1] = 1, F[2] = 1 and F[n] = F[n-1]+F[n-2] for 3 <= n . Here are the first 25 of them:

> restart;

> F[1] := 1: F[2] := 1: for k from 3 to 25 do
F[k] := F[k-1]+F[k-2] od:
seq(F[k], k = 1..25);

1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233, 377...
1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144, 233, 377...

>

It is an elementary exercise for a beginner to prove that the sequence F[2]/F[1], F[4]/F[3], F[6]/F[5], F[8]/F[7], `...`, ... is strictly monotonic in creasing (early terms followed by successive differences):

> seq(F[2*n]/F[2*n-1], n = 1..12);

1, 3/2, 8/5, 21/13, 55/34, 144/89, 377/233, 987/610...

> seq(F[2*n+2]/F[2*n+1] - F[2*n]/F[2*n-1], n = 1..11);

1/2, 1/10, 1/65, 1/442, 1/3026, 1/20737, 1/142130, ...

>

while the sequence F[3]/F[2], F[5]/F[4], F[7]/F[6], F[9]/F[8], `...`, ... is strictly monotonic de creasing:

> seq(F[2*n+1]/F[2*n], n = 1..12);

2, 5/3, 13/8, 34/21, 89/55, 233/144, 610/377, 1597/...

> seq(F[2*n+3]/F[2*n+2] - F[2*n+1]/F[2*n], n = 1..11);

-1/3, -1/24, -1/168, -1/1155, -1/7920, -1/54288, -1...

>

Also, the n- th term of the first sequence is less that the n- th term of the second sequence, and the difference between those tends to zero as proc (n) options operator, arrow; infinity end proc... :

> seq(F[2*n+1]/F[2*n] - F[2*n]/F[2*n-1], n = 1..11);

1, 1/6, 1/40, 1/273, 1/1870, 1/12816, 1/87841, 1/60...

>

Thus Limit(F[n+1]/F[n],n = infinity) exists, and is easily shown to be (sqrt(5)+1)/2 , namely ' phi '. In fact, every term of the sequence { F[n+1]/F[n] }, from the 12-th onwards, is greater than Brown's 1.618 'limit':

> seq(evalf(F[k+1]/F[k]), k = 1..24);

1., 2., 1.500000000, 1.666666667, 1.600000000, 1.62...
1., 2., 1.500000000, 1.666666667, 1.600000000, 1.62...
1., 2., 1.500000000, 1.666666667, 1.600000000, 1.62...
1., 2., 1.500000000, 1.666666667, 1.600000000, 1.62...

>

In short, just as Pi is not 3.14 (nor 22/7 ), phi is not 1.618.

Contact details 

After August 31st 2007 please use the following Gmail address: jbcosgrave at gmail.com


This page was last updated 18 February 2005 15:09:55 -0000