confronto fattorizzazioni - Università degli studi di Cagliari.

for n = 1:10;A = rand(n);A = A'*A;sol = ones(n,1);b = A*sol;%Comando "chol" di matlabR = chol(A);y1 = R'\b;x1 = R\y1;% Errore usando usando la funzione "chol" di Matlaberrore1(n) = norm(x1 – sol);%Comando "cholesky" implementato attraverso il mio algoritmo personalizzatoR=cholesky(A);y2=R'\b;x2=R\y2;% Errore usando usando la funzione "chol" di Matlaberrore2(n) = norm(x2 - sol);end7