Binet's Formula

Nov 5 2022

My last post talked a lot about the Golden Ratio and how we can see and derive it from many different sources. This post will focus on a special formula for finding the nth Fibonacci number that involves the Golden Ratio, named Binet's formula for Jacques Philippe Marie Binet, who derived the formula (although sources indicate that Abraham de Moivre and Daniel Bernoulli had known about it beforehand).

To start, let's express the formula itself:

F(n)=φnφˉnφφˉ=φnφˉn5F(n) = \frac{\varphi^n - \bar{\varphi}^n}{\varphi - \bar{\varphi}} = \frac{\varphi^n - \bar{\varphi}^n}{\sqrt{5}}

where φ\varphi is the Golden Ratio, and φˉ\bar{\varphi} is its conjugate. They are both solutions to the equation

x2x1=0x^2 - x - 1 = 0

To make things clearer, I want to write the polynomial in a different way:

x2x1x0=0x^2 - x^1 - x^0 = 0

We can multiply any amount of xxs to the equation, and it will still preserve equality (0 times anything is 0, while the relative powers of the polynomial are preserved). This means that we can rewrite the quadratic as

x3x2x1=0x^3 - x^2 - x^1 = 0

x4x3x2=0x^4 - x^3 - x^2 = 0

\vdots

xnxn1xn2=0x^n - x^{n - 1} - x^{n - 2} = 0

or:

xn=xn1+xn2x^n = x^{n - 1} + x^{n - 2}

Since we haven't changed the validity of the equation, the solutions to the above polynomial are still φ\varphi and φˉ\bar{\varphi}. This means that

φn=φn1+φn2\varphi^n = \varphi^{n-1} + \varphi^{n-2}

φˉn=φˉn1+φˉn2\bar{\varphi}^n = \bar{\varphi}^{n-1} + \bar{\varphi}^{n-2}

This is meaningful in that the powers of φ\varphi and φˉ\bar{\varphi} themselves satisfy the Fibonacci recurrence. So for any sequence PP defined as

Pn=aφn+bφˉnP_n = a\varphi^n + b\bar{\varphi}^n

also satisfies the Fibonacci recurrence:

Pn=a(φn1+φn2)+b(φˉn1+φˉn2)P_n = a(\varphi^{n-1} + \varphi^{n-2}) + b(\bar{\varphi}^{n-1} + \bar{\varphi}^{n-2})

Pn=aφn1+bφˉn1+aφn2+bφˉn2=Pn1+Pn2P_n = a\varphi^{n-1} + b\bar{\varphi}^{n-1} + a\varphi^{n-2} + b\bar{\varphi}^{n-2} = P_{n-1} + P_{n-2}

So, we can choose some aa and bb such that it follows the Fibonacci sequence: namely,

a+b=0a + b = 0

aφ+bφˉ=1a\varphi + b\bar{\varphi} = 1

solutions of which make sure that aa and bb fulfill the first two terms of the sequence. With some manipulation, we find that

a=1φφˉ=15a=\frac{1}{\varphi - \bar{\varphi}} = \frac{1}{\sqrt{5}}

and

b=a=15b = -a = -\frac{1}{\sqrt{5}}

So, with our arbitrary sequence PnP_n, we can plug our a, ba,\ b back in to find a formula for the Fibonacci sequence itself:

Pn=aφn+bφˉn=a(φnφˉn)P_n = a\varphi^n + b\bar{\varphi}^n = a(\varphi^n - \bar{\varphi}^n)

Fn=φnφˉn5F_n = \frac{\varphi^n - \bar{\varphi}^n}{\sqrt{5}}

Notice how we can use any value for a, ba,\ b here to find the nth term of any sequence that behaves like the Fibonacci sequence but with different starting numbers. All we'd have to do is find an a, ba,\ b that satisfies the system of equations

a+b=P0a + b = P_0

aφ+bφˉ=P1a\varphi + b\bar{\varphi} = P_1

We do not even have to be limited to Fibonacci sequences. Since our understanding of this formula comes from the quadratic

x2x1=0x^2 - x - 1 = 0

why shouldn't there exist some formula that would tell us the nth Pell number, or for that matter, the nth term of a given Lucas sequence?

We can start with the Pell numbers, of which growth is analogous to the silver ratio as the Fibonacci numbers are to the golden ratio. Once again, we can start with our 'silver ratio' quadratic:

x22x1=0x^2 - 2x - 1 = 0.

Once again, simply multiplying by xx throughout yields

xn=2xn1+xn2x^n = 2x^{n-1} + x^{n-2}

For the silver ratio δS\delta_S and its conjugate δSˉ\bar{\delta_S}, then,

δSn=2δSn1+δSn2\delta_S^n = 2\delta_S^{n-1} + \delta_S^{n-2} δSˉn=2δSˉn1+δSˉn2\bar{\delta_S}^n = 2\bar{\delta_S}^{n-1} + \bar{\delta_S}^{n-2}

Once again, any sequence satisfied by the Pell relation will also satisfy

Pn=aδSn+bδSˉnP_n = a\delta_S^n + b\bar{\delta_S}^n

Thus, we need to solve the system of equations

a+b=0a + b = 0 aδS+bδSˉ=1a\delta_S + b\bar{\delta_S} = 1

Which yields

a=1δSδSˉ=18=ba = \frac{1}{\delta_S - \bar{\delta_S}} = \frac{1}{\sqrt{8}} = -b

and we are left with the formula

Pn=δSnδSˉn22P_n = \frac{\delta_S^n - \bar{\delta_S}^n}{2\sqrt{2}}

For any Lucas sequence of the form

Pn=aPn1+bPn2P_n = aP_{n-1} + bP_{n-2}

and a corresponding quadratic equation

x2axb=0x^2 - ax - b = 0

we can derive a similar formula for the nth term.