Rank 1 Update for LU factors

% This function updates the LU factorization (PA = LU) of a Matrix A following a rank % one update of the matrix A = A + alpha*y*z'

% This function updates the LU factorization (PA = LU) of a Matrix A following a rank

% one update of the matrix A = A + alpha*y*z'

% L = Lower unit matrix

% U = Upper triangular matrix

% P = Permutation matrix

% alpha, yy(Column vector) and zy(Column vector) are the rank one update

% matrix alpha*y*z'