LU

bool matrix::LU(
  matrix&  L,    // 下三角行列
  matrix&  U      // 上三角行列
  );

  matrix matrix_a={{1,2,3,4},
                   {5,2,6,7},
                   {8,9,3,10},
                   {11,12,14,4}};
  matrix matrix_l,matrix_u;
//--- LU分解
  matrix_a.LU(matrix_l,matrix_u);
  Print("matrix_l\n",matrix_l);
  Print("matrix_u\n",matrix_u);
//--- A = L * Uかどうかを確認する
  Print("check\n",matrix_l.MatMul(matrix_u));
 
 
  /*
  matrix_l
  [[1,0,0,0]
   [5,1,0,0]
   [8,0.875,1,0]
   [11,1.25,0.5904761904761905,1]]
  matrix_u
  [[1,2,3,4]
   [0,-8,-9,-13]
   [0,0,-13.125,-10.625]
   [0,0,0,-17.47619047619047]]
  check
  [[1,2,3,4]
   [5,2,6,7]
   [8,9,3,10]
   [11,12,14,4]]
  */