牛顿法解非线性方程组python

牛顿法是一种用于解非线性方程组的迭代方法。在Python中，可以使用scipy库的optimize模块中的root函数来实现牛顿法解非线性方程组。以下是一个示例代码：

from scipy.optimize import root

# 定义一个包含多个方程的非线性方程组
def equations(x):
    # 方程1: x + y - 3 = 0
    eq1 = x[0] + x[1] - 3
    # 方程2: x^2 + y^2 - 4 = 0
    eq2 = x[0]**2 + x[1]**2 - 4
    return [eq1, eq2]

# 使用牛顿法求解非线性方程组
solution = root(equations, [0, 0])
x = solution.x
# 打印解
print("x =", x[0])
print("y =", x[1])