Initial commit

figs/energy.fig/energy.gnuplot

@@ -0,0 +1,41 @@
+set title ""
+set ylabel "$\\displaystyle\\frac{e}{\\rho}$" norotate offset -1,0
+set xlabel "$\\rho$"
+
+set xtics 1e-6, 100, 100
+set xtics add ("$10^{-6}$" 0.000001, "$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0, "$10^2$" 100)
+unset mxtics
+set xrange [0.000001:100]
+
+set ytics 8, 1, 13
+set mytics 2
+set yrange [7.5:13]
+
+
+# default output canvas size: 12.5cm x 8.75cm
+set term lua tikz size 8,6 standalone
+
+set key bottom right box linetype rgbcolor"#999999" width 0.3 height 0.3 spacing 1.3
+
+
+# set linestyle
+set style line 1 linetype rgbcolor "#4169E1" linewidth 2 dashtype "."
+set style line 2 linetype rgbcolor "#DC143C" linewidth 2
+set style line 3 linetype rgbcolor "#32CD32" linewidth 2 dashtype "-"
+set style line 4 linetype rgbcolor "#4B0082" linewidth 2
+set style line 5 linetype rgbcolor "#DAA520" linewidth 2
+set style line 6 linetype rgbcolor "#999999" linewidth 1.5 dashtype (1,5)
+
+set pointsize 1
+
+set logscale x
+
+a=1.2543564105910647
+
+plot \
+ "simpleq.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 1 title "simple" ,\
+ "medeq.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 3 title "med",\
+ "bigeq.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 5 title "big",\
+ "2020-10-15+energy.dat" using 1:($2/$1) with points ls 2 title "QMC" ,\
+ 2*pi*a*(1+128./(15*sqrt(pi))*sqrt(x*a**3)) ls 6 title "LHY"
+