三重対角行列を LU 分解する関数 trilu() と、 その結果を元に後退代入して連立1次方程式を解く関数 trisol() を 納めた trid-lu.c, trid-lu.h (とそれらの ``添字が から 始まる'' バージョン trid-lu1.c, trid-lu1.h) を公開している (http://www.math.meiji.ac.jp/~mk/program/ から入手できる)。