function b=binom(k, p)
% obliczanie i rysowanie rozk³adu dwumianowego dla danych k i p
for ii=1:k+1
    i=ii-1;
    b(ii)=p^i*(1-p)^(k-i)*factorial(k)/(factorial(i)*factorial(k-i));
end
ss=stem(0:k, b); set(ss(1), 'markersize', 2)
box off; set(gca,'xlim', [-.5 k+.5]);
title(sprintf('p=%.3f, n=%d', p, k));
