function [L, U] = nopivotlu(A)

[L, U, P] = lu(sparse(A), 0);
L = full(L);
U = full(U);

end