#!/usr/bin/env ruby
require("gsl")
n = 50
x = Vector.linspace(0.01, 14, n)
y05 = Sf::gamma_inc_P(0.5, x)
y1 = Sf::gamma_inc_P(1, x)
y3 = Sf::gamma_inc_P(3, x)
y10 = Sf::gamma_inc_P(10, x)
Vector.graph(x, y05, y1, y3, y10, "-T X -C -g 3 -y 0 1.1 -L 'Incomplete gamma functions'")