# Geometric series

In mathematics, a geometric series is the sum of an infinite number of terms that have a constant ratio between successive terms. For example, the series

The geometric series 1/4 + 1/16 + 1/64 + 1/256 + ... shown as areas of purple squares. Each of the purple squares has 1/4 of the area of the next larger square (1/2×1/2 = 1/4, 1/4×1/4 = 1/16, etc.). The sum of the areas of the purple squares is one third of the area of the large square.
Another geometric series (coefficient a = 4/9 and common ratio r = 1/9) shown as areas of purple squares. The total purple area is S = a / (1 - r) = (4/9) / (1 - (1/9)) = 1/2, which can be confirmed by observing that the unit square is partitioned into an infinite number of L-shaped areas each with four purple squares and four yellow squares, which is half purple.
${\displaystyle {\frac {1}{2}}\,+\,{\frac {1}{4}}\,+\,{\frac {1}{8}}\,+\,{\frac {1}{16}}\,+\,\cdots }$

is geometric, because each successive term can be obtained by multiplying the previous term by ${\displaystyle 1/2}$. In general, a geometric series is written as ${\displaystyle a+ar+ar^{2}+ar^{3}+...}$, where ${\displaystyle a}$ is the coefficient of each term and ${\displaystyle r}$ is the common ratio between adjacent terms. Geometric series are among the simplest examples of infinite series and can serve as an introduction to Taylor series and Fourier series. Geometric series had an important role in the early development of calculus, are used throughout mathematics, and have important applications in physics, engineering, biology, economics, computer science, queueing theory, and finance.

The name geometric series indicates each term is the geometric mean of its two neighboring terms, similar to how the name arithmetic series indicates each term is the arithmetic mean of its two neighboring terms. The sequence of geometric series terms (without any of the additions) is called a geometric sequence or, equivalently, a geometric progression.

## Coefficient a

The first nine terms of geometric series 1 + r + r2 + r3 ... drawn as functions (colored in the order red, green, blue, red, green, blue, ...) within the range |r| < 1. The closed form geometric series 1 / (1 - r) is the black dashed line.

The geometric series a + ar + ar2 + ar3 + ... is written in expanded form.[1] Every coefficient in the geometric series is the same. In contrast, the power series written as a0 + a1r + a2r2 + a3r3 + ... in expanded form has coefficients ai that can vary from term to term. In other words, the geometric series is a special case of the power series. The first term of a geometric series in expanded form is the coefficient a of that geometric series.

In addition to the expanded form of the geometric series, there is a generator form[1] of the geometric series written as

${\displaystyle \sum _{k=0}^{\infty }}$  ark

and a closed form of the geometric series written as

a / (1 - r) within the range |r| < 1.

The derivation of the closed form from the expanded form is shown in this article's Sum section. However even without that derivation, the result can be confirmed with long division: a divided by (1 - r) results in a + ar + ar2 + ar3 + ... , which is the expanded form of the geometric series.

It is often a convenience in notation to set the series equal to the sum s and work with the geometric series

s = a + ar + ar2 + ar3 + ar4 + ... in its normalized form
s / a = 1 + r + r2 + r3 + r4 + ... or in its normalized vector form
s / a = [1 1 1 1 1 ...][1 r r2 r3 r4 ...]T or in its normalized partial series form
sn / a = 1 + r + r2 + r3 + r4 + ... + rn, where n is the power (or degree) of the last term included in the partial sum sn.

Typically a geometric series is thought of as a sum of numbers a + ar + ar2 + ar3 + ... but can also be thought of as a sum of functions a + ar + ar2 + ar3 + ... that converges to the function a / (1 - r) within the range |r| < 1. The adjacent image shows the contribution each of the first nine terms (i.e., functions) make to the function a / (1 - r) within the range |r| < 1 when a = 1. Changing even one of the coefficients to something other than coefficient a would (in addition to changing the geometric series to a power series) change the resulting sum of functions to some function other than a / (1 - r) within the range |r| < 1. As an aside, a particularly useful change to the coefficients is defined by the Taylor series, which describes how to change the coefficients so that the sum of functions converges to any user selected, sufficiently smooth function within a range.

## Common ratio r

The convergence of the geometric series with r=1/2 and a=1/2

The convergence of the geometric series with r=1/2 and a=1

Close-up view of the cumulative sum of functions within the range -1 < r < -0.5 as the first 11 terms of the geometric series 1 + r + r2 + r3 + ... are added. The geometric series 1 / (1 - r) is the red dashed line.
Complex geometric series (coefficient a = 1 and common ratio r = 0.5 e0t) converging to a circle. In the animation, each term of the geometric series is drawn as a vector twice: once at the origin and again within the head-to-tail vector summation that converges to the circle. The circle intersects the real axis at 2 (= 1/(1-1/2) when θ = 0) and at 2/3 (= 1/(1-(-1/2)) when θ = 180 degrees).

The geometric series a + ar + ar2 + ar3 + ... is an infinite series defined by just two parameters: coefficient a and common ratio r. Common ratio r is the ratio of any term with the previous term in the series. Or equivalently, common ratio r is the term multiplier used to calculate the next term in the series. The following table shows several geometric series:

a r Example series
4 10 4 + 40 + 400 + 4000 + 40,000 + ···
3 1 3 + 3 + 3 + 3 + 3 + ···
1 2/3 1 + 2/3 + 4/9 + 8/27 + 16/81 + ···
1/2 1/2 1/2 + 1/4 + 1/8 + 1/16 + 1/32 + ···
9 1/3 9 + 3 + 1 + 1/3 + 1/9 + ···
7 1/10 7 + 0.7 + 0.07 + 0.007 + 0.0007 + ···
1 −1/2 1 − 1/2 + 1/4 − 1/8 + 1/16 − 1/32 + ···
3 −1 3 − 3 + 3 − 3 + 3 − ···

The convergence of the geometric series depends on the value of the common ratio r:

• If |r| < 1, the terms of the series approach zero in the limit (becoming smaller and smaller in magnitude), and the series converges to the sum a / (1 - r).
• If |r| = 1, the series does not converge. When r = 1, all of the terms of the series are the same and the series is infinite. When r = −1, the terms take two values alternately (for example, 2, −2, 2, −2, 2,... ). The sum of the terms oscillates between two values (for example, 2, 0, 2, 0, 2,... ). This is a different type of divergence. See for example Grandi's series: 1 − 1 + 1 − 1 + ···.
• If |r| > 1, the terms of the series become larger and larger in magnitude. The sum of the terms also gets larger and larger, and the series does not converge to a sum. (The series diverges.)

