N=10000;
dziewczynki=rand(4,N)>.5;
sum(sum(dziewczynki,1)>=3)