WOLFRAM

Wolfram言語 ＆ システムドキュメントセンター

アップグレード情報：

LinearAlgebra`Tridiagonal`

三重対角行列を生成するためには，Mathematica の組込み関数SparseArrayを使う．
三重対角行列のより高速なソルバとして，TridiagonalSolveの代りに Mathematica カーネルの組込みスパースソルバLinearSolveが使われるようになった．

Tridiagonalを置換する場合のヒント

バージョン5.2では，このようにパッケージをロードする：

Version 5.2 << LinearAlgebra`Tridiagonal`;

次で，行列の非零の対角要素を与えるリストを定義する：

Version 5.2 {a, b, c} = {{7, 1, 11}, {4, 8, 2, 12}, {5, 9, 3}};
Wolfram Language code: {a, b, c} = {{7, 1, 11}, {4, 8, 2, 12}, {5, 9, 3}};

リストから構築された行列mである：

Version 5.2 m = Table[ Switch[j - i, -1, a[[j]], 0, b[[j]], 1, c[[j - 1]], _, 0], {i, 4}, {j, 4}];
Wolfram Language code: m = SparseArray[{Band[{1, 2}] -> c, Band[{1, 1}] -> b, Band[{2, 1}] -> a}, 4]

MatrixFormを使って行列を表示する：

Version 5.2 MatrixForm[m];
Wolfram Language code: MatrixForm[Normal[m]]

TridiagonalSolveの代りに，Mathematica カーネルの組込みスパースソルバLinearSolveを使う：

Version 5.2 TridiagonalSolve[a, b, c, {2, 3, 4, 5}];
Wolfram Language code: LinearSolve[m, {2, 3, 4, 5}]