对一个线性方程组的求解（用sympy进行符号运算）

sympy是python的数学符号计算库，用它可以进行数学表达式的符号推导和演算。与一些成熟的专业符号运算软件相比，sympy功能以及速度比不上，但是对于一些计算工具来说，使用简单，体积小。

举个简单的例子

from sympy import *

a, b, c, x = symbols("a, b, c, x")
solve(a*x**2 + b*x + c, x)

求解一个一元二次方程

输出结果为 [-(b - sqrt(-4*a*c + b**2))/(2*a), -(b + sqrt(-4*a*c + b**2))/(2*a)]

可以用factor转换为标准的数学公式

正题

四元一次方程组

c1*x1 + c2*x2 - c5*x4 - D1=0
c4*x1 - c1*x2 - c5*x3 - D2=0
c3*x2 + c1*x3 + c6*x4 - D3=0
c3*x1 + c6*x3 - c1*x4 - D4=0

from sympy import *
init_printing()
c1, c2, c3, c4, c5, c6, D1, D2, D3, D4 = symbols('c1, c2, c3, c4, c5, c6, D1, D2, D3, D4')
x1, x2, x3, x4 = symbols('x1 x2 x3 x4 ')

result = solve([c1*x1 + c2*x2 - c5*x4 - D1,
           c4*x1 - c1*x2 - c5*x3 - D2,
           c3*x2 + c1*x3 + c6*x4 - D3,
           c3*x1 + c6*x3 - c1*x4 - D4],[x1, x2, x3, x4])
print result