Standard normal cumulative probability distribution.
num stdNormalCdf(num x) { return 0.5 + 0.5 * erf(x * invSqrt2); }