// Stabilite retrograde de la factorisation de Cholesky

W=triu(rand(50,50));
A=W'*W;

R=chol(A);
R1=R+1e-6*triu(rand(50,50));

// Ni R ni R1 ne sont proches de W
norm(W-R)/norm(W)
norm(W-R1)/norm(W)

// Mais R'*R est beaucoup plus proche de A que R1 * R1'
norm(A-R'*R)/norm(A)
norm(A-R1'*R1)/norm(A)