code: private static double contClassicLB(final double numSamplesF, final double theta,
        final double numSDev) {
      final double nHat = (numSamplesF - 0.5) / theta;
      final double b = numSDev * Math.sqrt((1.0 - theta) / theta);
      final double d  = 0.5 * b * Math.sqrt((b * b) + (4.0 * nHat));
      final double center = nHat + (0.5 * (b * b));
      return (center - d);
    }
docs: our "classic" bounds, but now with continuity correction