set terminal postscript eps enhanced color size 8.6cm, 6cm dl 3\
font "Helvetica, 16pt"
set output 'Vc_vs_nc.eps'

maj_f = 'Vc_nc.dat'
par_f = 'Vc_nc_par.dat'
set samples 10000
unset bars
# set key at 65,.8 reverse right
# set key reverse invert at 85,.2 Left
set key
set style line 1 lt 1 lc rgb "black" lw 2.5
set style line 2 lt 1 lc rgb "brown" lw 2.5
set style line 3 lt 1 lc rgb "orange" lw 2.5
set style line 4 lt 1 lc rgb "#7570B3" lw 2.5
set style line 5 lt 1 lc rgb "#6495ED" lw 2.5
set style line 6 lt 1 lc rgb "#034E7B" lw 2.5
set style line 7 lt 1 lc rgb "black" lw 2.5

set style line 11 lt 1 pt 6 ps 1.1 lc rgb "brown" lw 3
set style line 12 lt 1 pt 4 ps 1.1 lc rgb "orange" lw 3
set style line 13 lt 1 pt 10 ps 1.4 lc rgb "#7570B3" lw 3
set style line 14 lt 1 pt 12 ps 1.4 lc rgb "#6495ED" lw 3
set style line 15 lt 1 pt 14 ps 1.4 lc rgb "#034E7B" lw 3
set style line 16 lt 1 pt 16 ps 1.2 lc rgb "black" lw 1


Vc_maj(x) = 1 - 0.5807 / x + 0.3479 / x**2 - 0.06002 / x**3
Vc_par(x) = 1 - 0.1584 / x + 0.03987 / x**2 - 0.1743 / x**3
Vc_par_asym(x) = 1. - (1.-3.**(9./8.)/4.)/x
f_maj(x) = (x >= 1) ? (Vc_maj(x)) : 1/0
f_par(x) = (x >= 1) ? (Vc_par(x)) : 1/0
f_par_asym(x) = (x >= 1) ? (Vc_par_asym(x)) : 1/0

set ylabel 'V_c' offset 2,0
set xlabel 'n_c'
# set yrange [0.66:1.01]
# set yrange [:1.005]
set xrange [0:66]
set ytics .1
#set logscale x
#set logscale y
set key at graph 1., 0.22

set multiplot

# plot f_maj(x) ls 3 notitle, \
#     par_f u 1:($2) every 2 ls 15 notitle, \
#     f_par(x) ls 6 notitle, \
#     maj_f u 1:($2) ls 12 notitle

plot f_maj(x) ls 3 notitle, \
    par_f u 1:($2) every 2 ls 15 notitle, \
    f_par(x) ls 6 notitle, \
    maj_f u 1:($2) ls 12 notitle

unset yrange
set origin 0.29,0.19
set size 0.65,0.66
set xrange [20:66]
set ytics .01
set xtics 10
set key at graph 1.03, 0.35

# plot f_maj(x) ls 3 notitle, \
#     par_f u 1:($2) every 2 ls 15 title 'Parity', \
#     f_par(x) ls 6 notitle, \
#     maj_f u 1:($2) ls 12 title 'Majority'

plot f_maj(x) ls 3 notitle, \
    par_f u 1:($2) every 2 ls 15 title 'Parity', \
    f_par(x) ls 6 notitle, \
    maj_f u 1:($2) ls 12 title 'Majority'