
function pn = bernoulli( k , n , p)

%        pn = bernoulli( k , n , p )
%____________________________________________________________________
%
% Funzione che genera l'andamento della Distribuzione di Bernoulli 
% al variare del numero "n" di prove ripetute e del numero "k" di
% prove che hanno successo.
%
%    p = Probabilità che una prova ha successo
%
%
%               n!      k      n-k
%    Pn(k) = -------- p  (1-p)
%            k!(n-k)!
%
% N.B. 'k' può essere un vettore e Pn(k) il corrispondente vettore.
%____________________________________________________________________


for ind = 1 : length(k)
   
   
   BinCoeff(ind) = factorial(n)/factorial(k(ind))/factorial (n-k(ind)) ;
end

pn = BinCoeff .* p.^k .* (1-p).^(n-k) ;


