courseshub.world · Lives in Mathematics and Physics

Abraham de Moivre

Abraham de Moivre (born Abraham Moivre)

A course in two strands: the life of a Huguenot refugee who spent sixty-six years in London teaching mathematics by the hour, and the work in which he found the bell curve, the formula linking powers and angles, and the first mathematics of life insurance.

Born 26 May 1667, Vitry-le-François Died 27 November 1754, London Normal approximation November 1733 Level advanced undergraduate
Probabilities of k heads in 24 tosses of a fair coin, with the curve \(\frac{2}{\sqrt{2\pi n}}e^{-2(k-n/2)^2/n}\) that de Moivre found in 1733.
Lesson 00

How to use this course

Abraham de Moivre is remembered by students for one formula, \((\cos x+i\sin x)^n=\cos nx+i\sin nx\), and by statisticians for something larger: he was the first to see the normal curve emerge from the binomial distribution. He also wrote the standard English textbook of probability for half a century and founded the mathematics of life annuities. He did all of this without a university post, as a French Protestant exile who earned his living as a private tutor.

Prerequisites for Part II
  • Single-variable calculus, including Taylor series of \(\log(1+x)\) and \(e^x\).
  • Complex numbers in Cartesian and polar form.
  • Elementary combinatorics: binomial coefficients and the binomial theorem.
  • For Lessons 06–10, discrete probability at the level of a first course.

Each technical lesson begins with a paragraph for the general reader, then gives precise statements in modern notation, derivations, and exercises with folded solutions. De Moivre wrote in Latin and in English prose with little symbolism; where his notation differs from ours, the lesson says so. The order follows his career: algebra and series (1695–1707), probability (1711–1738), trigonometry and analysis (1722–1730), the 1733 limit theorem, and annuities (1725 onward).

Reference

Timeline

1667Born 26 May at Vitry-le-François in Champagne, son of Daniel Moivre, a Protestant surgeon.
c. 1678Sent to the Protestant Academy at Sedan; four years of Greek under Jacques du Rondel.
1682–84Studies logic at Saumur after Sedan is suppressed; reads Huygens on games of chance.
1684Moves to Paris; private lessons in mathematics from Jacques Ozanam.
1685Revocation of the Edict of Nantes. Detained at the priory of Saint-Martin-des-Champs for religious instruction.
168728 August: admitted with his brother to the French Protestant Savoy Church in London.
c. 1692Friendship with Edmond Halley; soon afterwards with Isaac Newton.
1695Halley presents his first paper, on fluxions, to the Royal Society.
1697The multinomial theorem. Elected Fellow of the Royal Society on 30 November.
1704Animadversiones against George Cheyne’s book on fluents.
1705Becomes a naturalized English subject.
1707Paper on equations of degrees 3, 5, 7, …: the germ of de Moivre’s formula.
1711De mensura sortis in the Philosophical Transactions.
1712Member of the Royal Society committee on the Newton–Leibniz priority dispute.
1718The Doctrine of Chances, first edition, dedicated to Newton.
1722De sectione anguli: the angle-multiplication formula.
1725Annuities upon Lives, first edition.
1730Miscellanea analytica with its Supplementum: recurring series, the central binomial term, the series for \(\log n!\).
1733November: the Latin pamphlet Approximatio ad summam terminorum binomii, the normal approximation.
1735Elected to the Berlin Academy of Sciences.
1738The Doctrine of Chances, second edition, with the 1733 theorem in English.
1743–52Revised editions of Annuities upon Lives (1743, 1750, 1752).
175427 June: foreign associate of the Paris Academy of Sciences. Dies in London on 27 November, aged 87.
1756Posthumous third edition of The Doctrine of Chances, with the annuities treatise bound in.
Part I · Formation
Lesson 01

Champagne, Sedan, Saumur, 1667–1684

Abraham Moivre was born on 26 May 1667 in Vitry-le-François, a small town on the Marne in Champagne. His father Daniel was a surgeon of modest but steady means, and the family was Protestant. The particle “de” appears only later, after he reached England. His first schooling was with the Catholic Christian Brothers in Vitry, a sign of how mixed provincial life still was before the persecutions of the 1680s.

At about eleven he was sent to the Protestant Academy at Sedan, one of the colleges the Huguenots had built under the protection of the Edict of Nantes. He studied Greek there for four years with Jacques du Rondel. In 1682 the crown suppressed the academy, part of a steady campaign against Protestant institutions, and he moved to the academy at Saumur on the Loire to study logic.

A book on dice

Mathematics was not part of his formal course. At Saumur he read on his own the Éléments des mathématiques of Jean Prestet and, more importantly, the short treatise De ratiociniis in ludo aleae (1657) by Christiaan Huygens, the first printed book on probability. Huygens computed fair shares in interrupted games and the odds of dice throws by the idea of expectation. The young reader later turned that slim tract into the largest book on chance of the century.

Lesson 02

Paris and the Revocation, 1684–1688

In 1684 Moivre went to Paris. Accounts differ on his formal enrolment; MacTutor places him at the Collège d’Harcourt for physics. What is clear is that he received his first real training in mathematics from Jacques Ozanam, a popular teacher and author of textbooks and books of mathematical recreations. Ozanam made his living by private lessons, the same trade his pupil would practise for the rest of his life.

