function [L, U] = mylu(A)
%MYLU   Perform LU factorisation without pivoting.
% Nick Hale, Aug 2019

[L, U] = lu(sparse(A), 0); % Sparse LU can be forced not to pivot
L = full(L); U = full(U);  % Convert output to full

end
