# Fibonacci number

(Redirected from Fibonacci Sequence)
A tiling with squares whose side lengths are successive Fibonacci numbers

In mathematics, the Fibonacci numbers or Fibonacci sequence are the numbers in the following integer sequence:[1][2]

$1,\;1,\;2,\;3,\;5,\;8,\;13,\;21,\;34,\;55,\;89,\;144,\; \ldots\;$

or (often, in modern usage):

$0,\;1,\;1,\;2,\;3,\;5,\;8,\;13,\;21,\;34,\;55,\;89,\;144,\; \ldots\;$ (sequence A000045 in OEIS).
The Fibonacci spiral: an approximation of the golden spiral created by drawing circular arcs connecting the opposite corners of squares in the Fibonacci tiling;[3] this one uses squares of sizes 1, 1, 2, 3, 5, 8, 13, 21, and 34.

By definition, the first two numbers in the Fibonacci sequence are 1 and 1, or 0 and 1, depending on the chosen starting point of the sequence, and each subsequent number is the sum of the previous two.

In mathematical terms, the sequence Fn of Fibonacci numbers is defined by the recurrence relation

$F_n = F_{n-1} + F_{n-2},\!\,$

with seed values[1][2]

$F_1 = 1,\; F_2 = 1$

or[4]

$F_0 = 0,\; F_1 = 1.$

The Fibonacci sequence is named after Fibonacci. His 1202 book Liber Abaci introduced the sequence to Western European mathematics,[5] although the sequence had been described earlier in Indian mathematics.[6][7][8] By modern convention, the sequence begins either with F0 = 0 or with F1 = 1. The Liber Abaci began the sequence with F1 = 1, without an initial 0.

Fibonacci numbers are closely related to Lucas numbers in that they are a complementary pair of Lucas sequences. They are intimately connected with the golden ratio; for example, the closest rational approximations to the ratio are 2/1, 3/2, 5/3, 8/5, ... . Applications include computer algorithms such as the Fibonacci search technique and the Fibonacci heap data structure, and graphs called Fibonacci cubes used for interconnecting parallel and distributed systems. They also appear in biological settings,[9] such as branching in trees, phyllotaxis (the arrangement of leaves on a stem), the fruit sprouts of a pineapple,[10] the flowering of an artichoke, an uncurling fern and the arrangement of a pine cone.[11]

## Origins

A page of Fibonacci's Liber Abaci from the Biblioteca Nazionale di Firenze showing (in box on right) the Fibonacci sequence with the position in the sequence labeled in Latin and Roman numerals and the value in Hindu-Arabic numerals.

The Fibonacci sequence appears in Indian mathematics, in connection with Sanskrit prosody.[7][12] In the Sanskrit tradition of prosody, there was interest in enumerating all patterns of long (L) syllables that are 2 units of duration, and short (S) syllables that are 1 unit of duration; counting the different patterns of L and S of a given duration results in the Fibonacci numbers: the number of patterns that are m short syllables long is the Fibonacci number Fm + 1.[8]

Susantha Goonatilake writes that the development of the Fibonacci sequence "is attributed in part to Pingala (200 BC), later being associated with Virahanka (c. 700 AD), Gopāla (c. 1135), and Hemachandra (c. 1150)".[6] Parmanand Singh cites Pingala's cryptic formula misrau cha ("the two are mixed") and cites scholars who interpret it in context as saying that the cases for m beats (Fm+1) is obtained by adding a [S] to Fm cases and [L] to the Fm−1 cases. He dates Pingala before 450 BC.[13]

However, the clearest exposition of the series arises in the work of Virahanka (c. 700 AD), whose own work is lost, but is available in a quotation by Gopala (c. 1135):

Variations of two earlier meters [is the variation]... For example, for [a meter of length] four, variations of meters of two [and] three being mixed, five happens. [works out examples 8, 13, 21]... In this way, the process should be followed in all mātrā-vṛttas [prosodic combinations].[14]

The series is also discussed by Gopala (before 1135 AD) and by the Jain scholar Hemachandra (c. 1150).

In the West, the Fibonacci sequence first appears in the book Liber Abaci (1202) by Leonardo of Pisa, known as Fibonacci.[5] Fibonacci considers the growth of an idealized (biologically unrealistic) rabbit population, assuming that: a newly born pair of rabbits, one male, one female, are put in a field; rabbits are able to mate at the age of one month so that at the end of its second month a female can produce another pair of rabbits; rabbits never die and a mating pair always produces one new pair (one male, one female) every month from the second month on. The puzzle that Fibonacci posed was: how many pairs will there be in one year?

• At the end of the first month, they mate, but there is still only 1 pair.
• At the end of the second month the female produces a new pair, so now there are 2 pairs of rabbits in the field.
• At the end of the third month, the original female produces a second pair, making 3 pairs in all in the field.
• At the end of the fourth month, the original female has produced yet another new pair, the female born two months ago produces her first pair also, making 5 pairs.

At the end of the nth month, the number of pairs of rabbits is equal to the number of new pairs (which is the number of pairs in month n − 2) plus the number of pairs alive last month (n − 1). This is the nth Fibonacci number.[15]

The name "Fibonacci sequence" was first used by the 19th-century number theorist Édouard Lucas.[16]

## List of Fibonacci numbers

The first 21 Fibonacci numbers Fn for n = 0, 1, 2, ..., 20 are:[17]

 F0 F1 F2 F3 F4 F5 F6 F7 F8 F9 F10 F11 F12 F13 F14 F15 F16 F17 F18 F19 F20 0 1 1 2 3 5 8 13 21 34 55 89 144 233 377 610 987 1597 2584 4181 6765

The sequence can also be extended to negative index n using the re-arranged recurrence relation

$F_{n-2} = F_n - F_{n-1},$

which yields the sequence of "negafibonacci" numbers[18] satisfying

$F_{-n} = (-1)^{n+1} F_n.$

Thus the bidirectional sequence is

 F−8 F−7 F−6 F−5 F−4 F−3 F−2 F−1 F0 F1 F2 F3 F4 F5 F6 F7 F8 −21 13 −8 5 −3 2 −1 1 0 1 1 2 3 5 8 13 21

## Use in mathematics

The Fibonacci numbers are the sums of the "shallow" diagonals (shown in red) of Pascal's triangle.

The Fibonacci numbers occur in the sums of "shallow" diagonals in Pascal's triangle (see Binomial coefficient).[19]

