% minst en sekser i fire kast med en terning
% minst en dobbel sekser i 24 kast med to terninger

hold on

ttt = cputime;

slutt = 20000;

stopp_her = [1 10 50 250 1000 1500 2500 4000 6000 10000 20000];

antall_gunstige_1 = 0;
antall_mulige_1 = slutt;

antall_gunstige_2 = 0;
antall_mulige_2 = slutt;

rf_1 = [0 0];
rf_2 = [0 0];

for i = 1:slutt,
  enkel_terning_fire_kast = ceil(6*rand(1,4));
  seksere_i_fire_kast = sum(enkel_terning_fire_kast == 6);
  if seksere_i_fire_kast >= 1,
    antall_gunstige_1 = antall_gunstige_1 + 1;
  end
  
  dobbel_terning_24_kast = ceil(6*rand(2,24));
  hjelp_dobbel_seks_i_24_kast = sum(dobbel_terning_24_kast == 6);
  dobbel_seks_i_24_kast = sum(hjelp_dobbel_seks_i_24_kast == 2);
  if dobbel_seks_i_24_kast >= 1,
    antall_gunstige_2 = antall_gunstige_2 + 1;
  end
  
  rf_1(2) = antall_gunstige_1/i;
  rf_2(2) = antall_gunstige_2/i;
  plot([i-1 i],rf_1,'r');
  plot([i-1 i],rf_2,'g');
  if sum(stopp_her == i) == 1,
    axis([0 i 0.4 0.6])
    plot([0 i],[.5177 .5177],'r--')
    plot([0 i],[.4914 .4914],'g--')
    pause
  end
  rf_1(1) = rf_1(2);
  rf_2(1) = rf_2(2);
end

cputime - ttt