code: function(distribution, coverageAlpha) {
          if (distribution.numTrials < 1000) {
              // don't even bother trying to optimize for small-ish sample sizes
              return [0, distribution.numTrials];
          } else {
              return [
                  Math.floor(distribution.inverseCdf(coverageAlpha / 2)),
                  Math.ceil(distribution.inverseSurvival(coverageAlpha / 2))
              ];
          }
      }
docs: /* For the given binomial distribution, compute an interval that covers at least
  (1 - coverageAlpha) of the total probability mass, centered at the expectation (unless we're
  at the boundary). Uses the normal approximation.