Sei Deine Matrix A, ok?
In der erste Spalte von A dividierst alle Spaltenwerte unterhalb der Diagonalen durch den negativen Wert der Diagonalen {1,1,6,-2}/(-1) , die berechneten Werte der Spalte schreibst Du in die erste Spalte der Einheitsmatrix, das sei L1
L1 : =⎝⎜⎜⎜⎛1−1−62010000100001⎠⎟⎟⎟⎞
L1⋅A=⎝⎜⎜⎜⎛100012−5−41−1−2213−3−6⎠⎟⎟⎟⎞
Jetzt wiederholt sich das Verfahren mit Spalte 2 und dem 2.Wert der Diagonalen von L1 A
{1,2,-5,-4}/(-2)
L2 : =⎝⎜⎜⎜⎛10000125200100001⎠⎟⎟⎟⎞
L2⋅L1⋅A=⎝⎜⎜⎜⎛100012001−1−29013290⎠⎟⎟⎟⎞
Jetzt hätten wir eine obere Dreiecksmatrix - genügt das so weit - wo willst Du hin?