N=10000;
na10=sum(rand(N,10)<0.6, 2);
ss=stem(hist(na10,1:10));
set(ss(1),'markersize',1.5); set(gca,'tickdir','out')
sum(na10>=9)
