
n = 100 ;
K = [ 0 : 100 ];
p = 1/2 ;

mk  = n*p ;
s2k = n*p*(1-p); 

Gauss     = exp( -(K-mk).^2 /(2*s2k) ) / sqrt(2*pi*s2k) ;
Bernoul   = bernoulli( K , n , p ) ;

plot( K , Gauss , 'or-.');

hold on ;

plot( K , Bernoul , 'sb-');

grid on ;

