SymPy对简单的线性系统给出了错误的答案？

Question

代码是

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在线shell，它也给出了相同的错误答案。

我在这里想念东西吗？

Answer 1

至少它是一致的。 x1（因为它未出现在方程组中）就像一个常数，并且如果x0和x2都等于该常数，则将满足第一个方程。来自solve的结果可能是首选：

>>> solve([x[0]-x[2], x[4]-x[3]], x)
{x3: x4, x0: x2}