LU 分解

由此，矩陣將被分解為上部三角矩陣和下三角矩陣。通常，它將用於提高 Gauß消除的效能和穩定性（如果通過排列完成）。

然而，這種方法通常不會或不好地工作，因為它不穩定。例如

A = [8 1 6
    3 5 7
    4 9 2];
[L,U] = lu(A);

新增置換矩陣就足以使 PA = LU：

[L,U,P]=lu(A);

在下文中，我們將繪製 lu 的執行時間，該執行時間取決於矩陣元素的平方根。 http://i.stack.imgur.com/cnhNK.jpg