function [L,U]=mylu0(A)
N=length(A);
for k=1:N-1
    for i=k+1:N
        A(i,k)=A(i,k)/A(k,k);
        for j=k+1:N
            A(i,j)=A(i,j)-A(i,k)*A(k,j);
        end
    end
end
U=triu(A);
L=eye(N)+tril(A,-1);