def bm_bn_estimate(nn): "Guess for bn constant defined by B+M" if np.iterable(nn): return np.array([bm_bn_estimate(n) for n in nn]) est = 0.87*nn-0.15 if nn < 0.2: est = 0.01*(nn/0.1)**4 return est