summaryrefslogtreecommitdiff
path: root/figs/scaling_func.gplot
diff options
context:
space:
mode:
authorJaron Kent-Dobias <jaron@kent-dobias.com>2017-05-26 10:20:59 -0400
committerJaron Kent-Dobias <jaron@kent-dobias.com>2017-05-26 10:20:59 -0400
commit7b4f4ee4813f8afc3e9949a368fe708351d3b89c (patch)
tree33b6d66b24cf57c193d3dabc8b891ca382940c92 /figs/scaling_func.gplot
downloadpaper-7b4f4ee4813f8afc3e9949a368fe708351d3b89c.tar.gz
paper-7b4f4ee4813f8afc3e9949a368fe708351d3b89c.tar.bz2
paper-7b4f4ee4813f8afc3e9949a368fe708351d3b89c.zip
changes
Diffstat (limited to 'figs/scaling_func.gplot')
-rw-r--r--figs/scaling_func.gplot16
1 files changed, 16 insertions, 0 deletions
diff --git a/figs/scaling_func.gplot b/figs/scaling_func.gplot
new file mode 100644
index 0000000..e25b340
--- /dev/null
+++ b/figs/scaling_func.gplot
@@ -0,0 +1,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)
+
+