Cumulative distribution function: probability of a value <= x.
x
double cdf(double x) => 0.5 * erfc(-(x - mean) / (std * math.sqrt2));