expect_numerator = (math.log(underlying / strike)
                            + (rate + vol * vol * 0.5) * maturity)