Octave で知っとくと便利な関数

＜N次元の疎な単位行列＞
speye(N)

＜N×Nの疎な3項帯行列＞
3つのN次元の列ベクトル La, Lb, Lc に要素を入れておき，
L = spdiags([Lc Lb La], -1:1, N, N)'

＜疎行列 S を通常の行列に変換＞
A=full(S)

＜通常の行列Aを疎行列に変換＞
S=sparse(A)

＜AとBをブロック対角要素に持つ行列＞
blkdiag(A, B)