eserc5

% Esercitazione 5
% risoluzione di un sistema lineare Ax=b
% (A matrice casuale nxn)
% mediante Gauss senza scambio di righe

clc, clear;

fprintf('\n Risoluzione del sistema Ax=b (A nxn)\n');

n=input(' Dai la dimensione del problema: ');

% creazione del problema
A=rand(n);
A(1,1)=0;
x_esatto=ones(n,1);
b=A*x_esatto;
A=[A,b];

% applicazione di Gauss;
for k=1:n-1
    for i=k+1:n
        m = A(i,k)/A(k,k);
        for j=k+1:n+1
            A(i,j)=A(i,j)-m*A(k,j);
        end
    end
end
disp(A);

%risoluzione del sistema triangolare superiore

x(n)=A(n,n+1)/A(n,n);
for i=(n-1):-1:1
    x(i)=A(i,n+1);
    for j=(i+1):n
        x(i)=x(i)-A(i,j)*x(j);
    end
    x(i)=x(i)/A(i,i);
end

% output

errore_relativo=norm(x_esatto-x')/norm(x_esatto);
fprintf('\n errore=%e\n\n ', errore_relativo);