In [ ]:
from pykinsol import solve
In [ ]:
def f(x, fout):
fout[0] = x[1] - 1 + 1/x[0]
fout[1] = x[1] + 1/(x[0] + 1)
def j(x, Jout, fx):
Jout[0, 0] = -1/(x[0]**2)
Jout[0, 1] = 1
Jout[1, 0] = -1/(x[0] + 1)**2
Jout[1, 1] = 1
In [ ]:
solve(f, j, [.5, -.5])