資源簡介
在scipy.optimize工具包中有兩個函數可以數值求解方程組，分別是root和solve，這兩個函數會找到方程組的一個近似解。
代碼片段和文件信息
#?方法1
#?==============================
from?scipy.optimize?import?root
def?f1（x）:
????return?[x[0]*（1j）+x[0]*x[1]+1x[0]+x[1]-1j]
print（root（f1[11]method=‘krylov‘）.x）
print（root（f1[11]method=‘krylov‘tol=1e-10）.x）#設置能夠允許的誤差為10的-10次方
#?==============================
from?scipy.optimize?import?fsolve
from?numpy?import?arraymat
def?f1（x）:
??return?[x[0]+x[0]*x[1]-2x[0]-x[1]-2]
def?jac1（x）:#方程組對應的雅可比矩陣
??return?mat（[[1+x[1]x[0]][1-1]]）
print（fsolve（f1[0-1]））#初始猜測值[0-1]
評論
共有 條評論
川公網安備 51152502000135號