docs: public static double quantile(double p, double alpha, double beta) {
      return (Double.isNaN(alpha) || Double.isNaN(beta) || Double.isNaN(p) || alpha < 0. || beta < 0. || p < 0 || p > 1) ? Double.NaN : //
          p == 0 ? 0. : p == 1 ? 1. : p > 0.5 //
              ? 1 - rawQuantile(1 - p, beta, alpha, logBeta(beta, alpha)) //
              : rawQuantile(p, alpha, beta, logBeta(alpha, beta));
    }