The rate of convergence also depends on the value of the common ratio r. Specifically, the rate of convergence gets slower as r approaches 1 or −1. For example, the geometric series with a = 1 is 1 + r + r2 + r3 + ... and converges to 1 / (1 - r) when |r| < 1. However, the number of terms needed to converge approaches infinity as r approaches 1 because a / (1 - r) approaches infinity and each term of the series is less than or equal to one. In contrast, as r approaches −1 the sum of the first several terms of the geometric series starts to converge to 1/2 but slightly flips up or down depending on whether the most recently added term has a power of r that is even or odd. That flipping behavior near r = −1 is illustrated in the adjacent image showing the first 11 terms of the geometric series with a = 1 and |r| < 1.

The common ratio r and the coefficient a also define the geometric progression, which is a list of the terms of the geometric series but without the additions. Therefore the geometric series a + ar + ar2 + ar3 + ... has the geometric progression (also called the geometric sequence) a, ar, ar2, ar3, ... The geometric progression - as simple as it is - models a surprising number of natural phenomena,

As an aside, the common ratio r can be a complex number such as |r|eiθ where |r| is the vector's magnitude (or length), θ is the vector's angle (or orientation) in the complex plane and i2 = -1. With a common ratio |r|eiθ, the expanded form of the geometric series is a + a|r|eiθ + a|r|2ei2θ + a|r|3ei3θ + ... Modeling the angle θ as linearly increasing over time at the rate of some angular frequency ω0 (in other words, making the substitution θ = ω0t), the expanded form of the geometric series becomes a + a|r|eiω0t + a|r|2ei2ω0t + a|r|3ei3ω0t + ... , where the first term is a vector of length a not rotating at all, and all the other terms are vectors of different lengths rotating at harmonics of the fundamental angular frequency ω0. The constraint |r|<1 is enough to coordinate this infinite number of vectors of different lengths all rotating at different speeds into tracing a circle, as shown in the adjacent video. Similar to how the Taylor series describes how to change the coefficients so the series converges to a user selected sufficiently smooth function within a range, the Fourier series describes how to change the coefficients (which can also be complex numbers in order to specify the initial angles of vectors) so the series converges to a user selected periodic function.

## Sum

### Closed-form formula

As an alternative to the algebraic derivation of the geometric series closed form formula, there is also the following geometric derivation. (TOP) Represent the first n+1 terms of sn/a as areas of overlapped similar triangles. For example, the area of the biggest overlapped (red) triangle is bh/2 = (2)(1)/2 = 1, which is the value of the first term of the geometric series. The area of the second biggest overlapped (green) triangle is bh/2 = (2r1/2)(r1/2)/2 = r, which is the value of the second term of the geometric series. Each progressively smaller triangle has its base and height scaled down by another factor of r1/2, resulting in a sequence of triangle areas 1, r, r2, r3, ... which is equal to the sequence of terms in the normalized geometric series. (MIDDLE) In the order from largest to smallest, remove each triangle's overlapped area, which is always a fraction r of its area, and scale the remaining 1−r of the triangle's non-overlapped area by 1/(1−r) so the area of the formerly overlapped triangle, now the area of a non-overlapped trapezoid, remains the same. (BOTTOM) Aggregate the resulting n+1 non-overlapped trapezoids into a single non-overlapped trapezoid and calculate its area. The area of that aggregated trapezoid represents the value of the partial series. That area is equal to the outermost triangle minus the empty triangle tip: sn/a = (1−rn+1) / (1−r), which simplifies to s/a = 1/(1−r) when n approaches infinity and |r| < 1.

For ${\displaystyle r\neq 1}$ , the sum of the first n+1 terms of a geometric series, up to and including the r n term, is

${\displaystyle s_{n}=a+ar+ar^{2}+ar^{3}+\cdots +ar^{n}=\sum _{k=0}^{n}ar^{k}=a\left({\frac {1-r^{n+1}}{1-r}}\right),}$

where r is the common ratio. One can derive that closed-form formula for the partial sum, sn, by subtracting out the many self-similar terms as follows:[2][3][4]

{\displaystyle {\begin{aligned}s_{n}=a\ +\ &ar\ +\ ar^{2}\ +\ ar^{3}\ +\ \cdots \ +\ ar^{n},\\rs_{n}=\ &ar\ +\ ar^{2}\ +\ ar^{3}\ +\ \cdots \ +\ ar^{n}\ +\ ar^{n+1},\\s_{n}-rs_{n}=\ &a\ -\ ar^{n+1},\\s_{n}(1-r)=\ &a(1-r^{n+1}),\\s_{n}=\ &a\left({\frac {1-r^{n+1}}{1-r}}\right)\quad {\text{(if }}r\neq 1{\text{)}}.\end{aligned}}}

As n approaches infinity, the absolute value of r must be less than one for the series to converge. The sum then becomes

${\displaystyle s=a+ar+ar^{2}+ar^{3}+ar^{4}+\cdots =\sum _{k=0}^{\infty }ar^{k}={\frac {a}{1-r}},{\text{ for }}|r|<1.}$

The formula also holds for complex r, with the corresponding restriction that the modulus of r is strictly less than one.

As an aside, the question of whether an infinite series converges is fundamentally a question about the distance between two values: given enough terms, does the value of the partial sum get arbitrarily close to the finite value it is approaching? In the above derivation of the closed form of the geometric series, the interpretation of the distance between two values is the distance between their locations on the number line. That is the most common interpretation of the distance between two values. However, the p-adic metric, which has become a critical notion in modern number theory, offers a definition of distance such that the geometric series 1 + 2 + 4 + 8 + ... with a = 1 and r = 2 actually does converge to a / (1 - r) = 1 / (1 - 2) = -1 even though r is outside the typical convergence range |r| < 1.

### Proof of convergence

We can prove that the geometric series converges using the sum formula for a geometric progression:

{\displaystyle {\begin{aligned}1+r+r^{2}+r^{3}+\cdots \ &=\lim _{n\rightarrow \infty }\left(1+r+r^{2}+\cdots +r^{n}\right)\\&=\lim _{n\rightarrow \infty }{\frac {1-r^{n+1}}{1-r}}.\end{aligned}}}

The second equality is true because if ${\displaystyle |r|<1,}$  then ${\displaystyle r^{n+1}\to 0}$  as ${\displaystyle n\to \infty }$  and

