行列式计算with python

问题源于同济大学《线性代数》（第五版），先看理论：

 

 问题是这样的：

教材上的算术求解方法如下：

接下来我们用python编程求解，代码如下：

import numpy as np

#定义矩阵
D=np.array([[2,1,-5,1],[1,-3,0,-6],[0,2,-1,2],[1,4,-7,6]])
#输出矩阵D行列式运算的值
d=np.linalg.det(D)
print "D=",d

D1=np.array([[8,1,-5,1],[9,-3,0,-6],[-5,2,-1,2],[0,4,-7,6]])
#输出矩阵D1行列式运算的值
d1=np.linalg.det(D1)
print "D1=",d1

D2=np.array([[2,8,-5,1],[1,9,0,-6],[0,-5,-1,2],[1,0,-7,6]])
#输出矩阵D2行列式运算的值
d2=np.linalg.det(D2)
print "D2=",d2

D3=np.array([[2,1,8,1],[1,-3,9,-6],[0,2,-5,2],[1,4,0,6]])
#输出矩阵D3行列式运算的值
d3=np.linalg.det(D3)
print "D2=",d3

D4=np.array([[2,1,-5,8],[1,-3,0,9],[0,2,-1,-5],[1,4,-7,0]])
#输出矩阵D4行列式运算的值
d4=np.linalg.det(D4)
print "D2=",d4

#求解x1,x2,x3,x4
x1=d1/d
x2=d2/d
x3=d3/d
x4=d4/d
print "x1=",x1
print "x2=",x2
print "x3=",x3
print "x4=",x4