Pompiao con MatLab

Bueno, aqui voy otra vez con mis cosas de geek... Resulta que este semestre me ha tocado llevar un curso de analisis numerico que me encanta. En este curso estamos utilizando MatLab como herramienta de estudio y tengo que decir WoW!!! Estoy super pompiao con el programita, esta super! Bueno, asumo que de esto se habran dado cuenta ya porque ultimamente he publicado casi todo el codigo que he escrito para esa clase y algun otro poco que he hecho por mi cuenta. Bien pues en este caso, les presento mi implementación de la factorizacion LU de una matriz A. No voy a entrar en los detalles de las ventajas y/o desventajas del mismo, ni en los detalles de la matemática detrás, pero si a alguien le interesa saber más sobre el tema, pueden dejar un comentario y puedo añadir más detalles, mas por el momento lo dejo hasta aquí. Espero les agrade.

%*************************************************************************
%Este programa halla la factorizacion LU de una matriz A si es *
%que esta existe o sea, el programa halla una matriz triangular *
%inferior y una matriz triangular superior tal que LU = A. *
% *
% *
%Input: La matriz A a ser factorizada *
% *
%Output: Las matrices L y U y un flag indicando si fue exitosa o no la *
% factorización. *
%*************************************************************************

function [flag,L,U] = LUFac(A)

flag =1; %Bandera indicando que no ha ocurrido ningun error
col=1; %Columna pivote
row=1; %fila pivote
nextRow = row+1;
[A_rowCant,A_colCant] = size(A); %tamano de A
L = eye(A_rowCant,A_colCant);

tic; %Comenzar a contar el tiempo de ejecución
for col=1:(A_colCant)
for rowIndex=1:(A_rowCant-row)
if abs(A(row,col)) < eps
flag = 0;
'Division by zero error'
return;
end
L(nextRow,col) = -(-A(nextRow,col)/A(row,col));
A(nextRow,:) = A(row,:)*(-A(nextRow,col)/A(row,col))+A(nextRow,:);
nextRow = nextRow + 1;
end
row = row+1;
col = col+1;
nextRow = row+1;
end

U=A;
toc;

return