{\displaystyle {\begin{aligned}(1+r+r^{2}+\cdots +r^{n})(1-r)&=((1-r)+(r-r^{2})+(r^{2}-r^{3})+...+(r^{n}-r^{n+1}))\\&=(1+(-r+r)+(-r^{2}+r^{2})+...+(-r^{n}+r^{n})-r^{n+1})\\&=1-r^{n+1}.\end{aligned}}}

Alternatively, a geometric interpretation of the convergence is shown in the adjacent diagram. The area of the white triangle is the series remainder = s - sn = arn+1 / (1 - r). Each additional term in the partial series reduces the area of that white triangle remainder by the area of the trapezoid representing the added term. The trapezoid areas (i.e., the values of the terms) get progressively thinner and shorter and closer to the origin. In the limit, as the number of trapezoids approaches infinity, the white triangle remainder vanishes as it is filled by trapezoids and therefore sn converges to s, provided |r|<1. In contrast, if |r|>1, the trapezoid areas representing the terms of the series instead get progressively wider and taller and farther from the origin, not converging to the origin and not converging as a series.

### Rate of convergence

The cumulative sum of functions within the range -1 < r < -0.5 as the first 51 terms of the geometric series 1 + r + r2 + r3 + ... are added. The magnitude of slope at r = -1 increases very gradually with each added term. The geometric series 1 / (1 - r) is the red dashed line.

As shown in the above proofs, the closed form of the geometric series partial sum up to and including the n-th power of r is a(1 - rn+1) / (1 - r) for any value of r, and the closed form of the geometric series is the full sum a / (1 - r) within the range |r| < 1.

If the common ratio is within the range 0 < r < 1, then the partial sum a(1 - rn+1) / (1 - r) increases with each added term and eventually gets within some small error, E, ratio of the full sum a / (1 - r). Solving for n at that error threshold,

{\displaystyle {\begin{aligned}a\left({\frac {1-r^{n+1}}{1-r}}\right)&\geq a\left({\frac {1-E}{1-r}}\right),\\1-r^{n+1}&\geq 1-E,\\r^{n+1}&\leq E,\\\ln(r^{n+1})&\leq \ln(E),\\(n+1)\ln(r)&\leq \ln(E),\\n+1&\geq \left\lceil {\frac {\ln(E)}{\ln(r)}}\right\rceil ,\end{aligned}}}

where 0 < r < 1, the ceiling operation ${\displaystyle \lceil \rceil }$  constrains n to integers, and dividing both sides by the natural log of r flips the inequality because it is negative. The result n+1 is the number of partial sum terms needed to get within aE / (1 - r) of the full sum a / (1 - r). For example to get within 1% of the full sum a / (1 - r) at r=0.1, only 2 (= ln(E) / ln(r) = ln(0.01) / ln(0.1)) terms of the partial sum are needed. However at r=0.9, 44 (= ln(0.01) / ln(0.9)) terms of the partial sum are needed to get within 1% of the full sum a / (1 - r).

If the common ratio is within the range -1 < r < 0, then the geometric series is an alternating series but can be converted into the form of a non-alternating geometric series by combining pairs of terms and then analyzing the rate of convergence using the same approach as shown for the common ratio range 0 < r < 1. Specifically, the partial sum

s = a + ar + ar2 + ar3 + ar4 + ar5 + ... + arn-1 + arn within the range -1 < r < 0 is equivalent to
s = a - ap + ap2 - ap3 + ap4 - ap5 + ... + apn-1 - apn with an n that is odd, with the substitution of p = -r, and within the range 0 < p < 1,
s = (a - ap) + (ap2 - ap3) + (ap4 - ap5) + ... + (apn-1 - apn) with adjacent and differently signed terms paired together,
s = a(1 - p) + a(1 - p)p2 + a(1 - p)p4 + ... + a(1 - p)p2(n-1)/2 with a(1 - p) factored out of each term,
s = a(1 - p) + a(1 - p)p2 + a(1 - p)p4 + ... + a(1 - p)p2m with the substitution m = (n - 1) / 2 which is an integer given the constraint that n is odd,

which is now in the form of the first m terms of a geometric series with coefficient a(1 - p) and with common ratio p2. Therefore the closed form of the partial sum is a(1 - p)(1 - p2(m+1)) / (1 - p2) which increases with each added term and eventually gets within some small error, E, ratio of the full sum a(1 - p) / (1 - p2). As before, solving for m at that error threshold,

{\displaystyle {\begin{aligned}a(1-p)\left({\frac {1-p^{2m+2}}{1-p^{2}}}\right)&\geq a(1-p)\left({\frac {1-E}{1-p^{2}}}\right),\\1-p^{2m+2}&\geq 1-E,\\p^{2m+2}&\leq E,\\\ln(p^{2m+2})&\leq \ln(E),\\(2m+2)\ln(p)&\leq \ln(E),\\m+1&\geq \left\lceil {\frac {\ln(E)}{2\ln(p)}}\right\rceil ,\\m+1&\geq \left\lceil {\frac {\ln(E)}{2\ln(-r)}}\right\rceil ,\end{aligned}}}

