def normal_lower_bound(probability, mu=0, sigma=1):
 return inverse_normal_cdf((1 - probability), mu, sigma)

def normal_lower_bound ( probability , mu = 0 , sigma = 1 ) : return inverse_normal_cdf ( ( 1 - probability ) , mu , sigma )
returns the z for which p = probability .