LU factorization is a method for decomposing an N × N {\displaystyle N\times N} matrix A {\displaystyle \mathbf {A} } into the product of a lower triangular matrix L {\displaystyle \mathbf {L} } and an upper triangular matrix U {\displaystyle \mathbf {U} } ,
