http://www2.ee.knct.ac.jp/el/E4/H15-E406/lu1.html/* LU分解による連立一次方程式の解法を用いたプログラム例 */ /* Lの対角要素は１とする。 */ #include <iostream.h> #include <iomanip.h> #include <math.h> const int N=3; /* 変数の数 */ const double Eps=1e-10; /* Pivotの最低値を規定する小さい数 */ /* ユーザー定義の型Matrixの定義 */ typedef double Matrix[N+1][N+1]; /* ユーザー定義の型Vectorの定義 */ typedef double Vector[N+1]; /* LU分解を行う関数の定義 */ int lu_decomposition( Matrix a,int n ) { int i,j,t; for ( t=1; t <=
