R Output for Discrete Probability Example 4

> dbinom(0:3, 4, 13/52)

0.31640625  0.421875  0.2109375  0.046875

> (39/52) * (38/51) * (37/50) * (36/49)

0.303817527010804

> choose(4, 1) * (13/52) * (39/51) * (38/50) * (37/49)

0.438847539015606

> choose(4, 2) * (13/52) * (12/51) * (39/50) * (38/49)

0.213493397358944

> choose(4, 3) * (13/52) * (12/51) * (11/50) * (39/49)

0.0412004801920768