class MoreMath::NormalDistribution

for the value +x+.
Returns the cumulative probability (p-value) of the NormalDistribution

def probability(x)
  0.5 * (1 + erf((x - @mu) / (@sigma * ROOT2)))
end