25.03.2015, 12:23, Camille Chambon kirjoitti: [clip] > Both work fine. But using timeit(...), the solution with > scipy.optimize.root is twenty times slower than the one with the while > loop. Is it normal? The default accuracy in optimize.root is 10^{-8}. Maybe your epsilon was much bigger?