function [L,U] = LUDeco2(A)

Na = size(A,2);

E = eye(Na,Na);

L =  eye(Na,Na);

for j=1:Na-1
    for i=j+1:Na
        if A(i,j)~=0
            E(i,j)=-A(i,j)/A(j,j);
            A=E*A;
            L=L*inv(E);
E = eye(Na,Na);
        else
            error('sdjfghskjg');
        end
    end
end
U=A(:,1:Na);