#!/usr/bin/gnuplot

reset

#set terminal postscript enhanced
#set output "slide-zavislost-vdrift-E-x2-model.ps"
set terminal postscript eps enhanced 24
set output "slide-zavislost-vdrift-E-x2-model.eps"

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

set lmargin 9
set rmargin 2


set mytics 2
set mxtics 2

set format y "%2.2f"

set grid

# ------------------------------------------------------------------------
#     +-----+  +-----+
#     |     |  |  X  |
#     +-----+  +-----+
#     +-----+  +-----+
#     |     |  |     |
#     +-----+  +-----+

set pointsize 2
set key top
#set ytics 0.005
set xlabel "{/Symbol a}"
set ylabel ""
set nolabel
set label 1 "v_{drift} [m.s^{-1}]" at graph -0.2,graph 0.5 center rotate

set logscale xy
set key bottom
set xtics ("10^{-4}" 1e-4, "10^{-3}" 1e-3, "10^{-2}" 1e-2, "10^{-1}" 1e-1, "1" 1, "10^1" 10, "10^2" 1e2, "10^3" 1e3)
set ytics (	"10^{-3}" 1e-3,	"10^{-2}" 1e-2,	"10^{-1}" 1e-1,	\
			"10^{0 }" 1,	"10^{1 }" 10,	"10^{2 }" 1e2,	\
			"10^{3 }" 1e3,	"10^{4 }" 1e4,	"10^{5 }" 1e5	)

k = 891.053540967816
#a               = 1.194            +/- 0.007405
a = 1.194

plot [1e-3:1e3] [1e-1:1e5] \
	"zavislost-vdrift-E_T-300-x2-model.dat" title "MC" with points 1, \
	k * sqrt(x) title "~E^{0.50}" with lines 2, \
	k * x		title "~E^{1.00}" with lines 2, \
	k * x**a title "~E^{1.19}" with lines 1


set nologscale xy
set format y "%2.2f"