blob: e25b34031103ce121161aec5120db634cddf9eb6 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
|
set terminal latex rotate
set ylabel "$\\chi B/At^{-\\gamma}$"
set xlabel "$h/Bt^{\\beta\\delta}$"
set key off
gamma = 0.57721566490153286060651209008240243104215933593992
Ei(m, x) = gamma + log(abs(x)) + sum[k = 1 : m] x**k / (k * k!)
Xi(m, A, B, x) = -A/(x * pi) * (1-B/x-(B/x)**2 * exp(B/x)* Ei(100, -B/x))
set yrange [0.001:0.6]
#set logscale xy
plot [0.001:5] Xi(200, -1, 1, x)
|
