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])

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

我使用了一些打印语句，并在调试器中运行了几次，但无济于事。但是，对于简单的2 x 2矩阵，代码在此高斯消去部分上工作正常。