def _cdf(z): """Cumulative density function (CDF) of the standard normal distribution.""" return .5 * (1. + erf(z/np.sqrt(2.)))