reset
resetsession
out=2
outputpath="./"
outputfilename="betadistribution"
if(out==1){
setterminalqtsize1000,750font",14"
}
if(out==2){
setterminalpdfcairosize15cm,12cmfont",20"linewidth2
setoutputoutputpath.outputfilename.".pdf"
}
if(out==3){
setterminalpngcairosize800,600font",14"
setoutputoutputpath.outputfilename.".png"
}
dbeta(x,a,b)=exp((a-1)*log(x)+(b-1)*log(1-x)-lgamma(a)-lgamma(b)+lgamma(a+b))
pbeta(x,a,b)=ibeta(a,b,x)
f(x,y,z)=dbeta(x,y,z)
setsamples500

setxrange[0:1]
setxtics0.2
setyrange[0:2.5]
setytics0.5
#set xlabel "x"
setylabel"PDF"
setkeyat.78,2.5
setgridxticsyticsmyticslw1lcrgb"#888888"dashtype3lt1
setborder3backlw2lcrgb"#222222"
setticsinnomirrorscale0.75
lww=2
setls1lwlwwlcrgb"#e41a1c"
setls2lwlwwlcrgb"#377eb8"
setls3lwlwwlcrgb"#4daf4a"
setls4lwlwwlcrgb"#984ea3"
setls5lwlwwlcrgb"#ff7f00"
plot\
f(x,0.5,0.5)wlls1title"{/Symbol a} = {/Symbol b} = 0.5",\
f(x,5,1)wlls2title"{/Symbol a} = 5, {/Symbol b} = 1",\
f(x,1,3)wlls3title"{/Symbol a} = 1, {/Symbol b} = 3",\
f(x,2,2)wlls4title"{/Symbol a} = 2, {/Symbol b} = 2",\
f(x,2,5)wlls5title"{/Symbol a} = 2, {/Symbol b} = 5"
if (out>1){
setoutput
}
if (out==1){
pause-1

from matplotlib.pyplot import *
from numpy import linspace
from scipy.stats import beta
x = linspace(0,1,75)
fig = figure()
ax = fig.add_subplot(111)
ax.plot(x,beta.pdf(x,0.5,0.5),label=r"$\alpha=\beta=0.5$")
ax.plot(x,beta.pdf(x,5,1),label=r"$\alpha=5, \beta=1$")
ax.plot(x,beta.pdf(x,1,3),label=r"$\alpha=1, \beta=3$")
ax.plot(x,beta.pdf(x,2,2),label=r"$\alpha=2, \beta=2$")
ax.plot(x,beta.pdf(x,2,5),label=r"$\alpha=2, \beta=5$")
ax.grid(True)
ax.minorticks_on()
ax.legend(loc=9)
setp(ax.get_legend().get_texts(),fontsize='small')
ax.set_ylim(0,2.6)
ax.set_xlabel("x")
ax.set_ylabel("PDF")
fig.savefig("Beta_distribution_pdf.svg",bbox_inches="tight",\
	pad_inches=.15)