docs: function rgeom(prob) {
        var scale;
  
        if (prob <= 0 || prob > 1) { return function() { return NaN; }; }
  
        scale = prob / (1 - prob);
  
        return function() {
           return rpois(rexp(scale)())();
        };
     }