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

    for i in range(n):
        L[i][0] = A[i][0]

    for j in range(n):
        U[0][j] = A[0][j] / L[0][0]

    for k in range(0, n):

        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