def bs_delta_to_strike(fwd, delta, vol, texp): ys = norm.ppf(delta) return fwd/math.exp((ys - 0.5 * vol * math.sqrt(texp)) * vol * math.sqrt(texp))