function [L,U] = my_lu2(eps)

%   [ 1   1;
%     eps -1 ]
    
L(1,1) = 1; L(2,2) = 1;
L(2,1) = single(eps);
    
U(1,1) = 1;
U(1,2) = 1;
U(2,2) = -1-single(eps);

end


% try
% eps = 10^7; A = [1 1; eps -1]; [L,U] = my_lu2(eps); norm(L*U-A)
% eps = 10^8; A = [1 1; eps -1]; [L,U] = my_lu2(eps); norm(L*U-A)

