factor2 = (math.sqrt(2.0 * math.pi / maturity)
               / (underlying + discount_strike))