# generic commands for twographbox output size load "twographbox.in" set data style lines set xlabel "Parameter {/Symbol q}" set ylabel "{/Times-Italic P}({/Symbol Q} = {/Symbol q})" set xrange [0:1] set nokey set label "[3,1]" at 0.25,0.7 set label "[6,2]" at 0.5,1.7 set label "[30,10]" at 0.57,5.0 a1=3 b1=1 alpha1 = 3 a2=6 b2=2 alpha2 = 42 a3=30 b3=10 alpha3 = 635.0e7 set output "beta-distributions2.eps" plot alpha1*(x**(a1-1))*((1-x)**(b1-1)) lw 2, alpha2*(x**(a2-1))*((1-x)**(b2-1)) lw 2, alpha3*(x**(a3-1))*((1-x)**(b3-1)) lw 2