docs: function rnorm(mu, sigma) {
        return function() {
           var v1, v2, rsq;
  
           do {
              v1 = 2 * rgen.random() - 1;
              v2 = 2 * rgen.random() - 1;
              rsq = v1 * v1 + v2 * v2;
           } while (rsq >= 1);
           return mu + sigma * v1 * Math.sqrt(-2 * Math.log(rsq) / rsq);
        };
     }