Newton-Pepys Problem

DOWNLOAD Mathematica Notebook NewtonPepysProblem

Samuel Pepys wrote Isaac Newton a long letter asking him to determine the probabilities for a set of dice rolls related to a wager he planned to make. Pepys asked which was more likely,

1. At least one six when six dice are rolled,

2. At least two sixes when 12 dice are rolled, or

3. At least three sixes when 18 dice are rolled.

Newton-PepysProblem

The probability of obtaining n or more sixes in a roll of 6n dice is given by

P_n=sum_(x=n)^(6n)(6n; x)(1/6)^x(5/6)^(6n-x)
(1)
=1-sum_(x=0)^(n-1)(6n; x)(1/6)^x(5/6)^(6n-x)
(2)
=(5^(5n))/(6^(6n))(6n; n)_2F^~_1(1,-5n;1+n;-1/5),
(3)

where _2F^~_1(a,b;c;x) is a regularized hypergeometric function. Values for n=1, 2, and 3 gives

P_1=(31031)/(46656) approx 0.6651
(4)
P_2=(1346704211)/(2176782336) approx 0.6187
(5)
P_3=(15166600495229)/(25389989167104) approx 0.5973
(6)

(OEIS A143162 and A143163). Therefore, P_1 is the most likely, with P_n asymptotically approaching 1/2 as n->infty.

Wolfram Web Resources

Mathematica »

The #1 tool for creating Demonstrations and anything technical.

Wolfram|Alpha »

Explore anything with the first computational knowledge engine.

Wolfram Demonstrations Project »

Explore thousands of free applications across science, mathematics, engineering, technology, business, art, finance, social sciences, and more.

Computerbasedmath.org »

Join the initiative for modernizing math education.

Online Integral Calculator »

Solve integrals with Wolfram|Alpha.

Step-by-step Solutions »

Walk through homework problems step-by-step from beginning to end. Hints help you try the next step on your own.

Wolfram Problem Generator »

Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet.

Wolfram Education Portal »

Collection of teaching and learning tools built by Wolfram education experts: dynamic textbook, lesson plans, widgets, interactive Demonstrations, and more.

Wolfram Language »

Knowledge-based programming for everyone.