用Python进行高斯消元法

Question

所以我一直在研究一个自定义函数，它使用LU分解来求解线性方程组。该函数假定等式和未知数相等。当我使用

    for k in range(0,n-1,1):
        for i in range(k+1,n,1):
            L[i,k] = A[i,k] / A[k,k] #Lower matrix, comprehension aid
            for j in range(0,n,1):
                A[i,j] = A[i,j] - (L[i,k] * A[k,j])
                d[i] = b[i] - (L[i,k] * b[j])

在一个3x3系数矩阵（a）where A = np.array([[8,4,-1],[-2,5,1],[2,-1,6]]) the A[1,2]值在我期望的0.75时变为0。不知道发生了什么。在

我使用了一些print语句，并在调试器中运行了几次，但都没有用。然而，对于一个简单的2x2矩阵，代码在这个高斯消去部分工作得很好。在