Initial commitHEADv0.0master

7 files changed, 273 insertions, 0 deletions

diff --git a/figs/action-angle.fig/Makefile b/figs/action-angle.fig/Makefile
new file mode 100644
index 0000000..0f08a68
--- /dev/null
+++ b/figs/action-angle.fig/Makefile
@@ -0,0 +1,33 @@
+PROJECTNAME=action-angle
+
+SIMPLEQ=~/Work/Research/2018+bose_gas/cmp/simpleq
+
+PDFS=$(addsuffix .pdf, $(PROJECTNAME))
+TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME))
+DATS=etas.dat
+
+all: $(PDFS)
+
+$(PDFS): $(DATS)
+	gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@)
+	pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@)
+
+etas.dat:
+	python etas.py > $@
+
+
+install: $(PDFS)
+	cp $^ $(INSTALLDIR)/
+
+clean-aux:
+	rm -f $(addsuffix .aux, $(PROJECTNAME))
+	rm -f $(addsuffix .log, $(PROJECTNAME))
+	rm -f $(addsuffix .tikz.tex, $(PROJECTNAME))
+
+clean-dat:
+	rm -f $(DATS)
+
+clean-tex:
+	rm -f $(PDFS)
+
+clean: clean-aux clean-tex
diff --git a/figs/action-angle.fig/action-angle.gnuplot b/figs/action-angle.fig/action-angle.gnuplot
new file mode 100644
index 0000000..2913583
--- /dev/null
+++ b/figs/action-angle.fig/action-angle.gnuplot
@@ -0,0 +1,25 @@
+set xlabel "$n$"
+set ylabel "$\\Delta_2\\gamma$" norotate
+
+set ytics 0.6660909, 0.0000001
+
+# 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 2
+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 2
+set style line 5 linetype rgbcolor "#DAA520" linewidth 2
+set style line 6 linetype rgbcolor "#555500" linewidth 2
+
+set pointsize 1
+
+plot \
+ "etas.dat" using ($0*2):1 ls 1 ,\
+ "etas.dat" using ($0*2+1):2 ls 2
+
diff --git a/figs/action-angle.fig/etas.dat b/figs/action-angle.fig/etas.dat
new file mode 100644
index 0000000..911f526
--- /dev/null
+++ b/figs/action-angle.fig/etas.dat
@@ -0,0 +1,49 @@
+0.6660909569512501 0.6660909569517841
+0.6660909569512068 0.6660909568572682
+0.666090881906312 0.6660908743475806
+0.6660909492953624 0.6660909569530551
+0.6660909569551257 0.666090956950617
+0.6660909569498561 0.6660909569513911
+0.6660909569508098 0.666090957737135
+0.6660908737928897 0.6660908734660769
+0.6660909574099513 0.6660909569513342
+0.6660909569510589 0.6660909569494657
+0.666090956949009 0.6660909569517685
+0.6660909569506912 0.6660909541449044
+0.6660908740414464 0.6660908770567131
+0.6660909571597132 0.6660909569500939
+0.666090956949887 0.666090956948812
+0.6660909569483953 0.666090956957927
+0.6660909270003836 0.6660909609769751
+0.6660909042589079 0.6660908702167302
+0.666090956890927 0.666090956948973
+0.6660909569482991 0.666090956948177
+0.666090956947643 0.6660909569220005
+0.666090854453731 0.6660909700231973
+0.6660909767366049 0.6660908612031291
+0.6660909569565692 0.666090956947536
+0.6660909569466511 0.6660909569475739
+0.6660909569469959 0.6660909569534238
+0.6660908799755063 0.6660908742427215
+0.6660909512195086 0.6660909569483087
+0.666090956948585 0.6660909569461111
+0.6660909569454749 0.6660909569473561
+0.6660909569467641 0.6660909579964863
+0.6660908731876072 0.6660908731982366
+0.666090958006647 0.666090956946789
+0.666090956946648 0.666090956945196
+0.666090956944789 0.6660909569490353
+0.666090956946594 0.6660909496555434
+0.6660908743084164 0.6660908815363281
+0.6660909568839747 0.6660909569457019
+0.6660909569453359 0.6660909569445963
+0.6660909569441398 0.6660909569542799
+0.6660908386690396 0.6660909880746821
+0.6660909925361512 0.6660908431096189
+0.6660909569319037 0.666090956944476
+0.6660909569436897 0.6660909569439311
+0.6660909569433824 0.6660909568669846
+0.666090876905717 0.6660908743223477
+0.6660909542783731 0.666090956945574
+0.6660909569495623 0.6660909569429838
+0.6660909569421549 0.6660909569434601
diff --git a/figs/ellipse.fig/Makefile b/figs/ellipse.fig/Makefile
new file mode 100644
index 0000000..a15ccff
--- /dev/null
+++ b/figs/ellipse.fig/Makefile
@@ -0,0 +1,24 @@
+PROJECTNAME=ellipse
+
+SIMPLEQ=~/Work/Research/2018+bose_gas/cmp/simpleq
+
+PDFS=$(addsuffix .pdf, $(PROJECTNAME))
+TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME))
+
+all: $(PDFS)
+
+$(PDFS): $(TEXS)
+	pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@)
+
+
+install: $(PDFS)
+	cp $^ $(INSTALLDIR)/
+
+clean-aux:
+	rm -f $(addsuffix .aux, $(PROJECTNAME))
+	rm -f $(addsuffix .log, $(PROJECTNAME))
+
+clean-tex:
+	rm -f $(PDFS)
+
+clean: clean-aux clean-tex
diff --git a/figs/ellipse.fig/ellipse.tikz.tex b/figs/ellipse.fig/ellipse.tikz.tex
new file mode 100644
index 0000000..96b6bff
--- /dev/null
+++ b/figs/ellipse.fig/ellipse.tikz.tex
@@ -0,0 +1,86 @@
+\documentclass{standalone}
+
+\usepackage{xcolor}
+
+\definecolor{darkgreen}{HTML}{329D32}
+
+\usepackage{tikz}
+
+\begin{document}
+\begin{tikzpicture}[scale=2.5]
+
+\draw[color=blue](0,0)ellipse(1.5 and 1.0);
+
+% angle of L
+\def\tt{15}
+% foci ae=sqrt(1-b^2/a^2)
+\def\ae{1.118}
+% collision point
+\def\Px{0.5}
+% Py=sqrt(1-Px^2/a^2)
+\def\Py{0.943}
+% slope Px/(a^2*sqrt(1-Px^2/a^2))
+\def\slope{0.236}
+% phi: atan(\Py/(\ae+\Px))
+\def\ph{15.2}
+% phi': pi-atan(\Py/(\ae-\Px))
+\def\php{108.2}
+% lambda: pi-atan(\slope)
+\def\lam{151.7}
+% psi: pi+phi-lambda
+\def\ps{43.5}
+% h=sin(\ph)*sqrt((\Px+\ae)^2+Py^2)
+\def\h{0.491}
+
+% foci
+\path(-\ae,0)coordinate(F1);
+\path(\ae,0)coordinate(F2);
+% collision point
+\path(\Px,\Py)coordinate(P);
+%Q
+\path(F1)++(90+\tt:\h)coordinate(Q);
+
+\draw[densely dotted](F1)++(\tt:3)--++(180+\tt:3.5);
+\draw[densely dotted](F2)++(\tt:0.7)--++(180+\tt:3.5);
+\draw[line width=1](P)++(\tt:1)--++(180+\tt:3.5);
+
+\draw[color=gray](F1)--(P)--(F2);
+
+\draw[color=darkgreen,dashed](P)--++(1,-\slope);
+\draw[color=darkgreen,dashed](P)--++(-1,\slope);
+
+\draw[rotate=\tt](F1)++(0.4,0)arc(0:\ph:0.4);
+\draw(F1)++(\tt+\ph/2:0.5)node{$\varphi$};
+\draw[rotate=\tt](F2)++(0.15,0)arc(0:\php:0.15);
+\draw(F2)++(\tt+\php/2:0.25)node{$\varphi'$};
+
+\draw[rotate=\tt](P)++(0.1,0)arc(0:\lam:0.1);
+\draw(P)++(\tt+\lam/2:0.2)node{$\lambda$};
+
+\draw[rotate=\tt+\lam](P)++(0.2,0)arc(0:\ps:0.2);
+\draw(P)++(\tt+\lam+\ps/2-7:0.3)node{$\psi$};
+\draw[rotate=180+\tt+\lam](P)++(0.2,0)arc(0:-\ps:0.2);
+\draw(P)++(\tt+180+\lam-\ps/2:0.3)node{$\psi$};
+
+
+\fill[color=red](F1)circle(0.03);
+\draw(F1)++(0,-0.15)node{$O$};
+\fill(F2)circle(0.03);
+\draw(F2)++(0,-0.15)node{$O'$};
+
+\fill(P)circle(0.03);
+\draw(P)++(\lam+\ps+90-\ps:0.15)node{$P$};
+
+\fill(0,0)circle(0.03);
+\draw(0,-0.15)node{$C$};
+
+\fill(Q)circle(0.03);
+\draw(Q)++(0,-0.15)node{$Q$};
+
+\draw[dotted,<->](F1)++(\tt:-0.5)--++(90+\tt:\h);
+\draw(F1)++(\tt:-0.6)++(90+\tt:\h/2)node{$h$};
+
+
+\end{tikzpicture}
+\end{document}
+
diff --git a/figs/trajectory.fig/Makefile b/figs/trajectory.fig/Makefile
new file mode 100644
index 0000000..3af70f0
--- /dev/null
+++ b/figs/trajectory.fig/Makefile
@@ -0,0 +1,30 @@
+PROJECTNAME=trajectory
+
+SIMPLEQ=~/Work/Research/2018+bose_gas/cmp/simpleq
+
+PDFS=$(addsuffix .pdf, $(PROJECTNAME))
+TEXS=trajectory.tikz.tex
+
+all: $(PDFS)
+
+$(PDFS): $(TEXS)
+	pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@)
+
+trajectory.tikz.tex:
+	python trajectory.py > $@
+
+
+install: $(PDFS)
+	cp $^ $(INSTALLDIR)/
+
+clean-aux:
+	rm -f $(addsuffix .aux, $(PROJECTNAME))
+	rm -f $(addsuffix .log, $(PROJECTNAME))
+
+clean-dat:
+	rm -f $(TEXS)
+
+clean-tex:
+	rm -f $(PDFS)
+
+clean: clean-aux clean-tex
diff --git a/figs/trajectory.fig/trajectory.tikz.tex b/figs/trajectory.fig/trajectory.tikz.tex
new file mode 100644
index 0000000..8854a5c
--- /dev/null
+++ b/figs/trajectory.fig/trajectory.tikz.tex
@@ -0,0 +1,26 @@
+\documentclass{standalone}
+\usepackage{tikz}
+\begin{document}
+\begin{tikzpicture}
+\draw[rotate=263.6206297915572,line width=1,dashed](-4.47213595499958,0) ellipse (5.0 and 2.2360679774997902);
+\begin{scope}
+\clip(-5,2)--++(15,0)--++(0,8)--++(-15,0);
+\draw[rotate=263.6206297915572,line width=2, color=blue](-4.47213595499958,0) ellipse (5.0 and 2.2360679774997902);
+\end{scope}
+\draw[rotate=231.45599405039505,line width=1,dashed](-3.725076664340498,0) ellipse (5.0 and 3.3352366999638683);
+\begin{scope}
+\clip(-5,2)--++(15,0)--++(0,8)--++(-15,0);
+\draw[rotate=231.45599405039505,line width=2, color=blue](-3.725076664340498,0) ellipse (5.0 and 3.3352366999638683);
+\end{scope}
+\draw[rotate=172.42531416885245,line width=1,dashed](-1.2502080963147824,0) ellipse (5.0 and 4.841175447751193);
+\begin{scope}
+\clip(-5,2)--++(15,0)--++(0,8)--++(-15,0);
+\draw[rotate=172.42531416885245,line width=2, color=blue](-1.2502080963147824,0) ellipse (5.0 and 4.841175447751193);
+\end{scope}
+\draw[line width=2](-5,2)--++(15,0);
+\fill[color=red](0,0)circle(0.2);
+\draw(0.5,0)node{\huge$O$};
+\draw[densely dotted,<->,line width=1](0,0)--++(0,2);
+\draw(0.3,1)node{\huge$h$};
+\end{tikzpicture}
+\end{document}