Jump to content

Simple continued fraction

From Wikipedia, the free encyclopedia

This is an old revision of this page, as edited by Luciotato (talk | contribs) at 10:16, 7 December 2014. The present address (URL) is a permanent link to this revision, which may differ significantly from the current revision.

A finite continued fraction, where n is a non-negative integer, a0 is an integer, and ai is a positive integer, for i=1,…,n.

In mathematics, a continued fraction is an expression obtained through an iterative process of representing a number as the sum of its integer part and the reciprocal of another number, then writing this other number as the sum of its integer part and another reciprocal, and so on.[1] In a finite continued fraction (or terminated continued fraction), the iteration/recursion is terminated after finitely many steps by using an integer in lieu of another continued fraction. In contrast, an infinite continued fraction is an infinite expression. In either case, all integers in the sequence, other than the first, must be positive. The integers ai are called the coefficients or terms of the continued fraction.[2]

Continued fractions have a number of remarkable properties related to the Euclidean algorithm for integers or real numbers. Every rational number p/q has two closely related expressions as a finite continued fraction, whose coefficients ai can be determined by applying the Euclidean algorithm to (p, q). The numerical value of an infinite continued fraction will be irrational; it is defined from its infinite sequence of integers as the limit of a sequence of values for finite continued fractions. Each finite continued fraction of the sequence is obtained by using a finite prefix of the infinite continued fraction's defining sequence of integers. Moreover, every irrational number α is the value of a unique infinite continued fraction, whose coefficients can be found using the non-terminating version of the Euclidean algorithm applied to the incommensurable values α and 1. This way of expressing real numbers (rational and irrational) is called their continued fraction representation.

It is generally assumed that the numerator of all of the fractions is 1. If arbitrary values and/or functions are used in place of one or more of the numerators or the integers in the denominators, the resulting expression is a generalized continued fraction. When it is necessary to distinguish the first form from generalized continued fractions, the former may be called a simple or regular continued fraction, or said to be in canonical form.

The term continued fraction may also refer to representations of rational functions, arising in their analytic theory. For this use of the term see Padé approximation and Chebyshev rational functions.

Motivation and notation

Consider a typical rational number 415/93, which is around 4.4624.

As a first approximation, start with 4, which is the integer part; 415/93 = 4 + 43/93.

Note that the fractional part is the reciprocal of 93/43 which is about 2.1628. Use the integer part, 2, as an approximation for the reciprocal, to get a second approximation of 4 + 1/2 = 4.5; 93/43 = 2 + 7/43.

The fractional part of 93/43 is the reciprocal of 43/7, and 43/7 is around 6.1429. Use 6 as an approximation for this to get 2 + 1/6 as an approximation for 93/43 and 4 + 1/2 + 1/6, about 4.4615, as the third approximation; 43/7 = 6 + 1/7 .

Finally, the fractional part of 43/7 is the reciprocal of 7, so its approximation in this scheme, 7, is exact (7/1 = 7 + 0/1) and produces the exact expression 4 + 1/2 + 1/6 + (1 / 7) for 415/93.

This expression is called the continued fraction representation of the number. Dropping some of the less essential parts of the expression 4 + 1/2 + 1/6 + (1 / 7) gives the abbreviated notation 415/93=[4;2,6,7]. Note that it is customary to replace only the first comma by a semicolon. Some older textbooks use all commas in the (n+1)-tuple, e.g. [4,2,6,7].[3][4]

If the starting number is rational then this process exactly parallels the Euclidean algorithm. In particular, it must terminate and produce a finite continued fraction representation of the number. If the starting number is irrational then the process continues indefinitely. This produces a sequence of approximations, all of which are rational numbers, and these converge to the starting number as a limit. This is the (infinite) continued fraction representation of the number. Examples of continued fraction representations of irrational numbers are:

  • 19 = [4;2,1,3,1,2,8,2,1,3,1,2,8,…] (sequence A010124 in the OEIS). The pattern repeats indefinitely with a period of 6.
  • e = [2;1,2,1,1,4,1,1,6,1,1,8,…] (sequence A003417 in the OEIS). The pattern repeats indefinitely with a period of 3 except that 2 is added to one of the terms in each cycle.
  • π = [3;7,15,1,292,1,1,1,2,1,3,1,…] (sequence A001203 in the OEIS). The terms in this representation are apparently random.
  • ϕ = [1;1,1,1,1,1,1,1,1,1,1,1,…] (sequence A000012 in the OEIS). The golden ratio, the most difficult irrational number to approximate rationally. See: A property of the golden ratio φ.

Continued fractions are, in some ways, more "mathematically natural" representations of a real number than other representations such as decimal representations, and they have several desirable properties:

  • The continued fraction representation for a rational number is finite and only rational numbers have finite representations. In contrast, the decimal representation of a rational number may be finite, for example 137/1600 = 0.085625, or infinite with a repeating cycle, for example 4/27 = 0.148148148148….
  • Every rational number has an essentially unique continued fraction representation. Each rational can be represented in exactly two ways, since [a0;a1,… an−1,an] = [a0;a1,… an−1,(an−1),1]. Usually the first, shorter one is chosen as the canonical representation.
  • The continued fraction representation of an irrational number is unique.
  • The real numbers whose continued fraction eventually repeats are precisely the quadratic irrationals.[5] For example, the repeating continued fraction [1;1,1,1,…] is the golden ratio, and the repeating continued fraction [1;2,2,2,…] is the square root of 2. In contrast, the decimal representations of quadratic irrationals are apparently random. The square roots of all (positive) integers, that are not perfect squares, are quadratic irrationals, hence are unique periodic continued fractions.
  • The successive approximations generated in finding the continued fraction representation of a number, i.e. by truncating the continued fraction representation, are in a certain sense (described below) the "best possible".

Basic formula

A continued fraction is an expression of the form

where ai, and bi are either rational numbers, real numbers, or complex numbers. If bi = 1 for all i the expression is called a simple continued fraction. If the expression contains a finite number of terms it is called a finite continued fraction. If the expression contains an infinite number of terms it is called an infinite continued fraction. [6]

Thus, all of the following illustrate valid finite simple continued fractions:

Examples of finite simple continued fractions
Formula Numeric Remarks
All integers are a degenerate case
Simplest possible fractional form
First integer may be negative
First integer may be zero

Calculating continued fraction representations

Consider a real number r. Let i be the integer part and f the fractional part of r. Then the continued fraction representation of r is [i;a1,a2,…], where [a1;a2,…] is the continued fraction representation of 1/f.

