python - 使用 SymPy 求解代数方程组

Question

我是 python 的 SymPy 的新手，正在尝试求解一个简单的方程组。我可以成功地计算变量“y”，但是当我尝试将这个值替换回“x”时，我无法让它简单地计算表达式并输出一个值。我的代码如下:

x = Symbol('x')
y = Symbol('y')
A = 1
B = 3
C = 6

x = solve(A*(1-x-y) + B*x + C*y - 4, x)
x = x[0]
y = solve(A*(1-x-y) + B*x - (8-C)*y, y)
y = y[0]

print x
print x.subs(y,0.5)

到目前为止，我所做的每一次尝试都只返回 -5*y/2 + 3/2。我试过将其转换为 float ，尝试以不同方式使用求解表达式，并将其转换为字符串，然后使用简化将其转换回要求解的表达式。

我所做的一切都行不通，我意识到这应该是一项简单的任务，但我无法弄清楚。任何建议都会有所帮助，谢谢!

Answer 1

不要忘记，您也可以使用 solve 一次解决这两个问题:

>>> x,y=symbols('x y')
>>> A = 1
>>> B = 3
>>> C = 6
>>> sol = solve((
...     A*(1-x-y) + B*x + C*y - 4,
...     A*(1-x-y) + B*x - (8-C)*y, ))
>>> sol
{x: 1/4, y: 1/2}
>>> print "x = {x} and y = {y}".format(x=sol[x],y=sol[y])
x = 1/4 and y = 1/2