code: public double[][] getU() {
      double[][] U = new double[n][n];
      for(int i = 0; i < n; i++) {
        System.arraycopy(LU[i], i, U[i], i, n - i);
      }
      return U;
    }
docs: Return upper triangular factor
  
  @return U