set term png
set output "pin_gumbel_pdf.png"

set title "gumbel pdf(a, b, x)=1/b*exp(-(x-a)/b)*exp(-exp(-(x-a)/b))"
set xlabel "x"
set xlabel "pdf(x)"
set xrange [-5:20]

f(a, b, x)=1/b*exp(-(x-a)/b)*exp(-exp(-(x-a)/b))
plot f(0.5, 2.0, x) w l lw 2 title "a=0.5, b=2.0"
replot f(1.0, 2.0, x) w l lw 2 title "a=1.0, b=2.0"
replot f(1.5, 3.0, x) w l lw 2 title "a=1.5, b=3.0"
replot f(3.0, 4.0, x) w l lw 2 title "a=3.0, b=4.0"
replot f(3.0, 6.0, x) w l lw 2 title "a=3.0, b=6.0"
set term png
set output "pin_gumbel_pdf.png"
replot

gnupot plot_gumble_pdf.plt