A method for Matrix Factorization and solving for

We then set , giving us . Since is lower triangular, and U is upper triangular, we an easily solve the two linear systems by substitution.

Using the matrix inverse is inefficient and numerically unstable.