# see also codes/EN/CPS
#                   QPS
#               and NORMAL
#
# for files that set standard options for display

	set term X11
# bhattychariya parameter for BSC
b(f) = 2.0 * sqrt(f*(1-f))
# error probability associated with A(w) words of weight w is
#  A(w) b(f)**w
Ab(A,w,f) = A*b(f)**w
# more exact expression
Abe(A,w,f) = A* (w!/((w/2)!)**2) * ( sqrt(f*(1-f)) )**w
Abeo(A,w,f) = A* (w!/( ((w-1)/2)! * ((w+1)/2)! ) ) * (  f**((w+1)/2) * (1-f)**((w-1)/2) ) 

H(x) = ( x*log(1/x) + (1-x)*log(1/(1-x)) ) / log(2)
gau(x,m,s) = (1.0/(s*sqrt(2.0*3.14159)) )*exp(-(((x-m)/s)**2)*0.5)
set size 0.8,0.8
set origin 0.1,0.1

# This overrides the normal margin calculation.
# Motivation: means that axes will always come in a known location
# even if labels on tics get large
       set bmargin 0
       set lmargin 0 
       set rmargin 0 
       set tmargin 0 
 
# this makes contours all the same colour
set noclabel

# this makes a thin crisp border
set linestyle 8 lt 0 lw 0.3
set linestyle 1 lt 1 lw 0.3
set linestyle 10 lt 1 lw 2
set linestyle 100 lt 1 lw 5
set linestyle 2 lt 2 lw 0.3
set linestyle 3 lt 3 lw 0.3
set linestyle 4 lt 4 lw 0.3
set linestyle 5 lt 5 lw 0.3
set linestyle 6 lt 6 lw 0.3
set linestyle 7 lt 7 lw 0.3
set linestyle 9 lt 9 lw 0.3
set ytics nomirror
set xtics nomirror
set border 3 ls 8

myerf(x) = (1-erf( x/sqrt(2.0) ) )/2.0