%LU razcep - demonstracija algoritma 2.4.1

n = 4;
A = rand(n)
b = rand(n,1)
Aorig = A;
borig = b;

alg241

L = tril(A,-1) + eye(n);
U = triu(A);
alg221
alg222
ostanek = Aorig*b - borig
nost = ostanek/eps