% Lecture 0: Monte Carlo 5-Card Stud
% Probability of 3 of a kind
tic
Pair3 = 0;

for i0 = 1:1e5
    
X = rand(1,52);
[a,Deck] = sort(X);
Hand = Deck(1:5);
Value = mod(Hand,13) + 1;
Suit = floor(Hand/13) + 1;

N = zeros(1,13);
for n=1:13
    N(n) = sum(Value == n);
end

[N,a] = sort(N, 'descend');

if ((N(1) == 3)*(N(2) < 2)) Pair3 = Pair3 + 1; end

end

[Pair3]
toc