a = zeros(50, 1);
i=1;
while i < 51
    if i < 51;
    y = zeros(20,1);
    x = rand(20,1);
    y(x<0.05) = 1;
    a(i) = sum(y(:)==1);
    i = i + 1;
    else
    end
end

sum(a(:)==0)/50
sum(a(:)==1)/50
sum(a(:)==2)/50