A = [1,2,1,-3,2; 3,7,5,-7,10; -2,-2,3,13,0; -1,2,5,6,24; 4,10,9,-4,13];
b = [-8; -24; 9; 28; -35];

A = [A b];

U1 = [A b];

#zamiana 1 i 5 wiersza
U1([5 1], :) = U1([1 5], :);

#zerowanie kolumny 1
U1(2, :) = U1(2, :) - (U1(2, 1) / U1(1, 1))*U1(1, :);
U1(3, :) = U1(3, :) - (U1(3, 1) / U1(1, 1))*U1(1, :);
U1(4, :) = U1(4, :) - (U1(4, 1) / U1(1, 1))*U1(1, :);
U1(5, :) = U1(5, :) - (U1(5, 1) / U1(1, 1))*U1(1, :);

#wartość maksymalna w 2 kolumnie (od wartości bezwzględnej)
maks = max(abs(U1(2:5, 2)));
#zamiana 2 i 4 wiersza
U1([2 4], :) = U1([4 2], :);

#zerowanie 2 kolumny:
U1(3, :) = U1(3, :) - (U1(3, 2) / U1(2, 2))*U1(2, :);
U1(4, :) = U1(4, :) - (U1(4, 2) / U1(2, 2))*U1(2, :);
U1(5, :) = U1(5, :) - (U1(5, 2) / U1(2, 2))*U1(2, :);

#wartość maksymalna w 3 kolumnie (od wartości bezwzględnej)
maks = max(abs(U1(3:5, 3)));
#zerowanie kolumny
U1(4, :) = U1(4, :) - (U1(4, 3) / U1(3, 3))*U1(3, :);
U1(5, :) = U1(5, :) - (U1(5, 3) / U1(3, 3))*U1(3, :);

#wartość maksymalna w 4 kolumnie (od wartości bezwzględnej)
maks = max(abs(U1(4:5, 4)));
#zerowanie 4 kolumny:
U1(5, :) = U1(5, :) - (U1(5, 4) / U1(4, 4))*U1(4, :);

b1 = U1(:, 6);

#majac macierz wspolczynnikow, mozna znalexc macierz L1
#deklaracja macierzy jednostkowej
L1 = eye(5, 5);

##A1 = [1  2  1 -3  2
##      3  7  5 -7 10
##     -2 -2  3 13  0
##     -1  2  5  6 24
##      4 10  9 -4 13];