#!/usr/bin/gnuplot

reset

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

set multiplot
set size 1.0,1.0

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

set format y "%2.2f"

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

set key top samplen 1.0

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  |  |     |
#     +-----+  +-----+
# distribucna funkcia pre T = 0K
set origin 0.0,0.0
set size 0.5,0.58
set lmargin 7
set rmargin 1

set grid
set xtics ( "10^3" 1e3, "10^4" 1e4, "10^5" 1e5, "10^6" 1e6, "10^7" 1e7, "10^8" 1e8)
set ytics ("10^{-4}" 1e-4, "10^{-3}" 1e-3, "10^{-2}" 1e-2, "10^{-1}" 1e-1)
set xlabel "N\n a)"
set ylabel "\n\n{/Symbol e}"

set logscale xy

plot [1e3:1e8][] "analytic-chyba.dat" title "MC" with points, \
	0.7 / sqrt(x) notitle


set nologscale xy


#
# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+
#     +-----+  +-----+
#     |     |  |  X  |
#     +-----+  +-----+
# presnost metody Monte Carlo
# zavislost chyby od poctu simulovanych ionov
set origin 0.5,0.0
set size 0.5,0.58
set lmargin 8
set rmargin 0

set nogrid
set xtics autofreq
set ytics 0.2
set xlabel "u\n b)"
set ylabel "{/Symbol j}(u, {/Symbol t})"

 set label 1 "{/Symbol t} = 0.3" at v_norm(18),0.85
#set label 2 "{/Symbol t} = 0.7" at v_norm(38),0.80
 set label 3 "{/Symbol t} = 1.1" at v_norm(62),0.65
#set label 4 "{/Symbol t} = 1.5" at v_norm(82),0.70
 set label 5 "{/Symbol t} = 1.9" at v_norm(105),0.45
#set label 6 "{/Symbol t} = 2.1" at v_norm(115),0.65
#set label 7 " {/Symbol t} = 2.5" at v_norm(136),0.60
#set label 8 "{/Symbol t} = 2.7" at v_norm(125),0.55
 set label 8 "{/Symbol t} = 2.7" at v_norm(118),0.25
#set label 9 "{/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