# File lib/backports/1.9.2/stdlib/matrix/lup_decomposition.rb, line 35
    def u
      Matrix.build(@col_size, @col_size) do |i, j|
        if (i <= j)
          @lu[i][j]
        else
          0
        end
      end
    end