Python中的二维二分法？

我正在寻找Python中的

optimize.bisect

函数的扩展，该函数允许我解两个耦合自洽方程组，即

y=f（x，y）

和

x=g（x，y）

。我在寻找一种二维二等分方法，但如果它运行良好，我不介意使用其他可能性。

据我所知，您希望找到以下函数的根

def fun(x):
    return [f(x[0], x[1]) - x[1],
            g(x[0], x[1]) - x[0]]

其中

x[0]

对应于定义中的

x

，而

x[1]

对应于

y

。 您可以使用以下方法解决此问题：

---

这里

[0，0]

是根的初始猜测。然而，选择合适的算法在很大程度上取决于您的问题。

您能举个例子吗

from scipy import optimize
solution = optimize.root(fun, [0, 0])