我有一些算法可以找到矩阵的行列式并将其预转换为上三角矩阵形式：

clear all;
A = rand([3 3])*100;
B = A;
N = size(A);
for K = 1 : N-1
    for I = K+1 : N
        for J = K:N
            if J == K
                aik = A(I,K);
            end
            A(I,J) = (A(I,J) - aik/A(K,K)*A(K,J));
        end
    end
end
det = 1;
for I = 1:N
    det = det*A(I,I);
end

这段代码运行得很好。现在我想开发一个类似的算法，但对于下三角矩阵的情况。我怎么做呢？重要的是，我需要使用MATLAB，而不需要调用任何外部函数。