Hi, I am trying to solve the follow equation: solve(x * ((1.0 / sqrt(2 * pi) * x * sigma) * exp(-0.5 * (log(x) - mu)**2 / sigma**2)) + 0.5 * erfc((mu - log(x)) / (sigma * sqrt(2))) - 1, x) I am not sure Scipy can do it and how it can do it ? Many thanks Xiong