function [P,L,A] = LU(A)
    n = length(A);
    P = eye(n);
    L = eye(n);
    % Your code here
end