Part of a series of articles about |

Calculus |
---|

In
mathematics, the **harmonic series** is the
infinite series formed by summing all positive
unit fractions:

The first terms of the series sum to approximately , where is the natural logarithm and is the Euler–Mascheroni constant. Because the logarithm has arbitrarily large values, the harmonic series does not have a finite limit: it is a divergent series. Its divergence was proven in the 14th century by Nicole Oresme using a precursor to the Cauchy condensation test for the convergence of infinite series. It can also be proven to diverge by comparing the sum to an integral, according to the integral test for convergence.

Applications of the harmonic series and its partial sums include Euler's proof that there are infinitely many prime numbers, the analysis of the coupon collector's problem on how many random trials are needed to provide a complete range of responses, the connected components of random graphs, the block-stacking problem on how far over the edge of a table a stack of blocks can be cantilevered, and the average case analysis of the quicksort algorithm.

The name of the harmonic series derives from the concept of
overtones or harmonics
in music: the
wavelengths of the overtones of a vibrating string are , , , etc., of the string's
fundamental wavelength.^{
[1]}^{
[2]} Every term of the harmonic series after the first is the
harmonic mean of the neighboring terms, so the terms form a
harmonic progression; the phrases *harmonic mean* and *harmonic progression* likewise derive from music.^{
[2]}
Beyond music, harmonic sequences have also had a certain popularity with architects. This was so particularly in the
Baroque period, when architects used them to establish the
proportions of
floor plans, of
elevations, and to establish harmonic relationships between both interior and exterior architectural details of churches and palaces.^{
[3]}

The divergence of the harmonic series was first proven in 1350 by
Nicole Oresme.^{
[2]}^{
[4]} Oresme's work, and the contemporaneous work of
Richard Swineshead on a different series, marked the first appearance of infinite series other than the
geometric series in mathematics.^{
[5]} However, this achievement fell into obscurity.^{
[6]} Additional proofs were published in the 17th century by
Pietro Mengoli^{
[2]}^{
[7]} and by
Jacob Bernoulli.^{
[8]}^{
[9]}^{
[10]} Bernoulli credited his brother
Johann Bernoulli for finding the proof,^{
[10]} and it was later included in Johann Bernoulli's collected works.^{
[11]}

The partial sums of the harmonic series were named
harmonic numbers, and given their usual notation , in 1968 by
Donald Knuth.^{
[12]}

The harmonic series is the infinite series

in which the terms are all of the positive
unit fractions. It is a
divergent series: as more terms of the series are included in
partial sums of the series, the values of these partial sums grow arbitrarily large, beyond any finite limit. Because it is a divergent series, it should be interpreted as a formal sum, an abstract mathematical expression combining the unit fractions, rather than as something that can be evaluated to a numeric value. There are many different proofs of the divergence of the harmonic series, surveyed in a 2006 paper by S. J. Kifowit and T. A. Stamps.

One way to prove divergence is to compare the harmonic series with another divergent series, where each denominator is replaced with the next-largest power of two:

Grouping equal terms shows that the second series diverges (because every grouping of convergent series is only convergent):

Because each term of the harmonic series is greater than or equal to the corresponding term of the second series (and the terms are all positive), it follows (by the
comparison test) that the harmonic series diverges as well. The same argument proves more strongly that, for every
positive
integer ,

This is the original proof given by
Nicole Oresme in around 1350.

It is possible to prove that the harmonic series diverges by comparing its sum with an improper integral. Specifically, consider the arrangement of rectangles shown in the figure to the right. Each rectangle is 1 unit wide and units high, so if the harmonic series converged then the total area of the rectangles would be the sum of the harmonic series. The curve stays entirely below the upper boundary of the rectangles, so the area under the curve (in the range of from one to infinity that is covered by rectangles) would be less than the area of the union of the rectangles. However, the area under the curve is given by a divergent improper integral,

Because this integral does not converge, the sum cannot converge either.

Replacing each rectangle by the next one in the sequence would produce a sequence of rectangles whose boundary lies below the curve rather than above it. This shows that the partial sums of the harmonic series differ from the integral by an amount that is bounded above and below by the unit area of the first rectangle:

Generalizing this argument, any infinite sum of values of a monotone decreasing positive function of (like the harmonic series) has partial sums that are within a bounded distance of the values of the corresponding integrals. Therefore, the sum converges if and only if the integral over the same range of the same function converges. When this equivalence is used to check the convergence of a sum by replacing it with an easier integral, it is known as the
integral test for convergence.

