ilosc_gier=100; powt=10000;
wyniki=zeros(40,1);
for i=1:powt
 a= sum(sum(rand(ilosc_gier,5)>5/6, 2)>=2);
 wyniki(round(a))=wyniki(round(a))+1;
end
ss=stem(wyniki); set(gca,'tickdir','out')
sum(wyniki(21:40))/powt
