return x - loc / (np.sqrt(c * c - 2 * logsf) - c)