function test_metodo_LU

warning off; % non scrive "warnings"
for n=5:5:20
    A=gallery('chebvand',n);b=rand(n,1); condA=condest(A);
    % risoluzione con "metodo_LU"
    tic; x_LU=metodo_LU(A,b);  t_metodo_LU=toc; % tempo impiegato
    % risoluzione con "backslash"
    tic; x_backslash=A\b; t_backslash=toc; % tempo impiegato
    % errore norma 2 soluzione
    err=norm(x_LU-x_backslash)/norm(x_backslash);
    % errore relativo soluzione backslash
    residuo_backslash=norm(b-A*x_backslash)/norm(b);
    % errore fattorizzazione LU
    [L,U,P]=lu(A); errLU=norm(P*A-L*U);
    
    fprintf('\n n: %3.0f cond: %1.3e err LU: %1.3e',n,condA,errLU);
    fprintf('\n err. LU vs backs.: %1.3e res. relativo: %1.3e',...
        err,residuo_backslash);
    fprintf('\n  tempo impiegato: LU: %1.3e backslash: %1.3e \n',...
        t_metodo_LU,t_backslash);
end
fprintf('\n')