#!/usr/bin/gnuplot -persist
#
# G N U P L O T
# Version 4.2 patchlevel 2
# last modified 31 Aug 2007
# System: Linux 2.6.24.3-34.fc8
set terminal svg
set output "Gumbel-Density.svg"
set key 15, 0.16, 0
set xrange [ -5.00000 : 20.0000 ]
set format y "%4.2f"
z(x,mu,beta) = exp(-(x-mu)/beta)
f(x,mu,beta) = z(x,mu,beta) * exp( -1 * z(x,mu,beta) ) / beta
F(x,mu,beta) = exp( -1*exp( -1*(x-mu)/beta ) )
plot \
f(x,0.5,2) ti sprintf("f(x, µ=%3.1f, β=%3.1f)", 0.5, 2) lt 1, \
f(x,1,2)   ti sprintf("f(x, µ=%3.1f, β=%3.1f)", 1, 2) lt 2, \
f(x,1.5,3) ti sprintf("f(x, µ=%3.1f, β=%3.1f)", 1.5,3) lt 3, \
f(x,3,4)   ti sprintf("f(x, µ=%3.1f, β=%3.1f)", 3, 4) lt 4