Partial sum of the harmonic series, | ||||
---|---|---|---|---|

expressed as a fraction | decimal | relative size | ||

1 | 1 | 1 | ||

2 | 3 | /2 | 1.5 | |

3 | 11 | /6 | ~1.83333 | |

4 | 25 | /12 | ~2.08333 | |

5 | 137 | /60 | ~2.28333 | |

6 | 49 | /20 | 2.45 | |

7 | 363 | /140 | ~2.59286 | |

8 | 761 | /280 | ~2.71786 | |

9 | 7129 | /2520 | ~2.82897 | |

10 | 7381 | /2520 | ~2.92897 | |

11 | 83711 | /27720 | ~3.01988 | |

12 | 86021 | /27720 | ~3.10321 | |

13 | 1145993 | /360360 | ~3.18013 | |

14 | 1171733 | /360360 | ~3.25156 | |

15 | 1195757 | /360360 | ~3.31823 | |

16 | 2436559 | /720720 | ~3.38073 | |

17 | 42142223 | /12252240 | ~3.43955 | |

18 | 14274301 | /4084080 | ~3.49511 | |

19 | 275295799 | /77597520 | ~3.54774 | |

20 | 55835135 | /15519504 | ~3.59774 |

Adding the first terms of the harmonic series produces a
partial sum, called a
harmonic number and denoted :^{
[12]}

These numbers grow very slowly, with
logarithmic growth, as can be seen from the integral test.^{
[15]} More precisely, by the
Euler–Maclaurin formula,

where is the
Euler–Mascheroni constant and which approaches 0 as goes to infinity.

No harmonic numbers are integers, except for .^{
[17]}^{
[18]} One way to prove that is not an integer is to consider the highest
power of two in the range from 1 to . If is the
least common multiple of the numbers from 1 to , then
can be rewritten as a sum of fractions with equal denominators

in which only one of the numerators, , is odd and the rest are even, and (when ) is itself even. Therefore, the result is a fraction with an odd numerator and an even denominator, which cannot be an integer.

Another proof that the harmonic numbers are not integers observes that the denominator of must be divisible by
all
prime numbers greater than , and uses
Bertrand's postulate to prove that this set of primes is non-empty. The same argument implies more strongly that, except for , , and , no harmonic number can have a
terminating decimal representation.^{
[17]} It has been conjectured that every prime number divides the numerators of only a finite subset of the harmonic numbers, but this remains unproven.^{
[19]}

The digamma function is defined as the logarithmic derivative of the gamma function

Just as the gamma function provides a continuous
interpolation of the
factorials, the digamma function provides a continuous interpolation of the harmonic numbers, in the sense that .

Many well-known mathematical problems have solutions involving the harmonic series and its partial sums.

The
jeep problem or desert-crossing problem is included in a 9th-century problem collection by
Alcuin, *
Propositiones ad Acuendos Juvenes* (formulated in terms of camels rather than jeeps), but with an incorrect solution.^{
[22]} The problem asks how far into the desert a jeep can travel and return, starting from a base with loads of fuel, by carrying some of the fuel into the desert and leaving it in depots. The optimal solution involves placing depots spaced at distances from the starting point and each other, where is the range of distance that the jeep can travel with a single load of fuel. On each trip out and back from the base, the jeep places one more depot, refueling at the other depots along the way, and placing as much fuel as it can in the newly placed depot while still leaving enough for itself to return to the previous depots and the base. Therefore, the total distance reached on the th trip is

where is the th harmonic number. The divergence of the harmonic series implies that crossings of any length are possible with enough fuel.

For instance, for Alcuin's version of the problem, : a camel can carry 30 measures of grain and can travel one leuca while eating a single measure, where a leuca is a unit of distance roughly equal to 2.3 kilometres (1.4 mi). The problem has : there are 90 measures of grain, enough to supply three trips. For the standard formulation of the desert-crossing problem, it would be possible for the camel to travel leucas and return, by placing a grain storage depot 5 leucas from the base on the first trip and 12.5 leucas from the base on the second trip. However, Alcuin instead asks a slightly different question, how much grain can be transported a distance of 30 leucas without a final return trip, and either strands some camels in the desert or fails to account for the amount of grain consumed by a camel on its return trips.^{
[22]}