In October 1685 Louis XIV revoked the Edict of Nantes by the Edict of Fontainebleau. Protestant worship became illegal and Protestant children could be taken for Catholic instruction. Moivre was confined in the priory of Saint-Martin-des-Champs in Paris, which then served as a school for the reconversion of young Protestants. The sources disagree on how long he stayed and on what terms he left: the priory’s records place his departure in 1688, while the register of the French Savoy Church in London records that he and his brother were admitted on 28 August 1687.

Whatever the exact sequence, by his early twenties he was in London, one of perhaps two hundred thousand Huguenots who left France in these years. He never returned. The Paris Academy elected him a foreign associate sixty-seven years later, five months before his death.

DiscussionThe two dates for de Moivre’s departure from France come from different kinds of records (a school register and a church register). How would a historian weigh them, and what possibilities reconcile them?

Lesson 03

A tutor in London, 1688–1710

In London de Moivre lived by teaching mathematics to the sons of the gentry, walking from house to house and meeting pupils in coffee houses. A well-known story says that at the house of the Earl of Devonshire he first saw Isaac Newton’s Principia (1687), bought a copy, and tore it into sheets that he carried in his pocket and read between lessons.

By about 1692 he had come to know Edmond Halley, then assistant secretary of the Royal Society, and through him Newton. In March 1695 Halley presented de Moivre’s first paper, on Newton’s method of fluxions. Two years later came the multinomial theorem (Lesson 04), and on 30 November 1697 he was elected a Fellow of the Royal Society. He was thirty.

Newton’s circle

De Moivre became one of the small group around Newton. He corrected the Latin translation of the Opticks, published in 1706. In 1704 he wrote a sharp Latin critique of George Cheyne’s treatise on the inverse method of fluxions, defending Newton’s methods against a careless popularizer. Newton’s regard for him is attested by a story, possibly embellished, that he would send questioners about the Principia to de Moivre as the man who knew it better than he did.

Friendship with the powerful did not produce a post. He became a naturalized English subject in 1705, but he hoped in vain for a university chair. In 1710 Johann Bernoulli, at his request, wrote to Gottfried Wilhelm Leibniz asking for help toward a professorship. Nothing came of it. Foreign birth counted against him at Cambridge and Oxford, and he remained a tutor for life.

Part II · The mathematics
Lesson 04

The multinomial theorem

For the general reader. Newton had shown how to expand \((a+b)^n\). De Moivre showed how to expand a sum of any number of terms raised to a power, and even an infinite series raised to a power. The coefficients count arrangements, and the result is the first general tool of combinatorial algebra.

Goals
  • State and prove the multinomial theorem.
  • Compute the coefficients of a power of a power series.
  • See the counting interpretation used throughout de Moivre’s probability.
Theorem (multinomial theorem, de Moivre 1697)

For a positive integer \(n\) and commuting quantities \(x_1,\dots,x_m\),

\[ (x_1+x_2+\cdots+x_m)^n=\sum_{k_1+\cdots+k_m=n}\binom{n}{k_1,\dots,k_m}x_1^{k_1}\cdots x_m^{k_m},\qquad \binom{n}{k_1,\dots,k_m}=\frac{n!}{k_1!\,k_2!\cdots k_m!}, \]

the sum running over all \(m\)-tuples of non-negative integers with sum \(n\).

Proof

Expanding the product of \(n\) factors, each term arises by choosing one \(x_j\) from each factor. The monomial \(x_1^{k_1}\cdots x_m^{k_m}\) arises once for each word of length \(n\) containing \(k_1\) letters \(x_1\), …, \(k_m\) letters \(x_m\). There are \(n!\) orderings of \(n\) labelled letters, and permuting equal letters among themselves gives the same word, so the number of words is \(n!/(k_1!\cdots k_m!)\).

Powers of an infinite series

De Moivre’s 1697 paper was titled A method of raising an infinite multinomial to any given power, or extracting any given root of the same. In modern terms, if \(f(z)=a_0+a_1z+a_2z^2+\cdots\) then the coefficient of \(z^N\) in \(f(z)^n\) is

\[ [z^N]\,f(z)^n=\sum\binom{n}{k_0,k_1,k_2,\dots}\,a_0^{k_0}a_1^{k_1}a_2^{k_2}\cdots,\qquad \sum_j k_j=n,\ \ \sum_j j\,k_j=N. \]

So the coefficient is a sum over partitions of \(N\) into at most \(n\) parts. This is the bookkeeping behind every later computation of the chance of throwing a given total with several dice.

Exercise 4.1Find the coefficient of \(x^2y^3z\) in \((x+y+z)^6\) and of \(x^2y^2z^2\) in the same expansion.

Show solution
\(\frac{6!}{2!\,3!\,1!}=60\) and \(\frac{6!}{2!\,2!\,2!}=90\).

Exercise 4.2Three fair dice are thrown. Use the coefficient of \(z^{10}\) in \((z+z^2+\cdots+z^6)^3\) to find the probability that the total is 10.

Show solution
Write \((z+\cdots+z^6)^3=z^3(1-z^6)^3(1-z)^{-3}\). We need the coefficient of \(z^7\) in \((1-3z^6+\cdots)\sum_j\binom{j+2}{2}z^j\), which is \(\binom{9}{2}-3\binom{3}{2}=36-9=27\). The probability is \(27/216=1/8\).

