
function test_metodo_LU3

n=100;

% definizione matrice A
for s=1:n
    for t=1:n
        if s <= t
            A(s,t)=s/t;
        else
            A(t,s)=t/s;
        end
    end
end

% definizione termine noto
for k=1:n
        if rem(k,2) == 0 %indice pari
            b(k,1)=+1;
        else
            b(k,1)=-1;
        end
end

% soluzione sistema lineare
x=metodo_LU(A,b);

% valutazione residuo relativo
err=norm(b-A*x)/norm(b);
fprintf('\n \t residuo relativo: %1.1e \n',err);
