function x = Eliminacion_Gaussiana(A, b)

N = length(b);

x = zeros(N,1);
    
for k =1:1:N-1
    for n = k+1:1:N
        b(n) = b(n) - A(n,k)*b(k)/A(k,k);
        for m=N:-1:k
            A(n,m) = A(n,m) - A(n,k)*A(k,m)/A(k,k);
        end;
    end;
disp(A)
disp(b)
end;

disp(A)
disp(b)

for k=N:-1:1
    suma = 0;
    for m=k+1:1:N
        suma = suma +A(k,m)*x(m);
    end;
    x(k) = (b(k)-suma)/A(k,k);
end;