To calculate a continued fraction representation of a number r, write down the integer part (technically the floor) of r. Subtract this integer part from r. If the difference is 0, stop; otherwise find the reciprocal of the difference and repeat. The procedure will halt if and only if r is rational. This process can be efficiently implemented using the Euclidean algorithm when the number is rational.

Find the continued fraction for 3.245 (= 349/200)
Step Real Number Integer part Fractional part Simplified Reciprocal of f Simplified
1 r = 349/200 i = 3 f = 349/200 − 3 = 49/200 1/f = 200/49 = 44/49
2 r = 44/49 i = 4 f = 44/49 − 4 = 4/49 1/f = 49/4 = 121/4
3 r = 121/4 i = 12 f = 121/4 − 12 = 1/4 1/f = 4/1 = 4
4 r = 4 i = 4 f = 4 − 4 = 0 STOP
Continued fraction form for 3.245 or 349/200 is [3; 4, 12, 4].
349/200 = 3 + 1/4 + 1/12 + 1/4

The number 3.245 can also be represented by the continued fraction expansion [3;4,12,3,1]; refer to Finite continued fractions below.

Notations for continued fractions

The integers a0, a1 etc., are called the coefficients or terms of the continued fraction.[2] One can abbreviate the continued fraction

in the notation of Carl Friedrich Gauss

or as

,

or in the notation of Pringsheim as

or in another related notation as

Sometimes angle brackets are used, like this:

The semicolon in the square and angle bracket notations is sometimes replaced by a comma.[3][4]

One may also define infinite simple continued fractions as limits:

This limit exists for any choice of a0 and positive integers a1a2, ... .

Finite continued fractions

Every finite continued fraction represents a rational number, and every rational number can be represented in precisely two different ways as a finite continued fraction, with the conditions that the first coefficient is an integer and other coefficients being positive integers. These two representations agree except in their final terms. In the longer representation the final term in the continued fraction is 1; the shorter representation drops the final 1, but increases the new final term by 1. The final element in the short representation is therefore always greater than 1, if present. In symbols:

[a0; a1, a2, …, an − 1, an, 1] = [a0; a1, a2, …, an − 1, an + 1].
[a0; 1] = [a0 + 1].

For example,

2.25 = 2 + 1/4 = [2; 4] = 2 + 1/3 + 1/1 = [2; 3, 1]
−4.2 = −5 + 4/5 = −5 + 1/1 + 1/4 = [−5; 1, 4] = −5 + 1/1 + 1/3 + 1/1 = [−5; 1, 3, 1].

Continued fractions of reciprocals

The continued fraction representations of a positive rational number and its reciprocal are identical except for a shift one place left or right depending on whether the number is less than or greater than one respectively. In other words, the numbers represented by [a0; a1, a2, …, an] and [0; a0, a1, …, an] are reciprocals. This is because if a is an integer then if x < 1 then x = 0 + 1/a + 1/b and 1/x = a + 1/b and if x > 1 then x = a + 1/b and 1/x = 0 + 1/a + 1/b with the last number that generates the remainder of the continued fraction being the same for both x and its reciprocal.

For example,

2.25 = 9/4 = [2; 4],
1/2.25 = 4/9 = [0; 2, 4].

Infinite continued fractions

Every infinite continued fraction is irrational, and every irrational number can be represented in precisely one way as an infinite continued fraction.

An infinite continued fraction representation for an irrational number is useful because its initial segments provide rational approximations to the number. These rational numbers are called the convergents of the continued fraction. The larger a term is in the continued fraction, the closer the corresponding convergent is to the irrational number being approximated. Numbers like π have occasional large terms in their continued fraction, which makes them easy to approximate with rational numbers. Other numbers like e have only small terms early in their continued fraction, which makes them more difficult to approximate rationally. The golden ratio ϕ has terms equal to 1 everywhere—the smallest values possible—which makes ϕ the most difficult number to approximate rationally. In this sense, therefore, it is the "most irrational" of all irrational numbers. Even-numbered convergents are smaller than the original number, while odd-numbered ones are larger.

For a continued fraction [a0; a1, a2, …], the first four convergents (numbered 0 through 3) are

a0/1, a1a0 + 1/a1, a2(a1a0 + 1) + a0/a2a1 + 1, a3(a2(a1a0 + 1) + a0) + (a1a0 + 1)/ a3(a2a1 + 1) + a1

In words, the numerator of the third convergent is formed by multiplying the numerator of the second convergent by the third quotient, and adding the numerator of the first convergent. The denominators are formed similarly. Therefore, each convergent can be expressed explicitly in terms of the continued fraction as the ratio of certain multivariate polynomials called continuants.

If successive convergents are found, with numerators h1, h2, … and denominators k1, k2, … then the relevant recursive relation is:

hn = anhn − 1 + hn − 2,
kn = ankn − 1 + kn − 2.

The successive convergents are given by the formula

hn/kn = anhn − 1 + hn − 2/ankn − 1 + kn − 2

Thus to incorporate a new term into a rational approximation, only the two previous convergents are necessary. The initial "convergents" (required for the first two terms) are 01 and 10. For example, here are the convergents for [0;1,5,2,2].

n −2 −1 0 1 2 3 4
an     0 1 5 2 2
hn 0 1 0 1 5 11 27
kn 1 0 1 1 6 13 32

When using the Babylonian method to generate successive approximations to the square root of an integer, if one starts with the lowest integer as first approximant, the rationals generated all appear in the list of convergents for the continued fraction. Specifically, the approximants will appear on the convergents list in positions 0, 1, 3, 7, 15, … , 2k−1, ... For example, the continued fraction expansion for 3 is [1;1,2,1,2,1,2,1,2,…]. Comparing the convergents with the approximants derived from the Babylonian method:

n −2 −1 0 1 2 3 4 5 6 7
an     1 1 2 1 2 1 2 1
hn 0 1 1 2 5 7 19 26 71 97
kn 1 0 1 1 3 4 11 15 41 56
x0 = 1 = 1/1
x1 = 1/2(1 + 3/1) = 2/1 = 2
x2 = 1/2(2 + 3/2) = 7/4
x3 = 1/2(7/4 + 3/7/4) = 97/56

Some useful theorems

If a0, a1, a2, … is an infinite sequence of positive integers, define the sequences hn and kn recursively:

Theorem 1. For any positive real number z

Theorem 2. The convergents of [a0; a1, a2, …] are given by

Theorem 3. If the nth convergent to a continued fraction is hn/kn, then

Corollary 1: Each convergent is in its lowest terms (for if hn and kn had a nontrivial common divisor it would divide knhn−1kn−1hn, which is impossible).

