ガウスの消去法をJavaで解かせるためのライブラリです。
下記リンク内にEclipse用のプロジェクトがあるので、参考にして下さい。
https://docs.google.com/a/sangi.jp/file/d/0B5hI-TWePsTWa3RDSjY4NHlHM1U/edit?usp=sharing
以下ソースコードです。
package Gauss_lib;
//1.Gauss_libクラスのオブジェクト生成
// Gauss_lib gauss_lib = new Gauss_lib(n);//次元数nをコンストラクタに設定
//2.係数行列の引き渡し
//double T[][] = {{2,4,2,8},{4,10,3,17},{3,7,1,11}};
//gauss_lib.setT(T);//係数行列の引き渡し
//3．計算用エンジン呼び出し
//gauss_lib.gauss();//ガウスエンジン呼び出し
//4.解を得る
//x = gauss_lib.getx();//解を取り出す
public class Gauss_lib {
//フィールド
int n;
//配列の初期化
double [] x; //解
double [][] T ;//係数行列格納用配列
//コンストラクタ
public Gauss_lib(int n){
this.n = n;
x = new double[n];//解を代入する配列
T = new double[n][n + 1];//係数行列
}
public int getN(){//数値を出力
return n;
}
public double[] getx(){//解の受け渡し
return x;
}
public void setT(double[][] n){//係数行列の受け渡し
T = n;
}
public void gauss() {
int i,j,k,l,pivot;
double p,q,m;
double b[] = new double [n + 1];
for(i=0;i<n;i++)
{
m=0;
pivot=i;
for(l=i;l<n;l++)
{
if(Math.abs(T[l][i])>m) //i列の中で一番値が大きい行を選ぶ
{
m=Math.abs(T[l][i]);
pivot=l;
}
}
if(pivot!=i) //pivotがiと違えば、行の入れ替え
{
for(j=0;j<n+1;j++)
{
b[j]=T[i][j];
T[i][j]=T[pivot][j];
T[pivot][j]=b[j];
}
}
}
for(k=0;k<n;k++)
{
p=T[k][k]; //対角要素を保存
T[k][k]=1; //対角要素は１になることがわかっているから
for(j=k+1;j<n+1;j++)
{
T[k][j]/=p;
}
for(i=k+1;i<n;i++)
{
q=T[i][k];
for(j=k+1;j<n+1;j++)
{
T[i][j]-=q*T[k][j];
}
T[i][k]=0; //0となることがわかっているところ
}
}
/*解の計算*/
for(i=n-1;i>=0;i--)
{
x[i]=T[i][n];
for(j=n-1;j>i;j--)
{
x[i]-=T[i][j]*x[j];
}
}
}
}
