On Fri, Jan 13, 2017 at 1:54 PM, Дмитрий Фролов <kain9...@gmail.com> wrote:
> Hellow! I have problem: I need to solve equation with one variable, but
> solution must be complex number. "find_root()" and "root()" in
> scipy.optimize can't do it. Help me please.

root() from scipy.optimize should work:

F.ex. to find the root of f(z) = z^2 + 1, write z = x + iy
so that f(x,y) = (x^2-y^2+1, 2*xy).
Following the example in:

import numpy as np
import scipy.optimize as opti

def func(x):
    f = [x[0]**2 - x[1]**2 + 1, 2*x[0]*x[1]]
    # The Jacobian:
    df = np.array([[2*x[0], -2*x[1]], [2*x[1], 2*x[0]]])
    return f, df

sol = opti.root(func, [0.5, 0.5], jac=True, method='lm')
print sol.x

Vegard Lima

You received this message because you are subscribed to the Google Groups 
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to sage-support+unsubscr...@googlegroups.com.
To post to this group, send email to sage-support@googlegroups.com.
Visit this group at https://groups.google.com/group/sage-support.
For more options, visit https://groups.google.com/d/optout.

Reply via email to