python解不等式方程_在SymPy中一次解决非线性方程组和不等式系统

我是SymPy和python的新手，我遇到了一个问题 . 我正在尝试解决系统'kunSys'：

>>> kunSys

[-w0 + w1 - 8*x1 + 20,

-2*w0 + w2 - 8*x2 + 4,

w0*(-x1 - 2*x2 + 2),

w1*x1,

w2*x2,

w0 >= 0,

w1 >= 0,

w2 >= 0]

使用变量列表'lagVars'：

>>> lagVars

(x1, x2, w0, w1, w2)

如您所见，我的系统包含方程和不等式 .

Trying:

>>> solve(kunSys,lagVars)

Get:

NotImplementedError:

inequality has more than one symbol of interest

但是在分别解决方程和不等式时它可以正常工作：

>>> kunSys[:5]

[-w0 + w1 - 8*x1 + 20,

-2*w0 + w2 - 8*x2 + 4,

w0*(-x1 - 2*x2 + 2),

w1*x1,

w2*x2]

>>> solve(kunSys[:5],lagVars)

[(0, 0, 0, -20, -4),

(0, 1/2, 0, -20, 0),

(0, 1, -2,