#!/usr/bin/env ruby
require("gsl")

# Singular at x = 0
f = Function.alloc { |x| 1.0/Math::sqrt(x) }

puts("QAG")
p f.qag(0, 1)
p f.qag([0, 1])

p Integration.qag(f, 0, 1)
p Integration.qag(f, [0, 1])
syntax highlighted by Code2HTML, v. 0.9.1