return np.log(c) - x - (c+1.0)*sc.log1p(np.exp(-x))