where 0 < p < 1 or equivalently -1 < r < 0, and the m+1 result is the number of partial sum pairs of terms needed to get within a(1 - p)E / (1 - p2) of the full sum a(1 - p) / (1 - p2). For example to get within 1% of the full sum a(1 - p) / (1 - p2) at p=0.1 or equivalently r=-0.1, only 1 (= ln(E) / (2 ln(p)) = ln(0.01) / (2 ln(0.1)) pair of terms of the partial sum are needed. However at p=0.9 or equivalently r=-0.9, 22 (= ln(0.01) / (2 ln(0.9))) pairs of terms of the partial sum are needed to get within 1% of the full sum a(1 - p) / (1 - p2). Comparing the rate of convergence for positive and negative values of r, n + 1 (the number of terms required to reach the error threshold for some positive r) is always twice as large as m + 1 (the number of term pairs required to reach the error threshold for the negative of that r) but the m + 1 refers to term pairs instead of single terms. Therefore, the rate of convergence is symmetric about r = 0, which can be a surprise given the asymmetry of a / (1 - r). One perspective that helps explain this rate of convergence symmetry is that on the r > 0 side each added term of the partial sum makes a finite contribution to the infinite sum at r = 1 while on the r < 0 side each added term makes a finite contribution to the infinite slope at r = -1.

As an aside, this type of rate of convergence analysis is particularly useful when calculating the number of Taylor series terms needed to adequately approximate some user-selected sufficiently-smooth function or when calculating the number of Fourier series terms needed to adequately approximate some user-selected periodic function.

## Historic insights

### Zeno of Elea (c.495 – c.430 BC)

2,500 years ago, Greek mathematicians had a problem with walking from one place to another. Physically, they were able to walk as well as we do today, perhaps better. Logically, however, they thought[5] that an infinitely long list of numbers greater than zero summed to infinity. Therefore, it was a paradox when Zeno of Elea pointed out that in order to walk from one place to another, you first have to walk half the distance, and then you have to walk half the remaining distance, and then you have to walk half of that remaining distance, and you continue halving the remaining distances an infinite number of times because no matter how small the remaining distance is you still have to walk the first half of it. Thus, Zeno of Elea transformed a short distance into an infinitely long list of halved remaining distances, all of which are greater than zero. And that was the problem: how can a distance be short when measured directly and also infinite when summed over its infinite list of halved remainders? The paradox revealed something was wrong with the assumption that an infinitely long list of numbers greater than zero summed to infinity.

### Euclid of Alexandria (c.300 BC)

Elements of Geometry, Book IX, Proposition 35. "If there is any multitude whatsoever of continually proportional numbers, and equal to the first is subtracted from the second and the last, then as the excess of the second to the first, so the excess of the last will be to all those before it."

For the same case of common ratio r>1, there are other geometric proofs of the geometric series partial sum closed form. One that represents terms as areas rather than lengths of line segments is summarized, as shown in the above figure, in three steps: (TOP) Represent the terms of a geometric series as the areas of overlapped similar triangles. (MIDDLE) From the largest to the smallest triangle, remove the overlapped left area portion (1/r) from the non-overlapped right area portion (1-1/r = (r-1)/r) and scale that non-overlapped trapezoid by r/(r-1) so its area is the same as the area of the original overlapped triangle. (BOTTOM) Notice that the area of the aggregate trapezoid is the area of a large triangle less the area of an empty small triangle at the large triangle's left tip. The large triangle is exactly the largest overlapped triangle scaled by r/(r-1). The empty small triangle started as a but that area was transformed into a non-overlapped scaled trapezoid leaving an empty left area portion (1/r). However, that empty triangle of area a/r must also be scaled by r/(r-1) so its slope matches the slope of all the non-overlapped scaled trapezoids. Therefore, Sn = area of large triangle - area of empty small triangle = arn+1/(r-1) - a/(r-1) = a(rn+1-1)/(r-1).

Euclid's Elements of Geometry[6] Book IX, Proposition 35, proof (of the proposition in adjacent diagram's caption):

Let AA', BC, DD', EF be any multitude whatsoever of continuously proportional numbers, beginning from the least AA'. And let BG and FH, each equal to AA', have been subtracted from BC and EF. I say that as GC is to AA', so EH is to AA', BC, DD'.

For let FK be made equal to BC, and FL to DD'. And since FK is equal to BC, of which FH is equal to BG, the remainder HK is thus equal to the remainder GC. And since as EF is to DD', so DD' to BC, and BC to AA' [Prop. 7.13], and DD' equal to FL, and BC to FK, and AA' to FH, thus as EF is to FL, so LF to FK, and FK to FH. By separation, as EL to LF, so LK to FK, and KH to FH [Props. 7.11, 7.13]. And thus as one of the leading is to one of the following, so (the sum of) all of the leading to (the sum of) all of the following [Prop. 7.12]. Thus, as KH is to FH, so EL, LK, KH to LF, FK, HF. And KH equal to CG, and FH to AA', and LF, FK, HF to DD', BC, AA'. Thus, as CG is to AA', so EH to DD', BC, AA'. Thus, as the excess of the second is to the first, so is the excess of the last is to all those before it. The very thing it was required to show.

The terseness of Euclid's propositions and proofs may have been a necessity. As is, the Elements of Geometry is over 500 pages of propositions and proofs. Making copies of this popular textbook was labor intensive given that the printing press was not invented until 1440. And the book's popularity lasted a long time: as stated in the cited introduction to an English translation, Elements of Geometry "has the distinction of being the world's oldest continuously used mathematical textbook." So being very terse was being very practical. The proof of Proposition 35 in Book IX could have been even more compact if Euclid could have somehow avoided explicitly equating lengths of specific line segments from different terms in the series. For example, the contemporary notation for geometric series (i.e., a + ar + ar2 + ar3 + ... + arn) does not label specific portions of terms that are equal to each other.

Also in the cited introduction the editor comments,

Most of the theorems appearing in the Elements were not discovered by Euclid himself, but were the work of earlier Greek mathematicians such as Pythagoras (and his school), Hippocrates of Chios, Theaetetus of Athens, and Eudoxus of Cnidos. However, Euclid is generally credited with arranging these theorems in a logical manner, so as to demonstrate (admittedly, not always with the rigour demanded by modern mathematics) that they necessarily follow from five simple axioms. Euclid is also credited with devising a number of particularly ingenious proofs of previously discovered theorems (e.g., Theorem 48 in Book 1).

To help translate the proposition and proof into a form that uses current notation, a couple modifications are in the diagram. First, the four horizontal line lengths representing the values of the first four terms of a geometric series are now labeled a, ar, ar2, ar3 in the diagram's left margin. Second, new labels A' and D' are now on the first and third lines so that all the diagram's line segment names consistently specify the segment's starting point and ending point.

Here is a phrase by phrase interpretation of the proposition:

Proposition in contemporary notation
"If there is any multitude whatsoever of continually proportional numbers" Taking the first n+1 terms of a geometric series Sn = a + ar + ar2 + ar3 + ... + arn
"and equal to the first is subtracted from the second and the last" and subtracting a from ar and arn
"then as the excess of the second to the first, so the excess of the last will be to all those before it." then (ar-a) / a = (arn-a) / (a + ar + ar2 + ar3 + ... + arn-1) = (arn-a) / Sn-1, which can be rearranged to the more familiar form Sn-1 = a(rn-1) / (r-1).

Similarly, here is a sentence by sentence interpretation of the proof:

Proof in contemporary notation
"Let AA', BC, DD', EF be any multitude whatsoever of continuously proportional numbers, beginning from the least AA'." Consider the first n+1 terms of a geometric series Sn = a + ar + ar2 + ar3 + ... + arn for the case r>1 and n=3.
"And let BG and FH, each equal to AA', have been subtracted from BC and EF." Subtract a from ar and ar3.
"I say that as GC is to AA', so EH is to AA', BC, DD'." I say that (ar-a) / a = (ar3-a) / (a + ar + ar2).
"For let FK be made equal to BC, and FL to DD'."
"And since FK is equal to BC, of which FH is equal to BG, the remainder HK is thus equal to the remainder GC."
"And since as EF is to DD', so DD' to BC, and BC to AA' [Prop. 7.13], and DD' equal to FL, and BC to FK, and AA' to FH, thus as EF is to FL, so LF to FK, and FK to FH."
"By separation, as EL to LF, so LK to FK, and KH to FH [Props. 7.11, 7.13]." By separation, (ar3-ar2) / ar2 = (ar2-ar) / ar = (ar-a) / a = r-1.
"And thus as one of the leading is to one of the following, so (the sum of) all of the leading to (the sum of) all of the following [Prop. 7.12]." The sum of those numerators and the sum of those denominators form the same proportion: ((ar3-ar2) + (ar2-ar) + (ar-a)) / (ar2 + ar + a) = r-1.
"And thus as one of the leading is to one of the following, so (the sum of) all of the leading to (the sum of) all of the following [Prop. 7.12]." And this sum of equal proportions can be extended beyond (ar3-ar2) / ar2 to include all the proportions up to (arn-arn-1) / arn-1.
"Thus, as KH is to FH, so EL, LK, KH to LF, FK, HF."
"And KH equal to CG, and FH to AA', and LF, FK, HF to DD', BC, AA'."
"Thus, as CG is to AA', so EH to DD', BC, AA'."
"Thus, as the excess of the second is to the first, so is the excess of the last is to all those before it." Thus, (ar-a) / a = (ar3-a) / S2. Or more generally, (ar-a) / a = (arn-a) / Sn-1, which can be rearranged in the more common form Sn-1 = a(rn-1) / (r-1).
"The very thing it was required to show." Q.E.D.

### Archimedes of Syracuse (c.287 – c.212 BC)

Archimedes' dissection of a parabolic segment into infinitely many triangles

Archimedes used the sum of a geometric series to compute the area enclosed by a parabola and a straight line. His method was to dissect the area into an infinite number of triangles.

Archimedes' Theorem states that the total area under the parabola is 4/3 of the area of the blue triangle.

Archimedes determined that each green triangle has 1/8 the area of the blue triangle, each yellow triangle has 1/8 the area of a green triangle, and so forth.

Assuming that the blue triangle has area 1, the total area is an infinite sum:

${\displaystyle 1\,+\,2\left({\frac {1}{8}}\right)\,+\,4\left({\frac {1}{8}}\right)^{2}\,+\,8\left({\frac {1}{8}}\right)^{3}\,+\,\cdots .}$

The first term represents the area of the blue triangle, the second term the areas of the two green triangles, the third term the areas of the four yellow triangles, and so on. Simplifying the fractions gives

${\displaystyle 1\,+\,{\frac {1}{4}}\,+\,{\frac {1}{16}}\,+\,{\frac {1}{64}}\,+\,\cdots .}$

This is a geometric series with common ratio 1/4 and the fractional part is equal to

${\displaystyle \sum _{n=0}^{\infty }4^{-n}=1+4^{-1}+4^{-2}+4^{-3}+\cdots ={4 \over 3}.}$

The sum is

${\displaystyle {\frac {1}{1-r}}\;=\;{\frac {1}{1-{\frac {1}{4}}}}\;=\;{\frac {4}{3}}.}$

This computation uses the method of exhaustion, an early version of integration. Using calculus, the same area could be found by a definite integral.

### Nicole Oresme (c.1323 – 1382)

A two dimensional geometric series diagram Nicole Oresme used to determine that the infinite series 1/2 + 2/4 + 3/8 + 4/16 + 5/32 + 6/64 + 7/128 + ... converges to 2.

Among his insights into infinite series, in addition to his elegantly simple proof of the divergence of the harmonic series, Nicole Oresme[7] proved that the series 1/2 + 2/4 + 3/8 + 4/16 + 5/32 + 6/64 + 7/128 + ... converges to 2. His diagram for his geometric proof, similar to the adjacent diagram, shows a two dimensional geometric series. The first dimension is horizontal, in the bottom row showing the geometric series S = 1/2 + 1/4 + 1/8 + 1/16 + ... , which is the geometric series with coefficient a = 1/2 and common ratio r = 1/2 that converges to S = a / (1-r) = (1/2) / (1-1/2) = 1. The second dimension is vertical, where the bottom row is a new coefficient aT equal to S and each subsequent row above it is scaled by the same common ratio r = 1/2, making another geometric series T = 1 + 1/2 + 1/4 + 1/8 + ... , which is the geometric series with coefficient aT = S = 1 and common ratio r = 1/2 that converges to T = aT / (1-r) = S / (1-r) = a / (1-r) / (1-r) = (1/2) / (1-1/2) / (1-1/2) = 2.

Although difficult to visualize beyond three dimensions, Oresme's insight generalizes to any dimension d. Using the sum of the d−1 dimension of the geometric series as the coefficient a in the d dimension of the geometric series results in a d-dimensional geometric series converging to Sd / a = 1 / (1-r)d within the range |r|<1. Pascal's triangle and long division reveals the coefficients of these multi-dimensional geometric series, where the closed form is valid only within the range |r|<1.

${\displaystyle {\begin{matrix}{\text{ 1}}\\{\text{ 1}}\quad {\text{ 1}}\\{\text{ 1}}\quad {\text{ 2}}\quad {\text{ 1}}\\{\text{ 1}}\quad {\text{ 3}}\quad {\text{ 3}}\quad {\text{ 1}}\\{\text{ 1}}\quad {\text{ 4}}\quad {\text{ 6}}\quad {\text{ 4}}\quad {\text{ 1}}\end{matrix}}}$

{\displaystyle {\begin{aligned}&d\quad S^{d}/a\ {\text{(closed form)}}\quad &&S^{d}/a\ {\text{(expanded form)}}\\&1\quad 1/(1-r)\quad &&1+r+r^{2}+r^{3}+r^{4}+\cdots \\&2\quad 1/(1-r)^{2}\quad &&1+2r+3r^{2}+4r^{3}+5r^{4}+\cdots \\&3\quad 1/(1-r)^{3}\quad &&1+3r+6r^{2}+10r^{3}+15r^{4}+\cdots \\&4\quad 1/(1-r)^{4}\quad &&1+4r+10r^{2}+20r^{3}+35r^{4}+\cdots \\\end{aligned}}}

Note that as an alternative to long division, it is also possible to calculate the coefficients of the d-dimensional geometric series by integrating the coefficients of dimension d−1. This mapping from division by 1-r in the power series sum domain to integration in the power series coefficient domain is a discrete form of the mapping performed by the Laplace transform. MIT Professor Arthur Mattuck shows how to derive the Laplace transform from the power series in this lecture video,[8] where the power series is a mapping between discrete coefficients and a sum and the Laplace transform is a mapping between continuous weights and an integral.