Exercise 4.3Find the first three terms of \((1+z+z^2+z^3+\cdots)^{3}\) by de Moivre’s rule, and check against \((1-z)^{-3}\).

Show solution
All \(a_j=1\). For \(N=1\): only \(k_0=2,k_1=1\), coefficient \(3\). For \(N=2\): \((k_0,k_1,k_2)=(2,0,1)\) gives 3 and \((1,2,0)\) gives 3, total 6. So \(1+3z+6z^2+\cdots\), matching \(\binom{N+2}{2}\).
Lesson 05

Powers and angles: de Moivre’s formula

For the general reader. Multiplying an angle by \(n\) and raising a quantity to the \(n\)-th power turn out to be the same operation in disguise. De Moivre used this to solve certain equations of high degree with trigonometric tables. The modern statement, with \(i=\sqrt{-1}\), packs the idea into one line.

Goals
  • State and prove de Moivre’s formula and derive multiple-angle identities.
  • Understand his 1707 solution of the equations \(T_n(x)=y\).
  • Prove the Cotes–de Moivre factorization of \(x^{2n}-2x^n\cos n\theta+1\).

The 1707 equations

In 1707 de Moivre considered the equations of odd degree that express \(\cos n\phi\) as a polynomial in \(\cos\phi\). In modern notation these are \(T_n(x)=y\), where \(T_n\) is the Chebyshev polynomial defined by \(T_n(\cos\phi)=\cos n\phi\); for instance \(T_3(x)=4x^3-3x\). He gave the root in the form

\[ x=\tfrac12\Big(\sqrt[n]{\,y+\sqrt{y^2-1}\,}+\sqrt[n]{\,y-\sqrt{y^2-1}\,}\Big), \]

a generalization of Cardano’s formula for the cubic. When \(|y|\le1\) the quantities under the \(n\)-th roots are complex, and the root is real only because of the identity below. In De sectione anguli (1722) he made the link with angle multiplication explicit, though always in terms of cosines and versed sines, never with \(i\) written as we do.

Theorem (de Moivre’s formula)

For every real \(x\) and every integer \(n\),

\[ (\cos x+i\sin x)^n=\cos nx+i\sin nx. \]

Proof. For \(n\ge0\) use induction and the addition formulas: \((\cos x+i\sin x)(\cos nx+i\sin nx)=\cos(n+1)x+i\sin(n+1)x\). For negative \(n\) note that \((\cos x+i\sin x)^{-1}=\cos x-i\sin x=\cos(-x)+i\sin(-x)\). Leonhard Euler stated the formula in this form and extended it to real exponents, and with \(e^{ix}=\cos x+i\sin x\) it reduces to \((e^{ix})^n=e^{inx}\).

Roots and factorization

The \(n\) solutions of \(z^n=\cos\theta+i\sin\theta\) are \(z_k=\cos\frac{\theta+2\pi k}{n}+i\sin\frac{\theta+2\pi k}{n}\), \(k=0,\dots,n-1\). Roger Cotes had found the factorization of \(x^{2n}\pm1\) into real quadratic factors. De Moivre, in Miscellanea analytica (1730), gave the general version:

Theorem (Cotes–de Moivre factorization) \[ x^{2n}-2x^n\cos n\theta+1=\prod_{k=0}^{n-1}\Big(x^2-2x\cos\big(\theta+\tfrac{2\pi k}{n}\big)+1\Big). \]

Proof. The left side equals \((x^n-e^{in\theta})(x^n-e^{-in\theta})\). The roots of \(x^n=e^{in\theta}\) are \(e^{i(\theta+2\pi k/n)}\), those of \(x^n=e^{-in\theta}\) their conjugates. Pairing each root with its conjugate gives \((x-e^{i\alpha})(x-e^{-i\alpha})=x^2-2x\cos\alpha+1\).

Exercise 5.1Use de Moivre’s formula with \(n=3\) to express \(\cos3x\) and \(\sin3x\) in terms of \(\cos x\) and \(\sin x\).

Show solution
Expand \((c+is)^3=c^3+3ic^2s-3cs^2-is^3\). Real part: \(\cos3x=c^3-3cs^2=4c^3-3c\). Imaginary part: \(\sin3x=3c^2s-s^3=3s-4s^3\).

Exercise 5.2Solve \(4x^3-3x=\tfrac12\) using the trigonometric method, and check that de Moivre’s 1707 expression gives one of the roots.

Show solution
Set \(x=\cos\phi\): \(\cos3\phi=\tfrac12\), so \(3\phi=\pm\frac{\pi}{3}+2\pi k\), giving \(x=\cos20^\circ,\ \cos100^\circ,\ \cos140^\circ\). In the formula, \(y\pm\sqrt{y^2-1}=\tfrac12\pm i\tfrac{\sqrt3}{2}=e^{\pm i\pi/3}\); taking principal cube roots \(e^{\pm i\pi/9}\) gives \(x=\tfrac12(e^{i\pi/9}+e^{-i\pi/9})=\cos20^\circ\).

Exercise 5.3Taking \(\theta=0\) and \(x=1\) in the factorization, prove \(\prod_{k=1}^{n-1}2\sin\frac{\pi k}{n}=n\).

