x = 0:6; n = 100, p = 0.01;  % number of events and individual probability
prob  = binopdf(x,n,p)       % generate the pdf
cprob = binocdf(x,n,p)       % generate the cdf
stem(x,prob,'b');    hold on  % plot both
stairs(x,cprob,'k'); hold off
xlabel('Anzahl fehlerhafter Schrauben in 100er Schachtel')
ylabel('Wahrscheinlichkeit')
legend('pdf','cdf','location','northwest')