$F_{n}=\sum_{k=0}^{\lfloor\frac{n-1}{2}\rfloor} \tbinom {n-k-1} k$

These numbers also give the solution to certain enumerative problems.[20] The most common such problem is that of counting the number of compositions of 1s and 2s that sum to a given total n: there are Fn+1 ways to do this. For example F6 = 8 counts the eight compositions:

1+1+1+1+1 = 1+1+1+2 = 1+1+2+1 = 1+2+1+1 = 2+1+1+1 = 2+2+1 = 2+1+2 = 1+2+2,

all of which sum to 6−1 = 5.

The Fibonacci numbers can be found in different ways among the set of binary strings, or equivalently, among the subsets of a given set.

• The number of binary strings of length n without consecutive 1s is the Fibonacci number Fn+2. For example, out of the 16 binary strings of length 4, there are F6 = 8 without consecutive 1s – they are 0000, 0001, 0010, 0100, 0101, 1000, 1001 and 1010. By symmetry, the number of strings of length n without consecutive 0s is also Fn+2. Equivalently, Fn+2 is the number of subsets S ⊂ {1,...,n} without consecutive integers: {i, i+1} ⊄ S for every i. The symmetric statement is: Fn+2 is the number of subsets S ⊂ {1,...,n} without two consecutive skipped integers: that is, S = {a1 < ... < ak} with ai+1 ≤ ai + 2.
• The number of binary strings of length n without an odd number of consecutive 1s is the Fibonacci number Fn+1. For example, out of the 16 binary strings of length 4, there are F5 = 5 without an odd number of consecutive 1s – they are 0000, 0011, 0110, 1100, 1111. Equivalently, the number of subsets S ⊂ {1,...,n} without an odd number of consecutive integers is Fn+1.
• The number of binary strings of length n without an even number of consecutive 0s or 1s is 2Fn. For example, out of the 16 binary strings of length 4, there are 2F4 = 6 without an even number of consecutive 0s or 1s – they are 0001, 0111, 0101, 1000, 1010, 1110. There is an equivalent statement about subsets.

## Relation to the golden ratio

### Closed-form expression

Like every sequence defined by a linear recurrence with constant coefficients, the Fibonacci numbers have a closed-form solution. It has become known as Binet's formula, even though it was already known by Abraham de Moivre:[21]

$F_n = \frac{\varphi^n-\psi^n}{\varphi-\psi} = \frac{\varphi^n-\psi^n}{\sqrt 5}$

where

$\varphi = \frac{1 + \sqrt{5}}{2} \approx 1.61803\,39887\cdots\,$

is the golden ratio (sequence A001622 in OEIS), and

$\psi = \frac{1 - \sqrt{5}}{2} = 1 - \varphi = - {1 \over \varphi} \approx -0.61803\,39887\cdots$[22]

Since $\psi = -\frac{1}{\varphi}$, this formula can also be written as

$F_n = \frac{\varphi^n-(-\varphi)^{-n}}{\sqrt 5}$

To see this,[23] note that φ and ψ are both solutions of the equations

$x^2 = x + 1,\, x^n = x^{n-1} + x^{n-2},\,$

so the powers of φ and ψ satisfy the Fibonacci recursion. In other words

$\varphi^n = \varphi^{n-1} + \varphi^{n-2}\,$

and

$\psi^n = \psi^{n-1} + \psi^{n-2}\, .$

It follows that for any values a and b, the sequence defined by

$U_n=a \varphi^n + b \psi^n\,$

satisfies the same recurrence

$U_n=a \varphi^{n-1} + b \psi^{n-1} + a \varphi^{n-2} + b \psi^{n-2} = U_{n-1} + U_{n-2}.\,$

If a and b are chosen so that U0 = 0 and U1 = 1 then the resulting sequence Un must be the Fibonacci sequence. This is the same as requiring a and b satisfy the system of equations:

$\left\{\begin{array}{l} a + b = 0\\ \varphi a + \psi b = 1\end{array}\right.$

which has solution

$a = \frac{1}{\varphi-\psi} = \frac{1}{\sqrt 5},\, b = -a$

producing the required formula.

### Computation by rounding

Since

$\frac{|\psi|^n}{\sqrt 5} < \frac{1}{2}$

for all n ≥ 0, the number Fn is the closest integer to

$\frac{\varphi^n}{\sqrt 5}\, .$

Therefore it can be found by rounding, or in terms of the floor function:

$F_n=\bigg\lfloor\frac{\varphi^n}{\sqrt 5} + \frac{1}{2}\bigg\rfloor,\ n \geq 0.$

Or the nearest integer function:

$F_n=\bigg[\frac{\varphi^n}{\sqrt 5}\bigg],\ n \geq 0.$

Similarly, if we already know that the number F > 1 is a Fibonacci number, we can determine its index within the sequence by

$n(F) = \bigg\lfloor \log_\varphi \left(F\cdot\sqrt{5} + \frac{1}{2}\right) \bigg\rfloor$

### Limit of consecutive quotients

Johannes Kepler observed that the ratio of consecutive Fibonacci numbers converges. He wrote that "as 5 is to 8 so is 8 to 13, practically, and as 8 is to 13, so is 13 to 21 almost", and concluded that the limit approaches the golden ratio $\varphi$.[24][25]

$\lim_{n\to\infty}\frac{F_{n+1}}{F_n}=\varphi$

This convergence does not depend on the starting values chosen, excluding 0, 0. For example, the initial values 19 and 31 generate the sequence 19, 31, 50, 81, 131, 212, 343, 555 ... etc. The ratio of consecutive terms in this sequence shows the same convergence towards the golden ratio.

In fact this holds for any sequence that satisfies the Fibonacci recurrence other than a sequence of 0s. This can be derived from Binet's formula.

Another consequence is that the limit of the ratio of two Fibonacci numbers offset by a particular finite deviation in index corresponds to the golden ratio raised by that deviation. Or, in other words:

$\lim_{n\to\infty}\frac{F_{n+\alpha}}{F_n}=\varphi^\alpha$
Animated GIF file showing successive tilings of the plane, and a graph of approximations to the Golden Ratio calculated by dividing successive pairs of Fibonacci numbers, one by the other. Uses the Fibonacci numbers 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, 144

### Decomposition of powers of the golden ratio

Since the golden ratio satisfies the equation

$\varphi^2 = \varphi + 1,\,$

this expression can be used to decompose higher powers $\varphi^n$ as a linear function of lower powers, which in turn can be decomposed all the way down to a linear combination of $\varphi$ and 1. The resulting recurrence relationships yield Fibonacci numbers as the linear coefficients:

$\varphi^n = F_n\varphi + F_{n-1}.$

This equation can be proved by induction on n.

This expression is also true for n < 1 if the Fibonacci sequence Fn is extended to negative integers using the Fibonacci rule $F_n = F_{n-1} + F_{n-2}.$

## Matrix form

A 2-dimensional system of linear difference equations that describes the Fibonacci sequence is

\begin{align} {F_{k+2} \choose F_{k+1}} &= \begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix} {F_{k+1} \choose F_{k}} \\ \vec F_{k+1} &= \mathbf{A} \vec F_{k} ~. \end{align}

