First somewhat working version of a binomial sampling function (for numpy and weave)