/*
 * <日付を入力>
 * 例題6.5
 * Gauss-Jordan法メソッドを使う実行クラス
 */


public class LinearEquationCalculation2{

  // メインメソッド
  public static void main(String [] args) {

    int degree = 3;  // ベクトル，行列の次数
    double [][] coefficient_matrix = { {  3.0, 10.0, -4.0}, 
                                       {  1.0,  2.0,  3.0}, 
                                       { -2.0, -4.0,  1.0} };  // 各要素の値を代入
    double [] unknown_vector = {0.0, 0.0, 0.0};  // 未知ベクトル．初期値は0
    double [] known_vector = {11.0, 14.0, -7.0 };  // 定数ベクトルの値を代入

    // 連立方程式の表示
    System.out.println(degree + "元連立一次方程式\n");
    for(int i=0; i<degree; i++) {
      for (int j=0; j<degree; j++) {
        System.out.print(" +(" + coefficient_matrix[i][j] + ")･x[" + j + "]");
      }
      System.out.println(" = " + known_vector[i]);
    }

    // Gauss-Jordan法のメソッドを持つクラスのオブジェクトを生成
    LinearEquationSolver2 les = new LinearEquationSolver2(degree, coefficient_matrix, unknown_vector, known_vector);

    // Gauss-Jordan法メソッドの呼び出し．戻り値が解ベクトルに代入される．
    unknown_vector = les.calculateGJE();
    
    // 解の表示
    System.out.println("解: ");
    for (int i=0; i<degree; i++) {
      System.out.println("x[" + i + "] = " + unknown_vector[i]);
    }

  }  // main

}  // class LinearEquationCalculation