1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
#
# Synchrotron function F in units of critical frequency v_c
#
save_encoding = GPVAL_ENCODING
set encoding utf8
set title "Synchrotron function F(x)"
set key rmargin reverse Left top title "\n\n" spacing 1.5
set log
set xrange [0.005 : 12.]
set yrange [0.005 : 1.9]
set xtics format "%g{/:Italic v_c}"
set xtics add ("{/:Italic v_c}" 1)
set tics nomirror
F(x) = SynchrotronF(x)
Power(x) = log(10.)*x*F(x)
plot F(x) lt 3 lw 2 title "F(x) = {/:Italic x{/*2 ∫@_{/*0.5 x}^{/*.5 ∞}}K_{5/3}(ν) dν}", \
Power(x) lt 4 lw 2 title "Power per unit log(ν)"
pause -1 "<cr> to continue"
reset
set encoding save_encoding
|
