function x = sol(a,b)

%SOL(A,b) Trova la soluzione di un sistema lineare 
%   A*x = b utilizzando la fattorizzazione di Gauss 
%   con pivoting per righe.

[l,u,p] = lu(a);
n=length(b);
b=p*b;
y(1)=b(1);for i=2:n;y=[y;b(i)-l(i,1:i-1)*y(1:i-1)];end;
x=y(n)/u(n,n);for i=n-1:-1:1;x=[(y(i)-u(i,i+1:n)*x)/u(i,i);x];end;
end