% 4-Card Stud Pair3 = 0; for i0 = 1:1e4 X = rand(1,52); [a,Deck] = sort(X); Hand = Deck(1:4); 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) Pair3 = Pair3 + 1; end end [Pair3]