docs: @Reference(title = "Fast and accurate computation of binomial probabilities", //
        authors = "C. Loader", booktitle = "", //
        url = "http://projects.scipy.org/scipy/raw-attachment/ticket/620/loader2000Fast.pdf", //
        bibkey = "web/Loader00")
    private static double stirlingError(int n) {
      // Try to use a table value:
      if(n < 16) {
        return STIRLING_EXACT_ERROR[n << 1];
      }
      final double nn = n * n;
      // Use the appropriate number of terms
      if(n > 500) {
        return (S0 - S1 / nn) / n;
      }
      if(n > 80) {
        return ((S0 - (S1 - S2 / nn)) / nn) / n;
      }
      if(n > 35) {
        return ((S0 - (S1 - (S2 - S3 / nn) / nn) / nn) / n);
      }
      return ((S0 - (S1 - (S2 - (S3 - S4 / nn) / nn) / nn) / nn) / n);
    }