In the
block-stacking problem, one must place a pile of identical rectangular blocks, one per layer, so that they hang as far as possible over the edge of a table without falling. The top block can be placed with of its length extending beyond the next lower block. If it is placed in this way, the next block down needs to be placed with at most of its length extending beyond the next lower block, so that the
center of mass of the top two block is supported and they do not topple. The third block needs to be placed with at most of its length extending beyond the next lower block, and so on. In this way, it is possible to place the blocks in such a way that they extend lengths beyond the table, where is the th harmonic number.^{
[24]}^{
[25]} The divergence of the harmonic series implies that there is no limit on how far beyond the table the block stack can extend.^{
[25]} For stacks with one block per layer, no better solution is possible, but significantly more overhang can be achieved using stacks with more than one block per layer.^{
[26]}

In 1737, Leonhard Euler observed that, as a formal sum, the harmonic series is equal to an Euler product in which each term comes from a prime number:

where denotes the set of prime numbers. The left equality comes from applying the
distributive law to the product and recognizing the resulting terms as the
prime factorizations of the terms in the harmonic series, and the right equality uses the standard formula for a
geometric series. The product is divergent, just like the sum, but if it converged one could take logarithms and obtain

Here, each logarithm is replaced by its
Taylor series, and the constant on the right is the evaluation of the convergent series of terms with exponent greater than one. It follows from these manipulations that the sum of reciprocals of primes, on the right hand of this equality, must diverge, for if it converged these steps could be reversed to show that the harmonic series also converges, which it does not. An immediate corollary is that
there are infinitely many prime numbers, because a finite sum cannot diverge.

Another problem in number theory closely related to the harmonic series concerns the average number of divisors of the numbers in a range from 1 to , formalized as the average order of the divisor function,

The operation of rounding each term in the harmonic series to the next smaller integer multiple of causes this average to differ from the harmonic numbers by a small constant, and
Peter Gustav Lejeune Dirichlet showed more precisely that the average number of divisors is (expressed in
big O notation). Bounding the final error term more precisely remains an open problem, known as
Dirichlet's divisor problem.

Several common games or recreations involve repeating a random selection from a set of items until all possible choices have been selected; these include the collection of
trading cards^{
[31]}^{
[32]} and the completion of
parkrun bingo, in which the goal is to obtain all 60 possible numbers of seconds in the times from a sequence of running events.^{
[33]} More serious applications of this problem include sampling all variations of a manufactured product for its
quality control,^{
[34]} and the
connectivity of
random graphs.^{
[35]} In situations of this form, once there are items remaining to be collected out of a total of equally-likely items, the probability of collecting a new item in a single random choice is and the expected number of random choices needed until a new item is collected is . Summing over all values of from down to 1 shows that the total expected number of random choices needed to collect all items is , where is the th harmonic number.^{
[36]}

The
quicksort algorithm for sorting a set of items can be analyzed using the harmonic numbers. The algorithm operates by choosing one item as a "pivot", comparing it to all the others, and recursively sorting the two subsets of items whose comparison places them before the pivot and after the pivot. In either its
average-case complexity (with the assumption that all input permutations are equally likely) or in its
expected time analysis of worst-case inputs with a random choice of pivot, all of the items are equally likely to be chosen as the pivot. For such cases, one can compute the probability that two items are ever compared with each other, throughout the recursion, as a function of the number of other items that separate them in the final sorted order. If items and are separated by other items, then the algorithm will make a comparison between and only when, as the recursion progresses, it picks or as a pivot before picking any of the other items between them. Because each of these items is equally likely to be chosen first, this happens with probability . The total expected number of comparisons, which controls the total running time of the algorithm, can then be calculated by summing these probabilities over all pairs, giving^{
[37]}

The divergence of the harmonic series corresponds in this application to the fact that, in the
comparison model of sorting used for quicksort, it is not possible to sort in
linear time.

The series

is known as the

Explicitly, the asymptotic expansion of the series is

Using alternating signs with only odd unit fractions produces a related series, the
Leibniz formula for π^{
[40]}

The Riemann zeta function is defined for real by the convergent series

which for would be the harmonic series. It can be extended by
analytic continuation to a
holomorphic function on all
complex numbers except , where the extended function has a
simple pole. Other important values of the zeta function include , the solution to the
Basel problem,
Apéry's constant , proved by
Roger Apéry to be an
irrational number, and the "critical line" of complex numbers with real part , conjectured by the
Riemann hypothesis to be the only values other than negative integers where the function can be zero.

