%% Uvodni hlavicka
 clear     %vymaze vsechny promenne z pameti
 close all %zavre vsechna okna
 clc

%% Definovani promennych
 S = 100000; %pocet opakovani hry (pocet simulaci)
 r = 0;      %pocitadlo stavu 5:3
 E_Bwin = 0; %stredni hodnota vyhry Boba pri stavu 5:3
 
%% Hra
 for s = 1:S
  p = rand;   
   %zahrajem 8 her  
     y = rand(8,1);
   %overime stav 5:3 => alice vyhraje prave 5x
    if (sum(y<p)==5)
       r = r + 1; %pocitadlo zaznamena stav 5:3
       E_Bwin = E_Bwin+(1-p)^3; %scitame pro kazdy stav 5:3 P(bob win)
    end
 end
 E_Bwin = E_Bwin/r; %stredni hodnota pravdepodobnosti vyhry Boba
 E_Awin = 1-E_Bwin; %stredni hodnota pravdepodobnosti vyhry Alice
 
 %% Textovy vypis na obrazovku
 disp('Prumerna pravdepodobnost vyhry Boba pri stavu 5:3')
 disp(E_Bwin)
 disp('Prumerna pravdepodobnost vyhry Alice pri stavu 5:3')
 disp(E_Awin)
 disp('Ferovy podil sanci')
 disp(E_Awin/E_Bwin)
 disp('Pocet stavu 5:3')
 disp(r)
 disp('Relativni pocet stavu 5:3 (z poctu simulaci)')
 disp(r/S)
 
%% graficke zobrazeni
 figure
 plot(1,E_Bwin,'*r')
 hold on
 plot(2,E_Awin,'+g') 
  
 
 
 