File: [Platon] / doc / diplomova-praca-rajo / images / distrib_T-0.gpi (download)
Revision 1.2, Sun Jan 16 15:01:30 2005 UTC (19 years, 3 months ago) by rajo
Changes since 1.1: +2 -2
lines
Fix: gnuplot 4.0 has different syntax (incompatible with gnuplot 3.7.x)
Change lines from:
set linestyle 1 linetype 7 linewidth 1
to:
set line style 1 linetype 7 linewidth 1
|
#!/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
Platon Group <platon@platon.org> http://platon.org/
|