...81-100 blue balls
npop = 100000 #population size
nsamp = 200 #sample size
niter = ma.floor(npop/nsamp)
population = ar.array("i",npop*[0]) #population
red = ar.array("i",niter*[0])
blue = ar.array("i",niter*[0])
pcred = ar.array("d",niter*[0])
cd = ar.array("d",niter*[0])
for i in...