A=[7 1 3 1 2; 2 8 6 1 5; 2 4 8 1 2; 3 4 1 7 1; 2 3 5 1 3];
[L,U] = lu_rozklad(A)