% Razcep Choleskega - demonstracija algoritma 2.5.1

n=5;
A = rand(n);
Rorig = triu(A);
A = Rorig'*Rorig; Aorig = A;

alg261

R = triu(A)
ost = R'*R-Aorig
nost = ost/eps