N = 5; p = 1/20; % number of events and the individual probability
x = [0:5];       % values for which the pdf is to be computed
prob = binopdf(x,N,p)  % values of the binomial distribution
stem(x,prob)            % generate the stem plot
xlabel('Anzahl der 20'); ylabel('Wahrscheinlichkeit')
xlim([-0.5 5.5])
