Matlab Factorizacion LU (sin pivotaje)

[skinner]

Buenas a todos. Estoy intentando programar la Factorización LU en Matlab. Intenté el siguiente código pero sorprendentemente no me da el resultado esperado:

function B=milu_sept(A)
dim=size(A); n=dim(1); %n filas=columnas

for col=1:(n-1)
    for fil=(col+1):n
        A(fil,col:n)=A(fil,col:n)-(A(fil,col)/A(col,col))*A(fil,col:n);
    end
end
B=A;
end

Según lo que yo creo, este código debería efectivamente devolverme la factorización LU de una matriz A cualquiera... ¿alguien puede ofrecerme ayuda?

Un saludo y muchas gracias!

[skinner]

Re: Factorizacion LU (sin pivotaje)

¡Ya lo he solucionado! Por si a alguien le interesa, el código correcto es (simplemente cambiar un "fil" por un "col"):

function B=milu_sept(A)
dim=size(A); n=dim(1); %n filas=columnas

for col=1:(n-1)
    for fil=(col+1):n
        A(fil,col:n)=A(fil,col:n)-(A(fil,col)/A(col,col))*A(col,col:n);
    end
end
B=A;
end

Un saludo!