The closed forms of Sd/a are related to but not equal to the derivatives of S = f(r) = 1 / (1-r). As shown in the following table, the relationship is Sk+1 = f(k)(r) / k!, where f(k)(r) denotes the kth derivative of f(r) = 1 / (1-r) and the closed form is valid only within the range |r| < 1.

{\displaystyle {\begin{aligned}&k\quad f^{(k)}(r)/a\ {\text{(closed form)}}\quad &&f^{(k)}(r)/a\ {\text{(generator form)}}\\&0\quad 1/(1-r)\quad &&\sum _{j=0}^{\infty }r^{j}\\&1\quad 1!/(1-r)^{2}\quad &&\sum _{j=1}^{\infty }jr^{j-1}\\&2\quad 2!/(1-r)^{3}\quad &&\sum _{j=2}^{\infty }j(j-1)r^{j-2}\\&3\quad 3!/(1-r)^{4}\quad &&\sum _{j=3}^{\infty }j(j-1)(j-2)r^{j-3}\\&4\quad 4!/(1-r)^{5}\quad &&\sum _{j=4}^{\infty }j(j-1)(j-2)(j-3)r^{j-4}\\\end{aligned}}}

## Applications

### Economics

In economics, geometric series are used to represent the present value of an annuity (a sum of money to be paid in regular intervals).

