set terminal postscript eps color enhanced solid font "Helvetica,22" size 12cm,10cm
set colorsequence classic
set output 'filterfunc.eps'

i = {0.0,1.0}

set label 1 at 300, 0.57
set label 1 "{/Symbol t}=1ms" tc lt 7
set key samplen 2

set xlabel '{/Symbol w}{/Symbol t} (rad.)'
set ylabel 'Filter function g_N({/Symbol w}, {/Symbol t})'
set xrange [0:100]
set yrange [0:0.5]
set xtics ("5{/Symbol p}" 15.70796, "15{/Symbol p}" 47.12389, "25{/Symbol p}" 78.53982)
set ytics 0.2
set samples 1000
fac_mod = 0
f_pdd1(x) = abs(1+exp(i*x)-2*exp(i*0.5*x)*cos(0.5*x*fac_mod))**2/x**2
f_pdd5(x) = abs(1+exp(i*x)-2*(exp(i*0.167*x)-exp(i*0.333*x)+exp(i*0.5*x)-exp(i*0.667*x)+exp(i*0.833*x))*cos(0.5*x*fac_mod))**2/x**2
f_pdd9(x) = abs(1+exp(i*x)-2*(exp(i*0.1*x)-exp(i*0.2*x)+exp(i*0.3*x)-exp(i*0.4*x)+exp(i*0.5*x)-exp(i*0.6*x)+exp(i*0.7*x)-exp(i*0.8*x)+exp(i*0.9*x))*cos(0.5*x*fac_mod))**2/x**2
f_udd5(x) = abs(1+exp(i*x)-2*(exp(i*0.067*x)-exp(i*0.25*x)+exp(i*0.5*x)-exp(i*0.75*x)+exp(i*0.933*x))*cos(0.5*x*fac_mod))**2/x**2
f_udd9(x) = abs(1+exp(i*x)-2*(exp(i*0.0245*x)-exp(i*0.0955*x)+exp(i*0.2061*x)-exp(i*0.3455*x)+exp(i*0.5*x)-exp(i*0.6545*x)+exp(i*0.7939*x)-exp(i*0.9045*x)+exp(i*0.9755*x))*cos(0.5*x*fac_mod))**2/x**2
f_pdd15(x) = abs(1+exp(i*x)-2*(exp(i*1/16*x)-exp(i*2/16*x)+exp(i*3/16*x)-exp(i*4/16*x)+exp(i*5/16*x)-exp(i*6/16*x)+exp(i*7/16*x)-exp(i*8/16*x)+exp(i*9/16*x)-exp(i*10/16*x)+exp(i*11/16*x)-exp(i*12/16*x)+exp(i*13/16*x)-exp(i*14/16*x)+exp(i*15/16*x))*cos(0.5*x*fac_mod))**2/x**2

plot f_pdd5(x) lw 2 dt 1 lc 1 t '5{/Symbol p} PDD',\
f_pdd9(x) lw 2 dt 1 lc 3 t '9{/Symbol p} PDD'


pause -1