The eigenvalues of the matrix A are $\scriptstyle \varphi=\frac12(1+\sqrt5)\,\!$ and $\scriptstyle 1-\varphi=\frac12(1-\sqrt5)$, for the respective eigenvectors $\scriptstyle {\varphi \choose 1}$ and $\scriptstyle {1-\varphi \choose 1}$ .

Since $\scriptstyle {\varphi \choose 1} - \scriptstyle {1-\varphi \choose 1}= \sqrt{5} ~\scriptstyle {1\choose 0}$, and $\vec F_{n} = \mathbf{A}^n \vec F_{0},$ the above closed-form expression for the nth element in the Fibonacci series as an analytic function of n is now read off directly,

$F_{n} = \cfrac{1}{\sqrt{5}}\cdot\left(\cfrac{1+\sqrt{5}}{2}\right)^n-\cfrac{1}{\sqrt{5}}\cdot\left(\cfrac{1-\sqrt{5}}{2}\right)^n~.$

The matrix has a determinant of −1, and thus it is a 2×2 unimodular matrix.

This property can be understood in terms of the continued fraction representation for the golden ratio:

$\varphi = 1 + \cfrac{1}{1 + \cfrac{1}{1 + \cfrac{1}{1 + \;\;\ddots\,}}}$

The Fibonacci numbers occur as the ratio of successive convergents of the continued fraction for φ, and the matrix formed from successive convergents of any continued fraction has a determinant of +1 or −1. The matrix representation gives the following closed expression for the Fibonacci numbers:

$\begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix}^n = \begin{pmatrix} F_{n+1} & F_n \\ F_n & F_{n-1} \end{pmatrix}.$

Taking the determinant of both sides of this equation yields Cassini's identity,

$(-1)^n = F_{n+1}F_{n-1} - F_n^2\,.$

Moreover, since An Am = An+m for any square matrix A, the following identities can be derived,

\begin{align} {F_m}{F_n} + {F_{m-1}}{F_{n-1}} &= F_{m+n-1}\\ F_{m} F_{n+1} + F_{m-1} F_n &= F_{m+n} ~ . \end{align}

In particular, with m = n,

\begin{align} F_{2n-1} &= F_n^2 + F_{n-1}^2\\ F_{2n} &= (F_{n-1}+F_{n+1})F_n\\ &= (2F_{n-1}+F_n)F_n ~ . \end{align}

These last two identities provide a way to compute Fibonacci numbers recursively in O(log(n)) arithmetic operations and in time O(M(n) log(n)), where M(n) is the time for the multplication of two numbers of n digits. This matches the time for computing the nth Fibonacci number from the closed-form matrix formula, but with fewer redundant steps if one avoids recomputing an already computed Fibonacci number (recursion with memoization).[26]

## Recognizing Fibonacci numbers

The question may arise whether a positive integer x is a Fibonacci number. This is true if and only if one or both of $5x^2+4$ or $5x^2-4$ is a perfect square.[27] This is because Binet's formula above can be rearranged to give

$n = \log_{\varphi}\left(\frac{F_n\sqrt{5} + \sqrt{5F_n^2 \pm 4}}{2}\right)$,

which allows one to find the position in the sequence of a given Fibonacci number.

This formula must return an integer for all n, so the expression under the radical must be an integer (otherwise the logarithm does not even return a rational number).

## Combinatorial identities

Most identities involving Fibonacci numbers can be proven using combinatorial arguments using the fact that Fn can be interpreted as the number of sequences of 1s and 2s that sum to n − 1. This can be taken as the definition of Fn, with the convention that F0 = 0, meaning no sum adds up to −1, and that F1 = 1, meaning the empty sum "adds up" to 0. Here, the order of the summand matters. For example, 1 + 2 and 2 + 1 are considered two different sums.

For example, the recurrence relation

$F_{n} = F_{n-1} + F_{n-2},\,$

or in words, the nth Fibonacci number is the sum of the previous two Fibonacci numbers, may be shown by dividing the Fn sums of 1s and 2s that add to n − 1 into two non-overlapping groups. One group contains those sums whose first term is 1 and the other those sums whose first term is 2. In the first group the remaining terms add to n − 2, so it has F(n − 1) sums, and in the second group the remaining terms add to n − 3, so there are Fn−2 sums. So there are a total of Fn−1 + Fn−2 sums altogether, showing this is equal to Fn.

Similarly, it may be shown that the sum of the first Fibonacci numbers up to the nth is equal to the (n + 2)-nd Fibonacci number minus 1.[28] In symbols:

$\sum_{i=1}^n F_i = F_{n+2} - 1$

This is done by dividing the sums adding to n + 1 in a different way, this time by the location of the first 2. Specifically, the first group consists of those sums that start with 2, the second group those that start 1 + 2, the third 1 + 1 + 2, and so on, until the last group, which consists of the single sum where only 1's are used. The number of sums in the first group is F(n), F(n − 1) in the second group, and so on, with 1 sum in the last group. So the total number of sums is F(n) + F(n − 1) + ... + F(1) + 1 and therefore this quantity is equal to F(n + 2).

A similar argument, grouping the sums by the position of the first 1 rather than the first 2, gives two more identities:

$\sum_{i=0}^{n-1} F_{2i+1} = F_{2n}$

and

$\sum_{i=1}^{n} F_{2i} = F_{2n+1}-1.$

In words, the sum of the first Fibonacci numbers with odd index up to F2n−1 is the (2n)th Fibonacci number, and the sum of the first Fibonacci numbers with even index up to F2n is the (2n + 1)th Fibonacci number minus 1.[29]

