#!/usr/bin/ruby
require 'gsl'
include GSL
include Math

f459 = Function.new { |x|
  1.0/(5.0*x*x*x + 6.0)
}

exp_result = -8.994400695837000137E-02
exp_abserr = 1.185290176227023727E-06

result = f459.qawc([-1.0, 5.0], 0, [0.0, 1e-3])
p result
puts("exp_result: #{exp_result}")
puts("exp_abserr: #{exp_abserr}")

p Integration.qawc(f459, [-1.0, 5.0], 0)
