#!/usr/bin/gnuplot

reset

#set terminal postscript enhanced
#set output "distrib-beta-model.ps"
set terminal postscript eps enhanced 16
set output "distrib-beta-model.eps"

set line style 1 linetype 7 linewidth 1
set line style 2 linetype 1 linewidth 3

set multiplot
set size 1.0,1.0
set lmargin 8
set rmargin 1

set key samplen 0.7

set format y "%2.3f"

lambda = 3e-7
acc    = 1e10
Nt     = 121
Nx     = 100
f_max_x          = 5.000000
f_max_v          = 5.000000
V         = sqrt(lambda * acc)
fi_norm   =   Nx / (f_max_v * sqrt(2/pi))
v_norm(x) = x / V

# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |  X  |  |     |
#     +-----+  +-----+
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
set origin 0.0,0.5
set size 0.5,0.5
set key bottom
set ytics 0.2

set xlabel "{/Symbol t}\n a)"
set ylabel "u_d({/Symbol t})"


plot [][] \
	"distrib-N-5e5_T-300_alpha-0.1_beta-0.3_speed-0.dat" \
	using (v_norm($6)):(fi_norm * $7 ) title "MC" with l 1

# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
#     +-----+  +-----+
#     |  X  |  |     |
#     +-----+  +-----+
set origin 0.0,0.0
set size 0.5,0.5
set key bottom
set ytics 0.001
set xlabel "{/Symbol t}\n c)"
set ylabel "u_d({/Symbol t})"

plot [][] \
	"distrib-N-5e5_T-300_alpha-0.1_beta-0.3_speed-0.dat" \
	using (v_norm($6)):(fi_norm * $7 ) title "MC" with l 1


# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |     |  |  X  |
#     +-----+  +-----+
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
set origin 0.5,0.5
set size 0.5,0.5
set lmargin 9
set rmargin 0

plot [][] \
	"distrib-N-5e5_T-300_alpha-0.1_beta-0.3_speed-0.dat" \
	using (v_norm($6)):(fi_norm * $7 ) title "MC" with l 1

# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
#     +-----+  +-----+
#     |     |  |  X  |
#     +-----+  +-----+
# Entropia:
set origin 0.5,0.0
set size 0.5,0.5
set lmargin 9
set rmargin 0
set key bottom

plot [][] \
	"distrib-N-5e5_T-300_alpha-0.1_beta-0.3_speed-0.dat" \
	using (v_norm($6)):(fi_norm * $7 ) title "MC" with l 1

set nomultiplot