/*
 * 日付を入力
 * 演習6.9
 * Gaussの消去法を使う実行クラス
 */


public class LinearEquationCalculation{

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

    try {
      int degree = 4;  // ベクトル，行列の次数
      double [][] coefficient_matrix = { {  1.0, -2.0,  1.0, -3.0}, 
                                         {  2.0,  1.0,  3.0,  2.0}, 
                                         {  1.0,  3.0, -2.0,  1.0},
                                         {  3.0,  2.0,  1.0,  1.0} };  // 各要素の値を代入
      double [] unknown_vector = {0.0, 0.0, 0.0, 0.0};  // 未知ベクトル．初期値は0
      double [] known_vector = {-12.0, 21.0, 5.0, 14.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]);
      }

      // ガウスの消去法を行うクラスのオブジェクトを生成
      LinearEquationSolver les = new LinearEquationSolver(degree, coefficient_matrix, unknown_vector, known_vector);

      // ガウスの消去法メソッドの呼び出し．戻り値が解ベクトルに代入される．
      unknown_vector = les.calculateGE();
    
      // 解の表示
      System.out.println("解: ");
      for (int i=0; i<degree; i++) {
        System.out.println("x[" + i + "] = " + unknown_vector[i]);
      }
    } catch (Exception e) {
      System.out.println(e);
    }
  }  // main

}  // class LinearEquationCalculation