#numeric stability

Illustrating the LU algorithm with pivots by example

[mathjax]

[Click here for a PDF of this post with nicer formatting]

Two previous examples of LU factorizations were given. I found one more to be the key to understanding how to implement this as a matlab algorithm, required for our problem set.

A matrix that contains both pivots and elementary matrix operations is

\begin{equation}\label{eqn:luAlgorithm:20} M= \begin{bmatrix} 0 & 0 & 2 & 1 \\ & 0…