docs: @Reference(authors = "G. Marsaglia", //
        title = "Evaluating the Normal Distribution", //
        booktitle = "Journal of Statistical Software 11(4)", //
        url = "https://doi.org/10.18637/jss.v011.i04", //
        bibkey = "doi:10.18637/jss.v011.i04")
    public static double cdf(double x, double mu, double sigma) {
      x = (x - mu) / sigma;
      if(x >= 8.22) {
        return 1.;
      }
      if(x <= -8.22) {
        return 0.;
      }
      if(x != x) {
        return Double.NaN;
      }
      double s = x, t = 0, b = x, q = x * x, i = 1;
      while(s != t && i < 1000) {
        t = s;
        s += (b *= q / (i += 2));
      }
      // Constant is 0.5 * log(2*pi)
      return .5 + s * FastMath.exp(-.5 * q - .91893853320467274178);
    }