def submatrix(A,i,j):
#矩阵A第i行第j列元素的余矩阵
p=len(A)#矩阵的行数
q=len(A[0])#矩阵的列数
C=[[A[x][y] for y in range(q) if y!=j] for x in range(p) if x!=i]#列表推导式
return C
def det(A):
#按第一行展开递归求矩阵的行列式
p=len(A)#矩阵的行数
q=len(A[0])#矩阵的列数
if(p==1 and q==1):
return A[0][0]
else:
value=0
for j in range(q):
value+=((-1)**(j+2