program='gaussian_uniform'
ps=2
x=findgen(10000)/1000.-5
y=1./sqrt(2*!pi)*exp(-0.5*x^2)

psdirect,program,ps

plot,x,y,xtitle='x/'+textoidl('\sigma'),ytitle='f',yr=[0,1]

L_uni=sqrt(12.)
ind=where(abs(x) le l_uni/2)
oplot,x(ind),x(ind)*0+1./L_uni,lines=2
oplot,-L_uni/2*[1,1],[0,1./L_uni],lines=2
oplot, L_uni/2*[1,1],[0,1./L_uni],lines=2

label_data,0.1,.9,['Gaussian','Uniform distribution with the same variance'],lines=[0.1,2]

psdirect,program,ps,/stop

end