def Crout(A):
    n = len(A)
    U = np.zeros((n, n))
    L = np.zeros((n, n))

    for k in range(0, n):
        U[k][k] = 1

        for j in range(k, n):
            ss = 0
            for s in range(0, k):
                ss += U[s][k] * L[j][s]
                L[j][k] = A[j][k] - ss

        for j in range(k + 1, n):
            ss = 0
            for s in range(0, k):
                ss += L[k][s] * U[s][j]
                U[k][j] = (A[k][j] - ss) / L[k][k]

    return L, U