aez-notes

Pepys wrote Newton to ask which of the tree evens is more likely: that a person get (a) at least 1 six when 6 dice are rolled, (b) at least 2 sixes when 12 dice are rolled, or (c) at least 3 sixes when 18 dice are rolled. What is the answer?

Straight computation shows that 6 dice is the best.

load(distrib);

f(n) := 1 - cdf_binomial(n - 1, n * 6, 1 / 6);
map(f, [1,2,3]), numer;

/* approx 0.665 for 6 dice */