#!/usr/bin/gnuplot

reset

#set terminal postscript enhanced
#set output "distrib_T-0.ps"
set terminal postscript eps enhanced 24
set output "slide-distrib_T-0.eps"

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

set xtics -4,0.5
set ytics 0,0.2
set mxtics 2
set mytics 4

set xlabel "u"
#set ylabel "{/Symbol t}"
set ylabel "{/Symbol j}(u, {/Symbol t})"

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

set label "{/Symbol t} = 0.3" at v_norm(18),0.85
#set label "{/Symbol t} = 0.7" at v_norm(38),0.80
set label "{/Symbol t} = 1.1" at v_norm(62),0.75
#set label "{/Symbol t} = 1.5" at v_norm(82),0.70
set label "{/Symbol t} = 1.9" at v_norm(105),0.65
#set label "{/Symbol t} = 2.1" at v_norm(115),0.65
#set label "{/Symbol t} = 2.5" at v_norm(136),0.60
set label "{/Symbol t} = 2.7" at v_norm(125),0.55

#set label "{/Symbol j}(u, {/Symbol \245)=\326(2/p)} e^{-0.5 u^2}" at v_norm(165),0.3
#set arrow 1 from 3,0.27 to 2,0.108

plot [0:3][0:0.9] \
	"distrib-N-1e6_lambda-3e-7_a-1e10_T-0.dat" \
	using (v_norm($6)):(fi_norm * $7 ) title "MC" with l 1, \
	sqrt(2/pi) * exp (- 0.5 * x * x ) title "{/Symbol j}(u, {/Symbol \245)"
	#"" using (-v_norm($6)):8 notitle with lines 1, \

# vim600: fdm=marker fdc=3