A different trick may be used to prove

$\sum_{i=1}^n {F_i}^2 = F_{n} F_{n+1},$

or in words, the sum of the squares of the first Fibonacci numbers up to Fn is the product of the nth and (n + 1)th Fibonacci numbers. In this case note that Fibonacci rectangle of size Fn by F(n + 1) can be decomposed into squares of size Fn, Fn−1, and so on to F1 = 1, from which the identity follows by comparing areas.

## Other identities

Numerous other identities can be derived using various methods. Some of the most noteworthy are:[30]

Catalan's Identity:

$F_n^2 - F_{n+r}F_{n-r} = (-1)^{n-r}F_r^2$

Cassini's Identity:

$F_n^2 - F_{n+1}F_{n-1} = (-1)^{n-1}$

d'Ocagne's identity:

$F_m F_{n+1} - F_{m+1} F_n = (-1)^n F_{m-n}$
$F_{2n} = F_{n+1}^2 - F_{n-1}^2 = F_n \left (F_{n+1}+F_{n-1} \right ) = F_nL_n$

where Ln is the n'th Lucas number. The last is an identity for doubling n; other identities of this type are

$F_{3n} = 2F_n^3 + 3F_n F_{n+1} F_{n-1} = 5F_n^3 + 3 (-1)^n F_n$

by Cassini's identity.

\begin{align} F_{3n+1} &= F_{n+1}^3 + 3 F_{n+1}F_n^2 - F_n^3 \\ F_{3n+2} &= F_{n+1}^3 + 3 F_{n+1}^2F_n + F_n^3 \\ F_{4n} &= 4F_nF_{n+1} \left (F_{n+1}^2 + 2F_n^2 \right ) - 3F_n^2 \left (F_n^2 + 2F_{n+1}^2 \right ) \end{align}

These can be found experimentally using lattice reduction, and are useful in setting up the special number field sieve to factorize a Fibonacci number.

More generally,[30]

$F_{kn+c} = \sum_{i=0}^k {k\choose i} F_{c-i} F_n^i F_{n+1}^{k-i}.$

Putting k = 2 in this formula, one gets again the formulas of the end of above section Matrix form.

## Power series

The generating function of the Fibonacci sequence is the power series

$s(x)=\sum_{k=0}^{\infty} F_k x^k.$

This series is convergent for $|x| < \frac{1}{\varphi},$ and its sum has a simple closed-form:[31]

$s(x)=\frac{x}{1-x-x^2}$

This can be proven by using the Fibonacci recurrence to expand each coefficient in the infinite sum:

\begin{align} s(x) &= \sum_{k=0}^{\infty} F_k x^k \\ &= F_0 + F_1x + \sum_{k=2}^{\infty} \left( F_{k-1} + F_{k-2} \right) x^k \\ &= x + \sum_{k=2}^{\infty} F_{k-1} x^k + \sum_{k=2}^{\infty} F_{k-2} x^k \\ &= x + x\sum_{k=0}^{\infty} F_k x^k + x^2\sum_{k=0}^{\infty} F_k x^k \\ &= x + x s(x) + x^2 s(x). \end{align}

Solving the equation

$s(x)=x+xs(x)+x^2s(x)$

for s(x) results in the above closed form.

If x is the inverse of an integer, the closed form of the series becomes

$\sum_{n=0}^\infty\,\frac{F_n}{k^{n}}\,=\,\frac{k}{k^{2}-k-1}.$

In particular,

$\sum_{n = 1}^{\infty}{\frac {F_n}{10^{(k + 1)(n + 1)}}} = \frac {1}{10^{2k + 2} - 10^{k + 1} - 1}$

for all non-negative integers k.

Some math puzzle-books present as curious the particular value $\frac{s(\frac{1}{10})}{10}=\frac{1}{89}$.[32]

## Reciprocal sums

Infinite sums over reciprocal Fibonacci numbers can sometimes be evaluated in terms of theta functions. For example, we can write the sum of every odd-indexed reciprocal Fibonacci number as

$\sum_{k=0}^\infty \frac{1}{F_{2k+1}} = \frac{\sqrt{5}}{4}\vartheta_2^2 \left(0, \frac{3-\sqrt 5}{2}\right) ,$

and the sum of squared reciprocal Fibonacci numbers as

$\sum_{k=1}^\infty \frac{1}{F_k^2} = \frac{5}{24} \left(\vartheta_2^4\left(0, \frac{3-\sqrt 5}{2}\right) - \vartheta_4^4\left(0, \frac{3-\sqrt 5}{2}\right) + 1 \right).$

If we add 1 to each Fibonacci number in the first sum, there is also the closed form

$\sum_{k=0}^\infty \frac{1}{1+F_{2k+1}} = \frac{\sqrt{5}}{2},$

and there is a nested sum of squared Fibonacci numbers giving the reciprocal of the golden ratio,

$\sum_{k=1}^\infty \frac{(-1)^{k+1}}{\sum_{j=1}^k {F_{j}}^2} = \frac{\sqrt{5}-1}{2}.$

No closed formula for the reciprocal Fibonacci constant

$\psi = \sum_{k=1}^{\infty} \frac{1}{F_k} = 3.359885666243 \dots$

is known, but the number has been proved irrational by Richard André-Jeannin.

Millin series gives the identity[33]

$\sum_{n=0}^{\infty} \frac{1}{F_{2^n}} = \frac{7 - \sqrt{5}}{2}\,,$

which follows from the closed form for its partial sums as N tends to infinity:

$\sum_{n=0}^N \frac{1}{F_{2^n}} = 3 - \frac{F_{2^N-1}}{F_{2^N}}.$

## Primes and divisibility

### Divisibility properties

Every 3rd number of the sequence is even and more generally, every kth number of the sequence is a multiple of Fk. Thus the Fibonacci sequence is an example of a divisibility sequence. In fact, the Fibonacci sequence satisfies the stronger divisibility property[34][35]

$\gcd(F_m,F_n) = F_{\gcd(m,n)}.$

Any three consecutive Fibonacci numbers are pairwise coprime, which means that, for every n,

gcd(Fn, Fn+1) = gcd(Fn, Fn+2) = gcd(Fn+1, Fn+2) = 1.

