function y=gen_binom(N,p,n)
y=zeros(1,n);
for i=1:n
testy=rand(1,N);
y(i)=sum(testy<=p);
end