Corollary 2: The difference between successive convergents is a fraction whose numerator is unity:

Corollary 3: The continued fraction is equivalent to a series of alternating terms:

Corollary 4: The matrix

has determinant plus or minus one, and thus belongs to the group of 2×2 unimodular matrices GL(2, Z).

Theorem 4. Each (s-th) convergent is nearer to a subsequent (n-th) convergent than any preceding (r-th) convergent is. In symbols, if the n-th convergent is taken to be [a0a1, ..., an] = xn, then

for all r < s < n.

Corollary 1: The even convergents (before the nth) continually increase, but are always less than xn.

Corollary 2: The odd convergents (before the nth) continually decrease, but are always greater than xn.

Theorem 5.

Corollary 1: Any convergent is nearer to the continued fraction than any other fraction whose denominator is less than that of the convergent

Corollary 2: Any convergent which immediately precedes a large quotient is a near approximation to the continued fraction.

Semiconvergents

If

hn − 1/kn − 1, hn/kn

are successive convergents, then any fraction of the form

hn − 1 + ahn/kn − 1 + akn

where a is a nonnegative integer and the numerators and denominators are between the n and n + 1 terms inclusive are called semiconvergents, secondary convergents, or intermediate fractions. Often the term is taken to mean that being a semiconvergent excludes the possibility of being a convergent, rather than that a convergent is a kind of semiconvergent.

The semiconvergents to the continued fraction expansion of a real number x include all the rational approximations which are better than any approximation with a smaller denominator. Another useful property is that consecutive semiconvergents a/b and c/d are such that adbc = ±1.

Best rational approximations

A best rational approximation to a real number x is a rational number n/d, d > 0, that is closer to x than any approximation with a smaller or equal denominator. The simple continued fraction for x generates all of the best rational approximations for x according to three rules:

  1. Truncate the continued fraction, and possibly decrement its last term.
  2. The decremented term cannot have less than half its original value.
  3. If the final term is even, half its value is admissible only if the corresponding semiconvergent is better than the previous convergent. (See below.)

For example, 0.84375 has continued fraction [0;1,5,2,2]. Here are all of its best rational approximations.

 [0;1]   [0;1,3]   [0;1,4]   [0;1,5]   [0;1,5,2]   [0;1,5,2,1]   [0;1,5,2,2] 
1 3/4 4/5 5/6 11/13 16/19 27/32

The strictly monotonic increase in the denominators as additional terms are included permits an algorithm to impose a limit, either on size of denominator or closeness of approximation.

The "half rule" mentioned above is that when ak is even, the halved term ak/2 is admissible if and only if |x − [a0 ; a1, …, ak − 1]| > |x − [a0 ; a1, …, ak − 1, ak/2]| [7] This is equivalent[7] to:[8]

[ak; ak − 1, …, a1] > [ak; ak + 1, …].

The convergents to x are best approximations in an even stronger sense: n/d is a convergent for x if and only if |dxn| is the least relative error among all approximations m/c with cd; that is, we have |dxn| < |cxm| so long as c < d. (Note also that |dkxnk| → 0 as k → ∞.)

Best rational within an interval

A rational that falls within the interval (x, y), for 0 < x < y, can be found with the continued fractions for x and y. When both x and y are irrational and

x = [a0; a1, a2, …, ak − 1, ak, ak + 1, …]
y = [a0; a1, a2, …, ak − 1, bk, bk + 1, …]

where x and y have identical continued fraction expansions up through ak−1, a rational that falls within the interval (x, y) is given by the finite continued fraction,

z(x,y) = [a0; a1, a2, …, ak − 1, min(ak, bk) + 1]

This rational will be best in that no other rational in (x, y) will have a smaller numerator or a smaller denominator.

If x is rational, it will have two continued fraction representations that are finite, x1 and x2, and similarly a rational y will have two representations, y1 and y2. The coefficients beyond the last in any of these representations should be interpreted as +∞; and the best rational will be one of z(x1, y1), z(x1, y2), z(x2, y1), or z(x2, y2).

For example, the decimal representation 3.1416 could be rounded from any number in the interval [3.14155, 3.14165]. The continued fraction representations of 3.14155 and 3.14165 are

3.14155 = [3; 7, 15, 2, 7, 1, 4, 1, 1] = [3; 7, 15, 2, 7, 1, 4, 2]
3.14165 = [3; 7, 16, 1, 3, 4, 2, 3, 1] = [3; 7, 16, 1, 3, 4, 2, 4]

and the best rational between these two is

[3; 7, 16] = 355/113 = 3.1415929....

Thus, in some sense, 355/113 is the best rational number corresponding to the rounded decimal number 3.1416.

Interval for a convergent

A rational number, which can be expressed as finite continued fraction in two ways,

z = [a0; a1, …, ak − 1, ak, 1] = [a0; a1, …, ak − 1, ak + 1]

will be one of the convergents for the continued fraction expansion of a number, if and only if the number is strictly between

x = [a0; a1, …, ak − 1, ak, 2] and
y = [a0; a1, …, ak − 1, ak + 2]

Note that the numbers x and y are formed by incrementing the last coefficient in the two representations for z, and that x < y when k is even, and x > y when k is odd.

For example, the number 355/113 has the continued fraction representations

355/113 = [3; 7, 15, 1] = [3; 7, 16]

and thus 355/113 is a convergent of any number strictly between

[3; 7, 15, 2] = 688/219 ≈ 3.1415525
[3; 7, 17] = 377/120 ≈ 3.1416667

Comparison of continued fractions

Consider x = [a0; a1, …] and y = [b0; b1, …]. If k is the smallest index for which ak is unequal to bk then x < y if (−1)k(akbk) < 0 and y < x otherwise.

If there is no such k, but one expansion is shorter than the other, say x = [a0; a1, …, an] and y = [b0; b1, …, bn, bn + 1, …] with ai = bi for 0 ≤ in, then x < y if n is even and y < x if n is odd.

Continued fraction expansions of π

To calculate the convergents of π we may set a0 = ⌊π⌋ = 3, define u1 = 1/π − 3 ≈ 7.0625 and a1 = ⌊u1⌋ = 7, u2 = 1/u1 − 7 ≈ 15.9665 and a2 = ⌊u2⌋ = 15, u3 = 1/u2 − 15 ≈ 1.003. Continuing like this, one can determine the infinite continued fraction of π as

[3;7,15,1,292,1,1,…] (sequence A001203 in the OEIS).

The fourth convergent of π is [3;7,15,1] = 355/113 = 3.14159292035..., sometimes called Milü, which is fairly close to the true value of π.

