複数の引数を用いたscipy.optimize.rootの使い方が分かりません

Question

例えば、
x - y = z + w
x + y = z - w
という連立方程式を考えるとして、(z, w) = (2, 1)と設定すれば解は(x, y) = (2, -1)となります。

これと同様のものをpythonライブラリのscipy.optimize.rootを使って解きたいのですが、全く違う解が出力されてしまいます。

def func(x, y):
    return [x[0]-x[1] -y[0] - y[1],
           x[0]-x[1] -y[0] + y[1]]
    
z=[0, 0]
z=optimize.root(func,(0, 0), args=(2, 1), method="lm").x
print(z)

[51. 49.]

コードで間違えているところがあれば教えてほしいです。
恐らくargsの使い方が誤っているかと思います。
宜しくお願いいたします。

Answer 1

import numpy as np
from scipy import optimize

def func(x, y):
    return [x[0]-x[1] -y[0] - y[1],
           x[0]-x[1] -y[0] + y[1]] # x[0] + x[1] では?

z=[0, 0]
z=optimize.root(func,[0, 0], args=[2, 1], method="lm").x
print(z)