Show solution
With \(\theta=0\) the left side is \((x^n-1)^2\). Divide both sides by \((x-1)^2\), the \(k=0\) factor, and let \(x\to1\): the left becomes \(n^2\). Each remaining factor at \(x=1\) is \(2-2\cos\frac{2\pi k}{n}=4\sin^2\frac{\pi k}{n}\). So \(\prod_{k=1}^{n-1}4\sin^2\frac{\pi k}{n}=n^2\), and all the sines are positive, so take square roots.
Lesson 06

The Doctrine of Chances

For the general reader. Gamblers wanted to know what a bet was worth. De Moivre turned their questions into a systematic theory, with rules for combining probabilities, a clear definition of independent events, and general methods in place of clever tricks. His book was the probability textbook of the English-speaking world for most of the eighteenth century.

Goals
  • Know the publication history and the main problems of the book.
  • Solve the gambler’s ruin problem by de Moivre’s method.
  • Derive his rule of thumb \(n\approx0.693/p\) for small probabilities.

The book

De Moivre’s Latin memoir De mensura sortis (“On the measurement of chance”) appeared in the Philosophical Transactions for 1711. The English book The Doctrine of Chances: or, a Method of Calculating the Probabilities of Events in Play followed in 1718, dedicated to Newton, with enlarged editions in 1738 and 1756. Its introduction states the multiplication rule for independent events in words: the probability that several independent events all happen is the product of their probabilities. It treats the problem of points, the dice and card games of the day, the matching problem (jeu de rencontre), and the duration of play.

The gambler’s ruin

Two players A and B play repeated independent rounds. A wins a round with probability \(p\), B with \(q=1-p\), and the loser pays the winner one counter. A starts with \(a\) counters, B with \(b\).

Proposition (de Moivre, 1711)

If \(p\ne q\), the probability that A ends with all \(a+b\) counters is

\[ P_A=\frac{1-(q/p)^{a}}{1-(q/p)^{a+b}}. \]

If \(p=q=\tfrac12\), then \(P_A=\dfrac{a}{a+b}\).

De Moivre’s argument. Give A’s counters the values \(q/p,\ (q/p)^2,\dots,(q/p)^a\) and B’s the values \((q/p)^{a+1},\dots,(q/p)^{a+b}\), stacked so that the counter that changes hands is always the one at the boundary. A round then becomes a fair game in value: A wins the counter worth \(r^{x+1}\) with probability \(p\) or loses the one worth \(r^{x}\) with probability \(q\), where \(r=q/p\) and \(x\) is A’s current number of counters, and \(p\,r^{x+1}=q\,r^{x}\). Since every round is fair, A’s expected final wealth equals his initial wealth. Equating \(P_A\sum_{j=1}^{a+b}r^j\) with \(\sum_{j=1}^{a}r^j\) and summing the geometric series gives the formula. In modern terms he had found the martingale \(r^{X_t}\).

How many trials for an even chance?

An event has probability \(p\) per trial. The chance it occurs at least once in \(n\) trials is \(1-(1-p)^n\). Setting this equal to \(\tfrac12\) gives \(n=\log2/(-\log(1-p))\). In De mensura sortis de Moivre observed that for small \(p\) this is close to \(n\approx(\log 2)/p\approx0.7/p\). The same computation, read the other way, gives the Poisson approximation \((1-p)^n\approx e^{-np}\); for this reason some historians argue that the Poisson distribution first appears in his work.

Exercise 6.1The Chevalier de Méré’s problem: how many throws of two dice are needed for an even chance of at least one double six? Compare the exact answer with de Moivre’s rule.

Show solution
\(p=1/36\). Exactly, \(1-(35/36)^{24}\approx0.491\) and \(1-(35/36)^{25}\approx0.506\), so 25 throws. De Moivre’s rule gives \(0.693\times36\approx24.95\), also 25.

Exercise 6.2A starts with 3 counters and B with 7, and A wins each round with probability \(0.6\). Find the probability that A ruins B. What is it for a fair game?

Show solution
\(r=0.4/0.6=2/3\). \(P_A=\frac{1-(2/3)^3}{1-(2/3)^{10}}=\frac{0.7037}{0.9827}\approx0.716\). For a fair game \(P_A=3/10\). A modest edge more than overcomes a large disadvantage in capital.

Exercise 6.3The matching problem, treated in the 1718 and 1738 editions: \(n\) cards numbered \(1,\dots,n\) are shuffled. Show by inclusion–exclusion that the probability that no card is in its own place is \(\sum_{k=0}^{n}(-1)^k/k!\), and evaluate it for \(n=4\).

Show solution
Let \(A_i\) be the event that card \(i\) is in place. Then \(P(A_{i_1}\cap\cdots\cap A_{i_k})=(n-k)!/n!\), and there are \(\binom nk\) such intersections, so \(P(\bigcup A_i)=\sum_{k\ge1}(-1)^{k+1}/k!\). The complement is \(\sum_{k=0}^{n}(-1)^k/k!\). For \(n=4\): \(1-1+\tfrac12-\tfrac16+\tfrac1{24}=\tfrac{9}{24}=0.375\), already close to \(1/e\approx0.368\).
Lesson 07