For example, suppose that a payment of $100 will be made to the owner of the annuity once per year (at the end of the year) in perpetuity. Receiving$100 a year from now is worth less than an immediate $100, because one cannot invest the money until one receives it. In particular, the present value of$100 one year in the future is $100 / (1 + ${\displaystyle I}$ ), where ${\displaystyle I}$ is the yearly interest rate. Similarly, a payment of$100 two years in the future has a present value of $100 / (1 + ${\displaystyle I}$ )2 (squared because two years' worth of interest is lost by not receiving the money right now). Therefore, the present value of receiving$100 per year in perpetuity is

${\displaystyle \sum _{n=1}^{\infty }{\frac {\100}{(1+I)^{n}}},}$

which is the infinite series:

${\displaystyle {\frac {\100}{(1+I)}}\,+\,{\frac {\100}{(1+I)^{2}}}\,+\,{\frac {\100}{(1+I)^{3}}}\,+\,{\frac {\100}{(1+I)^{4}}}\,+\,\cdots .}$

This is a geometric series with common ratio 1 / (1 + ${\displaystyle I}$  ). The sum is the first term divided by (one minus the common ratio):

${\displaystyle {\frac {\100/(1+I)}{1-1/(1+I)}}\;=\;{\frac {\100}{I}}.}$

For example, if the yearly interest rate is 10% (${\displaystyle I}$  = 0.10), then the entire annuity has a present value of $100 / 0.10 =$1000.

This sort of calculation is used to compute the APR of a loan (such as a mortgage loan). It can also be used to estimate the present value of expected stock dividends, or the terminal value of a financial asset assuming a stable growth rate.

### Fractal geometry

The interior of the Koch snowflake is a union of infinitely many triangles.

In the study of fractals, geometric series often arise as the perimeter, area, or volume of a self-similar figure.

For example, the area inside the Koch snowflake can be described as the union of infinitely many equilateral triangles (see figure). Each side of the green triangle is exactly 1/3 the size of a side of the large blue triangle, and therefore has exactly 1/9 the area. Similarly, each yellow triangle has 1/9 the area of a green triangle, and so forth. Taking the blue triangle as a unit of area, the total area of the snowflake is

${\displaystyle 1\,+\,3\left({\frac {1}{9}}\right)\,+\,12\left({\frac {1}{9}}\right)^{2}\,+\,48\left({\frac {1}{9}}\right)^{3}\,+\,\cdots .}$

The first term of this series represents the area of the blue triangle, the second term the total area of the three green triangles, the third term the total area of the twelve yellow triangles, and so forth. Excluding the initial 1, this series is geometric with constant ratio r = 4/9. The first term of the geometric series is a = 3(1/9) = 1/3, so the sum is

${\displaystyle 1\,+\,{\frac {a}{1-r}}\;=\;1\,+\,{\frac {\frac {1}{3}}{1-{\frac {4}{9}}}}\;=\;{\frac {8}{5}}.}$

Thus the Koch snowflake has 8/5 of the area of the base triangle.

### Geometric power series

The formula for a geometric series

${\displaystyle {\frac {1}{1-x}}=1+x+x^{2}+x^{3}+x^{4}+\cdots }$

can be interpreted as a power series in the Taylor's theorem sense, converging where ${\displaystyle |x|<1}$ . From this, one can extrapolate to obtain other power series. For example,

{\displaystyle {\begin{aligned}\tan ^{-1}(x)&=\int {\frac {dx}{1+x^{2}}}\\&=\int {\frac {dx}{1-(-x^{2})}}\\&=\int \left(1+\left(-x^{2}\right)+\left(-x^{2}\right)^{2}+\left(-x^{2}\right)^{3}+\cdots \right)dx\\&=\int \left(1-x^{2}+x^{4}-x^{6}+\cdots \right)dx\\&=x-{\frac {x^{3}}{3}}+{\frac {x^{5}}{5}}-{\frac {x^{7}}{7}}+\cdots \\&=\sum _{n=0}^{\infty }{\frac {(-1)^{n}}{2n+1}}x^{2n+1}.\end{aligned}}}

### Encoded numbers

Zeno of Elea's geometric series with coefficient a=1/2 and common ratio r=1/2 is the foundation of binary encoded approximations of fractions in digital computers. Concretely, the geometric series written in its normalized vector form is s/a = [1 1 1 1 1 …][1 r r2 r3 r4 …]T. Keeping the column vector of basis functions [1 r r2 r3 r4 …]T the same but generalizing the row vector [1 1 1 1 1 …] so that each entry can be either a 0 or a 1 allows for an approximate encoding of any fraction. For example, the value v = 0.34375 is encoded as  v/a = [0 1 0 1 1 0 …][1 r r2 r3 r4 …]T where coefficient a = 1/2 and common ratio r = 1/2. Typically, the row vector is written in the more compact binary form v = 0.010110 which is 0.34375 in decimal.

Similarly, the geometric series with coefficient a=1 and common ratio r=2 is the foundation for binary encoded integers in digital computers. Again, the geometric series written in its normalized vector form is s/a = [1 1 1 1 1 …][1 r r2 r3 r4 …]T. Keeping the column vector of basis functions [1 r r2 r3 r4 …]T the same but generalizing the row vector [1 1 1 1 1 …] so that each entry can be either a 0 or a 1 allows for an encoding of any integer. For example, the value v = 151 is encoded as  v/a = [1 1 1 0 1 0 0 1 0 …][1 r r2 r3 r4 r5 r6 r7 r8 …]T where coefficient a = 1 and common ratio r = 2. Typically, the row vector is written in reverse order (so that the most significant bit is first) in the more compact binary form v = …010010111 = 10010111 which is 151 in decimal.

Bit fields for encoding a 32-bit floating point number according to IEEE 754 standard.

As shown in the adjacent figure, the standard binary encoding of a 32-bit floating point number is a combination of a binary encoded integer and a binary encoded fraction, beginning at the most significant bit with

• the sign bit, followed by
• an 8-bit integer exponent field with an assumed offset of 127 (so a value of 127 represents an exponent value of 0) and with a base of 2 meaning that the exponent value specifies a bit shift of the fraction field, followed by
• a 23-bit fraction field with an assumed but not encoded 1 serving as the fraction's most significant nonzero bit which would be in bit position 23 if it were encoded.

Building upon the previous example of 0.34375 having binary encoding of 0.010110, a floating point encoding (according to the IEEE 754 standard) of 0.34375 is

• the sign bit which is 0 because the number is not negative,
• an 8-bit integer exponent field which must specify a shift that counters the 2 bit left shift to get the original binary encoding from 0.010110 to 1.0110, and that counter shift to recover the original binary encoding is a right shift of 2 bits which is specified by an exponent value of 125 (because 125 − 127 = -2 which is a right shift of 2 bits) which in binary is 0111 1101,
• a 23-bit fraction field: .0110 0000 0000 0000 0000 000.

Although encoding floating point numbers by hand like this is possible, letting a computer do it is easier and less error prone. The following Julia code confirms the hand calculated floating point encoding of the number 0.34375:

 julia> bitstring(Float32(0.34375)) "00111110101100000000000000000000" 

### Repeated decimals

Decimal fractions that have repeated patterns that continue forever (e.g., 0.7777... or 0.9999... or 0.123412341234...) can only be approximated when encoded as floating point numbers, but they can always be defined exactly as the ratio of two integers and those two integers can be calculated using the geometric series. For example, the repeated decimal fraction 0.7777... can be written as the geometric series

${\displaystyle 0.7777\ldots \;=\;{\frac {7}{10}}\,+\,{\frac {7}{10}}{\frac {1}{10}}\,+\,{\frac {7}{10}}{\frac {1}{10^{2}}}\,+\,{\frac {7}{10}}{\frac {1}{10^{3}}}\,+\,\cdots }$

where coefficient a = 7/10 and common ratio r = 1/10. The geometric series closed form reveals the two integers that specify the repeated pattern:

${\displaystyle 0.7777\ldots \;=\;{\frac {a}{1-r}}\;=\;{\frac {7/10}{1-1/10}}\;=\;{\frac {7/10}{9/10}}\;=\;{\frac {7}{9}}.}$

This approach extends beyond base-ten numbers. In fact, any fraction that has an infinitely repeated pattern in base-ten numbers also has an infinitely repeated pattern in numbers written in any other base. For example, looking at the floating point encoding for the number 0.7777...

 julia> bitstring(Float32(0.77777777777777777777)) "00111111010001110001110001110010" 

reveals the binary fraction 0.110001110001110001... where the binary pattern 0b110001 repeats indefinitely and can be written in mostly (except for the powers) binary numbers as

${\displaystyle 0.110001110001110001\ldots \;=\;{\frac {110001}{1000000}}\,+\,{\frac {110001}{1000000}}{\frac {1}{1000000}}\,+\,{\frac {110001}{1000000}}{\frac {1}{1000000^{2}}}\,+\,{\frac {110001}{1000000}}{\frac {1}{1000000^{3}}}\,+\,\cdots }$

where coefficient a = 0b110001 / 0b1000000 = 49 / 64 and common ratio r = 1 / 0b1000000 = 1 / 64. Using the geometric series closed form as before

${\displaystyle 0.7777\ldots \;=\;0b0.110001110001110001\ldots \;=\;{\frac {a}{1-r}}\;=\;{\frac {49/64}{1-1/64}}\;=\;{\frac {49/64}{63/64}}\;=\;{\frac {49}{63}}\;=\;{\frac {7}{9}}.}$

You may have noticed that the floating point encoding does not capture the 0b110001 repeat pattern in the last couple (least significant) bits. This is because floating point encoding rounds the remainder instead of truncating it. Therefore, if the most significant bit of the remainder is 1, the least significant bit of the encoded fraction gets incremented and that will cause a carry if the least significant bit of the fraction is already 1, which can cause another carry if that bit of the fraction is already a 1, which can cause another carry, etc. This floating point rounding and the subsequent carry propagation explains why the floating point encoding for 0.99999... is exactly the same as the floating point encoding for 1.

 julia> bitstring(Float32(0.99999999999999999999)) "00111111100000000000000000000000"julia> bitstring(Float32(1.0)) "00111111100000000000000000000000" 

As an example that has four digits in the repeated pattern, 0.123412341234... can be written as the geometric series

${\displaystyle 0.123412341234\ldots \;=\;{\frac {1234}{10000}}\,+\,{\frac {1234}{10000}}{\frac {1}{10000}}\,+\,{\frac {1234}{10000}}{\frac {1}{10000^{2}}}\,+\,{\frac {1234}{10000}}{\frac {1}{10000^{3}}}\,+\,\cdots }$

where coefficient a = 1234/10000 and common ratio r = 1/10000. The geometric series closed form reveals the two integers that specify the repeated pattern:

${\displaystyle 0.123412341234\ldots \;=\;{\frac {a}{1-r}}\;=\;{\frac {1234/10000}{1-1/10000}}\;=\;{\frac {1234/10000}{9999/10000}}\;=\;{\frac {1234}{9999}}.}$

### Complex Fourier series

The front and back planes show the sum of the first +/-n terms of the complex Fourier series that has coefficients set to converge to a tracing of the letter 'e' for exponential. (The Julia source code that generates the frames and audio of this animation is here[9] in Appendix B.)

As an example of the ability of the complex Fourier series to trace any 2D closed figure, in the adjacent animation a complex Fourier series traces the letter 'e' (for exponential). Given the intricate coordination of motions shown in the animation, a definition of the complex Fourier series can be surprisingly compact in just two equations:

{\displaystyle {\begin{aligned}s(t)&=\sum _{n=-\infty }^{\infty }c_{n}e^{2\pi int}\\c_{n}&=\int _{0}^{1}s(t)e^{-2\pi int}dt\\\end{aligned}}}

where parameterized function s(t) traces some 2D closed figure in the complex plane as the parameter t progresses through the period from 0 to 1.

To help make sense of these compact equations defining the complex Fourier series, note that the complex Fourier series summation looks similar to the complex geometric series except that the complex Fourier series is basically two complex geometric series (one set of terms rotating in the positive direction and another set of terms rotating in the negative direction), and the coefficients of the complex Fourier series are complex constants that can vary from term to term. By allowing terms to rotate in either direction, the series becomes capable of tracing any 2D closed figure. In contrast, the complex geometric series has all the terms rotating in the same direction and it can trace only circles. Allowing the coefficients of the complex geometric series to vary from term to term would expand upon the shapes it can trace but all the possible shapes would still be limited to being puffy and cloud-like, not able to trace the shape of a simple line segment, for example going back and forth between 1 + i0 and -1 + i0. However, Euler's formula shows that the addition of just two terms rotating in opposite directions can trace that line segment between 1 + i0 and -1 + i0:

{\displaystyle {\begin{aligned}e^{i\theta }&=\cos \theta +i\sin \theta \\e^{-i\theta }&=\cos \theta -i\sin \theta \\\cos \theta &={\frac {e^{i\theta }+e^{-i\theta }}{2}}.\\\end{aligned}}}

Concerning the complex Fourier series second equation defining how to calculate the coefficients, the coefficient of the non-rotating term c0 can be calculated by integrating the complex Fourier series first equation over the range of one period from 0 to 1. Over that range, all the rotating terms integrate to zero, leaving just c0. Similarly, any of the terms in the complex Fourier series first equation can be made to be a non-rotating term by multiplying both sides of the equation by ${\displaystyle e^{-2\pi int}}$  before integrating to calculate cn, and that is the complex Fourier series second equation.

## Notes

1. ^ a b Riddle, Douglas F. Calculus and Analytic Geometry, Second Edition Belmont, California, Wadsworth Publishing, p. 566, 1970.
2. ^ Abramowitz & Stegun (1972, p. 10)
3. ^ Moise (1967, p. 48)
4. ^ Protter & Morrey (1970, pp. 639–640)
5. ^ Riddle, Douglas E (1974). Calculus and Analytic Geometry (2nd ed.). Wadsworth Publishing. p. 556. ISBN 053400301-X.
6. ^ Euclid; J.L. Heiberg (2007). Euclid's Elements of Geometry (PDF). Translated by Richard Fitzpatrick. Richard Fitzpatrick. ISBN 978-0615179841.
7. ^ Babb, J (2003). "Mathematical Concepts and Proofs from Nicole Oresme: Using the History of Calculus to Teach Mathematics" (PDF). Department of Mathematics & Statistics, University of Winnipeg, 515 Portage Avenue: The Seventh International History, Philosophy and Science Teaching conference. pp. 11–12, 21.{{cite web}}: CS1 maint: location (link)
8. ^ Mattuck, Arthur. "Lecture 19, MIT 18.03 Differential Equations, Spring 2006". MIT OpenCourseWare. Archived from the original on 2021-11-11.
9. ^ Sepesi, G (13 February 2022). "Zeno's Enduring Example". Towards Data Science. pp. Appendix B.

## References

• Abramowitz, M. and Stegun, I. A. (Eds.). Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 9th printing. New York: Dover, p. 10, 1972.
• Andrews, George E. (1998). "The geometric series in calculus". The American Mathematical Monthly. Mathematical Association of America. 105 (1): 36–40. doi:10.2307/2589524. JSTOR 2589524.
• Arfken, G. Mathematical Methods for Physicists, 3rd ed. Orlando, FL: Academic Press, pp. 278–279, 1985.
• Beyer, W. H. CRC Standard Mathematical Tables, 28th ed. Boca Raton, FL: CRC Press, p. 8, 1987.
• Courant, R. and Robbins, H. "The Geometric Progression." §1.2.3 in What Is Mathematics?: An Elementary Approach to Ideas and Methods, 2nd ed. Oxford, England: Oxford University Press, pp. 13–14, 1996.
• James Stewart (2002). Calculus, 5th ed., Brooks Cole. ISBN 978-0-534-39339-7
• Larson, Hostetler, and Edwards (2005). Calculus with Analytic Geometry, 8th ed., Houghton Mifflin Company. ISBN 978-0-618-50298-1
• Pappas, T. "Perimeter, Area & the Infinite Series." The Joy of Mathematics. San Carlos, CA: Wide World Publ./Tetra, pp. 134–135, 1989.
• Protter, Murray H.; Morrey, Charles B. Jr. (1970), College Calculus with Analytic Geometry (2nd ed.), Reading: Addison-Wesley, LCCN 76087042
• Roger B. Nelsen (1997). Proofs without Words: Exercises in Visual Thinking, The Mathematical Association of America. ISBN 978-0-88385-700-7

### History and philosophy

• C. H. Edwards Jr. (1994). The Historical Development of the Calculus, 3rd ed., Springer. ISBN 978-0-387-94313-8.
• Swain, Gordon and Thomas Dence (April 1998). "Archimedes' Quadrature of the Parabola Revisited". Mathematics Magazine. 71 (2): 123–30. doi:10.2307/2691014. JSTOR 2691014.
• Eli Maor (1991). To Infinity and Beyond: A Cultural History of the Infinite, Princeton University Press. ISBN 978-0-691-02511-7
• Morr Lazerowitz (2000). The Structure of Metaphysics (International Library of Philosophy), Routledge. ISBN 978-0-415-22526-7

### Economics

• Carl P. Simon and Lawrence Blume (1994). Mathematics for Economists, W. W. Norton & Company. ISBN 978-0-393-95733-4
• Mike Rosser (2003). Basic Mathematics for Economists, 2nd ed., Routledge. ISBN 978-0-415-26784-7

### Biology

• Edward Batschelet (1992). Introduction to Mathematics for Life Scientists, 3rd ed., Springer. ISBN 978-0-387-09648-3
• Richard F. Burton (1998). Biology by Numbers: An Encouragement to Quantitative Thinking, Cambridge University Press. ISBN 978-0-521-57698-7

### Computer science

• John Rast Hubbard (2000). Schaum's Outline of Theory and Problems of Data Structures With Java, McGraw-Hill. ISBN 978-0-07-137870-3