要使用SciPy解決非線性方程組,可以使用scipy.optimize
模塊中的root
函數。以下是一個簡單的示例代碼,演示如何使用SciPy解決非線性方程組:
import numpy as np
from scipy.optimize import root
# 定義一個包含兩個未知數的非線性方程組
def func(x):
return [x[0] + 2*x[1] - 2, x[0]**2 + 4*x[1]**2 - 4]
# 初始猜測值
x0 = np.array([0, 0])
# 使用root函數求解非線性方程組
sol = root(func, x0)
# 輸出結果
print(sol.x)
在這個示例中,func
函數定義了一個包含兩個未知數的非線性方程組。我們使用root
函數來求解該方程組,并提供一個初始猜測值x0
。最后,打印出求解得到的結果sol.x
,即方程組的解。