1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
|
set ylabel norotate "$g(r_*)$"
set xlabel "$\\rho$"
set xtics 1e-4, 100, 10
set xtics add ("$10^{-4}$" 0.0001, "$10^{-3}$" 0.001, "$10^{-2}$" 0.01, "$10^{-1}$" 0.1, "$1$" 1, "$10$" 10)
unset mxtics
set xrange [0.0001:10]
#set ytics 0, 1e-5, 2e-5
#set ytics add ("$10^{-5}$" 1e-5)
#set mytics
set yrange [0.999:]
# default output canvas size: 12.5cm x 8.75cm
set term lua tikz size 8,6 standalone
set key top right box linetype rgbcolor"#999999" width 2.5 height 0.3 spacing 1.3
set style line 1 linetype rgbcolor "#4169E1" linewidth 2 dashtype "."
set style line 2 linetype rgbcolor "#DC143C" linewidth 2
set style line 3 linetype rgbcolor "#32CD32" linewidth 2
set style line 4 linetype rgbcolor "#4B0082" linewidth 2
set style line 5 linetype rgbcolor "#DAA520" linewidth 2
set style line 6 linetype rgbcolor "#000000" linewidth 1
se logscale x
plot \
1 ls 6 notitle,\
(x<=1.46e-3 ? 1 : sqrt(-1)) with lines ls 1 notitle,\
"medeq.dat" using 1:($3/$1**2>=1 ? $3/$1**2 : 1) with lines ls 1 title "med",\
(x<=7.23e-4 ? 1 : sqrt(-1)) with lines ls 2 notitle,\
"bigeq.dat" using 1:($3/$1**2>=1 ? $3/$1**2 : 1) with lines ls 2 title "big"
|
