6.2 Die LU-Zerlegung

Sei

$A_{n\times n}$ - vollbesetzt, zeilenweise gespeichert,
$L$ - untere Dreiecksmatrix mit Normierung $\ell_{i,i}=1$, spaltenweise gespeichert,
$U$ - obere Dreiecksmatrix, zeilenweise gespeichert.



und es soll gelten $\sum\limits_{j=1}^n \ell_{i,j} \cdot u_{j,k} \,=\, a_{i,k}$ .