Every prime number p divides a Fibonacci number that can be determined by the value of p modulo 5. If p is congruent to 1 or 4 (mod 5), then p divides Fp − 1, and if p is congruent to 2 or 3 (mod 5), then, p divides Fp + 1. The remaining case is that p = 5, and in this case p divides Fp. These cases can be combined into a single formula, using the Legendre symbol:[36]

$p \mid F_{p-\left(\frac{5}{p}\right)}.$

### Primality Testing

The above formula can be used as a primality test in the sense that if

$n \mid F_{n-\left(\frac{5}{n}\right)}$, where the Legendre symbol has been replaced by the Jacobi symbol, then this is evidence that n is a prime, and if it fails to hold, then n is definitely not a prime. If n is composite and satisfies the formula, then n is a Fibonacci pseudoprime.

When m is large--say a 500-bit number--then we can calculate Fm (mod n) efficiently using the matrix form. Thus

$\begin{pmatrix} F_{m+1} & F_m \\ F_m & F_{m-1} \end{pmatrix}$$\begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix}^m$ (mod m).

Here the matrix power Am is calculated using Modular exponentiation, which can be adapted to matrices--modular exponentiation for matrices[37]

### Fibonacci primes

Main article: Fibonacci prime

A Fibonacci prime is a Fibonacci number that is prime. The first few are:

2, 3, 5, 13, 89, 233, 1597, 28657, 514229, ... (sequence A005478 in OEIS).

Fibonacci primes with thousands of digits have been found, but it is not known whether there are infinitely many.[38]

Fkn is divisible by Fn, so, apart from F4 = 3, any Fibonacci prime must have a prime index. As there are arbitrarily long runs of composite numbers, there are therefore also arbitrarily long runs of composite Fibonacci numbers.

No Fibonacci number greater than F6 = 8 is one greater or one less than a prime number.[39]

The only nontrivial square Fibonacci number is 144.[40] Attila Pethő proved in 2001 that there is only a finite number of perfect power Fibonacci numbers.[41] In 2006, Y. Bugeaud, M. Mignotte, and S. Siksek proved that 8 and 144 are the only such non-trivial perfect powers.[42]

### Prime divisors of Fibonacci numbers

