ports//math/ruby-gsl/work/rb-gsl-1.9.2/samples/pdf/rayleigh.rb

#!/usr/bin/env ruby
require("gsl")
x = Vector.linspace(0, 5, 100)
y1 = Ran::rayleigh_pdf(x, 1)
y2 = Ran::rayleigh_pdf(x, 2)
Vector.graph(x, y1, y2, "-T X -g 3 -C -X x -Y 'p(x)' -x 0 5 -L 'Rayleigh distribution, a = 1' --toggle-rotate-y-label")

syntax highlighted by Code2HTML, v. 0.9.1