% ECE 341 Random Proecsses
% Lecture #4: Farkle Odds
% Enueration

disp('------------------');
disp('Farkle Odds using enumeration');

tic

Pair6 = 0;
Pair5 = 0;
Pair4 = 0;
Pair3 = 0;
Pair222 = 0;
Pair33 = 0;
Pair42 = 0;

for d1=1:6
  for d2=1:6
    for d3=1:6
      for d4=1:6
        for d5=1:6
          for d6 = 1:6
            Dice = [d1,d2,d3,d4,d5,d6];
            Dice = sort(Dice);
     
            N = zeros(1,6);
            for i=1:6
                N(i) = sum(Dice == i);
            end
            [N,b] = sort(N, 'descend');

            if (N(1) == 6) Pair6 = Pair6 + 1;
            elseif (N(1) == 5) Pair5 = Pair5 + 1;
            elseif ((N(1)==4)*(N(2)==2)) Pair42 = Pair42 + 1;
            elseif (N(1)==4) Pair4 = Pair4 + 1;
            elseif ((N(1)==3)*(N(2)==3)) Pair33 = Pair33 + 1;
            elseif (N(1)==3) Pair3 = Pair3 + 1;
            elseif ((N(1)==2)*(N(2)==2)*(N(3)==2)) Pair222 = Pair222+1;
            end
         end
        end
      end
    end
  end
end

disp('  6   5x    4xy   3xyz   222   42  33')
[Pair6,Pair5,Pair4,Pair3,Pair222,Pair42,Pair33]

toc