value = optimize.newton(func, x0, tol=1e-10)