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