Consider dealing a five-card hand from a deck of playing cards. A deck of cards contains 52 cards, 13 of each suit. Suppose that we wish to determine the probability, when dealing a five-card hand, that the hand contains a 'flush' (i.e., all five cards belong to the same suit).

a) Describe how you could estimate this probability using Monte Carlo simulation. Be clear and specific about how your simulation would operate, how the probability would be estimated, and what underlying principles would make this a good estimate.

b) Is a monte carlo simulation, or a direct analytical solution, a better choice here to determine the probability?