

function [  ] = GaussianaSimple ()

% Se introducen la matriz y el vector de constantes 
A=matrizA;
b=vectorb;
[m n]=size(A);
% Se halla la matriz aumentada 
G=[A b];

if m=n 

if det(A)~=0

	for i=1:m-1 	% El primer contador recorre las filas 
		for k=i+1:m		% El segundo contador recorre las filas 
		aux=-G(k,i)/G(i,i);	% Número que convierte en ceros las posiciones debajo del pivote
			for j=1:n+1 		% El tercer contador recorre las columnas 
				G(k,j)=G(k,j)+(aux*G(i,j)); % i es constante mientras k y j recorren todas las posiciones  
			endfor	
		endfor
		G
	endfor

	% Se define una matriz X vacía 
	for t=1:m
		X(t,1)=0; 
	endfor

	for h=n:-1:1	% Se hace una sustitución regresiva 
		suma=0;		% Esta suma es un elemento de ayuda para el despeje 
		for p=1:n
			suma=suma+G(h,p)*X(p,1);
		endfor
		X(h,1)=(G(h,m+1)-suma)/G(h,h); % Despeje de la variable en la posición (h,1)
	endfor
	disp('');
	disp('MATRIZ AUMENTADA G');
	disp('');
	disp([A b]);
	disp('');
	disp('VECTOR X ');
	disp('');
	disp(X);
	disp('');
	else
	disp('el sistema no es cuadrado y no tiene unica solucion')
	end

	else
		disp('el sistema no es cuadrado y no tiene unica solucion')
	end


endfunction
