代码是
from sympy import *
x = symbols('x0:5')
linsolve([x[0]-x[2], x[4]-x[3]], x)
这应该给出类似的东西
(x2, x1, x2, x4, x4)
作为答案，但它给出了
(x1, x1, x2, x4, x4)
我还尝试了 live.sympy.org 在线外壳，它也给出了同样的错误答案。
我在这里错过了什么吗？
至少它是一致的。x1，因为它没有出现在方程集中，所以就像一个常数，如果x0和x2都等于该常数，那么第一个方程将得到满足。solve的结果可能是首选：
>>> solve([x[0]-x[2], x[4]-x[3]], x)
{x3: x4, x0: x2}