#a = randi (6004799503160661, 1e6, 1);
#a = randi (9,1e6,1);
#a = randi_rej (6004799503160661);
#a = randi_rej (9);
#a = randi_new (6004799503160661);
a = randi_new (9);

b = a(a>median(a));
c = hist(mod(b,2),[0 1])
c ./ sum (c)