Let us suppose that the quotients found are, as above, [3;7,15,1]. The following is a rule by which we can write down at once the convergent fractions which result from these quotients without developing the continued fraction.

The first quotient, supposed divided by unity, will give the first fraction, which will be too small, namely, 3/1. Then, multiplying the numerator and denominator of this fraction by the second quotient and adding unity to the numerator, we shall have the second fraction, 22/7, which will be too large. Multiplying in like manner the numerator and denominator of this fraction by the third quotient, and adding to the numerator the numerator of the preceding fraction, and to the denominator the denominator of the preceding fraction, we shall have the third fraction, which will be too small. Thus, the third quotient being 15, we have for our numerator Template:J, and for our denominator, Template:J. The third convergent, therefore, is 333/106. We proceed in the same manner for the fourth convergent. The fourth quotient being 1, we say 333 times 1 is 333, and this plus 22, the numerator of the fraction preceding, is 355; similarly, 106 times 1 is 106, and this plus 7 is 113.

In this manner, by employing the four quotients [3;7,15,1], we obtain the four fractions:

3/1, 22/7, 333/106, 355/113, ….

These convergents are alternately smaller and larger than the true value of π, and approach nearer and nearer to π. The difference between a given convergent and π is less than the reciprocal of the product of the denominators of that convergent and the next convergent. For example, the fraction 22/7 is greater than π, but 22/7π is less than 1/7 × 106 = 1/742 (in fact, 22/7π is just less than 1/790).

The demonstration of the foregoing properties is deduced from the fact that if we seek the difference between one of the convergent fractions and the next adjacent to it we shall obtain a fraction of which the numerator is always unity and the denominator the product of the two denominators. Thus the difference between 22/7 and 3/1 is 1/7, in excess; between 333/106 and 22/7, 1/742, in deficit; between 355/113 and 333/106, 1/11978, in excess; and so on. The result being, that by employing this series of differences we can express in another and very simple manner the fractions with which we are here concerned, by means of a second series of fractions of which the numerators are all unity and the denominators successively be the product of every two adjacent denominators. Instead of the fractions written above, we have thus the series:

3/1 + 1/1 × 71/7 × 106 + 1/106 × 113 − …

The first term, as we see, is the first fraction; the first and second together give the second fraction, 22/7; the first, the second and the third give the third fraction 333/106, and so on with the rest; the result being that the series entire is equivalent to the original value.

Generalized continued fraction

A generalized continued fraction is an expression of the form

where the an (n > 0) are the partial numerators, the bn are the partial denominators, and the leading term b0 is called the integer part of the continued fraction.

To illustrate the use of generalized continued fractions, consider the following example. The sequence of partial denominators of the simple continued fraction of π does not show any obvious pattern:

or

However, several generalized continued fractions for π have a perfectly regular structure, such as:

The first two of these are special cases of the arctangent function with π = 4 arctan (1).

Other continued fraction expansions

Periodic continued fractions

The numbers with periodic continued fraction expansion are precisely the irrational solutions of quadratic equations with rational coefficients (rational solutions have finite continued fraction expansions as previously stated). The simplest examples are the golden ratio φ = [1;1,1,1,1,1,…] and 2 = [1;2,2,2,2,…]; while 14 = [3;1,2,1,6,1,2,1,6…] and 42 = [6;2,12,2,12,2,12…]. All irrational square roots of integers have a special form for the period; a symmetrical string, like the empty string (for 2) or 1,2,1 (for 14), followed by the double of the leading integer.

A property of the golden ratio φ

Because the continued fraction expansion for φ doesn't use any integers greater than 1, φ is one of the most "difficult" real numbers to approximate with rational numbers. Hurwitz's theorem[9] states that any real number k can be approximated by infinitely many rational m/n with

While virtually all real numbers k will eventually have infinitely many convergents m/n whose distance from k is significantly smaller than this limit, the convergents for φ (i.e., the numbers 5/3, 8/5, 13/8, 21/13, etc.) consistently "toe the boundary", keeping a distance of almost exactly away from φ, thus never producing an approximation nearly as impressive as, for example, 355/113 for π. It can also be shown that every real number of the form a + bφ/c + dφ, where a, b, c, and d are integers such that adbc = ±1, shares this property with the golden ratio φ; and that all other real numbers can be more closely approximated.

Regular patterns in continued fractions

While there is no discernable pattern in the simple continued fraction expansion of π, there is one for e, the base of the natural logarithm:

which is a special case of this general expression for positive integer n:

Another, more complex pattern appears in this continued fraction expansion for positive odd n:

with a special case for n = 1:

Other continued fractions of this sort are

where n is a positive integer; also, for integral n:

with a special case for n = 1:

If In(x) is the modified, or hyperbolic, Bessel function of the first kind, we may define a function on the rationals p/q by

which is defined for all rational numbers, with p and q in lowest terms. Then for all nonnegative rationals, we have

with similar formulas for negative rationals; in particular we have

Many of the formulas can be proved using Gauss's continued fraction.

Typical continued fractions

