I wonder if there is any good-working solver for equations like norm.cdf(x,0,1)-k in Python? Here we need to find x and k is known. I tried using fsolve, but it throws runti