def BD_ratio(params, zp): n = params['nB'] bnb = S.get_b_n(n) bnd = S.get_b_n(params['nD']) pre = n * gamma(2*n) * np.exp(bnb) / (bnb ** (2 * n)) h = params['ReD'] / ((bnd) ** params['nD']) Ie = 10 ** ((params['MeD'] - zp) / -2.5) I0 = Ie * np.exp(bnd) return pre * ((params['ReB'] / h)**2.) * (Ie / I0)