Recurring series and generating functions

For the general reader. Many probability problems lead to sequences where each term is built from the few before it. De Moivre called these recurring series. He showed that the whole infinite sequence can be packed into a single fraction, and that unpacking the fraction gives a closed formula for every term. This is the method of generating functions.

Goals
  • Show that linear recurrences with constant coefficients have rational generating functions.
  • Use partial fractions to find closed forms.
  • Derive the closed form for the Fibonacci numbers.

De Moivre called the coefficients of the recurrence the scale of relation. Suppose a sequence satisfies

\[ u_{n}=c_1u_{n-1}+c_2u_{n-2}+\cdots+c_ku_{n-k}\qquad(n\ge k). \]
Theorem (de Moivre)

The series \(U(z)=\sum_{n\ge0}u_nz^n\) of a sequence with scale of relation \(c_1,\dots,c_k\) is a rational function

\[ U(z)=\frac{P(z)}{1-c_1z-c_2z^2-\cdots-c_kz^k}, \]

with \(\deg P<k\) determined by \(u_0,\dots,u_{k-1}\). Conversely, every such rational function expands as a recurring series with that scale.

Proof. Multiply \(U(z)\) by \(Q(z)=1-c_1z-\cdots-c_kz^k\). For \(n\ge k\) the coefficient of \(z^n\) in \(Q(z)U(z)\) is \(u_n-c_1u_{n-1}-\cdots-c_ku_{n-k}=0\), so \(QU\) is a polynomial of degree less than \(k\). Conversely, \(P/Q\) has a power series whose coefficients satisfy the same relation for the same reason.

If \(Q(z)=\prod_j(1-\alpha_jz)\) with distinct \(\alpha_j\), partial fractions give \(U(z)=\sum_j\frac{A_j}{1-\alpha_jz}\), hence \(u_n=\sum_jA_j\alpha_j^n\). De Moivre used this to solve the duration-of-play problem, in which the probability that a game ends at a given round obeys a linear recurrence.

Exercise 7.1The Fibonacci numbers satisfy \(F_0=0,\ F_1=1,\ F_n=F_{n-1}+F_{n-2}\). Find their generating function and derive \(F_n=\frac{\varphi^n-\psi^n}{\sqrt5}\), with \(\varphi=\frac{1+\sqrt5}2,\ \psi=\frac{1-\sqrt5}2\). (This closed form was known to de Moivre more than a century before Jacques Binet published it in 1843.)

Show solution
\(Q=1-z-z^2\) and \(P=F_0+(F_1-F_0)z=z\), so \(F(z)=\frac{z}{1-z-z^2}=\frac{z}{(1-\varphi z)(1-\psi z)}\). Partial fractions: \(\frac{z}{(1-\varphi z)(1-\psi z)}=\frac{1}{\varphi-\psi}\big(\frac{1}{1-\varphi z}-\frac{1}{1-\psi z}\big)\), and \(\varphi-\psi=\sqrt5\).

Exercise 7.2A sequence satisfies \(u_n=5u_{n-1}-6u_{n-2}\) with \(u_0=1,\ u_1=4\). Find its generating function and \(u_n\).

Show solution
\(Q=1-5z+6z^2=(1-2z)(1-3z)\), \(P=u_0+(u_1-5u_0)z=1-z\). Then \(\frac{1-z}{(1-2z)(1-3z)}=\frac{A}{1-2z}+\frac{B}{1-3z}\) with \(A+B=1,\ 3A+2B=1\), so \(A=-1,\ B=2\): \(u_n=2\cdot3^n-2^n\). Check: \(u_1=6-2=4\).
Lesson 08

Factorials and the central term

For the general reader. To find the chance of exactly 5000 heads in 10,000 tosses, one must handle numbers like \(10000!\), far beyond hand calculation. De Moivre found an accurate formula for such factorials. His friend James Stirling supplied the missing constant, \(\sqrt{2\pi}\), and the formula has carried Stirling’s name ever since.

Goals
  • Derive \(\log n!\approx(n+\tfrac12)\log n-n+C\) by comparing sums and integrals.
  • Identify \(C=\log\sqrt{2\pi}\) using Wallis’s product.
  • Approximate the central binomial term.

The history

In Miscellanea analytica (1730) de Moivre studied the ratio of the middle term of \((1+1)^n\) to the sum \(2^n\). In the Supplementum he gave a series for \(\log n!\) of the form below, with the constant expressed as a series \(1-\frac1{12}+\frac1{360}-\frac1{1260}+\cdots\) which he evaluated numerically. James Stirling, who corresponded with him in 1729, showed in his Methodus differentialis (1730) that the constant is \(\log\sqrt{2\pi}\). De Moivre acknowledged this openly. Thomas Bayes later showed that de Moivre’s series for the constant does not converge, although its first terms give an excellent approximation.

Theorem (de Moivre–Stirling formula) \[ n!\sim\sqrt{2\pi n}\,\Big(\frac ne\Big)^{n},\qquad \log n!=\Big(n+\tfrac12\Big)\log n-n+\log\sqrt{2\pi}+\frac1{12n}-\frac1{360n^3}+O(n^{-5}). \]

Derivation

Let \(d_n=\log n!-(n+\tfrac12)\log n+n\). Then

