Wir testen den Algorithmus mit einer 3·3-Matrix:

A = Matrix.fromString("5,4,2; 1,9,7; 3,0,6")
L, U = luRecursive(A)
print(L, U)

Zur Kontrolle vergleichen wir nun die Originalmatrix mit dem Produkt der beiden Dreiecksmatrizen:

print(A, L*U)

Fortsetzung: 3. Umformung in eine Iteration