function A=lugauss(A)
%LUGAUSS LU factorizacion sin pivoteo.
% A = LUGAUSS(A) almacena una matriz triangular
% superior en la parte triangular superior de A
% y una matriz triangular inferior en la parte
% estrictamente inferior de A (los elementos
% diagonales de L son 1).
[n,m]= size(A); 
if n ~= m; error('A no es una matriz cuadrada '); else
for k = 1:n-1
for i = k+1:n
A(i,k) = A(i,k)/A(k,k);
if A(k,k) == 0 , error('Elemento diagonal nulo '); end
j = [k+1:n]; A(i,j) = A(i,j) - A(i,k)*A(k,j);
end
end
end
return
