function [ L,U ] = h2oLU( A )
U=A;
[rA cA]=size(A);
L=eye(rA);
for j=1:rA-1
    for k=j+1:rA
        L(k,j)= U(k,j)/U(j,j);
        U(k,:)=U(k,:)-U(k,j)/U(j,j)*U(j,:);
    end
end
end

