%muchy
wektor=zeros(20,1);
wektor(1:10)=1;
N=1000;
liczba_1=zeros(N,1);
wektor_nowy=zeros(20,1);
for i=1:N
wektor_nowy=randsample(wektor,20,true);    
liczba_1(i)=sum(wektor_nowy);
end
[n,x]=hist(liczba_1,1:20);
bar(x,n/N)
prawd_symul=(sum(liczba_1<=6)+sum(liczba_1>=14))/N
prawd_dwum=sum(binopdf([0 1 2 3 4 5 6 14 15 16 17 18 19 20],20,0.5))