为了让矩阵可以更容易求逆，我们将其化为上三角矩阵

# 首先，将矩阵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列右侧所有元