The random harmonic series is

where the values are
independent and identically distributed random variables that take the two values and with equal probability . It converges
with probability 1, as can be seen by using the
Kolmogorov three-series theorem or of the closely related
Kolmogorov maximal inequality. The sum of the series is a
random variable whose
probability density function is close to for values between and , and decreases to near-zero for values greater than or less than . Intermediate between these ranges, at the values , the probability density is for a nonzero but very small value .

The depleted harmonic series where all of the terms in which the digit 9 appears anywhere in the denominator are removed can be shown to converge to the value 22.92067661926415034816....^{
[44]} In fact, when all the terms containing any particular string of digits (in any
base) are removed, the series converges.^{
[45]}

- ^
^{a}^{b}Rice, Adrian (2011). "The harmonic series: A primer". In Jardine, Dick; Shell-Gellasch, Amy (eds.).*Mathematical Time Capsules: Historical Modules for the Mathematics Classroom*. MAA Notes. Vol. 77. Washington, DC: Mathematical Association of America. pp. 269–276. ISBN 978-0-88385-984-1. - ^
^{a}^{b}^{c}^{d}Kullman, David E. (May 2001). "What's harmonic about the harmonic series?".*The College Mathematics Journal*.**32**(3): 201–203. doi: 10.2307/2687471. JSTOR 2687471. **^**Hersey, George L. (2001).*Architecture and Geometry in the Age of the Baroque*. University of Chicago Press. pp. 11–12, 37–51. ISBN 978-0-226-32783-9.**^**Oresme, Nicole (c. 1360).*Quaestiones super Geometriam Euclidis*[*Questions concerning Euclid's Geometry*] (in Latin).**^**Stillwell, John (2010).*Mathematics and its History*. Undergraduate Texts in Mathematics (3rd ed.). New York: Springer. p. 182. doi: 10.1007/978-1-4419-6053-5. ISBN 978-1-4419-6052-8. MR 2667826.**^**Derbyshire, John (2003).*Prime Obsession: Bernhard Riemann and the Greatest Unsolved Problem in Mathematics*. Washington, DC: Joseph Henry Press. p. 10. ISBN 0-309-08549-7. MR 1968857.**^**Mengoli, Pietro (1650). "Praefatio [Preface]".*Novae quadraturae arithmeticae, seu De additione fractionum*[*New arithmetic quadrature (i.e., integration), or On the addition of fractions*] (in Latin). Bologna: Giacomo Monti. Mengoli's proof is by contradiction: Let denote the sum of the series. Group the terms of the series in triplets: . Since for , , then , which is impossible for any finite . Therefore, the series diverges.**^**Bernoulli, Jacob (1689).*Propositiones arithmeticae de seriebus infinitis earumque summa finita*[*Arithmetical propositions about infinite series and their finite sums*]. Basel: J. Conrad.**^**Bernoulli, Jacob (1713).*Ars conjectandi, opus posthumum. Accedit Tractatus de seriebus infinitis*[*Theory of inference, posthumous work. With the Treatise on infinite series…*]. Basel: Thurneysen. pp. 250–251.

From p. 250, prop. 16:- "
*XVI. Summa serei infinita harmonicè progressionalium, &c. est infinita. Id primus deprehendit Frater:…*" - [16. The sum of an infinite series of harmonic progression, , is infinite. My brother first discovered this…]

- "
- ^
^{a}^{b}Dunham, William (January 1987). "The Bernoullis and the harmonic series".*The College Mathematics Journal*.**18**(1): 18–23. doi: 10.1080/07468342.1987.11973001. JSTOR 2686312. **^**Bernoulli, Johann (1742). "Corollary III of*De seriebus varia*".*Opera Omnia*. Lausanne & Basel: Marc-Michel Bousquet & Co. vol. 4, p. 8. Johann Bernoulli's proof is also by contradiction. It uses a telescopic sum to represent each term asChanging the order of summation in the corresponding double series gives, in modern notation.- ^
^{a}^{b}Knuth, Donald E. (1968). "1.2.7 Harmonic numbers".*The Art of Computer Programming, Volume I: Fundamental Algorithms*(1st ed.). Addison-Wesley. pp. 73–78. Knuth writes, of the partial sums of the harmonic series "This sum does not occur very frequently in classical mathematics, and there is no standard notation for it; but in the analysis of algorithms it pops up nearly every time we turn around, and we will consistently use the symbol ... The letter stands for "harmonic", and we call a "harmonic number" because [the infinite series] is customarily called the harmonic series." - ^
^{a}^{b}^{c}^{d}Kifowit, Steven J.; Stamps, Terra A. (Spring 2006). "The harmonic series diverges again and again" (PDF).*AMATYC Review*. American Mathematical Association of Two-Year Colleges.**27**(2): 31–43. See also unpublished addendum, " More proofs of divergence of the harmonic series" by Kifowit. **^**Roy, Ranjan (December 2007). "Review of*A Radical Approach to Real Analysis*by David M. Bressoud".*SIAM Review*.**49**(4): 717–719. JSTOR 20454048.One might point out that Cauchy's condensation test is merely the extension of Oresme's argument for the divergence of the harmonic series

- ^
^{a}^{b}Bressoud, David M. (2007).*A Radical Approach to Real Analysis*. Classroom Resource Materials Series (2nd ed.). Washington, DC: Mathematical Association of America. pp. 137–138. ISBN 978-0-88385-747-2. MR 2284828. **^**Boas, R. P. Jr.; Wrench, J. W. Jr. (1971). "Partial sums of the harmonic series".*The American Mathematical Monthly*.**78**(8): 864–870. doi: 10.1080/00029890.1971.11992881. JSTOR 2316476. MR 0289994.- ^
^{a}^{b}^{c}Havil, Julian (2003). "Chapter 2: The harmonic series".*Gamma: Exploring Euler's Constant*. Princeton University Press. pp. 21–25. ISBN 978-0-691-14133-6. - ^
^{a}^{b}Osler, Thomas J. (November 2012). "96.53 Partial sums of series that cannot be an integer".*The Mathematical Gazette*.**96**(537): 515–519. doi: 10.1017/S0025557200005167. JSTOR 24496876. S2CID 124359670. See in particular Theorem 1, p. 516. **^**Sanna, Carlo (2016). "On the -adic valuation of harmonic numbers".*Journal of Number Theory*.**166**: 41–46. doi: 10.1016/j.jnt.2016.02.020. MR 3486261.**^**Ross, Bertram (1978). "The psi function".*Mathematics Magazine*.**51**(3): 176–179. doi: 10.1080/0025570X.1978.11976704. JSTOR 2689999. MR 1572267.**^**Sofo, Anthony; Srivastava, H. M. (2015). "A family of shifted harmonic sums".*The Raman. J*.**37**: 89–108. doi: 10.1007/s11139-014-9600-9. S2CID 254990799.- ^
^{a}^{b}Hadley, John; Singmaster, David (March 1992). "Problems to sharpen the young: An annotated translation of*Propositiones ad acuendos juvenes*".*The Mathematical Gazette*.**76**(475): 102–126. doi: 10.2307/3620384. JSTOR 3620384. S2CID 125835186. See problem 52: De homine patrefamilias – A lord of the manor, pp. 124–125. **^**Gale, David (May 1970). "The jeep once more or jeeper by the dozen".*The American Mathematical Monthly*.**77**(5): 493–501. doi: 10.1080/00029890.1970.11992525. JSTOR 2317382.**^**Graham, Ronald; Knuth, Donald E.; Patashnik, Oren (1989). "6.3 Harmonic numbers".*Concrete Mathematics*(2e ed.). Addison-Wesley. pp. 272–278. ISBN 978-0-201-55802-9.- ^
^{a}^{b}Sharp, R. T. (1954). "Problem 52: Overhanging dominoes" (PDF).*Pi Mu Epsilon Journal*.**1**(10): 411–412. **^**Paterson, Mike; Peres, Yuval; Thorup, Mikkel; Winkler, Peter; Zwick, Uri (2009). "Maximum overhang".*The American Mathematical Monthly*.**116**(9): 763–787. doi: 10.4169/000298909X474855. MR 2572086. S2CID 1713091.**^**Euler, Leonhard (1737). "Variae observationes circa series infinitas" [Various observations concerning infinite series].*Commentarii Academiae Scientiarum Petropolitanae*(in Latin).**9**: 160–188.- ^
^{a}^{b}Rubinstein-Salzedo, Simon (2017). "Could Euler have conjectured the prime number theorem?".*Mathematics Magazine*.**90**(5): 355–359. arXiv: 1701.04718. doi: 10.4169/math.mag.90.5.355. JSTOR 10.4169/math.mag.90.5.355. MR 3738242. S2CID 119165483. **^**Pollack, Paul (2015). "Euler and the partial sums of the prime harmonic series".*Elemente der Mathematik*.**70**(1): 13–20. doi: 10.4171/EM/268. MR 3300350.**^**Tsang, Kai-Man (2010). "Recent progress on the Dirichlet divisor problem and the mean square of the Riemann zeta-function".*Science China*.**53**(9): 2561–2572. Bibcode: 2010ScChA..53.2561T. doi: 10.1007/s11425-010-4068-6. MR 2718848. S2CID 6168120.**^**Maunsell, F. G. (October 1938). "A problem in cartophily".*The Mathematical Gazette*.**22**(251): 328–331. doi: 10.2307/3607889. JSTOR 3607889. S2CID 126381029.**^**Gerke, Oke (April 2013). "How much is it going to cost me to complete a collection of football trading cards?".*Teaching Statistics*.**35**(2): 89–93. doi: 10.1111/test.12005. S2CID 119887116.**^**Parker, Matt (February 12, 2022). "The coupon collector's problem (with Geoff Marshall)".*Stand-up maths*. YouTube.**^**Luko, Stephen N. (March 2009). "The "coupon collector's problem" and quality control".*Quality Engineering*.**21**(2): 168–181. doi: 10.1080/08982110802642555. S2CID 109194745.**^**Frieze, Alan; Karoński, Michał (2016). "4.1 Connectivity".*Introduction to Random Graphs*. Cambridge University Press, Cambridge. pp. 64–68. doi: 10.1017/CBO9781316339831. ISBN 978-1-107-11850-8. MR 3675279.**^**Isaac, Richard (1995). "8.4 The coupon collector's problem solved".*The Pleasures of Probability*. Undergraduate Texts in Mathematics. New York: Springer-Verlag. pp. 80–82. doi: 10.1007/978-1-4612-0819-8. ISBN 0-387-94415-X. MR 1329545.**^**Cormen, Thomas H.; Leiserson, Charles E.; Rivest, Ronald L.; Stein, Clifford (2009) [1990]. "Chapter 7: Quicksort".*Introduction to Algorithms*(3rd ed.). MIT Press and McGraw-Hill. pp. 170–190. ISBN 0-262-03384-4.**^**Cormen et al. (2009), Section 8.1, "Lower bounds for sorting", pp. 191–193.**^**Freniche, Francisco J. (2010). "On Riemann's rearrangement theorem for the alternating harmonic series" (PDF).*The American Mathematical Monthly*.**117**(5): 442–448. doi: 10.4169/000298910X485969. JSTOR 10.4169/000298910x485969. MR 2663251. S2CID 20575373.**^**Soddy, F. (1943). "The three infinite harmonic series and their sums (with topical reference to the Newton and Leibniz series for )".*Proceedings of the Royal Society*.**182**(989): 113–129. Bibcode: 1943RSPSA.182..113S. doi: 10.1098/rspa.1943.0026. MR 0009207. S2CID 202575422.**^**Bombieri, E. (2010). "The classical theory of zeta and -functions".*Milan Journal of Mathematics*.**78**(1): 11–59. doi: 10.1007/s00032-010-0121-8. MR 2684771. S2CID 120058240.**^**Schmuland, Byron (May 2003). "Random harmonic series" (PDF).*The American Mathematical Monthly*.**110**(5): 407–416. doi: 10.2307/3647827. JSTOR 3647827.**^**Bettin, Sandro; Molteni, Giuseppe; Sanna, Carlo (2018). "Small values of signed harmonic sums".*Comptes Rendus Mathématique*.**356**(11–12): 1062–1074. arXiv: 1806.05402. doi: 10.1016/j.crma.2018.11.007. hdl: 2434/634047. MR 3907571. S2CID 119160796.**^**Baillie, Robert (May 1979). "Sums of reciprocals of integers missing a given digit".*The American Mathematical Monthly*.**86**(5): 372–374. doi: 10.1080/00029890.1979.11994810. JSTOR 2321096.**^**Schmelzer, Thomas; Baillie, Robert (June 2008). "Summing a curious, slowly convergent series".*The American Mathematical Monthly*.**115**(6): 545–540. doi: 10.1080/00029890.2008.11920559. JSTOR 27642532. S2CID 11461182.

Wikimedia Commons has media related to
Harmonic series (mathematics).