I've been told that the method scipy.optimize.newton() will solve complex functions so long as the first derivative is provided. I can't make it work. The documentation for newton() mentions nothing of complex functions. Could someone show me how one would find the roots of a function like f(z) = 1 + z^2 in SciPy? I need to solve something much more complex, but a simple example will help me immensely.