% GAUSS
%
%
% X=GAUSS(A,B)
%
% Devuelve en X la solucion del sistema AX=B 
% con el metodo de Gauss sin pivotaje
%
% Extraido del texto:
%  
%   "Matlab en cinco lecciones de numerico" 
%    por V. Dominguez y M.L. Rapun.
% 
% Más informacion en 
%
%   http://www.unavarra.es/personal/victor_dominguez
function x = gauss(a,b)

n=length(a);
% transformacion del sistema en uno triangular
for i=1:n-1
    for k=i+1:n
        l=a(k,i)/a(i,i);
        for j=i+1:n
            a(k,j)=a(k,j)-l*a(i,j);
        end
        b(k)=b(k)-l*b(i);
    end
end 

% resolucion del sistema triangular
x=zeros(n,1); % tambien vale x=b*0;
x(n)=b(n)/a(n,n);
for i=n-1:-1:1
    s=0;
    for j=i+1:n
        s=s+a(i,j)*x(j); % sumatorio
    end
    x(i)=(b(i)-s)/a(i,i);
end

return