\[ d_n-d_{n+1}=\Big(n+\tfrac12\Big)\log\frac{n+1}{n}-1=\frac{1}{2t}\log\frac{1+t}{1-t}-1=\frac{t^2}{3}+\frac{t^4}{5}+\cdots,\qquad t=\frac{1}{2n+1}. \]

So \(d_n\) decreases, \(0<d_n-d_{n+1}<\frac{t^2}{3(1-t^2)}=\frac{1}{12}\big(\frac1n-\frac1{n+1}\big)\), and \(d_n\) converges to a limit \(C\) with \(C<d_n<C+\frac1{12n}\). To identify \(C\), take Wallis’s product in the form \(\lim_n\frac{2^{4n}(n!)^4}{((2n)!)^2(2n+1)}=\frac{\pi}{2}\). Substituting \(n!\approx e^{C}n^{n+1/2}e^{-n}\) gives \(e^{2C}/4=\pi/2\), so \(e^{C}=\sqrt{2\pi}\).

Corollary (the central term)

For even \(n\), \(\displaystyle\binom{n}{n/2}2^{-n}\sim\frac{2}{\sqrt{2\pi n}}\).

Exercise 8.1Compare \(10!\) with \(\sqrt{20\pi}(10/e)^{10}\), and then include the factor \(e^{1/120}\).

Show solution
\(10!=3\,628\,800\). The leading term gives about \(3\,598\,696\), a relative error of 0.83%, close to \(1/(12\cdot10)\). Multiplying by \(e^{1/120}\) gives about \(3\,628\,810\), an error of 3 parts per million.

Exercise 8.2Use the corollary to estimate the probability of exactly 50 heads in 100 tosses of a fair coin. The exact value is \(0.07959\).

Show solution
\(2/\sqrt{200\pi}=2/25.066\approx0.07979\), an error of about 0.25%.

Exercise 8.3Prove the corollary from the theorem.

Show solution
With \(n=2m\), \(\binom{2m}{m}=\frac{(2m)!}{(m!)^2}\sim\frac{\sqrt{4\pi m}(2m/e)^{2m}}{2\pi m\,(m/e)^{2m}}=\frac{2^{2m}}{\sqrt{\pi m}}\). Dividing by \(2^{2m}\) gives \(1/\sqrt{\pi m}=2/\sqrt{2\pi n}\).
Lesson 09

The normal approximation, 1733

For the general reader. Toss a coin many times and draw a bar chart of how often each number of heads comes up. The bars trace a bell-shaped curve. De Moivre found the equation of that curve in 1733 and used it to show how tightly the proportion of heads clusters around one half. It was the first appearance of the normal distribution and the first case of the central limit theorem.

Goals
  • Derive the local approximation of binomial terms by \(e^{-2l^2/n}\).
  • State the de Moivre–Laplace theorem for general \(p\).
  • Use it to quantify Jacob Bernoulli’s law of large numbers.

In November 1733 de Moivre circulated privately a short Latin pamphlet, Approximatio ad summam terminorum binomii \((a+b)^n\) in seriem expansi. He wanted to improve Jacob Bernoulli’s law of large numbers (1713), which proved convergence of frequencies but gave hopeless bounds on the number of trials needed. An English version appeared in the 1738 and 1756 editions of The Doctrine of Chances.

The local step

Let \(n\) be even and \(M=\binom{n}{n/2}\). For the term at distance \(l\) from the middle,

\[ \log\frac{\binom{n}{n/2+l}}{M}=\sum_{j=1}^{l}\log\frac{n/2-j+1}{n/2+j}\approx\sum_{j=1}^{l}\Big(-\frac{4j-2}{n}\Big)=-\frac{2l^2}{n}, \]

using \(\log\frac{1-x}{1+y}\approx-x-y\) for small \(x,y\). Combined with Lesson 08, the term is approximately \(\frac{2}{\sqrt{2\pi n}}e^{-2l^2/n}\). This is the curve drawn at the top of the page.

Summing

Summing over \(|l|\le L\) and replacing the sum by an integral, de Moivre obtained the probability that the number of heads lies within \(L\) of \(n/2\):

\[ P\big(|S_n-\tfrac n2|\le L\big)\approx\frac{4}{\sqrt{2\pi n}}\int_0^{L}e^{-2t^2/n}\,dt. \]

He evaluated it by expanding the exponential in series, for \(L=\tfrac12\sqrt n\), \(\sqrt n\) and \(\tfrac32\sqrt n\). In modern terms these are 1, 2 and 3 standard deviations, and the answers are about 0.68, 0.95 and 0.997. He called \(\tfrac12\sqrt n\) the modulus, the unit by which deviations should be measured. He also indicated the extension to unequal chances \(p\ne q\).

Theorem (de Moivre 1733; Laplace 1812)

Let \(S_n\) be the number of successes in \(n\) independent trials with success probability \(p\in(0,1)\), \(q=1-p\). For all real \(a<b\),

\[ \lim_{n\to\infty}P\Big(a\le\frac{S_n-np}{\sqrt{npq}}\le b\Big)=\frac{1}{\sqrt{2\pi}}\int_a^b e^{-x^2/2}\,dx. \]

