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.