function [L U] = LUslow(A)

% Compute the LU decomposition for the provided matrix, A, in which L is unit diagonal.
% Assume that A can be reduced to row-echelon form without row interchanges.

return;
