Birthday problem with 3+ people?

I would like to use Matlab to figure out the probability that at least 3 people have the same birthday in a sample size of 50 people, but am having trouble figuring out how to do the counting. A series of loops perhaps?

Answer 1

[Jon]

Here are some ideas that you might use

To generate random birthdays for groups of 50 people - look at randi, ignoring leap years birthdays are numbers between 1 and 365

For each random sample group of 50 people you can use histcounts with 365 bins to count how many people have each birthday

with result of histcounts, you can determine if the sample has three or more people with the same birthday by taking the max and seeing if it is greater than or equal 3

do this N times and then estimate the probability as the number of times that you have 3 or more people with same birthday divided by N