初等变换
# 为了让矩阵可以更容易求逆，我们将其化为上三角矩阵
# 首先，将矩阵A化为上三角矩阵U：
U = np.copy(A)
for k in range(n-1):
for i in range(k+1, n):
factor = U[i, k] / U[k, k]
for j in range(k+1, n):
U[i, j] -= factor * U[k, j]
U[i, k] = 0
print("U=",U)
# 然后，对上三角矩阵U求逆得到A的逆矩阵invA：
invA = np.zeros((n, n))
for i in range(n-1,-1,-1): # 从最后一行开始倒序遍历整个上三角矩阵U。
invA[i][i]=float(1/U[i][i]) # 先将对角元素的倒数存储在invA中。
for j in range(0, i): # 接下来处理当前行的其余元素。
temp=0; # temp用来存储当前行j列所有右侧元素乘以对应位置的倒数之和。
for k in range (j+1 , i + 1): # 遍历当前行j列右侧所有元
免责声明：
以上内容除特别注明外均来源于网友提问，创作工场回答，未经许可，严谨转载。