function gamma = randbin(lambda)

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

if rand < lambda
    gamma = 1;
else
    gamma = 0;
end