Most irrational numbers do not have any periodic or regular behavior in their continued fraction expansion. Nevertheless Khinchin proved that for almost all real numbers x, the ai (for i = 1, 2, 3, …) have an astonishing property: their geometric mean is a constant (known as Khinchin's constant, K ≈ 2.6854520010…) independent of the value of x. Paul Lévy showed that the nth root of the denominator of the nth convergent of the continued fraction expansion of almost all real numbers approaches an asymptotic limit, approximately 3.27582, which is known as Lévy's constant. Lochs' theorem states that nth convergent of the continued fraction expansion of almost all real numbers determines the number to an average accuracy of just over n decimal places.

Generalized continued fraction for square roots

Continued fraction techniques are one method of computing square roots.

The identity

(1)

leads via recursion to the generalized continued fraction for any square root:[10]

(2)

Pell's equation

Continued fractions play an essential role in the solution of Pell's equation. For example, for positive integers p and q, p2 − 2q2 = ±1 if and only if p/q is a convergent of 2.

Continued fractions and dynamical systems

Continued fractions also play a role in the study of dynamical systems, where they tie together the Farey fractions which are seen in the Mandelbrot set with Minkowski's question mark function and the modular group Gamma.

The backwards shift operator for continued fractions is the map h(x) = 1/x − ⌊1/x called the Gauss map, which lops off digits of a continued fraction expansion: h([0; a1, a2, a3, …]) = [0; a2, a3, …]. The transfer operator of this map is called the Gauss–Kuzmin–Wirsing operator. The distribution of the digits in continued fractions is given by the zero'th eigenvector of this operator, and is called the Gauss–Kuzmin distribution.

Eigenvalues and eigenvectors

The Lanczos algorithm uses a continued fraction expansion to iteratively approximate the eigenvalues and eigenvectors of a large sparse matrix.[11]

Examples of rational and irrational numbers

An infinite continued fraction is defined by the sequences , for , with .

A continued fraction is a mathematical expression that can be writen as a fraction with a denominator that is a sum that contains another simple or continued fraction. Depending on whether this iteration terminates with a simple fraction or not, the continued fraction is finite or infinite.

Different fields of mathematics have different terminology and notation for continued fraction. In number theory the standard unqualified use of the term continued fraction refers to the special case where all numerators are 1, and is treated in the article Simple continued fraction. The present article treats the case where numerators and denominators are sequences of constants or functions. From the perspective of number theory, these are called generalized continued fraction. From the perspective of complex analysis or numerical analysis, however, they are just standard, and in the present article they will simply be called "continued fraction".

Formulation

A continued fraction is an expression of the form

where the an (n > 0) are the partial numerators, the bn are the partial denominators, and the leading term b0 is called the integer part of the continued fraction.

The successive convergents of the continued fraction are formed by applying the fundamental recurrence formulas:

where An is the numerator and Bn is the denominator, called continuants,[12][13] of the nth convergent. They are given by the three-term recurrence relation [14]

with initial values

If the sequence of convergents {xn} approaches a limit, the continued fraction is convergent and has a definite value. If the sequence of convergents never approaches a limit, the continued fraction is divergent. It may diverge by oscillation (for example, the odd and even convergents may approach two different limits), or it may produce an infinite number of zero denominators Bn.

History

The story of continued fractions begins with the Euclidean algorithm,[15] a procedure for finding the greatest common divisor of two natural numbers m and n. That algorithm introduced the idea of dividing to extract a new remainder – and then dividing by the new remainder repeatedly.

Nearly two thousand years passed before Bombelli (1579) devised a technique for approximating the roots of quadratic equations with continued fractions in the mid-sixteenth century. Now the pace of development quickened. Just 24 years later, in 1613, Pietro Cataldi introduced the first formal notation for the generalized continued fraction.[16] Cataldi represented a continued fraction as

with the dots indicating where the next fraction goes, and each & representing a modern plus sign.

Late in the seventeenth century John Wallis introduced the term "continued fraction" into mathematical literature.[17] New techniques for mathematical analysis (Newton's and Leibniz's calculus) had recently come onto the scene, and a generation of Wallis' contemporaries put the new phrase to use.

In 1748 Euler published a theorem showing that a particular kind of continued fraction is equivalent to a certain very general infinite series.[18] Euler's continued fraction formula is still the basis of many modern proofs of convergence of continued fractions.

In 1761, Johann Heinrich Lambert gave the first proof that π is irrational, by using the following continued fraction for tan x:[19]

Continued fractions can also be applied to problems in number theory, and are especially useful in the study of Diophantine equations. In the late eighteenth century Lagrange used continued fractions to construct the general solution of Pell's equation, thus answering a question that had fascinated mathematicians for more than a thousand years.[20] Lagrange's discovery implies that the canonical continued fraction expansion of the square root of every non-square integer is periodic and that, if the period is of length p > 1, it contains a palindromic string of length p − 1.

In 1813 Gauss derived from complex-valued hypergeometric functions what is now called Gauss's continued fractions.[21] They can be used to express many elementary functions and some more advanced functions (such as the Bessel functions), as continued fractions that are rapidly convergent almost everywhere in the complex plane.

Notation

The long continued fraction expression displayed in the introduction is easy for an unfamiliar reader to interpret. However, it takes up a lot of space and can be difficult to typeset. So mathematicians have devised several alternative notations. One convenient way to express a generalized continued fraction sets each nested fraction on the same line, indicating the nesting by dangling plus signs in the denominators:

Sometimes the plus signs are typeset to vertically align with the denominators but not under the fraction bars:

Pringsheim wrote a generalized continued fraction this way:

Carl Friedrich Gauss evoked the more familiar infinite product Π when he devised this notation:

Here the "K" stands for Kettenbruch, the German word for "continued fraction". This is probably the most compact and convenient way to express continued fractions; however, it is not widely used by English typesetters.

Some elementary considerations

Here are some elementary results that are of fundamental importance in the further development of the analytic theory of continued fractions.

Partial numerators and denominators

If one of the partial numerators an + 1 is zero, the infinite continued fraction

is really just a finite continued fraction with n fractional terms, and therefore a rational function of a1 to an and b0 to bn + 1. Such an object is of little interest from the point of view adopted in mathematical analysis, so it is usually assumed that all ai ≠ 0. There is no need to place this restriction on the partial denominators bi.

The determinant formula

When the nth convergent of a continued fraction

is expressed as a simple fraction xn = An/Bn we can use the determinant formula

(1)

to relate the numerators and denominators of successive convergents xn and xn − 1 to one another. The proof for this can be easily seen by induction.

Proof

Base case

The case n = 1 results from a very simple computation.

Inductive step

Assume that (1) holds for n − 1. Then we need to see the same relation holding true for n. Substituting the value of An and Bn in (1) we obtain:
which is true because of our induction hypothesis.
Specifically, if neither Bn nor Bn − 1 is zero (n > 0) we can express the difference between the (n − 1)th and nth convergents like this:

The equivalence transformation

If {ci} = {c1, c2, c3, ...} is any infinite sequence of non-zero complex numbers we can prove, by induction, that

where equality is understood as equivalence, which is to say that the successive convergents of the continued fraction on the left are exactly the same as the convergents of the fraction on the right.

The equivalence transformation is perfectly general, but two particular cases deserve special mention. First, if none of the ai are zero, a sequence {ci} can be chosen to make each partial numerator a 1:

where c1 = 1/a1, c2 = a1/a2, c3 = a2/a1a3, and in general cn + 1 = 1/an + 1cn.

Second, if none of the partial denominators bi are zero we can use a similar procedure to choose another sequence {di} to make each partial denominator a 1:

where d1 = 1/b1 and otherwise dn + 1 = 1/bnbn + 1.

These two special cases of the equivalence transformation are enormously useful when the general convergence problem is analyzed.

Notions of convergence

As mentioned in the introduction, the continued fraction

converges if the sequence of convergents {xn} tends to a finite limit. This notion of convergence is very natural, but it is sometimes too restrictive. It is therefore useful to introduce the notion of general convergence of a continued fraction. Roughly speaking, this consists in replacing the part of the fraction by wn, instead of by 0, to compute the convergents. The convergents thus obtained are called modified convergents. We say that the continued fraction converges generally if there exists a sequence such that the sequence of modified convergents converges for all sufficiently distinct from . The sequence is then called an exceptional sequence for the continued fraction. See Chapter 2 of Lorentzen & Waadeland (1992) for a rigorous definition.

There also exists a notion of absolute convergence for continued fractions, which is based on the notion of absolute convergence of a series: a continued fraction is said to be absolutely convergent when the series

where are the convergents of the continued fraction, converges absolutely.[22] The Śleszyński–Pringsheim theorem provides a sufficient condition for absolute convergence.

Finally, a continued fraction of one or more complex variables is uniformly convergent in an open neighborhood Ω when its convergents converge uniformly on Ω; that is, when for every ε > 0 there exists M such that for all n > M, for all ,

Even and odd convergents

It is sometimes necessary to separate a continued fraction into its even and odd parts. For example, if the continued fraction diverges by oscillation between two distinct limit points p and q, then the sequence {x0, x2, x4, ...} must converge to one of these, and {x1, x3, x5, ...} must converge to the other. In such a situation it may be convenient to express the original continued fraction as two different continued fractions, one of them converging to p, and the other converging to q.

The formulas for the even and odd parts of a continued fraction can be written most compactly if the fraction has already been transformed so that all its partial denominators are unity. Specifically, if

is a continued fraction, then the even part xeven and the odd part xodd are given by

and

respectively. More precisely, if the successive convergents of the continued fraction x are {x1, x2, x3, ...}, then the successive convergents of xeven as written above are {x2, x4, x6, ...}, and the successive convergents of xodd are {x1, x3, x5, ...}.[23]

Conditions for irrationality

If a1, a2,... and b1, b2,... are positive integers with akbk for all sufficiently large k, then

converges to an irrational limit.[24]

Fundamental recurrence formulas

The partial numerators and denominators of the fraction's successive convergents are related by the fundamental recurrence formulas:

The continued fraction's successive convergents are then given by

These recurrence relations are due to John Wallis (1616–1703) and Leonhard Euler (1707–1783).[25] These recurrence relations are simply a different notation for the relations obtained by Pietro Antonio Cataldi (1548-1626).

As an example, consider the regular continued fraction in canonical form that represents the golden ratio φ:

Applying the fundamental recurrence formulas we find that the successive numerators An are {1, 2, 3, 5, 8, 13, ...} and the successive denominators Bn are {1, 1, 2, 3, 5, 8, ...}, the Fibonacci numbers. Since all the partial numerators in this example are equal to one, the determinant formula assures us that the absolute value of the difference between successive convergents approaches zero quite rapidly.

Linear fractional transformations

A linear fractional transformation (LFT) is a complex function of the form

where z is a complex variable, and a, b, c, d are arbitrary complex constants such that c + dz ≠ 0. An additional restriction that adbc is customarily imposed, to rule out the cases in which w = f(z) is a constant. The linear fractional transformation, also known as a Möbius transformation, has many fascinating properties. Four of these are of primary importance in developing the analytic theory of continued fractions.

  • If d ≠ 0 the LFT has one or two fixed points. This can be seen by considering the equation
which is clearly a quadratic equation in z. The roots of this equation are the fixed points of f(z). If the discriminant (cb)2 + 4ad is zero the LFT fixes a single point; otherwise it has two fixed points.
such that f(g(z)) = g(f(z)) = z for every point z in the extended complex plane, and both f and g preserve angles and shapes at vanishingly small scales. From the form of z = g(w) we see that g is also an LFT.
  • The composition of two different LFTs for which adbc is itself an LFT for which adbc. In other words, the set of all LFTs for which adbc is closed under composition of functions. The collection of all such LFTs, together with the "group operation" composition of functions, is known as the automorphism group of the extended complex plane.
  • If b = 0 the LFT reduces to
which is a very simple meromorphic function of z with one simple pole (at c/d) and a residue equal to a/d. (See also Laurent series.)

The continued fraction as a composition of LFTs

Consider a sequence of simple linear fractional transformations

Here we use τ to represent each simple LFT, and we adopt the conventional circle notation for composition of functions. We also introduce a new symbol Τn to represent the composition of n + 1 transformations τi; that is,

and so forth. By direct substitution from the first set of expressions into the second we see that

and, in general,

where the last partial denominator in the finite continued fraction K is understood to be bn + z. And, since bn + 0 = bn, the image of the point z = 0 under the iterated LFT Τn is indeed the value of the finite continued fraction with n partial numerators:

A geometric interpretation

Defining a finite continued fraction as the image of a point under the iterated linear functional transformation Τn(z) leads to an intuitively appealing geometric interpretation of infinite continued fractions.

The relationship

can be understood by rewriting Τn(z) and Τn + 1(z) in terms of the fundamental recurrence formulas:

In the first of these equations the ratio tends toward An/Bn as z tends toward zero. In the second, the ratio tends toward An/Bn as z tends to infinity. This leads us to our first geometric interpretation. If the continued fraction converges, the successive convergents An/Bn are eventually arbitrarily close together. Since the linear fractional transformation Τn(z) is a continuous mapping, there must be a neighborhood of z = 0 that is mapped into an arbitrarily small neighborhood of Τn(0) = An/Bn. Similarly, there must be a neighborhood of the point at infinity which is mapped into an arbitrarily small neighborhood of Τn(∞) = An − 1/Bn − 1. So if the continued fraction converges the transformation Τn(z) maps both very small z and very large z into an arbitrarily small neighborhood of x, the value of the continued fraction, as n gets larger and larger.

For intermediate values of z, since the successive convergents are getting closer together we must have

where k is a constant, introduced for convenience. But then, by substituting in the expression for Τn(z) we obtain

so that even the intermediate values of z (except when z ≈ −k−1) are mapped into an arbitrarily small neighborhood of x, the value of the continued fraction, as n gets larger and larger. Intuitively, it is almost as if the convergent continued fraction maps the entire extended complex plane into a single point.[26]

Notice that the sequence {Τn} lies within the automorphism group of the extended complex plane, since each Τn is a linear fractional transformation for which abcd. And every member of that automorphism group maps the extended complex plane into itself: not one of the Τn can possibly map the plane into a single point. Yet in the limit the sequence {Τn} defines an infinite continued fraction which (if it converges) represents a single point in the complex plane.

When an infinite continued fraction converges, the corresponding sequence {Τn} of LFTs "focuses" the plane in the direction of x, the value of the continued fraction. At each stage of the process a larger and larger region of the plane is mapped into a neighborhood of x, and the smaller and smaller region of the plane that's left over is stretched out ever more thinly to cover everything outside that neighborhood.[27]

For divergent continued fractions, we can distinguish three cases:

  1. The two sequences {Τ2n − 1} and {Τ2n} might themselves define two convergent continued fractions that have two different values, xodd and xeven. In this case the continued fraction defined by the sequence {Τn} diverges by oscillation between two distinct limit points. And in fact this idea can be generalized: sequences {Τn} can be constructed that oscillate among three, or four, or indeed any number of limit points. Interesting instances of this case arise when the sequence {Τn} constitutes a subgroup of finite order within the group of automorphisms over the extended complex plane.
  2. The sequence {Τn} may produce an infinite number of zero denominators Bi while also producing a subsequence of finite convergents. These finite convergents may not repeat themselves or fall into a recognizable oscillating pattern. Or they may converge to a finite limit, or even oscillate among multiple finite limits. No matter how the finite convergents behave, the continued fraction defined by the sequence {Τn} diverges by oscillation with the point at infinity in this case.[28]
  3. The sequence {Τn} may produce no more than a finite number of zero denominators Bi. while the subsequence of finite convergents dances wildly around the plane in a pattern that never repeats itself and never approaches any finite limit either.

Interesting examples of cases 1 and 3 can be constructed by studying the simple continued fraction

where z is any real number such that z < −1/4.[29]

Euler's continued fraction formula

Euler proved the following identity:[18]

From this many other results can be derived, such as

and

Euler's formula connecting continued fractions and series is the motivation for the fundamental inequalities[link or clarification needed], and also the basis of elementary approaches to the convergence problem.

Examples

Transcendental functions and numbers

Here are two continued fractions that can be built via Euler's identity.

Here are additional generalized continued fractions:

This last is based on an algorithm derived by Aleksei Nikolaevich Khovansky in the 1970s.[30]

Example: the natural logarithm of 2 (= [0; 1, 2, 3, 1, 5, 2/3, 7, 1/2, 9, 2/5,..., 2k − 1, 2/k,...] ≈ 0.693147...):[31]

π

Here are three of π's best-known generalized continued fractions, the first and third of which are derived from their respective arctangent formulas above by setting x = y = 1 and multiplying by 4. The Leibniz formula for π:

converges too slowly, requiring roughly 3 × 10n terms to achieve n correct decimal places. The series derived by Nilakantha Somayaji:

is a much more obvious expression but still converges quite slowly, requiring nearly 50 terms for five decimals and nearly 120 for six. Both converge sublinearly to π. On the other hand:

converges linearly to π, adding at least three digits of precision per four terms, a pace slightly faster than the arcsine formula for π:

which adds at least three decimal digits per five terms.[32]

  • Note: this continued fraction's rate of convergence μ tends to 3 − 8 ≈ 0.1715729, hence 1/μ tends to 3 + 8 ≈ 5.828427, whose common logarithm is 0.7655... ≈ 13/17 > 3/4. The same 1/μ = 3 + 8 (the silver ratio squared) also is observed in the unfolded general continued fractions of both the natural logarithm of 2 and the nth root of 2 (which works for any integer n > 1) if calculated using 2 = 1 + 1. For the folded general continued fractions of both expressions, the rate convergence μ = (3 − 8)2 = 17 − 288 ≈ 0.02943725, hence 1/μ = (3 + 8)2 = 17 + 288 ≈ 33.97056, whose common logarithm is 1.531... ≈ 26/17 > 3/2, thus adding at least three digits per two terms. This is because the folded GCF folds each pair of fractions from the unfolded GCF into one fraction, thus doubling the convergence pace. The Manny Sardina reference further explains "folded" continued fractions.
  • Note: Using the continued fraction for arctan x/y cited above with the best-known Machin-like formula provides an even more rapidly, although still linearly, converging expression:

with u = 5 and v = 239.

Roots of positive numbers

The nth root of any positive number zm can be expressed by restating z = xn + y, resulting in

which can be simplified, by folding each pair of fractions into one fraction, to

The square root of z is a special case with m = 1 and n = 2:

which can be simplified by noting that 5/10 = 3/6 = 1/2:

The square root can also be expressed by a periodic continued fraction, but the above form converges more quickly with the proper x and y.

Example 1

The cube root of two (21/3 or 32 ≈ 1.259921...) can be calculated in two ways:

Firstly, "standard notation" of x = 1, y = 1, and 2zy = 3:

Secondly, a rapid convergence with x = 5, y = 3 and 2zy = 253:

Example 2

Pogson's ratio (1001/5 or 5100 ≈ 2.511886...), with x = 5, y = 75 and 2zy = 6325:

Example 3

The twelfth root of two (21/12 or 122 ≈ 1.059463...), using "standard notation":

Example 4

Equal temperament's perfect fifth (27/12 or 1227 ≈ 1.498307...), with m = 7:

With "standard notation":

A rapid convergence with x = 3, y = −7153, and 2zy = 219 + 312:

More details on this technique can be found in General Method for Extracting Roots using (Folded) Continued Fractions.

Higher dimensions

Another meaning for generalized continued fraction is a generalization to higher dimensions. For example, there is a close relationship between the simple continued fraction in canonical form for the irrational real number α, and the way lattice points in two dimensions lie to either side of the line y = αx. Generalizing this idea, one might ask about something related to lattice points in three or more dimensions. One reason to study this area is to quantify the mathematical coincidence idea; for example, for monomials in several real numbers, take the logarithmic form and consider how small it can be. Another reason is to find a possible solution to Hermite's problem.

There have been numerous attempts to construct a generalized theory. Notable efforts in this direction were made by Felix Klein (the Klein polyhedron), Georges Poitou and George Szekeres.

See also

Notes

  1. ^ http://www.britannica.com/EBchecked/topic/135043/continued-fraction
  2. ^ a b Pettofrezzo & Byrkit (1970, p. 150)
  3. ^ a b Long (1972, p. 173)
  4. ^ a b Pettofrezzo & Byrkit (1970, p. 152)
  5. ^ Weisstein, Eric W. "Periodic Continued Fraction". MathWorld.
  6. ^ Darren C. Collins, Continued Fractions, MIT Undergraduate Journal of Mathematics, [1]
  7. ^ a b M. Thill (2008), "A more precise rounding algorithm for rational numbers", Computing, 82: 189–198, doi:10.1007/s00607-008-0006-7
  8. ^ Shoemake, Ken (1995), "I.4: Rational Approximation", in Paeth, Alan W. (ed.), Graphic Gems V, San Diego, California: Academic Press, pp. 25–31, ISBN 0-12-543455-3
  9. ^ Theorem 193: Hardy, G.H.; Wright, E.M. (1979). An Introduction to the Theory of Numbers (Fifth ed.). Oxford.
  10. ^ Ben Thurston, "Estimating square roots, generalized continued fraction expression for every square root", The Ben Paul Thurston Blog
  11. ^ Martin, Richard M. (2004), Electronic Structure: Basic Theory and Practical Methods, Cambridge University Press, p. 557, ISBN 9781139643658.
  12. ^ Cusick & Flahive 1989.
  13. ^ Chrystal 1999.
  14. ^ Jones & Thron 1980, p. 20.
  15. ^ Euclid (2008) - The Euclidean algorithm generates a continued fraction as a by-product.
  16. ^ Cataldi 1613.
  17. ^ Wallis 1699.
  18. ^ a b Euler 1748, Chapter 18.
  19. ^ Havil 2012, pp. 104–105.
  20. ^ Brahmagupta (598–670) was the first mathematician to make a systematic study of Pell's equation.
  21. ^ Gauss 1813.
  22. ^ Lorentzen & Waadeland 1992.
  23. ^ Oskar Perron derives even more general extension and contraction formulas for continued fractions. See Perron (1977a), Perron (1977b).
  24. ^ Angell 2021.
  25. ^ Porubský 2008.
  26. ^ This intuitive interpretation is not rigorous because an infinite continued fraction is not a mapping: it is the limit of a sequence of mappings. This construction of an infinite continued fraction is roughly analogous to the construction of an irrational number as the limit of a Cauchy sequence of rational numbers.
  27. ^ Because of analogies like this one, the theory of conformal mapping is sometimes described as "rubber sheet geometry".
  28. ^ One approach to the convergence problem is to construct positive definite continued fractions, for which the denominators Bi are never zero.
  29. ^ This periodic fraction of period one is discussed more fully in the article convergence problem.
  30. ^ An alternative way to calculate log(x)
  31. ^ Borwein, Crandall & Fee 2004, p. 278, 280.
  32. ^ Beckmann 1971.

References

  • Angell, David (2021). Irrationality and Transcendence in Number Theory. Chapman and Hall/CRC. ISBN 9780367628376.
  • Chrystal, George (1999). Algebra, an Elementary Text-book for the Higher Classes of Secondary Schools and for Colleges: Pt. 1. American Mathematical Society. p. 500. ISBN 0-8218-1649-7.
  • Lorentzen, Lisa; Waadeland, Haakon (1992). Continued Fractions with Applications. Reading, MA: North Holland. ISBN 978-0-444-89265-2. (Covers primarily analytic theory and some arithmetic theory.)
  • Perron, Oskar (1977a) [1954]. Die Lehre von den Kettenbrüchen. Vol. Band I: Elementare Kettenbrüche (3 ed.). Vieweg + Teubner Verlag. ISBN 9783519020219.
  • Perron, Oskar (1977b) [1954]. Die Lehre von den Kettenbrüchen. Vol. Band II: Analytisch-funktionentheoretische Kettenbrüche (3 ed.). Vieweg + Teubner Verlag. ISBN 9783519020226.
  • Porubský, Štefan (2008). "Basic definitions for continued fractions". Interactive Information Portal for Algorithmic Mathematics. Prague, Czech Republic: Institute of Computer Science of the Czech Academy of Sciences. Retrieved 2 May 2022.
  • Szekeres, George (1970). "Multidimensional continued fractions". Ann. Univ. Sci. Budapest. Eötvös Sect. Math. 13: 113–140.
  • Wall, Hubert Stanley (1967). Analytic Theory of Continued Fractions (Reprint ed.). Chelsea Pub Co. ISBN 0-8284-0207-8. (This reprint of the D. Van Nostrand edition of 1948 covers both history and analytic theory.)
  • Wallis, John (1699). Opera mathematica [Mathematical Works].

History of continued fractions

  • 300 BC Euclid's Elements contains an algorithm for the greatest common divisor which generates a continued fraction as a by-product
  • 499 The Aryabhatiya contains the solution of indeterminate equations using continued fractions
  • 1579 Rafael Bombelli, L'Algebra Opera – method for the extraction of square roots which is related to continued fractions
  • 1613 Pietro Cataldi, Trattato del modo brevissimo di trovar la radice quadra delli numeri – first notation for continued fractions
Cataldi represented a continued fraction as a0 & n1/d1. & n2/d2. & n3/d3 with the dots indicating where the following fractions went.

See also

Notes

  1. ^ Sandifer, Ed (February 2006). "How Euler Did It: Who proved e is irrational?" (PDF). MAA Online.
  2. ^ "E101 – Introductio in analysin infinitorum, volume 1". Retrieved 2008-03-16.

References

  • Jones, William B.; Thron, W. J. (1980). Continued Fractions: Analytic Theory and Applications. Encyclopedia of Mathematics and its Applications. Vol. 11. Reading. Massachusetts: Addison-Wesley Publishing Company. ISBN 0-201-13510-8.
  • Khinchin, A. Ya. (1964) [Originally published in Russian, 1935]. Continued Fractions. University of Chicago Press. ISBN 0-486-69630-8.
  • Long, Calvin T. (1972), Elementary Introduction to Number Theory (2nd ed.), Lexington: D. C. Heath and Company, LCCN 77-171950
  • Oskar Perron, Die Lehre von den Kettenbrüchen, Chelsea Publishing Company, New York, NY 1950.
  • Pettofrezzo, Anthony J.; Byrkit, Donald R. (1970), Elements of Number Theory, Englewood Cliffs: Prentice Hall, LCCN 77-81766
  • Rockett, Andrew M.; Szüsz, Peter (1992). Continued Fractions. World Scientific Press. ISBN 981-02-1047-7.
  • H. S. Wall, Analytic Theory of Continued Fractions, D. Van Nostrand Company, Inc., 1948 ISBN 0-8284-0207-8
  • A. Cuyt, V. Brevik Petersen, B. Verdonk, H. Waadeland, W.B. Jones, Handbook of Continued fractions for Special functions, Springer Verlag, 2008 ISBN 978-1-4020-6948-2
  • Rieger, G. J. A new approach to the real numbers (motivated by continued fractions). Abh. Braunschweig.Wiss. Ges. 33 (1982), 205–217

Warning: Default sort key "Continued Fraction" overrides earlier default sort key "Generalized Continued Fraction".