function [L,U] = LUDecomp(A)
%LUDECOMP LU factorization of A.
%   [L,U]= LUDECOMP(A) returns the LU factorization of A.
    [L,U] = lu(A);
