# Submitted by M. Edward (Ed) Borasky
require 'hilbert' # also brings in mathn, matrix, rational and complex

def run_hilbert(dimension)
  m = hilbert(dimension)
  print "Hilbert matrix of dimension #{dimension} times its inverse = identity? "
  k = m * m.inv
  print "#{k==Matrix.I(dimension)}\n"
  m = nil # for the garbage collector
  k = nil
end

Bench.run [10, 20, 30, 40, 50, 60] do |n|
  run_hilbert(n)
end