［Octave］LU分解を行う

lu関数は引数に与えた行列をLU分解する。以下、実行例。

> mx = [2 -3 -1; 1 / 2 1 -1; 0 1 -1];
> mx
mx =
   2.00000  -3.00000  -1.00000
   0.50000   1.00000  -1.00000
   0.00000   1.00000  -1.00000
> [mxl, mxu] = lu(mx)
mxl =
   1.00000   0.00000   0.00000
   0.25000   1.00000   0.00000
   0.00000   0.57143   1.00000
mxu =
   2.00000  -3.00000  -1.00000
   0.00000   1.75000  -0.75000
   0.00000   0.00000  -0.57143