function [L,U, A] = MyLU(A)
%
%  LU factorization without partial pivoting
%  9/15/2006
%

[m,n] = size(A);
for j = 1:m-1
    for k = j+1:m
        alpha = A(k,j)/A(j,j);
        A(k,j) = alpha;
        for i = j+1:n
            A(k,i) = A(k,i)-alpha*A(j,i);
        end
    end
end

U = zeros(size(A));
for i=1:m
    for j = i:n
        U(i,j)=A(i,j);
    end
end

L = eye(m);
for i=2:m
    for j= 1:(i-1)
        L(i,j) = A(i,j);
    end
end