double cdf(double x) { final double inner = 1 + math.erf((x - mu) / (sigma * math.sqrt(2))); return inner.abs() / 2; }