% Este programa Ilustra Eliminación Gaussiana
% que consiste en:
%          i) eliminación hacia adelante y
%         ii) sustitución hacia atrás.
%
clear
a=[-0.04  0.04  0.12 3;
   0.56  -1.56  0.32 1;
   -0.24 1.24  -0.28 0]
%
x=[0 0 0];  %  x se inicializa como vector columna.
tempo= a(2,:); a(2,:)=a(1,:); a(1,:)=tempo;
%
a(2,:) = a(2,:) - a(1,:)*a(2,1)/a(1,1);
a(3,:) = a(3,:) - a(1,:)*a(3,1)/a(1,1); a
%
tempo = a(3,:); a(3,:) = a(2,:); a(2,:)=tempo ;a
%
a(3,:) = a(3,:) - a(2,:)*a(3,2)/a(2,2);a
%
x(3) = a(3,4)/a(3,3);
x(2) = (a(2,4) - a(2,3)*x(3))/a(2,2);
%x(1) = (a(1,4) - a(1,2:3)*x(2:3))/a(1,1);x