Initial commitv0.0

21 files changed, 828 insertions, 0 deletions

diff --git a/figs/condensate.fig/Makefile b/figs/condensate.fig/Makefile
new file mode 100644
index 0000000..626032e
--- /dev/null
+++ b/figs/condensate.fig/Makefile
@@ -0,0 +1,43 @@
+PROJECTNAME=condensate16 condensate05
+
+SIMPLEQ=simpleq
+
+DATS=simpleq16.dat mueq16.dat bigeq16.dat simpleq0.5.dat mueq0.5.dat bigeq0.5.dat
+PDFS=$(addsuffix .pdf, $(PROJECTNAME))
+TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME))
+
+all: $(PDFS)
+
+$(PDFS): $(DATS)
+	gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@)
+	pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@)
+
+simpleq16.dat:
+	julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16" -M simpleq condensate_fraction_rho > $@
+mueq16.dat:
+	julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16" -M mueq condensate_fraction_rho > $@
+bigeq16.dat:
+	julia -p 8 $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;v_a=16" -M medeq condensate_fraction_rho > $@
+simpleq0.5.dat:
+	julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=0.5" -M simpleq condensate_fraction_rho > $@
+mueq0.5.dat:
+	julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=0.5" -M mueq condensate_fraction_rho > $@
+bigeq0.5.dat:
+	julia -p 8 $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;v_a=0.5" -M medeq condensate_fraction_rho > $@
+
+
+install: $(PDFS)
+	cp $^ $(INSTALLDIR)/
+
+clean-aux:
+	rm -f $(addsuffix .tikz.tex, $(PROJECTNAME))
+	rm -f $(addsuffix .aux, $(PROJECTNAME))
+	rm -f $(addsuffix .log, $(PROJECTNAME))
+
+clean-dat:
+	rm -f $(DATS)
+
+clean-tex:
+	rm -f $(PDFS)
+
+clean: clean-aux clean-tex
diff --git a/figs/condensate.fig/condensate05.gnuplot b/figs/condensate.fig/condensate05.gnuplot
new file mode 100644
index 0000000..66a642e
--- /dev/null
+++ b/figs/condensate.fig/condensate05.gnuplot
@@ -0,0 +1,37 @@
+set ylabel "$\\tilde\\eta$" 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 ytics 0.6, 0.1
+set mytics 2
+
+set xrange [4e-7:100]
+
+set yrange [:0.06]
+
+# default output canvas size: 12.5cm x 8.75cm
+set term lua tikz size 8,6 standalone
+
+set key top left 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 pointsize 1
+
+set logscale x
+
+plot \
+ "simpleq0.5.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 1 title "simple",\
+ "mueq0.5.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 3 title "med",\
+ "bigeq0.5.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 5 title "big",\
+ "2020-08-27+0.5.dat" using 1:(1-$3) with points ls 2 title "QMC"
+
diff --git a/figs/condensate.fig/condensate16.gnuplot b/figs/condensate.fig/condensate16.gnuplot
new file mode 100644
index 0000000..b581781
--- /dev/null
+++ b/figs/condensate.fig/condensate16.gnuplot
@@ -0,0 +1,36 @@
+set ylabel "$\\tilde\\eta$" 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 ytics 0.6, 0.1
+set mytics 2
+
+set xrange [1e-7:100]
+
+
+# default output canvas size: 12.5cm x 8.75cm
+set term lua tikz size 8,6 standalone
+
+set key top left 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 pointsize 1
+
+set logscale x
+
+plot \
+ "simpleq16.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 1 title "simple",\
+ "mueq16.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 3 title "med",\
+ "bigeq16.dat" using 1:($3<1e-5 ? $2 : 1/0) with lines ls 5 title "big",\
+ "2020-08-27+16.dat" using 1:(1-$3) with points ls 2 title "QMC"
+
diff --git a/figs/correlation.fig/Makefile b/figs/correlation.fig/Makefile
new file mode 100644
index 0000000..4a87b58
--- /dev/null
+++ b/figs/correlation.fig/Makefile
@@ -0,0 +1,43 @@
+PROJECTNAME=correlation0001 correlation02
+
+SIMPLEQ=simpleq
+
+DATS=bigeq02.dat mueq02.dat simpleq02.dat bigeq0001.dat simpleq0001.dat mueq0001.dat
+PDFS=$(addsuffix .pdf, $(PROJECTNAME))
+TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME))
+
+all: $(PDFS)
+
+$(PDFS): $(DATS)
+	gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@)
+	pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@)
+
+simpleq0001.dat:
+	julia $(SIMPLEQ)/main.jl -p "rho=1e-4;minlrho=-6;nlrho=20;eq=simpleq;N=12;P=8;xmax=55" -M medeq 2pt > $@
+mueq0001.dat:
+	julia $(SIMPLEQ)/main.jl -p "rho=1e-4;minlrho=-6;nlrho=20;eq=mueq;N=12;P=8;xmax=55" -M medeq 2pt > $@
+bigeq0001.dat:
+	julia $(SIMPLEQ)/main.jl -p "rho=1e-4;minlrho=-6;nlrho=20;v_a=16;eq=bigeq;N=12;P=8;xmax=55" -M medeq 2pt > $@
+
+simpleq02.dat:
+	julia $(SIMPLEQ)/main.jl -p "rho=2e-2;minlrho=-6;nlrho=60;eq=simpleq;N=12;P=8;xmax=10" -M medeq 2pt > $@
+mueq02.dat:
+	julia $(SIMPLEQ)/main.jl -p "rho=2e-2;minlrho=-6;nlrho=60;eq=mueq;N=12;P=8;xmax=10" -M medeq 2pt > $@
+bigeq02.dat:
+	julia $(SIMPLEQ)/main.jl -p "rho=2e-2;minlrho=-6;nlrho=60;v_a=16;eq=bigeq;N=12;P=8;xmax=10" -M medeq 2pt > $@
+
+install: $(PDFS)
+	cp $^ $(INSTALLDIR)/
+
+clean-aux:
+	rm -f $(addsuffix .tikz.tex, $(PROJECTNAME))
+	rm -f $(addsuffix .aux, $(PROJECTNAME))
+	rm -f $(addsuffix .log, $(PROJECTNAME))
+
+clean-dat:
+	rm -f $(DATS)
+
+clean-tex:
+	rm -f $(PDFS)
+
+clean: clean-aux clean-tex
diff --git a/figs/correlation.fig/correlation0001.gnuplot b/figs/correlation.fig/correlation0001.gnuplot
new file mode 100644
index 0000000..24dd900
--- /dev/null
+++ b/figs/correlation.fig/correlation0001.gnuplot
@@ -0,0 +1,28 @@
+set title "$\\rho=0.0001$"
+set ylabel "$\\frac{\\tilde C_2}{\\rho^2}$" norotate offset -1,0
+set xlabel "$|\\mathbf x|$"
+
+# 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 1
+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 "#555500" linewidth 2
+
+set pointsize 0.5
+
+rho=1e-4
+
+plot \
+ "simpleq0001.dat" using 1:($2/rho**2) with lines ls 1 title "simple",\
+ "mueq0001.dat" using 1:($2/rho**2) with lines ls 3 title "med",\
+ "bigeq0001.dat" using 1:($2/rho**2) with lines ls 5 title "big",\
+ "2020-10-01+correlation16-0001.dat" using 1:2 ls 2 title "QMC"
+
diff --git a/figs/correlation.fig/correlation02.gnuplot b/figs/correlation.fig/correlation02.gnuplot
new file mode 100644
index 0000000..eefa865
--- /dev/null
+++ b/figs/correlation.fig/correlation02.gnuplot
@@ -0,0 +1,28 @@
+set title "$\\rho=0.02$"
+set ylabel "$\\frac{\\tilde C_2}{\\rho^2}$" norotate offset -1,0
+set xlabel "$|\\mathbf x|$"
+
+# 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 1
+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 "#555500" linewidth 2
+
+set pointsize 0.5
+
+rho=2e-2
+
+plot \
+ "simpleq02.dat" using 1:($2/rho**2) with lines ls 1 title "simple",\
+ "mueq02.dat" using 1:($2/rho**2) with lines ls 3 title "med",\
+ "bigeq02.dat" using 1:($2/rho**2) with lines ls 5 title "big",\
+ "2020-08-27+correlation.dat" using 1:2 ls 2 title "QMC"
+
diff --git a/figs/energy.fig/Makefile b/figs/energy.fig/Makefile
new file mode 100644
index 0000000..6ffa316
--- /dev/null
+++ b/figs/energy.fig/Makefile
@@ -0,0 +1,43 @@
+PROJECTNAME=energy16 energy1
+
+SIMPLEQ=simpleq
+
+DATS=simpleq16.dat mueq16.dat bigeq16.dat simpleq1.dat mueq1.dat bigeq1.dat
+PDFS=$(addsuffix .pdf, $(PROJECTNAME))
+TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME))
+
+all: $(PDFS)
+
+$(PDFS): $(DATS)
+	gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@)
+	pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@)
+
+simpleq16.dat:
+	julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16" -M simpleq energy_rho > $@
+mueq16.dat:
+	julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16" -M mueq energy_rho > $@
+bigeq16.dat:
+	julia -p 8 $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;v_a=16" -M medeq energy_rho > $@
+simpleq1.dat:
+	julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=1" -M simpleq energy_rho > $@
+mueq1.dat:
+	julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=1" -M mueq energy_rho > $@
+bigeq1.dat:
+	julia -p 8 $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;v_a=1" -M medeq energy_rho > $@
+
+
+install: $(PDFS)
+	cp $^ $(INSTALLDIR)/
+
+clean-aux:
+	rm -f $(addsuffix .tikz.tex, $(PROJECTNAME))
+	rm -f $(addsuffix .aux, $(PROJECTNAME))
+	rm -f $(addsuffix .log, $(PROJECTNAME))
+
+clean-dat:
+	rm -f $(DATS)
+
+clean-tex:
+	rm -f $(PDFS)
+
+clean: clean-aux clean-tex
diff --git a/figs/energy.fig/energy1.gnuplot b/figs/energy.fig/energy1.gnuplot
new file mode 100644
index 0000000..9667dba
--- /dev/null
+++ b/figs/energy.fig/energy1.gnuplot
@@ -0,0 +1,37 @@
+set title "$v(\\mathbf x)=e^{-|\\mathbf x|}$"
+set ylabel "$\\displaystyle\\frac{\\tilde 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 pointsize 1
+
+set logscale x
+
+plot \
+ "simpleq1.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 1 title "simple" ,\
+ "mueq1.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 3 title "med",\
+ "bigeq1.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"
+
diff --git a/figs/energy.fig/energy16.gnuplot b/figs/energy.fig/energy16.gnuplot
new file mode 100644
index 0000000..4fd8de6
--- /dev/null
+++ b/figs/energy.fig/energy16.gnuplot
@@ -0,0 +1,37 @@
+set title "$v(\\mathbf x)=16e^{-|\\mathbf x|}$"
+set ylabel "$\\displaystyle\\frac{\\tilde 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 20, 40, 220
+set mytics 2
+set yrange [20:220]
+
+
+# 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 pointsize 1
+
+set logscale x
+
+plot \
+ "simpleq16.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 1 title "simple",\
+ "mueq16.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 3 title "med",\
+ "bigeq16.dat" using 1:($3<1e-5 ? ($2/$1) : 1/0) with lines ls 5 title "big",\
+ "2020-08-27+16.dat" using 1:($2/$1) with points ls 2 title "QMC"
+
diff --git a/figs/energy_relative_error.fig/Makefile b/figs/energy_relative_error.fig/Makefile
new file mode 100644
index 0000000..8c75f0f
--- /dev/null
+++ b/figs/energy_relative_error.fig/Makefile
@@ -0,0 +1,71 @@
+PROJECTNAME=cmp_energy16 cmp_energy1
+
+SIMPLEQ=simpleq
+
+DATS= simpleq16.dat mueq16.dat bigeq16.dat simpleq1.dat mueq1.dat bigeq1.dat
+PDFS=$(addsuffix .pdf, $(PROJECTNAME))
+TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME))
+
+all: $(PDFS)
+
+$(PDFS): $(DATS)
+	gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@)
+	pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@)
+
+bigeq16.dat:
+	for rho in 1e-3 5e-3 1e-2 2e-2 1e-1 1e+0; do \
+	  echo -n $$rho "  " >> $@-tmp ;\
+	  julia $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;J=10;maxiter=100;v_a=16;rho=$$rho;minlrho=-6;nlrho=50" -M medeq energy >> $@-tmp ;\
+	done
+	paste $@-tmp 2020-08-27+16.dat > $@
+	rm $@-tmp
+simpleq16.dat:
+	for lrho in -3 -2.3010299956639813 -2 -1.6989700043360187 -1 0; do \
+	  julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16;minlrho=$$lrho;nlrho=1" -M simpleq energy_rho >> $@-tmp ;\
+	done
+	paste $@-tmp 2020-08-27+16.dat > $@
+	rm $@-tmp
+mueq16.dat:
+	for lrho in -3 -2.3010299956639813 -2 -1.6989700043360187 -1 0; do \
+	  julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=16;minlrho=$$lrho;nlrho=1" -M mueq energy_rho >> $@-tmp ;\
+	done
+	paste $@-tmp 2020-08-27+16.dat > $@
+	rm $@-tmp
+
+bigeq1.dat:
+	for rho in 1e-6 1e-4 1e-3 5e-3 1e-2 1.5e-2 2e-2 5e-2 1e-1 1.0 1e1 5e1; do \
+	  echo -n $$rho "  " >> $@-tmp ;\
+	  julia $(SIMPLEQ)/main.jl -p "N=12;P=8;J=10;maxiter=100;v_a=1;rho=$$rho;minlrho=-6;nlrho=50;eq=bigeq" -M medeq energy >> $@-tmp ;\
+	done
+	paste $@-tmp 2020-10-15+energy.dat > $@
+	rm $@-tmp
+simpleq1.dat:
+	for lrho in -6 -4 -3 -2.3010299956639813 -2 -1.8239087409443189 -1.6989700043360187 -1.3010299956639813 -1 0 1 1.6989700043360187; do\
+	  julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=1;minlrho=$$lrho;nlrho=1" -M simpleq energy_rho >> $@-tmp ;\
+	done
+	paste $@-tmp 2020-10-15+energy.dat > $@
+	rm $@-tmp
+mueq1.dat:
+	for lrho in -6 -4 -3 -2.3010299956639813 -2 -1.8239087409443189 -1.6989700043360187 -1.3010299956639813 -1 0 1 1.6989700043360187; do\
+	  julia $(SIMPLEQ)/main.jl -p "tolerance=1e-11;order=100;maxiter=100;v_a=1;minlrho=$$lrho;nlrho=1" -M mueq energy_rho >> $@-tmp ;\
+	done
+	paste $@-tmp 2020-10-15+energy.dat > $@
+	rm $@-tmp
+
+
+install: $(PDFS)
+	cp $^ $(INSTALLDIR)/
+
+clean-aux:
+	rm -f $(addsuffix .tikz.tex, $(PROJECTNAME))
+	rm -f $(addsuffix .aux, $(PROJECTNAME))
+	rm -f $(addsuffix .log, $(PROJECTNAME))
+
+clean-dat:
+	rm -f $(DATS)
+	rm -f Abar.dat
+
+clean-tex:
+	rm -f $(PDFS)
+
+clean: clean-aux clean-tex
diff --git a/figs/energy_relative_error.fig/cmp_energy1.gnuplot b/figs/energy_relative_error.fig/cmp_energy1.gnuplot
new file mode 100644
index 0000000..977d249
--- /dev/null
+++ b/figs/energy_relative_error.fig/cmp_energy1.gnuplot
@@ -0,0 +1,38 @@
+set title "$v(\\mathbf x)=e^{-|\\mathbf x|}$"
+set ylabel "energy relative error" offset +3
+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 [1e-7:1e2]
+
+set ytics 1e-5, 10, 1e-1
+set ytics add ("$10^{-5}$" 1e-5, "$10^{-4}$" 1e-4, "$10^{-3}$" 1e-3, "$10^{-2}$" 1e-2, "$10^{-1}$" 1e-1)
+set yrange [1e-6:1e-1]
+
+# default output canvas size: 12.5cm x 8.75cm
+set term lua tikz size 8,6 standalone
+
+set key bottom left 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 pointsize 1
+
+max(x,y)=(x>y ? x : y)
+
+set logscale xy
+
+plot \
+ "simpleq1.dat" using 1:($3<1e-5 ? max(abs($2-$5)/$5,abs($2-$6)/$6) : 1/0) every ::1 with lp ls 1 title "simple" ,\
+ "mueq1.dat" using 1:($3<1e-5 ? max(abs($2-$5)/$5,abs($2-$6)/$6) : 1/0) every ::1 with lp ls 3 title "med" ,\
+ "bigeq1.dat" using 1:($3<1e-5 ? max(abs($2-$5)/$5,abs($2-$6)/$6) : 1/0) every ::1 with lp ls 5 title "big" ,\
+ "2020-10-15+energy.dat" using 1:(max(max(abs($6-$2)/$2,abs($6-$3)/$3),max(abs($7-$2)/$2,abs($7-$3)/$3))) every ::1 with lp ls 2 title "BDJ"
+
diff --git a/figs/energy_relative_error.fig/cmp_energy16.gnuplot b/figs/energy_relative_error.fig/cmp_energy16.gnuplot
new file mode 100644
index 0000000..b2ffb4c
--- /dev/null
+++ b/figs/energy_relative_error.fig/cmp_energy16.gnuplot
@@ -0,0 +1,37 @@
+set title "$v(\\mathbf x)=16e^{-|\\mathbf x|}$"
+set ylabel "energy relative error" offset +3
+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 1e-5, 10, 1
+set ytics add ("$10^{-5}$" 1e-5, "$10^{-4}$" 1e-4, "$10^{-3}$" 1e-3, "$10^{-2}$" 1e-2, "$10^{-1}$" 1e-1, "$1$" 1)
+set yrange [1e-5:1]
+
+# default output canvas size: 12.5cm x 8.75cm
+set term lua tikz size 8,6 standalone
+
+set key bottom left 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 pointsize 1
+
+set logscale xy
+
+max(x,y)=(x>y ? x : y)
+
+plot \
+ "simpleq16.dat" using 1:($3<1e-5 ? max(abs($2-$5)/$5,abs($2-$6)/$6) : 1/0) with lp ls 1 title "simple" ,\
+ "mueq16.dat" using 1:($3<1e-5 ? max(abs($2-$5)/$5,abs($2-$6)/$6) : 1/0) with lp ls 3 title "med" ,\
+ "bigeq16.dat" using 1:($3<1e-5 ? max(abs($2-$5)/$5,abs($2-$6)/$6) : 1/0) with lp ls 5 title "big" ,\
+ "2020-08-27+16.dat" using 1:(max(max(abs($6-$2)/$2,abs($6-$3)/$3),max(abs($7-$2)/$2,abs($7-$3)/$3))) with lp ls 2 title "BDJ"
diff --git a/figs/hardcore.fig/Makefile b/figs/hardcore.fig/Makefile
new file mode 100644
index 0000000..1cf8737
--- /dev/null
+++ b/figs/hardcore.fig/Makefile
@@ -0,0 +1,47 @@
+PROJECTNAME=hardcore_energy hardcore_condensate
+SIMPLEQ=simpleq
+
+DATS=simpleq-energy.dat simpleq_soft-energy.dat mueq_soft-energy.dat bigeq_soft-energy.dat simpleq-condensate.dat simpleq_soft-condensate.dat mueq_soft-condensate.dat bigeq_soft-condensate.dat
+
+PDFS=$(addsuffix .pdf, $(PROJECTNAME))
+TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME))
+
+all: $(PDFS)
+
+$(PDFS): $(DATS)
+	gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@)
+	pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@)
+
+simpleq-energy.dat:
+	julia -p 8 $(SIMPLEQ)/main.jl -p "N=12;P=8;J=5;nlrho=24;maxlrho=0.15" -M simpleq-hardcore energy_rho > $@
+simpleq_soft-energy.dat:
+	julia $(SIMPLEQ)/main.jl -p "maxiter=100;maxlrho=0.15;v_a=512;v_b=1.2494804800401416" -U tent -M simpleq energy_rho > $@
+mueq_soft-energy.dat:
+	julia $(SIMPLEQ)/main.jl -p "maxiter=100;maxlrho=0.15;v_a=512;v_b=1.2494804800401416" -U tent -M mueq energy_rho > $@
+bigeq_soft-energy.dat:
+	julia -p 8 $(SIMPLEQ)/main.jl -p "N=12;P=8;eq=bigeq;maxiter=100;maxlrho=0.15;nlrho=64;v_a=512;v_b=1.2494804800401416" -U tent -M medeq energy_rho_init_prevrho > $@
+
+simpleq-condensate.dat:
+	julia -p 8 $(SIMPLEQ)/main.jl -p "N=12;P=8;J=5;maxlrho=0.15;nlrho=24" -M simpleq-hardcore condensate_fraction_rho > $@
+simpleq_soft-condensate.dat:
+	julia $(SIMPLEQ)/main.jl -p "maxiter=100;maxlrho=0.15;v_a=512;v_b=1.2494804800401416" -U tent -M simpleq condensate_fraction_rho > $@
+mueq_soft-condensate.dat:
+	julia $(SIMPLEQ)/main.jl -p "maxiter=100;maxlrho=0.15;v_a=512;v_b=1.2494804800401416" -U tent -M mueq condensate_fraction_rho > $@
+bigeq_soft-condensate.dat:
+	julia -p 8 $(SIMPLEQ)/main.jl -p "N=12;P=8;eq=bigeq;maxiter=100;maxlrho=0.15;nlrho=64;v_a=512;v_b=1.2494804800401416" -U tent -M medeq condensate_fraction_rho > $@
+
+install: $(PDFS)
+	cp $^ $(INSTALLDIR)/
+
+clean-aux:
+	rm -f $(addsuffix .tikz.tex, $(PROJECTNAME))
+	rm -f $(addsuffix .aux, $(PROJECTNAME))
+	rm -f $(addsuffix .log, $(PROJECTNAME))
+
+clean-dat:
+	rm -f $(DATS)
+
+clean-tex:
+	rm -f $(PDFS)
+
+clean: clean-aux clean-tex
diff --git a/figs/hardcore.fig/hardcore_condensate.gnuplot b/figs/hardcore.fig/hardcore_condensate.gnuplot
new file mode 100644
index 0000000..974c9fd
--- /dev/null
+++ b/figs/hardcore.fig/hardcore_condensate.gnuplot
@@ -0,0 +1,39 @@
+set ylabel "$\\tilde\\eta$" norotate
+set xlabel "$\\rho$"
+
+set xtics 1e-6, 100, 1
+set xtics add ("$10^{-6}$" 0.000001, "$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0)
+unset mxtics
+
+set xrange [:sqrt(2)]
+
+#set ytics 1e-6, 10, 1
+#set ytics add ("$10^{-6}$" 0.000001, "$10^{-5}$" 1e-5, "$10^{-4}$" 0.0001, "$10^{-3}$" 1e-3, "$10^{-2}$" 0.01, "$10^{-1}$" 1e-1, "$1$" 1.0)
+#unset mytics
+set yrange [0:1]
+
+# default output canvas size: 12.5cm x 8.75cm
+set term lua tikz size 8,6 standalone
+
+set key top left 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 2
+
+set pointsize 1
+
+set logscale x
+
+plot \
+ "simpleq-condensate.dat" using 1:($3<1e-2 ? $2 : 1/0) ls 6 title "simple hc",\
+ "simpleq_soft-condensate.dat" using 1:($3<1e-2 ? $2 : 1/0) with lines ls 1 title "simple" ,\
+ "mueq_soft-condensate.dat" using 1:($3<1e-2 ? $2 : 1/0) with lines ls 3 title "med",\
+ "bigeq_soft-condensate.dat" using 1:($3<1e-2 ? $2 : 1/0) with lines ls 5 title "big",\
+ "GBC99-condensate.dat" using 1:(1-$2) ls 2 title "QMC [15]"
+
diff --git a/figs/hardcore.fig/hardcore_energy.gnuplot b/figs/hardcore.fig/hardcore_energy.gnuplot
new file mode 100644
index 0000000..9967232
--- /dev/null
+++ b/figs/hardcore.fig/hardcore_energy.gnuplot
@@ -0,0 +1,38 @@
+set ylabel "$\\frac{\\tilde e}\\rho$" norotate
+set xlabel "$\\rho$"
+
+set xtics 1e-6, 100, 1
+set xtics add ("$10^{-6}$" 0.000001, "$10^{-4}$" 0.0001, "$10^{-2}$" 0.01, "$1$" 1.0)
+unset mxtics
+
+set xrange [:sqrt(2)]
+
+#set ytics 1e-6, 10, 1
+#set ytics add ("$10^{-6}$" 0.000001, "$10^{-5}$" 1e-5, "$10^{-4}$" 0.0001, "$10^{-3}$" 1e-3, "$10^{-2}$" 0.01, "$10^{-1}$" 1e-1, "$1$" 1.0)
+#unset mytics
+
+# default output canvas size: 12.5cm x 8.75cm
+set term lua tikz size 8,6 standalone
+
+set key top left 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 2
+
+set pointsize 1
+
+set logscale xy
+
+plot \
+ "simpleq-energy.dat" using 1:($3<1e-2 ? $2/$1 : 1/0) ls 6 title "simple hc",\
+ "simpleq_soft-energy.dat" using 1:($3<1e-2 ? $2/$1 : 1/0) with lines ls 1 title "simple",\
+ "mueq_soft-energy.dat" using 1:($3<1e-2 ? $2/$1 : 1/0) every ::::95 with lines ls 3 title "med",\
+ "bigeq_soft-energy.dat" using 1:($3<1e-2 ? $2/$1 : 1/0) with lines ls 5 title "big",\
+ "GBC99-energy.dat" using 1:($2/(2*$1)) ls 2 title "QMC [15]"
+
diff --git a/figs/tan.fig/Makefile b/figs/tan.fig/Makefile
new file mode 100644
index 0000000..02ea899
--- /dev/null
+++ b/figs/tan.fig/Makefile
@@ -0,0 +1,35 @@
+PROJECTNAME=tan7 tan4
+
+SIMPLEQ=simpleq
+
+DATS=bigeq7.dat bigeq4.dat
+PDFS=$(addsuffix .pdf, $(PROJECTNAME))
+TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME))
+
+all: $(PDFS)
+
+$(PDFS): $(DATS)
+	gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@)
+	pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@)
+
+bigeq7.dat:
+	julia $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;J=10;maxiter=100;v_a=1;rho=1e-7;minlrho=-7;nlrho=1" -M medeq momentum_distribution > $@
+bigeq4.dat:
+	julia $(SIMPLEQ)/main.jl -p "eq=bigeq;N=12;P=8;J=10;maxiter=100;v_a=1;rho=1e-4;minlrho=-4;nlrho=1" -M medeq momentum_distribution > $@
+
+install: $(PDFS)
+	cp $^ $(INSTALLDIR)/
+
+clean-aux:
+	rm -f $(addsuffix .tikz.tex, $(PROJECTNAME))
+	rm -f $(addsuffix .aux, $(PROJECTNAME))
+	rm -f $(addsuffix .log, $(PROJECTNAME))
+	rm -f fit.log
+
+clean-dat:
+	rm -f $(DATS)
+
+clean-tex:
+	rm -f $(PDFS)
+
+clean: clean-aux clean-tex
diff --git a/figs/tan.fig/tan4.gnuplot b/figs/tan.fig/tan4.gnuplot
new file mode 100644
index 0000000..f1424b6
--- /dev/null
+++ b/figs/tan.fig/tan4.gnuplot
@@ -0,0 +1,48 @@
+set title "$\\rho=10^{-4}$"
+set ylabel norotate "$\\tilde{\\mathcal M}$"
+set xlabel "$|\\mathbf k|$"
+
+set xtics 1e-3, 100, 1e3
+set xtics add ("$10^{-3}$" 0.001, "$10^{-1}$" 1e-1, "$10^{1}$" 10, "$10^{3}$" 1000)
+unset mxtics
+set xrange [1e-3:2e3]
+
+set ytics 1e-40, 1e10, 1
+set ytics add ("$10^{-40}$" 1e-40, "$10^{-30}$" 1e-30, "$10^{-20}$" 1e-20, "$10^{-10}$" 1e-10)
+unset mytics
+set yrange [1e-41:1e10]
+
+# default output canvas size: 12.5cm x 8.75cm
+set term lua tikz size 8,6 standalone
+
+set key off
+
+
+# set linestyle
+set style line 1 linetype rgbcolor "#4169E1" linewidth 3
+set style line 2 linetype rgbcolor "#DC143C" linewidth 2
+set style line 3 linetype rgbcolor "#32CD32" linewidth 2
+set style line 4 linetype rgbcolor "#4B0082" linewidth 3
+set style line 5 linetype rgbcolor "#DAA520" linewidth 3
+set style line 6 linetype rgbcolor "#660000" linewidth 3
+set style line 7 linetype rgbcolor "#003333" linewidth 3
+
+set pointsize 1
+
+set logscale xy
+
+f(x)=a-4*x
+g(x)=c-12*x
+
+fit [-2:-1] f(x) "bigeq4.dat" u (log10($1)):(log10(abs($2))) via a
+fit [1:] g(x) "bigeq4.dat" u (log10($1)):(log10(abs($2))) via c
+
+set label "$|\\mathbf k|^{-4}$" at 1e-2,10**(a) textcolor "#660000"
+set label "$|\\mathbf k|^{-12}$" at 1e2,10**(c-22) textcolor "#003333"
+
+plot \
+ "bigeq4.dat" using 1:(abs($2)) every::2::118 with lines ls 5 ,\
+ 10**(f(log10(x))) ls 6 dashtype "." ,\
+ 10**(g(log10(x))) ls 7 dashtype "." ,\
+
+
diff --git a/figs/tan.fig/tan7.gnuplot b/figs/tan.fig/tan7.gnuplot
new file mode 100644
index 0000000..6ccda8b
--- /dev/null
+++ b/figs/tan.fig/tan7.gnuplot
@@ -0,0 +1,48 @@
+set title "$\\rho=10^{-7}$"
+set ylabel norotate "$\\tilde{\\mathcal M}$"
+set xlabel "$|\\mathbf k|$"
+
+set xtics 1e-3, 1000, 1e6
+set xtics add ("$10^{-3}$" 0.001, "$10^{0}$" 1, "$10^{3}$" 1000, "$10^6$" 1000000)
+unset mxtics
+set xrange [1e-3:1000000]
+
+set ytics 1e-60, 1e20, 1
+set ytics add ("$10^{-60}$" 1e-60, "$10^{-40}$" 1e-40, "$10^{-20}$" 1e-20)
+unset mytics
+set yrange [1e-61:1e10]
+
+# default output canvas size: 12.5cm x 8.75cm
+set term lua tikz size 8,6 standalone
+
+set key off
+
+
+# set linestyle
+set style line 1 linetype rgbcolor "#4169E1" linewidth 3
+set style line 2 linetype rgbcolor "#DC143C" linewidth 2
+set style line 3 linetype rgbcolor "#32CD32" linewidth 2
+set style line 4 linetype rgbcolor "#4B0082" linewidth 3
+set style line 5 linetype rgbcolor "#DAA520" linewidth 3
+set style line 6 linetype rgbcolor "#660000" linewidth 3
+set style line 7 linetype rgbcolor "#003333" linewidth 3
+
+set pointsize 1
+
+set logscale xy
+
+f(x)=a-4*x
+g(x)=c-12*x
+
+fit [-2:-1] f(x) "bigeq7.dat" u (log10($1)):(log10(abs($2))) via a
+fit [1:] g(x) "bigeq7.dat" u (log10($1)):(log10(abs($2))) via c
+
+set label "$|\\mathbf k|^{-4}$" at 1e-2,10**(a-4) textcolor "#660000"
+set label "$|\\mathbf k|^{-12}$" at 1e3,10**(c-34) textcolor "#003333"
+
+plot \
+ "bigeq7.dat" using 1:(abs($2)) with lines ls 5 ,\
+ 10**(f(log10(x))) ls 6 dashtype "." ,\
+ 10**(g(log10(x))) ls 7 dashtype "." ,\
+
+
diff --git a/figs/ux.fig/Makefile b/figs/ux.fig/Makefile
new file mode 100644
index 0000000..bd20f97
--- /dev/null
+++ b/figs/ux.fig/Makefile
@@ -0,0 +1,43 @@
+PROJECTNAME=ux-02 ux-0001
+
+SIMPLEQ=simpleq
+
+DATS=simpleq16-0.02.dat mueq16-0.02.dat bigeq16-0.02.dat simpleq16-0.0001.dat mueq16-0.0001.dat bigeq16-0.0001.dat
+PDFS=$(addsuffix .pdf, $(PROJECTNAME))
+TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME))
+
+all: $(PDFS)
+
+$(PDFS): $(DATS)
+	gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@)
+	pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@)
+
+simpleq16-0.02.dat:
+	julia $(SIMPLEQ)/main.jl -p "v_a=16;xmax=10;nx=200;rho=2e-2" -M simpleq ux >$@
+mueq16-0.02.dat:
+	julia $(SIMPLEQ)/main.jl -p "v_a=16;xmax=10;nx=200;rho=2e-2" -M mueq ux > $@
+bigeq16-0.02.dat:
+	julia $(SIMPLEQ)/main.jl -p "rho=2e-2;minlrho=-6;nlrho=60;v_a=16;eq=bigeq;N=16;P=8;xmax=10;nx=200;N=16;P=8" -M medeq ux > $@
+
+simpleq16-0.0001.dat:
+	julia $(SIMPLEQ)/main.jl -p "v_a=16;xmax=40;nx=200;rho=1e-4" -M simpleq ux >$@
+mueq16-0.0001.dat:
+	julia $(SIMPLEQ)/main.jl -p "v_a=16;xmax=40;nx=200;rho=1e-4" -M mueq ux > $@
+bigeq16-0.0001.dat:
+	julia $(SIMPLEQ)/main.jl -p "rho=1e-4;minlrho=-6;nlrho=20;v_a=16;eq=bigeq;N=16;P=8;xmax=40;nx=200;N=16;P=8" -M medeq ux > $@
+
+install: $(PDFS)
+	cp $^ $(INSTALLDIR)/
+
+clean-aux:
+	rm -f $(addsuffix .tikz.tex, $(PROJECTNAME))
+	rm -f $(addsuffix .aux, $(PROJECTNAME))
+	rm -f $(addsuffix .log, $(PROJECTNAME))
+
+clean-dat:
+	rm -f $(DATS)
+
+clean-tex:
+	rm -f $(PDFS)
+
+clean: clean-aux clean-tex
diff --git a/figs/ux.fig/ux-0001.gnuplot b/figs/ux.fig/ux-0001.gnuplot
new file mode 100644
index 0000000..4cc6fd7
--- /dev/null
+++ b/figs/ux.fig/ux-0001.gnuplot
@@ -0,0 +1,26 @@
+set title "$\\rho=0.0001$"
+set ylabel "$\\tilde g_2$" norotate offset -1,0
+set xlabel "$|\\mathbf x|$"
+
+# 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 "#555500" linewidth 2
+
+set pointsize 1
+
+plot \
+ "simpleq16-0.0001.dat" using 1:(1-$2) with lines ls 1 title "simple",\
+ "mueq16-0.0001.dat" using 1:(1-$2) with lines ls 3 title "med",\
+ "bigeq16-0.0001.dat" using 1:(1-$2) with lines ls 5 title "big",\
+ "2020-10-01+u16-0001.dat" using 1:2 ls 2 title "QMC"
+
diff --git a/figs/ux.fig/ux-02.gnuplot b/figs/ux.fig/ux-02.gnuplot
new file mode 100644
index 0000000..6379449
--- /dev/null
+++ b/figs/ux.fig/ux-02.gnuplot
@@ -0,0 +1,26 @@
+set title "$\\rho=0.02$"
+set ylabel "$\\tilde g_2$" norotate offset -1,0
+set xlabel "$|\\mathbf x|$"
+
+# 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 "#555500" linewidth 2
+
+set pointsize 0.5
+
+plot \
+ "simpleq16-0.02.dat" using 1:(1-$2) with lines ls 1 title "simple",\
+ "mueq16-0.02.dat" using 1:(1-$2) with lines ls 3 title "med",\
+ "bigeq16-0.02.dat" using 1:(1-$2) with lines ls 5 title "big",\
+ "2020-08-27+u.dat" using 1:2  ls 2 title "QMC"
+