z1 = (x1 - mean) / std_dev z2 = (x2 - mean) / std_dev
probability = 0.6827 # approximately
p2 = 1 - normal_cdf(x2, mean, std)