De Moivre proved the symmetric case \(p=\tfrac12\) and sketched the general one. Pierre-Simon Laplace gave the general proof in his Théorie analytique des probabilités (1812), hence the name de Moivre–Laplace theorem. The curve was later attached to Carl Friedrich Gauss’s theory of errors and became the “Gaussian”, but de Moivre found it first.

De Moivre drew a philosophical conclusion: the regularity of frequencies in the long run reflects a stable underlying order, which he attributed to design. The square-root law he discovered, that fluctuations grow like \(\sqrt n\) while totals grow like \(n\), is the quantitative content.

Exercise 9.1A fair coin is tossed 10,000 times. Estimate the probability that the number of heads lies between 4950 and 5050 inclusive. The exact value is about 0.6875.

Show solution
\(\sqrt{npq}=50\), so the range is \(\pm1\) standard deviation: \(\approx0.6827\). With the continuity correction (use \(\pm50.5\)): \(\Phi(1.01)-\Phi(-1.01)\approx0.6875\).

Exercise 9.2How many tosses of a fair coin make the probability at least 0.95 that the proportion of heads lies within 0.01 of \(\tfrac12\)?

Show solution
We need \(0.01n\ge1.96\sqrt{n}/2\), so \(\sqrt n\ge98\), \(n\ge9604\). Bernoulli’s own 1713 bounds required far larger numbers of trials for comparable guarantees; this is the improvement de Moivre sought.

Exercise 9.3Using the local step, estimate \(\binom{100}{55}\big/\binom{100}{50}\) and compare with the exact value \(0.6091\).

Show solution
\(e^{-2\cdot25/100}=e^{-1/2}\approx0.6065\), within 0.5%.
Lesson 10

Annuities upon lives

For the general reader. Governments and towns in the seventeenth and eighteenth centuries raised money by selling life annuities: pay a lump sum now and receive a yearly income until death. Pricing them requires combining interest with the chances of survival. Halley had shown how to do it with a mortality table, but the arithmetic was enormous. De Moivre found a simple model of mortality that made the calculation short.

Goals
  • Define the present value of a life annuity from a life table.
  • Derive de Moivre’s closed formula under linear mortality.
  • Understand the joint-life problem.

Halley’s method

In 1693 Halley published a life table built from the burial and birth records of Breslau, with a method for valuing annuities. If \(l_x\) is the number of survivors at age \(x\), \(i\) the rate of interest and \(v=1/(1+i)\), then an annuity of 1 per year, paid at the end of each year while a life now aged \(x\) survives, is worth

\[ a_x=\sum_{t\ge1}v^t\,\frac{l_{x+t}}{l_x}. \]

De Moivre’s hypothesis

In Annuities upon Lives (1725) de Moivre observed that over most of adult life Halley’s numbers of survivors decrease roughly in arithmetic progression. He assumed that \(l_x\) is proportional to \(86-x\) from about age 12 onward, with 86 the limiting age. The remaining span \(n=86-x\) he called the complement of life.

Proposition (de Moivre, 1725)

If \(l_{x+t}/l_x=1-t/n\) for \(0\le t\le n\), then

\[ a_x=\frac{1}{i}\Big(1-\frac{1+i}{n}\,a_{\overline{n}|}\Big),\qquad a_{\overline{n}|}=\sum_{t=1}^{n}v^t=\frac{1-v^n}{i}. \]

Derivation. We need \(\sum_{t=1}^{n}v^t-\frac1n\sum_{t=1}^{n}t\,v^t\). With \((Ia)_{\overline{n}|}=\sum tv^t=\frac{(1+i)a_{\overline{n}|}-nv^n}{i}\) (sum by parts), the result is \(a_{\overline{n}|}-\frac{(1+i)a_{\overline{n}|}}{ni}+\frac{v^n}{i}=\frac{1-v^n}{i}+\frac{v^n}{i}-\frac{(1+i)a_{\overline{n}|}}{ni}\), which is the formula. One table of annuities certain then values every life annuity.

The same book treats annuities on two or more lives, reversions and tontines. De Moivre approximated the curve of survivors by straight segments where a single line fitted poorly, a piecewise-linear method that kept the arithmetic short. His later editions (1743, 1750, 1752) and the 1756 printing with The Doctrine of Chances made the work a standard reference for the new insurance societies.

Exercise 10.1Value an annuity of £1 per year on a life aged 56 at 5% interest under de Moivre’s hypothesis.

Show solution
\(n=30\), \(a_{\overline{30}|}=(1-1.05^{-30})/0.05\approx15.372\). Then \(a_{56}=\frac{1}{0.05}\big(1-\frac{1.05}{30}\times15.372\big)\approx20(1-0.5380)\approx9.24\). A direct sum \(\sum_{t=1}^{30}1.05^{-t}(1-t/30)\) gives the same.

Exercise 10.2Under de Moivre’s hypothesis, show that the expected remaining lifetime at age \(x\) is \(n/2\), and that the probability two independent lives with complements \(n\) and \(m\ge n\) both survive \(t\) years is \((1-t/n)(1-t/m)\).

Show solution
The time of death is uniform on \([0,n]\), so its mean is \(n/2\). For two independent lives the survival probabilities multiply: \(\frac{n-t}{n}\cdot\frac{m-t}{m}\) for \(t\le n\). The value of a joint-life annuity is \(\sum_t v^t(1-t/n)(1-t/m)\), which also has a closed form in annuities certain.
Part III · Life and legacy
Lesson 11

