function gamma = randbin(lambda,N)

% questa funzione genera una variabile binomiale con valori {0,1} con
% P[gamma=1]=lambda

gamma = (rand(1,N) < lambda);