With the exceptions of 1, 8 and 144 (F1 = F2, F6 and F12) every Fibonacci number has a prime factor that is not a factor of any smaller Fibonacci number (Carmichael's theorem).[43]

The divisibility of Fibonacci numbers by a prime p is related to the Legendre symbol $\left(\tfrac{p}{5}\right)$ which is evaluated as follows:

$\left(\frac{p}{5}\right) = \begin{cases} 0 & \textrm{if}\;p =5\\ 1 &\textrm{if}\;p \equiv \pm1 \pmod 5\\ -1 &\textrm{if}\;p \equiv \pm2 \pmod 5.\end{cases}$

If p is a prime number then

$F_p \equiv \left(\frac{p}{5}\right) \pmod p \quad \text{and}\quad F_{p-\left(\frac{p}{5}\right)} \equiv 0 \pmod p.$[44][45]

For example,

\begin{align} (\tfrac{2}{5}) &= -1, &F_3 &= 2, &F_2&=1, \\ (\tfrac{3}{5}) &= -1, &F_4 &= 3,&F_3&=2, \\ (\tfrac{5}{5}) &= 0, &F_5 &= 5, \\ (\tfrac{7}{5}) &= -1, &F_8 &= 21,&F_7&=13, \\ (\tfrac{11}{5})& = +1, &F_{10}& = 55, &F_{11}&=89. \end{align}

It is not known whether there exists a prime p such that

$F_{p-\left(\frac{p}{5}\right)} \equiv 0 \pmod{p^2}.$

Such primes (if there are any) would be called Wall–Sun–Sun primes.

Also, if p ≠ 5 is an odd prime number then:[45]

$5F^2_{\frac{p \pm 1}{2}} \equiv \begin{cases} \tfrac{1}{2} \left (5\left(\frac{p}{5}\right)\pm 5 \right ) \pmod p & \textrm{if}\;p \equiv 1 \pmod 4\\ \tfrac{1}{2} \left (5\left(\frac{p}{5}\right)\mp 3 \right ) \pmod p & \textrm{if}\;p \equiv 3 \pmod 4. \end{cases}$

Example 1. p = 7, in this case p ≡ 3 (mod 4) and we have:

$(\tfrac{7}{5}) = -1: \qquad \tfrac{1}{2}\left (5(\tfrac{7}{5})+3 \right ) =-1, \quad \tfrac{1}{2} \left (5(\tfrac{7}{5})-3 \right )=-4.$
$F_3=2 \text{ and } F_4=3.$
$5F_3^2=20\equiv -1 \pmod {7}\;\;\text{ and }\;\;5F_4^2=45\equiv -4 \pmod {7}$

Example 2. p = 11, in this case p ≡ 3 (mod 4) and we have:

$(\tfrac{11}{5}) = +1: \qquad \tfrac{1}{2}\left (5(\tfrac{11}{5})+3 \right )=4, \quad \tfrac{1}{2} \left (5(\tfrac{11}{5})- 3 \right )=1.$
$F_5=5 \text{ and } F_6=8.$
$5F_5^2=125\equiv 4 \pmod {11} \;\;\text{ and }\;\;5F_6^2=320\equiv 1 \pmod {11}$

Example 3. p = 13, in this case p ≡ 1 (mod 4) and we have:

$(\tfrac{13}{5}) = -1: \qquad \tfrac{1}{2}\left (5(\tfrac{13}{5})-5 \right ) =-5, \quad \tfrac{1}{2}\left (5(\tfrac{13}{5})+ 5 \right )=0.$
$F_6=8 \text{ and } F_7=13.$
$5F_6^2=320\equiv -5 \pmod {13} \;\;\text{ and }\;\;5F_7^2=845\equiv 0 \pmod {13}$

Example 4. p = 29, in this case p ≡ 1 (mod 4) and we have:

$(\tfrac{29}{5}) = +1: \qquad \tfrac{1}{2}\left (5(\tfrac{29}{5})-5 \right )=0, \quad \tfrac{1}{2}\left (5(\tfrac{29}{5})+5 \right )=5.$
$F_{14}=377 \text{ and } F_{15}=610.$
$5F_{14}^2=710645\equiv 0 \pmod {29} \;\;\text{ and }\;\;5F_{15}^2=1860500\equiv 5 \pmod {29}$

For odd n, all odd prime divisors of Fn are congruent to 1 modulo 4, implying that all odd divisors of Fn (as the products of odd prime divisors) are congruent to 1 modulo 4.[46]

For example,

$F_1 = 1, F_3 = 2, F_5 = 5, F_7 = 13, F_9 = 34 = 2 \cdot 17, F_{11} = 89, F_{13} = 233, F_{15} = 610 = 2 \cdot 5 \cdot 61.$

All known factors of Fibonacci numbers F(i) for all i < 50000 are collected at the relevant repositories.[47][48]

### Periodicity modulo n

Main article: Pisano period

It may be seen that if the members of the Fibonacci sequence are taken mod n, the resulting sequence must be periodic with period at most n2−1. The lengths of the periods for various n form the so-called Pisano periods (sequence A001175 in OEIS). Determining the Pisano periods in general is an open problem, although for any particular n it can be solved as an instance of cycle detection.

## Right triangles

Starting with 5, every second Fibonacci number is the length of the hypotenuse of a right triangle with integer sides, or in other words, the largest number in a Pythagorean triple. The length of the longer leg of this triangle is equal to the sum of the three sides of the preceding triangle in this series of triangles, and the shorter leg is equal to the difference between the preceding bypassed Fibonacci number and the shorter leg of the preceding triangle.

The first triangle in this series has sides of length 5, 4, and 3. Skipping 8, the next triangle has sides of length 13, 12 (5 + 4 + 3), and 5 (8 − 3). Skipping 21, the next triangle has sides of length 34, 30 (13 + 12 + 5), and 16 (21 − 5). This series continues indefinitely. The triangle sides a, b, c can be calculated directly:

$\displaystyle a_n = F_{2n-1}$
$\displaystyle b_n = 2 F_n F_{n-1}$
$\displaystyle c_n = F_n^2 - F_{n-1}^2.$

These formulas satisfy $a_n ^2 = b_n ^2 + c_n ^2$ for all n, but they only represent triangle sides when n > 2.

Any four consecutive Fibonacci numbers Fn, Fn+1, Fn+2 and Fn+3 can also be used to generate a Pythagorean triple in a different way:[49]

$a = F_n F_{n+3} \, ; \, b = 2 F_{n+1} F_{n+2} \, ; \, c = F_{n+1}^2 + F_{n+2}^2 \, ; \, a^2 + b^2 = c^2 \,.$

Example 1: let the Fibonacci numbers be 1, 2, 3 and 5. Then:

$\displaystyle a = 1 \times 5 = 5$
$\displaystyle b = 2 \times 2 \times 3 = 12$
$\displaystyle c = 2^2 + 3^2 = 13 \,$
$\displaystyle 5^2 + 12^2 = 13^2 \,.$

## Magnitude

Since Fn is asymptotic to $\varphi^n/\sqrt5$, the number of digits in Fn is asymptotic to $n\,\log_{10}\varphi\approx0.2090\,n$. As a consequence, for every integer d > 1 there are either 4 or 5 Fibonacci numbers with d decimal digits.

More generally, in the base b representation, the number of digits in Fn is asymptotic to $n\,\log_b\varphi$.

## Applications

The Fibonacci numbers are important in the computational run-time analysis of Euclid's algorithm to determine the greatest common divisor of two integers: the worst case input for this algorithm is a pair of consecutive Fibonacci numbers.[50]

Yuri Matiyasevich was able to show that the Fibonacci numbers can be defined by a Diophantine equation, which led to his original solution of Hilbert's tenth problem.

The Fibonacci numbers are also an example of a complete sequence. This means that every positive integer can be written as a sum of Fibonacci numbers, where any one number is used once at most.

Moreover, every positive integer can be written in a unique way as the sum of one or more distinct Fibonacci numbers in such a way that the sum does not include any two consecutive Fibonacci numbers. This is known as Zeckendorf's theorem, and a sum of Fibonacci numbers that satisfies these conditions is called a Zeckendorf representation. The Zeckendorf representation of a number can be used to derive its Fibonacci coding.

Fibonacci numbers are used by some pseudorandom number generators.

Fibonacci numbers are used in a polyphase version of the merge sort algorithm in which an unsorted list is divided into two lists whose lengths correspond to sequential Fibonacci numbers – by dividing the list so that the two parts have lengths in the approximate proportion φ. A tape-drive implementation of the polyphase merge sort was described in The Art of Computer Programming.

Fibonacci numbers arise in the analysis of the Fibonacci heap data structure.

The Fibonacci cube is an undirected graph with a Fibonacci number of nodes that has been proposed as a network topology for parallel computing.

A one-dimensional optimization method, called the Fibonacci search technique, uses Fibonacci numbers.[51]

The Fibonacci number series is used for optional lossy compression in the IFF 8SVX audio file format used on Amiga computers. The number series compands the original audio wave similar to logarithmic methods such as µ-law.[52][53]

Since the conversion factor 1.609344 for miles to kilometers is close to the golden ratio (denoted φ), the decomposition of distance in miles into a sum of Fibonacci numbers becomes nearly the kilometer sum when the Fibonacci numbers are replaced by their successors. This method amounts to a radix 2 number register in golden ratio base φ being shifted. To convert from kilometers to miles, shift the register down the Fibonacci sequence instead.[54]

## In nature

Further information: Patterns in nature and Phyllotaxis
Yellow Chamomile head showing the arrangement in 21 (blue) and 13 (aqua) spirals. Such arrangements involving consecutive Fibonacci numbers appear in a wide variety of plants.

Fibonacci sequences appear in biological settings,[9] in two consecutive Fibonacci numbers, such as branching in trees, arrangement of leaves on a stem, the fruitlets of a pineapple,[10] the flowering of artichoke, an uncurling fern and the arrangement of a pine cone,[11] and the family tree of honeybees.[55] However, numerous poorly substantiated claims of Fibonacci numbers or golden sections in nature are found in popular sources, e.g., relating to the breeding of rabbits in Fibonacci's own unrealistic example, the seeds on a sunflower, the spirals of shells, and the curve of waves.[56]

Przemysław Prusinkiewicz advanced the idea that real instances can in part be understood as the expression of certain algebraic constraints on free groups, specifically as certain Lindenmayer grammars.[57]

Illustration of Vogel's model for n=1 ... 500

A model for the pattern of florets in the head of a sunflower was proposed by H. Vogel in 1979.[58] This has the form

$\theta = \frac{2\pi}{\phi^2} n,\ r = c \sqrt{n}$

where n is the index number of the floret and c is a constant scaling factor; the florets thus lie on Fermat's spiral. The divergence angle, approximately 137.51°, is the golden angle, dividing the circle in the golden ratio. Because this ratio is irrational, no floret has a neighbor at exactly the same angle from the center, so the florets pack efficiently. Because the rational approximations to the golden ratio are of the form F(j):F(j + 1), the nearest neighbors of floret number n are those at n ± F(j) for some index j, which depends on r, the distance from the center. It is often said that sunflowers and similar arrangements have 55 spirals in one direction and 89 in the other (or some other pair of adjacent Fibonacci numbers), but this is true only of one range of radii, typically the outermost and thus most conspicuous.[59]

### The bee ancestry code

Fibonacci numbers also appear in the pedigrees of idealized honeybees, according to the following rules:

• If an egg is laid by an unmated female, it hatches a male or drone bee.
• If, however, an egg was fertilized by a male, it hatches a female.

Thus, a male bee always has one parent, and a female bee has two.

If one traces the pedigree of any male bee (1 bee), he has 1 parent (1 bee), 2 grandparents, 3 great-grandparents, 5 great-great-grandparents, and so on. This sequence of numbers of parents is the Fibonacci sequence. The number of ancestors at each level, Fn, is the number of female ancestors, which is Fn−1, plus the number of male ancestors, which is Fn−2.[60] This is under the unrealistic assumption that the ancestors at each level are otherwise unrelated.

## Generalizations

The Fibonacci sequence has been generalized in many ways. These include:

• Generalizing the index to negative integers to produce the negafibonacci numbers.
• Generalizing the index to real numbers using a modification of Binet's formula.[30]
• Starting with other integers. Lucas numbers have L1 = 1, L2 = 3, and Ln = Ln−1 + Ln−2. Primefree sequences use the Fibonacci recursion with other starting points to generate sequences in which all numbers are composite.
• Letting a number be a linear function (other than the sum) of the 2 preceding numbers. The Pell numbers have Pn = 2Pn − 1 + Pn − 2.
• Not adding the immediately preceding numbers. The Padovan sequence and Perrin numbers have P(n) = P(n − 2) + P(n − 3).
• Generating the next number by adding 3 numbers (tribonacci numbers), 4 numbers (tetranacci numbers), or more. The resulting sequences are known as n-Step Fibonacci numbers.[61]
• Adding other objects than integers, for example functions or strings – one essential example is Fibonacci polynomials.

## Notes

1. ^ a b
2. ^ a b Bona 2011, p. 180.
3. ^ John Hudson Tiner (200). Exploring the World of Mathematics: From Ancient Record Keeping to the Latest Advances in Computers. New Leaf Publishing Group. ISBN 9781614581550.
4. ^ Lucas 1891, p. 3.
5. ^ a b Pisano 2002, pp. 404–5.
6. ^ a b Goonatilake, Susantha (1998), Toward a Global Science, Indiana University Press, p. 126, ISBN 978-0-253-33388-9
7. ^ a b Singh, Parmanand (1985), The So-called Fibonacci numbers in ancient and medieval India, Historia Mathematica 12 (3): 229–44, doi:10.1016/0315-0860(85)90021-7
8. ^ a b Knuth, Donald (2006), The Art of Computer Programming, 4. Generating All Trees – History of Combinatorial Generation, Addison–Wesley, p. 50, ISBN 978-0-321-33570-8, "it was natural to consider the set of all sequences of [L] and [S] that have exactly m beats. ...there are exactly Fm+1 of them. For example the 21 sequences when m = 7 are: [gives list]. In this way Indian prosodists were led to discover the Fibonacci sequence, as we have observed in Section 1.2.8 (from v.1)"
9. ^ a b Douady, S; Couder, Y (1996), Phyllotaxis as a Dynamical Self Organizing Process (PDF), Journal of Theoretical Biology 178 (178): 255–74, doi:10.1006/jtbi.1996.0026
10. ^ a b Jones, Judy; Wilson, William (2006), "Science", An Incomplete Education, Ballantine Books, p. 544, ISBN 978-0-7394-7582-9
11. ^ a b Brousseau, A (1969), Fibonacci Statistics in Conifers, Fibonacci Quarterly (7): 525–32
12. ^ Knuth, Donald (1968), The Art of Computer Programming 1, Addison Wesley, ISBN 81-7758-754-4, "Before Fibonacci wrote his work, the sequence Fn had already been discussed by Indian scholars, who had long been interested in rhythmic patterns... both Gopala (before 1135 AD) and Hemachandra (c. 1150) mentioned the numbers 1,2,3,5,8,13,21 explicitly [see P. Singh Historia Math 12 (1985) 229–44]" p. 100 (3d ed)..."
13. ^ Agrawala, VS (1969), Pāṇinikālīna Bhāratavarṣa (Hn.). Varanasi-I: TheChowkhamba Vidyabhawan, "SadgurushiShya writes that Pingala was a younger brother of Pāṇini [Agrawala 1969, lb]. There is an alternative opinion that he was a maternal uncle of Pāṇini [Vinayasagar 1965, Preface, 121. ... Agrawala [1969, 463–76], after a careful investigation, in which he considered the views of earlier scholars, has concluded that Pāṇini lived between 480 and 410 BC"
14. ^ Velankar, HD (1962), ‘Vṛttajātisamuccaya’ of kavi Virahanka, Jodhpur: Rajasthan Oriental Research Institute, p. 101, ""For four, variations of meters of two [and] three being mixed, five happens. For five, variations of two earlier – three [and] four, being mixed, eight is obtained. In this way, for six, [variations] of four [and] of five being mixed, thirteen happens. And like that, variations of two earlier meters being mixed, seven morae [is] twenty-one. In this way, the process should be followed in all mātrā-vṛttas"
15. ^ Knott, Ron. "Fibonacci's Rabbits". University of Surrey Faculty of Engineering and Physical Sciences.
16. ^ Gardner, Martin (1996), Mathematical Circus, The Mathematical Association of America, p. 153, ISBN 0-88385-506-2, "It is ironic that Leonardo, who made valuable contributions to mathematics, is remembered today mainly because a 19th-century French number theorist, Édouard Lucas... attached the name Fibonacci to a number sequence that appears in a trivial problem in Liber abaci"
17. ^ Knott, R, "Fib table", Fibonacci, UK: Surrey has the first 300 Fn factored into primes and links to more extensive tables.
18. ^ Knuth, Donald (2008-12-11), "Negafibonacci Numbers and the Hyperbolic Plane", Annual meeting, The Fairmont Hotel, San Jose, CA: The Mathematical Association of America
19. ^ Lucas 1891, p. 7.
20. ^ Stanley, Richard (2011). Enumerative Combinatorics I (2nd ed.). Cambridge Univ. Press. p. "121, Ex 1.35". ISBN 978-1107602625.
21. ^
22. ^ Ball 2003, p. 156.
23. ^ Ball 2003, pp. 155–6.
24. ^ Kepler, Johannes (1966), A New Year Gift: On Hexagonal Snow, Oxford University Press, p. 92, ISBN 0-19-858120-3
25. ^ Strena seu de Nive Sexangula, 1611
26. ^
27. ^ Gessel, Ira (October 1972), Fibonacci is a Square (PDF), The Fibonacci Quarterly 10 (4): 417–19, retrieved April 11, 2012
28. ^ Lucas 1891, p. 4.
29. ^ Vorobiev, Nikolaĭ Nikolaevich; Martin, Mircea (2002), "Chapter 1", Fibonacci Numbers, Birkhäuser, pp. 5–6, ISBN 3-7643-6135-2
30. ^ a b c
31. ^ Glaister, P (1995), Fibonacci power series, The Mathematical Gazette 79 (486): 521, doi:10.2307/3618079
32. ^ Köhler, Günter (February 1985), Generating functions of Fibonacci-like sequences and decimal expansions of some fractions (PDF), The Fibonacci Quarterly 23 (1): 29–35, retrieved December 31, 2011
33. ^
34. ^ Ribenboim, Paulo (2000), My Numbers, My Friends, Springer-Verlag
35. ^ Su, Francis E (2000), "Fibonacci GCD's, please", Mudd Math Fun Facts, et al, HMC
36. ^ Williams, H. C. (1982), A note on the Fibonacci quotient $F_{p-\varepsilon}/p$, Canadian Mathematical Bulletin 25 (3): 366–370, doi:10.4153/CMB-1982-053-0, MR 668957. Williams calls this property "well known".
37. ^ Prime Numbers, Richard Crandall, Carl Pomerance, Springer, second edition, 2005, p.142.
38. ^
39. ^ Honsberger, Ross (1985), Mathematical Gems III, AMS Dolciani Mathematical Expositions (9): 133, ISBN 0-88385-318-3
40. ^ Cohn, JHE (1964), Square Fibonacci Numbers etc, Fibonacci Quarterly 2: 109–13
41. ^ Pethő, Attila (2001), Diophantine properties of linear recursive sequences II, Acta Math. Paedagogicae Nyíregyháziensis 17: 81–96
42. ^ Bugeaud, Y; Mignotte, M; Siksek, S (2006), Classical and modular approaches to exponential Diophantine equations. I. Fibonacci and Lucas perfect powers, Ann. Math. 2 (163): 969–1018, arXiv:math/0403046, Bibcode:2004math......3046B, doi:10.4007/annals.2006.163.969
43. ^ Knott, Ron, The Fibonacci numbers, UK: Surrey
44. ^ Ribenboim, Paulo (1996), The New Book of Prime Number Records, New York: Springer, p. 64, ISBN 0-387-94457-5
45. ^ a b Lemmermeyer 2000, pp. 73–4.
46. ^ Lemmermeyer 2000, p. 73.
47. ^ Fibonacci and Lucas factorizations, Mersennus collects all known factors of F(i) with i < 10000.
48. ^ Factors of Fibonacci and Lucas numbers, Red golpe collects all known factors of F(i) with 10000 < i < 50000.
49. ^ Koshy, Thomas (2007), Elementary number theory with applications, Academic Press, p. 581, ISBN 0-12-372487-2
50. ^ Knuth, Donald E (1997), The Art of Computer Programming, 1: Fundamental Algorithms (3rd ed.), Addison–Wesley, p. 343, ISBN 0-201-89683-4
51. ^ Avriel, M; Wilde, DJ (1966), Optimality of the Symmetric Fibonacci Search Technique, Fibonacci Quarterly (3): 265–9
52. ^ Amiga ROM Kernel Reference Manual, Addison–Wesley, 1991
53. ^ "IFF", Multimedia Wiki
54. ^ "Zeckendorf representation", Encyclopedia of Math
55. ^ "Marks for the da Vinci Code: B–". Maths. Computer Science For Fun: CS4FN.
56. ^ Simanek, D. "Fibonacci Flim-Flam". LHUP.
57. ^ Prusinkiewicz, Przemyslaw; Hanan, James (1989), Lindenmayer Systems, Fractals, and Plants (Lecture Notes in Biomathematics), Springer-Verlag, ISBN 0-387-97092-4
58. ^ Vogel, H (1979), A better way to construct the sunflower head, Mathematical Biosciences 44 (44): 179–89, doi:10.1016/0025-5564(79)90080-4
59. ^ Prusinkiewicz, Przemyslaw; Lindenmayer, Aristid (1990), The Algorithmic Beauty of Plants, Springer-Verlag, pp. 101–7, ISBN 978-0-387-97297-8
60. ^ The Fibonacci sequence as it appears in nature, The Fibonacci Quarterly 1 (1), 1963: 53–56
61. ^

## References

• Ball, Keith M (2003), "8: Fibonacci's Rabbits Revisited", Strange Curves, Counting Rabbits, and Other Mathematical Explorations, Princeton, NJ: Princeton University Press, ISBN 0-691-11321-1.
• Beck, Matthias; Geoghegan, Ross (2010), The Art of Proof: Basic Training for Deeper Mathematics, New York: Springer.
• Bóna, Miklós (2011), A Walk Through Combinatorics (3rd ed.), New Jersey: World Scientific.
• Lemmermeyer, Franz (2000), Reciprocity Laws, New York: Springer, ISBN 3-540-66957-4.
• Lucas, Édouard (1891), Théorie des nombres (in French) 1, Gauthier-Villars.
• Pisano, Leonardo (2002), Fibonacci's Liber Abaci: A Translation into Modern English of the Book of Calculation, Sources and Studies in the History of Mathematics and Physical Sciences, Sigler, Laurence E, trans, Springer, ISBN 0-387-95419-8