Gauss C++ (codigo)

% A es la matriz de coeficientes
% b es el vector de terminos independientes
% x es el vector solucion
% n es el orden del sistema

[n,m] = size(A);

% eliminacion
for i=1:n
  for k=i+1:n
    factor = (A(k,i)/A(i,i));
    A(k,i:n) = A(k,i:n) - factor*A(i,i:n);
    b(k) = b(k) - factor*b(i);
  end
end

%sustitucion inversa
x = zeros(1,n);
x(n) = b(n)/A(n,n);
for k=n-1:-1:1
  x(k) = (b(k) - sum(x(k+1:n).*A(k,k+1:n)))/A(k,k);
end
x