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


syntax highlighted by Code2HTML, v. 0.9.1