Algoritmo de Gauss - Jordan

Se denomina matriz escalonada a una matriz en la que las filas posteriores a una fila cuyos elementos son todos ceros, tienen todos sus elementos igual a cero, y el número de elementos nulos al comienzo de cada fila no nula es estrictamente menor que en la siguiente. El proceso de eliminación Gausseana viene establecido por el siguiente algoritmo:

A efecto de implementarlo en MatLab podemos establecer el método a través del siguiente Script cuyas entradas sería la Matriz aumentada del sistema de ecuaciones y como salida los valores de cada variable:

clear all ;
clc;

a=[2 1 -1 8; -3 -1 2 -11; -2 1 2 -3]
[f,c] = size(a)

for k=1:c-1
  a(k,:)=a(k,:)/a(k,k)
    for j=k+1:f
    a(j,:)=a(j,:)-a(k,:)*a(j,k)
    j=j+1;
    a
    pause
    end
  k=k+1;
  a
  pause
end

for k=f:-1:2
  for j=k-1:-1:1
    a(j,:)=a(j,:)-a(k,:)*a(j,k)
    j=j-1;
    a
    pause
  end
  k=k-1;
  a
  pause
end

fprintf('Result\n');

Fuente:

• Métodos Numéricos para Ingeniería
• Ing. Ricardo Seminario Vasquez
• Enlace
• Métodos de Resolución de Ecuaciones Iterativos
• Wladimiro Diaz Villanueva 
• Enlace