docs: private double getExpectedProbability(int events)
      {
          BigDecimal kFactorial = new BigDecimal(Maths.bigFactorial(events));
          double numerator = Math.exp(-mean) * Math.pow(mean, events);
          return new BigDecimal(numerator).divide(kFactorial, RoundingMode.HALF_UP).doubleValue();
      }