| 1 |
We wish to decompose the matrix A into the product LU, where L is a lower triangular matrix with 1's on the main diagonal, and U is an upper triangular matrix.
|
| 2 |
We assume A is a full M by M matrix.
|
| 3 |
In general, pivoting is necessary to ensure numerical stability.
|
| 4 |
LU decomposition is often used in the solution of systems of linear equations, Ax = b. The equations can be written as two triangular systems,
|
| 5 |
The first equation is solved for y by forward reduction, and the solution x is then obtained from the second equation by back substitution.
|
