求解非線性方程組在Python（scipy.optimize.fsolve）的系統

Question

我試圖解決非線性方程的以下簡單系統（Source(second example)）：

(I) y - x^2 = 7 - 5x 
(II) 4y - 8x = -21 

應僅具有一個溶液（x = 3.5，y = 1.75）。

我目前使用scipy協議棧的方法如下：

from scipy.optimize import fsolve 

def equations(p): 
    x, y = p 
    return (y - x**2 -7 + 5*x, 4*y - 8*x + 21) 

x, y = fsolve(equations, (5, 5)) 

print(equations((x, y))) 

，併產生以下（這是沒有結果的）：

(0.0, 0.0) 

我已經嘗試過不同的起始預期，但它沒有提供正確的解決方案。

我的方法有什麼問題？我錯過了什麼嗎？

提前致謝！

Answer 1

這工作完全正常：

In [1]: %paste 
from scipy.optimize import fsolve 

def equations(p): 
    x, y = p 
    return (y - x**2 -7 + 5*x, 4*y - 8*x + 21) 

x, y = fsolve(equations, (5, 5)) 

print(equations((x, y))) 

## -- End pasted text -- 
(0.0, 0.0) 

In [2]: x 
Out[2]: 3.5000000414181831 

In [3]: y 
Out[3]: 1.7500000828363667 

equations(x, y)是(0, 0)意味着這兩個y - x**2 -7 + 5*x和4*y - 8*x + 21爲0

也許你糊塗了，並打算print(x, y)而不是print(equations(x, y))？