LU decomposition for a matrix A allows us to solve the linear system
Ax=b using a combination of forward and back substitution with.
LU factorization with full pivoting involves both row and column permutations, PAQ=LU, where L and U, and P are defined as before, and Q is a permutation matrix to the columns of A