Priority, patrons and poverty, 1710–1740

In 1712 the Royal Society appointed de Moivre to the committee that examined the dispute between Newton and Leibniz over the invention of the calculus. The committee, drawn from Newton’s supporters, found for Newton. De Moivre’s presence there reflected his closeness to Newton rather than impartiality, and it sits awkwardly with the fact that Leibniz and Johann Bernoulli had tried to help him two years earlier.

His own priority quarrels concerned probability. Pierre Rémond de Montmort had published his Essay d’analyse sur les jeux de hazard in 1708, with a second edition in 1713 that included correspondence with Nicolaus Bernoulli. Montmort felt that De mensura sortis took insufficient account of his work, and de Moivre felt the reverse. The exchange was sharp in print but ended on civil terms. Later, de Moivre accused Thomas Simpson of taking his annuity methods in cheaper books, and complained of it in the preface to his 1743 edition.

The relationship with James Stirling was the opposite. Stirling wrote to him in June 1729 about the central binomial term, and de Moivre printed both his own result and Stirling’s improvement, crediting his worthy and learned friend. He was elected to the Berlin Academy in 1735.

None of this paid. He supported himself mainly by tutoring, and he was known to sit at Slaughter’s Coffee House in St Martin’s Lane, where he played chess for small stakes and answered questions. He never married. The second edition of The Doctrine of Chances (1738) and the revisions of Annuities were his main sources of reputation, if not of income.

DiscussionDe Moivre sat on a committee judging a dispute in which he had personal loyalties, and was himself party to priority disputes. Compare the standards of credit in 1712 with those of a modern journal. What has changed, and what has not?

Lesson 12

Last years and legacy, 1740–1754

De Moivre worked on new editions into his eighties. On 27 June 1754 the Paris Academy of Sciences elected him a foreign associate, the recognition his native country had withheld. He died in London on 27 November 1754, aged 87, and was buried at St Martin-in-the-Fields. A famous story says that he noticed he was sleeping a quarter of an hour longer each night, summed the progression, and predicted that he would die on the day he slept for 24 hours. The story circulated after his death and has no contemporary documentation. The third edition of The Doctrine of Chances appeared in 1756.

Legacy

  • Probability and statistics: the normal approximation, the square-root law and the idea of measuring deviations in a natural unit, which Laplace generalized and Gauss tied to the theory of errors. The central limit theorem in its modern forms descends from the 1733 pamphlet.
  • Method: generating functions and recurring series became the standard analytic tools of combinatorics and of the probability of Laplace.
  • Analysis: the de Moivre–Stirling formula for \(n!\), and the angle-multiplication formula through which every student first meets complex numbers in polar form.
  • Actuarial science: “de Moivre’s law” of uniform mortality remains the first example in actuarial textbooks, and his annuity methods served the life insurance offices founded in Britain in the eighteenth century.
  • Influence: Thomas Bayes read him closely; Laplace built on him; Euler gave his trigonometric results their modern form.

Final projectReproduce de Moivre’s 1733 computation. Starting from the local approximation, compute by series expansion the probability that 3600 tosses of a fair coin give between 1770 and 1830 heads. Compare with a direct numerical sum of binomial terms, and explain the source of the remaining error.

Reference

Connections

People linked to Abraham de Moivre across this collection. A name without a link has no page here yet. See Abraham de Moivre in the network of all 337 lives →

Teachers and advisors

  • Jacques Ozanam 1684-1685

Influenced by

Influenced

Colleagues and collaborators

Rivals

  • Pierre Rémond de Montmort 1711-1719
  • Thomas Simpson

Friends

Also mentioned in this life

Fields

Reference

Sources and further reading

Biography

  1. J. J. O’Connor and E. F. Robertson, Abraham de Moivre, MacTutor History of Mathematics.
  2. Abraham de Moivre, Wikipedia.
  3. R. J. Pulskamp, Abraham de Moivre: works and translations.
  4. I. Schneider, Abraham de Moivre, Encyclopedia of Mathematics.
  5. D. R. Bellhouse, Abraham de Moivre: Setting the Stage for Classical Probability and Its Applications, CRC Press, 2011.
  6. H. M. Walker, Abraham de Moivre, in the 1967 reprint of The Doctrine of Chances, Chelsea.

Mathematics

  1. A. de Moivre, A method of raising an infinite multinomial to any given power, Phil. Trans. 19 (1697).
  2. A. de Moivre, De mensura sortis, Phil. Trans. 27 (1711); English translation by B. McClintock, International Statistical Review 52 (1984).
  3. A. de Moivre, The Doctrine of Chances, London, 1718; 2nd ed. 1738; 3rd ed. 1756.
  4. A. de Moivre, Annuities upon Lives, London, 1725.
  5. A. de Moivre, Miscellanea analytica de seriebus et quadraturis, London, 1730.
  6. A. Hald, A History of Probability and Statistics and Their Applications before 1750, Wiley, 1990.
  7. S. M. Stigler, The History of Statistics, Harvard University Press, 1986, chapter 2.
  8. P.-S. Laplace, Théorie analytique des probabilités, Paris, 1812.