#!/usr/bin/env ruby
require("gsl")
x = Vector.linspace(0, 20, 100)
y = Sf::dawson(x)
y.graph(x, "-T X -C -g 3 -L 'Dawson integral'")