redo of the phase diagrams

1 files changed, 1609 insertions, 0 deletions

diff --git a/plots.nb b/plots.nb
new file mode 100644
index 0000000..2df6f03
--- /dev/null
+++ b/plots.nb
@@ -0,0 +1,1609 @@
+(* Content-type: application/vnd.wolfram.mathematica *)
+
+(*** Wolfram Notebook File ***)
+(* http://www.wolfram.com/nb *)
+
+(* CreatedBy='Mathematica 12.0' *)
+
+(*CacheID: 234*)
+(* Internal cache information:
+NotebookFileLineBreakTest
+NotebookFileLineBreakTest
+NotebookDataPosition[       158,          7]
+NotebookDataLength[     62035,       1601]
+NotebookOptionsPosition[     57527,       1519]
+NotebookOutlinePosition[     57859,       1534]
+CellTagsIndexPosition[     57816,       1531]
+WindowFrame->Normal*)
+
+(* Beginning of Notebook Content *)
+Notebook[{
+Cell[BoxData[
+ RowBox[{"<<", "\"\<VariationalMethods`\>\""}]], "Input",
+ CellChangeTimes->{{3.7705578793146677`*^9, 3.770557893242214*^9}},
+ CellLabel->"In[50]:=",ExpressionUUID->"11c1c939-5ea7-4e35-820b-196b7afc7932"],
+
+Cell[BoxData[{
+ RowBox[{
+  RowBox[{"fo", "=", 
+   RowBox[{
+    RowBox[{
+     FractionBox["1", "2"], 
+     RowBox[{"(", 
+      RowBox[{
+       RowBox[{"r", " ", 
+        SuperscriptBox[
+         RowBox[{"\[Eta]", "[", "x", "]"}], "2"]}], "+", 
+       RowBox[{"c", " ", 
+        SuperscriptBox[
+         RowBox[{
+          RowBox[{"\[Eta]", "'"}], "[", "x", "]"}], "2"]}], "+", 
+       RowBox[{"d", " ", 
+        SuperscriptBox[
+         RowBox[{
+          RowBox[{"\[Eta]", "''"}], "[", "x", "]"}], "2"]}]}], ")"}]}], "+", 
+    RowBox[{"u", " ", 
+     SuperscriptBox[
+      RowBox[{"\[Eta]", "[", "x", "]"}], "4"]}]}]}], 
+  ";"}], "\[IndentingNewLine]", 
+ RowBox[{
+  RowBox[{"fe", "=", 
+   RowBox[{
+    FractionBox["1", "2"], 
+    RowBox[{"(", 
+     RowBox[{
+      RowBox[{"\[Lambda]x", " ", 
+       SuperscriptBox[
+        RowBox[{"\[Epsilon]x", "[", "x", "]"}], "2"]}], "+", 
+      RowBox[{"\[Lambda]a1", " ", 
+       SuperscriptBox[
+        RowBox[{"\[Epsilon]a1", "[", "x", "]"}], "2"]}], "+", 
+      RowBox[{"\[Lambda]a2", " ", 
+       SuperscriptBox[
+        RowBox[{"\[Epsilon]a2", "[", "x", "]"}], "2"]}]}], ")"}]}]}], 
+  ";"}], "\[IndentingNewLine]", 
+ RowBox[{
+  RowBox[{"fi", "=", 
+   RowBox[{
+    FractionBox["1", "2"], 
+    RowBox[{"(", 
+     RowBox[{
+      RowBox[{"b", " ", 
+       RowBox[{"\[Epsilon]x", "[", "x", "]"}], " ", 
+       RowBox[{"\[Eta]", "[", "x", "]"}]}], "+", 
+      RowBox[{"e1", " ", 
+       SuperscriptBox[
+        RowBox[{"\[Epsilon]a1", "[", "x", "]"}], "2"], 
+       RowBox[{"\[Eta]", "[", "x", "]"}]}], "+", 
+      RowBox[{"e2", " ", 
+       SuperscriptBox[
+        RowBox[{"\[Epsilon]a2", "[", "x", "]"}], "2"], 
+       RowBox[{"\[Eta]", "[", "x", "]"}]}], "+", 
+      RowBox[{"h1", " ", 
+       RowBox[{"\[Epsilon]a1", "[", "x", "]"}], 
+       RowBox[{"\[Epsilon]x", "[", "x", "]"}], " ", 
+       RowBox[{"\[Eta]", "[", "x", "]"}]}], "+", 
+      RowBox[{"h2", " ", 
+       RowBox[{"\[Epsilon]a2", "[", "x", "]"}], 
+       RowBox[{"\[Epsilon]x", "[", "x", "]"}], " ", 
+       RowBox[{"\[Eta]", "[", "x", "]"}]}]}], ")"}]}]}], 
+  ";"}], "\[IndentingNewLine]", 
+ RowBox[{
+  RowBox[{"f", "=", 
+   RowBox[{
+    RowBox[{"fo", "+", "fe", "+", "fi"}], "/.", 
+    RowBox[{"{", 
+     RowBox[{
+      RowBox[{"\[Lambda]a2", "\[Rule]", "0"}], ",", 
+      RowBox[{"e2", "\[Rule]", "0"}], ",", 
+      RowBox[{"h2", "\[Rule]", "0"}]}], "}"}]}]}], ";"}]}], "Input",
+ CellChangeTimes->{{3.770556239149859*^9, 3.770556315944047*^9}, {
+  3.7705578472891817`*^9, 3.770557967428043*^9}, {3.770565309038463*^9, 
+  3.770565318491414*^9}},
+ CellLabel->"In[51]:=",ExpressionUUID->"98997921-fc2d-40a3-b750-05550f4b5c59"],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{"s\[Epsilon]", "=", 
+   RowBox[{"Solve", "[", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{
+        RowBox[{"VariationalD", "[", 
+         RowBox[{"f", ",", 
+          RowBox[{"\[Epsilon]x", "[", "x", "]"}], ",", "x"}], "]"}], 
+        "\[Equal]", "0"}], ",", 
+       RowBox[{
+        RowBox[{"VariationalD", "[", 
+         RowBox[{"f", ",", 
+          RowBox[{"\[Epsilon]a1", "[", "x", "]"}], ",", "x"}], "]"}], 
+        "\[Equal]", "0"}], ",", 
+       RowBox[{
+        RowBox[{"VariationalD", "[", 
+         RowBox[{"f", ",", 
+          RowBox[{"\[Epsilon]a2", "[", "x", "]"}], ",", "x"}], "]"}], 
+        "\[Equal]", "0"}]}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"\[Epsilon]x", "[", "x", "]"}], ",", 
+       RowBox[{"\[Epsilon]a1", "[", "x", "]"}], ",", 
+       RowBox[{"\[Epsilon]a2", "[", "x", "]"}]}], "}"}]}], "]"}]}], 
+  ";"}]], "Input",
+ CellChangeTimes->{{3.770554820338584*^9, 3.770554829738276*^9}, {
+   3.770554860325389*^9, 3.770554886102066*^9}, {3.770557979668975*^9, 
+   3.77055802659697*^9}, 3.770560903908511*^9},
+ CellLabel->"In[58]:=",ExpressionUUID->"2ebc8031-6231-481f-b049-707f66f906b6"],
+
+Cell[BoxData[
+ TemplateBox[{
+  "Solve","svars",
+   "\"Equations may not give solutions for all \\\"solve\\\" variables.\"",2,
+   58,17,33203065702852886569,"Local"},
+  "MessageTemplate"]], "Message", "MSG",
+ CellChangeTimes->{3.770565319391838*^9},
+ CellLabel->
+  "During evaluation of \
+In[58]:=",ExpressionUUID->"42656e0b-6b9d-4ade-a02a-1e1a2c506fba"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"s\[Epsilon]", "/.", 
+  RowBox[{"e1", "\[Rule]", "0"}]}]], "Input",
+ CellChangeTimes->{{3.770565327531459*^9, 3.77056532843566*^9}, {
+  3.77056537362363*^9, 3.770565376951768*^9}},
+ CellLabel->"In[60]:=",ExpressionUUID->"b1ffed5d-b41b-4c8d-9f41-b0dd65eaa34f"],
+
+Cell[BoxData[
+ RowBox[{"{", 
+  RowBox[{"{", 
+   RowBox[{
+    RowBox[{
+     RowBox[{"\[Epsilon]x", "[", "x", "]"}], "\[Rule]", 
+     FractionBox[
+      RowBox[{"2", " ", "b", " ", "\[Lambda]a1", " ", 
+       RowBox[{"\[Eta]", "[", "x", "]"}]}], 
+      RowBox[{
+       RowBox[{
+        RowBox[{"-", "4"}], " ", "\[Lambda]a1", " ", "\[Lambda]x"}], "+", 
+       RowBox[{
+        SuperscriptBox["h1", "2"], " ", 
+        SuperscriptBox[
+         RowBox[{"\[Eta]", "[", "x", "]"}], "2"]}]}]]}], ",", 
+    RowBox[{
+     RowBox[{"\[Epsilon]a1", "[", "x", "]"}], "\[Rule]", 
+     RowBox[{"-", 
+      FractionBox[
+       RowBox[{"b", " ", "h1", " ", 
+        SuperscriptBox[
+         RowBox[{"\[Eta]", "[", "x", "]"}], "2"]}], 
+       RowBox[{
+        RowBox[{
+         RowBox[{"-", "4"}], " ", "\[Lambda]a1", " ", "\[Lambda]x"}], "+", 
+        RowBox[{
+         SuperscriptBox["h1", "2"], " ", 
+         SuperscriptBox[
+          RowBox[{"\[Eta]", "[", "x", "]"}], "2"]}]}]]}]}]}], "}"}], 
+  "}"}]], "Output",
+ CellChangeTimes->{3.77056532898451*^9, 3.7705653771424007`*^9},
+ CellLabel->"Out[60]=",ExpressionUUID->"c2d3dc1c-b8ea-46f2-b457-b0dbb99b7d7d"]
+}, Open  ]],
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{"f\[Eta]eff", "=", 
+   RowBox[{"Expand", "[", 
+    RowBox[{"f", "/.", 
+     RowBox[{"s\[Epsilon]", "[", 
+      RowBox[{"[", "1", "]"}], "]"}]}], "]"}]}], ";"}]], "Input",
+ CellChangeTimes->{{3.770561243988637*^9, 3.770561256154991*^9}, {
+  3.770561337365841*^9, 3.7705613556017838`*^9}},
+ CellLabel->"In[33]:=",ExpressionUUID->"d2d5bc38-647e-41ae-9585-f6514cd2c6c7"],
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{"\[Chi]", "=", 
+   RowBox[{
+    RowBox[{"D", "[", 
+     RowBox[{
+      RowBox[{
+       RowBox[{"VariationalD", "[", 
+        RowBox[{
+         RowBox[{"f", "/.", 
+          RowBox[{"s\[Epsilon]", "[", 
+           RowBox[{"[", "1", "]"}], "]"}]}], ",", 
+         RowBox[{"\[Eta]", "[", "x", "]"}], ",", "x"}], "]"}], "/.", 
+       RowBox[{"{", 
+        RowBox[{
+         RowBox[{"\[Eta]", "[", "x", "]"}], "\[Rule]", 
+         RowBox[{
+          RowBox[{"\[Eta]", "[", "x", "]"}], "[", 
+          RowBox[{"\[Eta]", "[", "xp", "]"}], "]"}]}], "}"}]}], ",", 
+      RowBox[{"\[Eta]", "[", "xp", "]"}]}], "]"}], "/.", 
+    RowBox[{"{", 
+     RowBox[{
+      RowBox[{
+       RowBox[{
+        RowBox[{"\[Eta]", "[", "x_", "]"}], "[", 
+        RowBox[{"\[Eta]", "[", "xp_", "]"}], "]"}], "\[RuleDelayed]", 
+       RowBox[{"\[Eta]", "[", "x", "]"}]}], ",", 
+      RowBox[{
+       RowBox[{
+        RowBox[{
+         RowBox[{"\[Eta]", "[", "x_", "]"}], "'"}], "[", 
+        RowBox[{"\[Eta]", "[", "xp_", "]"}], "]"}], "\[RuleDelayed]", 
+       RowBox[{"DiracDelta", "[", 
+        RowBox[{"x", "-", "xp"}], "]"}]}], ",", 
+      RowBox[{
+       RowBox[{
+        RowBox[{
+         RowBox[{"Derivative", "[", "i_", "]"}], "[", "\[Eta]", "]"}], "[", 
+        "x", "]"}], "\[RuleDelayed]", 
+       RowBox[{
+        RowBox[{
+         RowBox[{"Derivative", "[", "i", "]"}], "[", "DiracDelta", "]"}], "[", 
+        RowBox[{"x", "-", "xp"}], "]"}]}]}], "}"}]}]}], ";"}]], "Input",
+ CellChangeTimes->{{3.7705609237352962`*^9, 3.7705609537314863`*^9}, {
+  3.770561084362817*^9, 3.770561153048244*^9}, {3.770561183376079*^9, 
+  3.770561205422473*^9}, {3.7705615497853127`*^9, 3.7705616904555693`*^9}, {
+  3.770561744959778*^9, 3.770561935829486*^9}},
+ CellLabel->"In[48]:=",ExpressionUUID->"0f544c33-f435-459c-a473-d005863f95f2"],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"FourierTransform", "[", 
+  RowBox[{"\[Chi]", ",", "x", ",", "q"}], "]"}]], "Input",
+ CellChangeTimes->{{3.77056193699454*^9, 3.770561944259122*^9}},
+ CellLabel->"In[49]:=",ExpressionUUID->"b4dcf7cc-1bf9-49f0-a9af-a4dcf6893896"],
+
+Cell[BoxData["$Aborted"], "Output",
+ CellChangeTimes->{3.7705619589957113`*^9},
+ CellLabel->"Out[49]=",ExpressionUUID->"639ccfaa-595b-4074-8a0f-26ae9b6470cb"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{"1", "/", 
+   RowBox[{"(", 
+    RowBox[{"\[Lambda]", "-", 
+     RowBox[{
+      RowBox[{"b", "/", "2"}], "/", 
+      RowBox[{"(", 
+       RowBox[{
+        RowBox[{"b", "/", 
+         RowBox[{"(", 
+          RowBox[{"2", " ", "\[Lambda]"}], ")"}]}], "+", 
+        RowBox[{
+         RowBox[{"2", "/", "b"}], " ", "/", "\[Chi]t"}]}], ")"}]}]}], ")"}]}],
+   "//", "Simplify"}]], "Input",
+ CellChangeTimes->{{3.770567125728854*^9, 3.770567177945568*^9}, 
+   3.7705674449229307`*^9},
+ CellLabel->"In[64]:=",ExpressionUUID->"5b3333e0-2528-421d-98d3-c7414514d461"],
+
+Cell[BoxData[
+ FractionBox[
+  RowBox[{
+   RowBox[{"4", " ", "\[Lambda]"}], "+", 
+   RowBox[{
+    SuperscriptBox["b", "2"], " ", "\[Chi]t"}]}], 
+  RowBox[{"4", " ", 
+   SuperscriptBox["\[Lambda]", "2"]}]]], "Output",
+ CellChangeTimes->{{3.770567175117193*^9, 3.770567178165471*^9}, 
+   3.770567445208769*^9},
+ CellLabel->"Out[64]=",ExpressionUUID->"81ba9537-513b-4bf7-9467-d533e9b63498"]
+}, Open  ]],
+
+Cell[BoxData[{
+ RowBox[{
+  RowBox[{"Fu", "=", 
+   RowBox[{
+    FractionBox["V", "2"], 
+    RowBox[{"(", 
+     RowBox[{
+      RowBox[{"r", " ", 
+       SuperscriptBox["\[Eta]", "2"]}], "+", 
+      RowBox[{"2", " ", "u", " ", 
+       SuperscriptBox["\[Eta]", "4"]}]}], ")"}]}]}], 
+  ";"}], "\[IndentingNewLine]", 
+ RowBox[{
+  RowBox[{"Fm", "=", 
+   RowBox[{
+    FractionBox["V", "2"], 
+    RowBox[{"(", 
+     RowBox[{
+      RowBox[{
+       FractionBox["1", "2"], 
+       RowBox[{"(", 
+        RowBox[{"r", "+", 
+         RowBox[{"c", " ", 
+          SuperscriptBox["q", "2"]}], "+", 
+         RowBox[{"d", " ", 
+          SuperscriptBox["q", "4"]}]}], ")"}], 
+       SuperscriptBox["\[Eta]", "2"]}], "+", 
+      RowBox[{
+       FractionBox["3", "4"], "u", " ", 
+       SuperscriptBox["\[Eta]", "4"]}]}], ")"}]}]}], ";"}]}], "Input",
+ CellChangeTimes->{{3.77057072935113*^9, 3.7705707715959597`*^9}, {
+  3.7705708106509323`*^9, 3.770570826028399*^9}},
+ CellLabel->"In[55]:=",ExpressionUUID->"19f1d159-2940-4211-a220-00ea16b951c8"],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"Fm", "/.", 
+  RowBox[{"{", 
+   RowBox[{"q", "\[Rule]", 
+    SqrtBox[
+     RowBox[{
+      RowBox[{"-", "c"}], "/", 
+      RowBox[{"(", 
+       RowBox[{"2", "d"}], ")"}]}]]}], "}"}]}]], "Input",
+ CellChangeTimes->{{3.7707358225505323`*^9, 3.7707358412169533`*^9}},
+ CellLabel->"In[59]:=",ExpressionUUID->"7886c291-33ec-4ace-a89c-ce01acce2fe0"],
+
+Cell[BoxData[
+ RowBox[{
+  FractionBox["1", "2"], " ", "V", " ", 
+  RowBox[{"(", 
+   RowBox[{
+    RowBox[{
+     FractionBox["1", "2"], " ", 
+     RowBox[{"(", 
+      RowBox[{
+       RowBox[{"-", 
+        FractionBox[
+         SuperscriptBox["c", "2"], 
+         RowBox[{"4", " ", "d"}]]}], "+", "r"}], ")"}], " ", 
+     SuperscriptBox["\[Eta]", "2"]}], "+", 
+    FractionBox[
+     RowBox[{"3", " ", "u", " ", 
+      SuperscriptBox["\[Eta]", "4"]}], "4"]}], ")"}]}]], "Output",
+ CellChangeTimes->{{3.770735834827647*^9, 3.7707358414325438`*^9}},
+ CellLabel->"Out[59]=",ExpressionUUID->"691b31d2-84e5-47dc-abab-1a896d8bf522"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"eq", "=", 
+  RowBox[{"FullSimplify", "[", 
+   RowBox[{
+    RowBox[{
+     RowBox[{
+      RowBox[{"(", 
+       RowBox[{"Fm", "/.", 
+        RowBox[{"{", 
+         RowBox[{
+          RowBox[{"q", "\[Rule]", 
+           RowBox[{"Sqrt", "[", 
+            RowBox[{
+             RowBox[{"-", "c"}], "/", 
+             RowBox[{"(", 
+              RowBox[{"2", "d"}], ")"}]}], "]"}]}], ",", 
+          RowBox[{"\[Eta]", "\[Rule]", 
+           RowBox[{"Sqrt", "[", 
+            RowBox[{
+             RowBox[{"(", 
+              RowBox[{
+               RowBox[{"c", "^", "2"}], "-", 
+               RowBox[{"4", " ", "d", " ", "r"}]}], ")"}], "/", 
+             RowBox[{"(", 
+              RowBox[{"12", " ", "d", " ", "u"}], ")"}]}], "]"}]}]}], "}"}]}],
+        ")"}], "/", 
+      RowBox[{"(", 
+       RowBox[{"Fu", "/.", 
+        RowBox[{"\[Eta]", "\[Rule]", 
+         RowBox[{"Sqrt", "[", 
+          RowBox[{
+           RowBox[{"(", 
+            RowBox[{
+             RowBox[{"c", "^", "2"}], "-", 
+             RowBox[{"4", " ", "d", " ", "r"}]}], ")"}], "/", 
+           RowBox[{"(", 
+            RowBox[{"16", " ", "d", " ", "u"}], ")"}]}], "]"}]}]}], ")"}]}], "//",
+      "Expand"}], ",", 
+    RowBox[{"Assumptions", "\[Rule]", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"r", "<", "0"}], ",", 
+       RowBox[{"c", "<", "0"}], ",", 
+       RowBox[{"d", ">", "0"}], ",", 
+       RowBox[{"u", ">", "0"}]}], "}"}]}]}], "]"}]}]], "Input",
+ CellChangeTimes->{{3.7705719684602833`*^9, 3.7705719897201157`*^9}, {
+  3.770572024010268*^9, 3.770572030983783*^9}, {3.770572061179989*^9, 
+  3.77057207487043*^9}},
+ CellLabel->
+  "In[103]:=",ExpressionUUID->"4da03829-80f7-4d3c-b1df-3061c9d0fa86"],
+
+Cell[BoxData[
+ RowBox[{
+  FractionBox["2", "3"], "-", 
+  FractionBox[
+   RowBox[{"4", " ", 
+    SuperscriptBox["c", "2"]}], 
+   RowBox[{"3", " ", 
+    RowBox[{"(", 
+     RowBox[{
+      SuperscriptBox["c", "2"], "+", 
+      RowBox[{"4", " ", "d", " ", "r"}]}], ")"}]}]]}]], "Output",
+ CellChangeTimes->{{3.770571971025888*^9, 3.77057199001871*^9}, {
+   3.770572027244319*^9, 3.770572031316922*^9}, 3.770572075525268*^9},
+ CellLabel->
+  "Out[103]=",ExpressionUUID->"de2415bd-ffa4-42c0-a20c-1aaf19c5b3cd"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{"-", "T"}], " ", 
+  RowBox[{"D", "[", 
+   RowBox[{
+    RowBox[{
+     RowBox[{"Fm", "/.", 
+      RowBox[{"\[Eta]", "\[Rule]", 
+       RowBox[{"Sqrt", "[", 
+        RowBox[{
+         RowBox[{"(", 
+          RowBox[{
+           RowBox[{"c", "^", "2"}], "-", 
+           RowBox[{"4", " ", "d", " ", "r"}]}], ")"}], "/", 
+         RowBox[{"(", 
+          RowBox[{"12", " ", "d", " ", "u"}], ")"}]}], "]"}]}]}], "/.", 
+     RowBox[{"r", "\[Rule]", 
+      RowBox[{"a", 
+       RowBox[{"(", 
+        RowBox[{"T", "-", "Tc"}], ")"}]}]}]}], ",", 
+    RowBox[{"{", 
+     RowBox[{"T", ",", "2"}], "}"}]}], "]"}]}]], "Input",
+ CellChangeTimes->{{3.770642179466674*^9, 3.770642237343441*^9}},
+ CellLabel->
+  "In[264]:=",ExpressionUUID->"deee272c-c109-48f7-a88d-1ee72a679856"],
+
+Cell[BoxData[
+ FractionBox[
+  RowBox[{
+   SuperscriptBox["a", "2"], " ", "T", " ", "V"}], 
+  RowBox[{"12", " ", "u"}]]], "Output",
+ CellChangeTimes->{{3.77064218035308*^9, 3.770642237571024*^9}},
+ CellLabel->
+  "Out[264]=",ExpressionUUID->"d46cf980-3684-47bc-a454-7e432e809dda"]
+}, Open  ]],
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{"ss", "=", 
+   RowBox[{"Solve", "[", 
+    RowBox[{
+     RowBox[{"eq", "\[Equal]", "1"}], ",", "c"}], "]"}]}], ";"}]], "Input",
+ CellChangeTimes->{{3.770571271302767*^9, 3.77057133636125*^9}, {
+  3.770572037947668*^9, 3.7705720404974947`*^9}},
+ CellLabel->
+  "In[104]:=",ExpressionUUID->"67413114-8b38-4eff-8a22-b1615b0ea9bb"],
+
+Cell[CellGroupData[{
+
+Cell[BoxData["ss"], "Input",
+ CellChangeTimes->{{3.770572089368773*^9, 3.770572089551257*^9}},
+ CellLabel->
+  "In[106]:=",ExpressionUUID->"415b3bc6-5bea-472c-b832-a17a204ee60e"],
+
+Cell[BoxData[
+ RowBox[{"{", 
+  RowBox[{
+   RowBox[{"{", 
+    RowBox[{"c", "\[Rule]", 
+     RowBox[{"-", 
+      FractionBox[
+       RowBox[{"2", " ", "\[ImaginaryI]", " ", 
+        SqrtBox["d"], " ", 
+        SqrtBox["r"]}], 
+       SqrtBox["5"]]}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{"c", "\[Rule]", 
+     FractionBox[
+      RowBox[{"2", " ", "\[ImaginaryI]", " ", 
+       SqrtBox["d"], " ", 
+       SqrtBox["r"]}], 
+      SqrtBox["5"]]}], "}"}]}], "}"}]], "Output",
+ CellChangeTimes->{3.770572089716961*^9},
+ CellLabel->
+  "Out[106]=",ExpressionUUID->"c4009c3e-79a6-4a49-aa2e-fc55f5e134e4"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"FullSimplify", "[", 
+  RowBox[{
+   RowBox[{"c", "/.", 
+    RowBox[{"ss", "[", 
+     RowBox[{"[", "1", "]"}], "]"}]}], ",", 
+   RowBox[{"Assumptions", "\[Rule]", 
+    RowBox[{"{", 
+     RowBox[{
+      RowBox[{"r", "<", "0"}], ",", 
+      RowBox[{"c", "<", "0"}], ",", 
+      RowBox[{"d", ">", "0"}], ",", 
+      RowBox[{"u", ">", "0"}]}], "}"}]}]}], "]"}]], "Input",
+ CellChangeTimes->{{3.770570828320014*^9, 3.770570953489955*^9}, {
+  3.7705709875289993`*^9, 3.770571005837789*^9}, {3.770571136462347*^9, 
+  3.770571146581094*^9}, {3.770571270289277*^9, 3.770571299891721*^9}, {
+  3.7705713595269117`*^9, 3.7705713613770313`*^9}},
+ CellLabel->
+  "In[105]:=",ExpressionUUID->"7dbbfd98-b060-4a55-86f6-4759247a6580"],
+
+Cell[BoxData[
+ FractionBox[
+  RowBox[{"2", " ", 
+   SqrtBox[
+    RowBox[{
+     RowBox[{"-", "d"}], " ", "r"}]]}], 
+  SqrtBox["5"]]], "Output",
+ CellChangeTimes->{
+  3.770570862632203*^9, {3.770570908731566*^9, 3.770570919950509*^9}, 
+   3.770570954410015*^9, 3.770570989007962*^9, {3.770571132154612*^9, 
+   3.770571142707281*^9}, 3.770571266388103*^9, 3.77057139945437*^9, 
+   3.770572044736827*^9, 3.770572078258441*^9},
+ CellLabel->
+  "Out[105]=",ExpressionUUID->"de6f3c84-c35e-4336-9744-7ae8779837e0"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"sp", "=", 
+  RowBox[{"FiniteGroupData", "[", 
+   RowBox[{
+    RowBox[{"{", 
+     RowBox[{"\"\<CrystallographicPointGroup\>\"", ",", " ", "15"}], "}"}], 
+    ",", "\"\<SpaceRepresentation\>\""}], "]"}]}]], "Input",
+ CellChangeTimes->{{3.770571376872231*^9, 3.7705714181475554`*^9}, {
+  3.770571458262583*^9, 3.7705714653390207`*^9}, {3.7705715375305443`*^9, 
+  3.7705716215777903`*^9}},
+ CellLabel->"In[90]:=",ExpressionUUID->"68ae3535-4b40-4f87-b02a-90936fba2195"],
+
+Cell[BoxData[
+ RowBox[{"{", 
+  RowBox[{
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{"0", ",", 
+       RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"-", "1"}], ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"-", "1"}], ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", 
+       RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", 
+       RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", 
+       RowBox[{"-", "1"}]}], "}"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"-", "1"}], ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", 
+       RowBox[{"-", "1"}]}], "}"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", 
+       RowBox[{"-", "1"}]}], "}"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{"0", ",", 
+       RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"-", "1"}], ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", 
+       RowBox[{"-", "1"}]}], "}"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"-", "1"}], ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", 
+       RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", 
+       RowBox[{"-", "1"}]}], "}"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{"0", ",", 
+       RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", 
+       RowBox[{"-", "1"}]}], "}"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"-", "1"}], ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", 
+       RowBox[{"-", "1"}]}], "}"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", 
+       RowBox[{"-", "1"}]}], "}"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"-", "1"}], ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", 
+       RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}], ",", 
+   RowBox[{"{", 
+    RowBox[{
+     RowBox[{"{", 
+      RowBox[{"0", ",", 
+       RowBox[{"-", "1"}], ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"-", "1"}], ",", "0", ",", "0"}], "}"}], ",", 
+     RowBox[{"{", 
+      RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}]}], "}"}]], "Output",
+ CellChangeTimes->{{3.770571400483303*^9, 3.77057141869128*^9}, {
+  3.770571538202324*^9, 3.770571621806196*^9}},
+ CellLabel->"Out[90]=",ExpressionUUID->"aa33fa6c-e061-4c62-a7cd-54f992e9b260"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{"(", 
+   RowBox[{"sp", ".", 
+    RowBox[{"{", 
+     RowBox[{"\[Delta]x", ",", "\[Delta]y", ",", "\[Delta]z"}], "}"}]}], 
+   ")"}], "[", 
+  RowBox[{"[", 
+   RowBox[{"All", ",", "3"}], "]"}], "]"}]], "Input",
+ CellChangeTimes->{{3.770571623707589*^9, 3.770571653620319*^9}},
+ CellLabel->"In[94]:=",ExpressionUUID->"32ce610b-fed3-46f7-bf71-317ac1393c49"],
+
+Cell[BoxData[
+ RowBox[{"{", 
+  RowBox[{
+  "\[Delta]z", ",", "\[Delta]z", ",", "\[Delta]z", ",", "\[Delta]z", ",", 
+   RowBox[{"-", "\[Delta]z"}], ",", 
+   RowBox[{"-", "\[Delta]z"}], ",", 
+   RowBox[{"-", "\[Delta]z"}], ",", 
+   RowBox[{"-", "\[Delta]z"}], ",", 
+   RowBox[{"-", "\[Delta]z"}], ",", 
+   RowBox[{"-", "\[Delta]z"}], ",", 
+   RowBox[{"-", "\[Delta]z"}], ",", 
+   RowBox[{"-", "\[Delta]z"}], ",", "\[Delta]z", ",", "\[Delta]z", ",", 
+   "\[Delta]z", ",", "\[Delta]z"}], "}"}]], "Output",
+ CellChangeTimes->{{3.7705716421912518`*^9, 3.770571653803969*^9}},
+ CellLabel->"Out[94]=",ExpressionUUID->"29074888-bd44-4674-a7b9-f356fb297e11"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"pp1", "=", 
+  RowBox[{"Plot", "[", 
+   RowBox[{
+    RowBox[{"{", 
+     RowBox[{
+      RowBox[{"Piecewise", "[", 
+       RowBox[{"{", 
+        RowBox[{
+         RowBox[{"{", 
+          RowBox[{"\[ImaginaryI]", ",", 
+           RowBox[{"c", "\[GreaterEqual]", "0"}]}], "}"}], ",", 
+         RowBox[{"{", 
+          RowBox[{
+           RowBox[{
+            RowBox[{"c", "^", "2"}], "/", "4"}], ",", 
+           RowBox[{"c", "<", "0"}]}], "}"}]}], "}"}], "]"}], ",", 
+      RowBox[{"Piecewise", "[", 
+       RowBox[{"{", 
+        RowBox[{
+         RowBox[{"{", 
+          RowBox[{"0", ",", 
+           RowBox[{"c", "\[GreaterEqual]", "0"}]}], "}"}], ",", 
+         RowBox[{"{", 
+          RowBox[{"\[ImaginaryI]", ",", 
+           RowBox[{"c", "<", "0"}]}], "}"}]}], "}"}], "]"}], ",", 
+      RowBox[{"Piecewise", "[", 
+       RowBox[{"{", 
+        RowBox[{
+         RowBox[{"{", 
+          RowBox[{"\[Iota]", ",", 
+           RowBox[{"c", "\[GreaterEqual]", "0"}]}], "}"}], ",", 
+         RowBox[{"{", 
+          RowBox[{
+           RowBox[{
+            RowBox[{"-", "5"}], 
+            RowBox[{
+             RowBox[{"c", "^", "2"}], "/", "4"}]}], ",", 
+           RowBox[{"c", "<", "0"}]}], "}"}]}], "}"}], "]"}]}], "}"}], ",", 
+    RowBox[{"{", 
+     RowBox[{"c", ",", 
+      RowBox[{"-", "1.2"}], ",", "1.2"}], "}"}], ",", 
+    StyleBox[
+     RowBox[{"PlotRange", "\[Rule]", 
+      RowBox[{"{", 
+       RowBox[{
+        RowBox[{"{", 
+         RowBox[{
+          RowBox[{"-", "1.1"}], ",", "1.1"}], "}"}], ",", 
+        RowBox[{"{", 
+         RowBox[{
+          RowBox[{"-", "0.3"}], ",", "0.3"}], "}"}]}], "}"}]}],
+     FontWeight->"Bold"], 
+    StyleBox[",",
+     FontWeight->"Bold"], 
+    RowBox[{"Frame", "\[Rule]", "True"}], ",", 
+    RowBox[{"PlotStyle", "\[Rule]", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{
+        RowBox[{"ColorData", "[", "97", "]"}], "[", "2", "]"}], ",", 
+       RowBox[{
+        RowBox[{"ColorData", "[", "97", "]"}], "[", "1", "]"}], ",", 
+       RowBox[{"{", 
+        RowBox[{
+         RowBox[{
+          RowBox[{"ColorData", "[", "97", "]"}], "[", "1", "]"}], ",", 
+         "Dashed"}], "}"}]}], "}"}]}], ",", 
+    RowBox[{"PlotRange", "\[Rule]", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"-", "0.3"}], ",", "0.3"}], "}"}]}], ",", 
+    RowBox[{"FrameLabel", "\[Rule]", 
+     RowBox[{"{", 
+      RowBox[{
+      "\"\<\!\(\*SubscriptBox[\(c\), \(\[UpTee]\)]\)\>\"", ",", 
+       "\"\<\!\(\*OverscriptBox[\(r\), \(~\)]\) \!\(\*SubscriptBox[\(D\), \(\
+\[UpTee]\)]\)\>\""}], "}"}]}], ",", 
+    RowBox[{"FrameTicks", "\[Rule]", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"{", 
+        RowBox[{"False", ",", "False"}], "}"}], ",", 
+       RowBox[{"{", 
+        RowBox[{"False", ",", "False"}], "}"}]}], "}"}]}], ",", 
+    RowBox[{"FrameTicks", "\[Rule]", "False"}], ",", 
+    RowBox[{"ImageSize", "\[Rule]", "150"}], ",", 
+    RowBox[{"Epilog", "\[Rule]", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"Inset", "[", 
+        RowBox[{"\"\<Unordered\>\"", ",", 
+         RowBox[{"{", 
+          RowBox[{"0.5", ",", "0.16"}], "}"}]}], "]"}], ",", 
+       RowBox[{"Inset", "[", 
+        RowBox[{"\"\<Uniform\>\"", ",", 
+         RowBox[{"{", 
+          RowBox[{"0.5", ",", 
+           RowBox[{"-", "0.16"}]}], "}"}]}], "]"}], ",", 
+       RowBox[{"Inset", "[", 
+        RowBox[{"\"\<Modulated\>\"", ",", 
+         RowBox[{"{", 
+          RowBox[{
+           RowBox[{"-", "0.67"}], ",", 
+           RowBox[{"-", "0.06"}]}], "}"}]}], "]"}], ",", 
+       RowBox[{"Inset", "[", 
+        RowBox[{"\"\<(b)\>\"", ",", 
+         RowBox[{"{", 
+          RowBox[{
+           RowBox[{"-", "0.9"}], ",", 
+           RowBox[{"-", "0.25"}]}], "}"}]}], "]"}]}], "}"}]}], ",", 
+    RowBox[{"LabelStyle", "\[Rule]", "Black"}], ",", 
+    RowBox[{"Filling", "\[Rule]", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"1", "\[Rule]", 
+        RowBox[{"{", "3", "}"}]}], ",", 
+       RowBox[{"2", "\[Rule]", "Bottom"}], ",", 
+       RowBox[{"3", "\[Rule]", "Bottom"}]}], "}"}]}]}], "]"}]}]], "Input",
+ CellChangeTimes->{{3.770582847772399*^9, 3.7705832223990307`*^9}, {
+   3.770583287222756*^9, 3.7705833378257504`*^9}, {3.7705834138103867`*^9, 
+   3.770583430268839*^9}, {3.770583514179763*^9, 3.770583514808462*^9}, {
+   3.770583827834421*^9, 3.770583883010873*^9}, {3.770583914572714*^9, 
+   3.770584064549211*^9}, {3.770584799621869*^9, 3.770584875260147*^9}, 
+   3.770585053921783*^9, {3.770585084166893*^9, 3.7705850850093603`*^9}, {
+   3.774042541938417*^9, 3.7740425632301273`*^9}, {3.774042733486842*^9, 
+   3.774042741015345*^9}, {3.774042828843865*^9, 3.774042831009488*^9}, {
+   3.774042888527259*^9, 3.774042888711533*^9}, {3.7740429724914083`*^9, 
+   3.77404298458764*^9}, {3.7740433130185966`*^9, 3.7740433133040457`*^9}, {
+   3.77404363278094*^9, 3.77404383606604*^9}, 3.774044155463279*^9, {
+   3.774044195194737*^9, 3.774044285426847*^9}, {3.7740443306428967`*^9, 
+   3.77404434744316*^9}},
+ CellLabel->
+  "In[395]:=",ExpressionUUID->"5b48adc3-c181-4ba9-8919-787aedf6881e"],
+
+Cell[BoxData[
+ GraphicsBox[{GraphicsComplexBox[CompressedData["
+1:eJxl2Hk01N//B/BBqWzZkn1LSKRQKZmLijZZsmQNUfYQhZDQ2Gls2UtCKsre
+xzIvRMnaCC1KStZKki1J38/vd9yZ7+n7/mfO89zlzH3Nmft6zEjYnTN0YCQQ
+CEcYCIT/exUyY5dWVv4Byv//vEAkb/V38V8YcIbgl4tUz/2PiTj75b/Nq343
+Bc6JNmsvPKUid/bCXftEpyDwAcO1kn3PUWF/Dml11iTYhhW+cT3cgZr5bNWP
+CXyFpbANipM+raj5wMvPwj6fofI32Uk3/in6Dj3G11+Nw1zf5HnpwSbkWqFp
+1WY6BsXs1RM5YY2oNcvBMrJyBNZZ9QhxiNWji4SfC4QdwxCYF5KrylWL9B/y
+GEiWDcEFn2JOHc1HiO3E8M0O3o/gOe5YOe5ejjZ6SMcpJQwCa3/cfOaWEsRO
+UB/RGngH9xd5GT8M3kP7q9VkvAveQKVc/epI5Xy00Ukgl//OS3iyTtWNLzAH
+lf2256Ce7YEsE00teeZ09JUUX75kTYWWHvIY29sE1P9ZL8lXrh3IkS9to1sj
+UORyVlF2RTPwXF7a1OkdiEZnTviZnAL4Wu2a/KLBHfWHbjIGvVJ4VTj289eM
+ERJOrvHsuZgKZgVb7C9XSKJ2pEAQMb+OCCvPyueDcP31tfw/nBnqhrhIaR+d
+YynoTUDVsqp3F2RViUX3TMQhpdm5m6xTLbAkMSwmeC4UBS1tmr1d2gj9a8Oa
+6l77INtatdGol49AkVKluHToNJJ0SvAsd7sDtt3pj5nztRDBeJj/Yr8mcFkV
+2ItNnicy3vw2n5H3AjgUchNXiaeioluzaTExz4Fd+GDvZxEyInnkCmrFtoLk
+sOXi4hoSCvttQ5VQaYLsiaizyyF+6PYrnYceQrUgxqJSK/zQEbk++hKbRiqG
+KT0GnhRRXdQEIXIPmSPhCqf2He8EJmReHvrsZFcHCMp+8uJsjkYDfkdkeg2f
+gpiGR/fYn2Dk92uEX3ZHA3S9Ht3L5ueFRG8qiuYlVsC+4m/vh/it0CsBVeMu
+o1uwplKNNYmwE0ltOKOj2Uai1fPNcbnQSuYe4Oaa5G/jTEMdSVoFTbPPwSLS
+1/3QIhltL/YzfWfRBtv7/fax24cj5qDLaHCmCYKtU+x6vC6hJyMdSafr6uBz
+r+qpvFwXtHaKjYP07CGEcWafnMs0QJf9g/anBCaArzXxtYv1BpSeKvKpJacT
+armEs6MMY1FElbi9h0QLDBtFhR9VCkFLl1aHj3xvAEJGhPAJcW+kUVzBmLqz
+CliJGtr2vTYoUNNhSNM2D14fHZDKJO9DPfCyN7LIh3aekNaGtkC/aogytb6/
+Z9QBPb0Gf6T978HpQXGBmes6iHcDV+ZZKX8wtWFpMbAYJp75tCHun2vZkPEk
+w3tihwJS3DxFHThFpu2nnf+0Zf/uHsgv4DMbi0pD79jHi2zkqKB73vFBoHMC
+Yu74NCv3uA1MM+q7lzdEoEOUuAvIuBl2CgS1xngGoL6vnGus4ihwYPjL07MK
+bmh6T2rJrZgS4Kv6PVCETqDONVnmFf7JYJIbuelZrzCqjnL+TmDtgs1u9pnG
+InHIqi+akRLRAmdu65k6dIUg9VpeykHDRji5VD3/eNEb6VcseJOZH8GsYEpE
+/zo7pN97Se7SaD7MPfDre2GvgVQ21EwLep2mnSeRuop48UgNFOkOXTHRPovS
+pfq0bQhF8LPt6pXiQ0eQTN7BJU+2UEjpP+qQWjRHZCZJcKR03oRDKsV+1wt2
+oAbFXJ/LzTG0/WzfwMlJQjyoEoa6O9ZxoLtFxzPX9udCtOPc5ESuKgpflBVg
+ogbT5peDwuW+U64gNFfv0Cf7jNhQO/ZHOywDVEg7080LZVBQKcf4b8sk2vwl
+gxSXKL0e4K2N9nVoSUMOR/aVFB+gguDsLJtgbgLysBRsps61QQmbUs33cxHI
+0ZqP+9zVZqi7wGKWOR2ArD2V6l/PUeC4KU/IoTk3NDJgxbRjsgSmnnwQYOQ3
+Qj7H5UcO/EwBQRM9H6MgcdTKxLWQibrgcnlZMiU0DrU4NFDta1vAwyrY8oto
+KLK2zDTzCWmEnUVf4IGTD3qlKriN49wjaDkf+4Dvjh3aI7am8/rdAhBzHy14
+GqyJhI2l53+U1oCIspfXp59nUai77UTA4yIYqxbPWyo+ii4o2fOQ15GAve1+
+kZH5MvFb1/2SUnIOrNdT9G09ooxeKX+rYv5zDfwl15pMsnOjpPnXaHP2bRAq
+v418yXvRxcUXblvkvOHd88D7fU6viUn6f7SDWbJgt/m7aVYLOSTdwN61eSkR
+Nqjwnj81x49ebZQZNlwfDDcMnS6Wi04RuXVzNLuYYqBmckb6Hmkt8o8qkzR3
+sYE/S9MxYa8fEl8V/K4tsuiB676id9BoGiocGy7tNaCCigDjvQhIQNJSFfJ5
+nO1QcphSuSUrAvlue/OelNUMsc1zc9eUA5Gk8mNh/q0AE3fT6h9ZuCObEg35
+BdFSYBjQMJ64aoRavv+4F/PoOoheNYz2SpZAhhr2vybMuiBFdGRKuSoOBSpY
+7lHta4GDRNczk9qhyF3eRiY2vREe5N3zT8jxQZmLd6OO5f5b/70BUioLdmhb
+U7uYNNsdCHLxbVfn1qL1b9mabXoWvDPqg3xVJ5RGa8ClXv7YPSNH5Msjs491
+czHc1W+ZvuFwDIn4/cgWjwiHTKORn991GdAnby8mzR85wCwgdc5ykwqSf1kZ
+eJxCBm3dLiNiEQ+KyLxBLvtzG7z4Hp7aHa+GnoRPhx2SuQiPWa763qkfJAam
+R01kdmSBZV1aXG/DVuRv35vkVpMEIQv+EZvfCaJg60qPj7FX4Ggry2kXrR/E
+6sTtHdYhscDDdWDeehcr2pg3pZ/qcAbSCtTnR82A2Pm5J5whOgyCSfJBajaL
+RFm3ChMWBU8onz7XqOLRTbS7pVhGkgiEreOdB0bDJojV77+M7f5oBFkl4Uos
+WhnEv/spx2THQFlnKS2X/sObwDNOoWXheepjMtMT+rhlwPPp8+20fNW1vKxy
+L5WWJ5pq19sE9NByiezI9hTSS1qWmCb0+155Q8utvUXJOnXvaNnlY10U1WaQ
+ltlHvxr6MH2kZRvP4sBTUUO0HFE8+uXGumFabl9wsjcpHKHlhpYA4gmFMVru
+znsl/oQ6TsuPdhbpxLt/pmW9o9zU7G9faJnRvd5OijxJy0rRl7f6sk7Rcqdq
+v05kKz0vuFo49qZ9p+XFoU3VByym/8cvOGNP4hyVzT+l9pQ+H3sAZy1NUWMN
+Mn38737e/Ys7evdv+jjuvzhnc+eEa12gj//dP7d1JNft+Ugfx/0S5wty107u
+5Ke/f9y/cD526ePNx4709X/3n7c7X/ZV99DHcb+hnSeR6lnJQt8f9xeccf/A
+uaXhOCrfTJ+fqxiqsMqOvv/YZSYz/zldWlaeuHj6Vwd9/JKO1/IcE3292IVv
+bTOi9OzFXlx6Xv6/zhv54vxta/r658b7gpl3ydKy9YGWwoxW+nh+TuOB6wT6
++iPPEXeyED1nnV1e4NlCz9plVq4vFOn5r98jQPO6sG5OjCsVsNftd27Xcalv
+A+x1VTbSWsu3zYC9LtnnWNi/vx6w1xcPV0pyJlYB9vqsi+7ebb+LAXvdtSuw
+j+B0C7DXe2ZOq0vtSQTs9eM9dgt7lfwBe12RS0X9vZ0O7XPBD/b59FYPQ7P3
+NYB93tj9o6/JsAywzy85DxQMHvn33l7xOTdr3k4F/QzAPq8pylNTE4oA7POf
+pCveEXpnAfvcWIWqXjCeSsE+N1yIfrRalgLY55ynU1rPOFQA9rmz6Myjn5H3
+APtcO8+NgTJ0A7DPC9g0P296GwfY53PUwbZ/ezjg72NOkmk310VZwD7nTOvU
+kbMoAexzGX22NUMoD7DP9z/POay38d8+t+LzkfZZlzjFEMA+j8xl4fhHw3zl
+PDtWzhNBwT43Sd7y/rcRAPZ56T2+qDe+lYB9vmSQGs2pXgTY55y/suqaZHIA
++7zsaqiFjjYZsM/zkponlVddAHw/8M3FdjFX7gXsc75o6q2UE6WAfS464eo0
+QMkH7PP6hrtDZ2XTAPs8Tve5aR3LVcA+d/zGcp2h3xaC/5jmKT3bCxqjlu7q
+DxIo2OeuX302rrWNBuzzAH4nse+BboDvG7cFtr3SpSyAfS5h+SFCCRmu1Edi
+pT6+FOzzawaW5259BcA+F7+mJSSyugqwzxUaPcKDhIoB+7yhK255/kkOYJ8f
+WT/z1mZfAmCfv5WqM0lb7wv4fuNj3lWZtVoLsM93N1GuGHeXAvZ5cOymImJo
+AWCfu9cUznrvSQfs89KcGmfDGBJgn88I8bevDbKHeiOyC6cGglb2jSZe7ckU
+7HN+CUvJLxtjAfs84MYPTjc3D8D3ZUB3dtDySWHAPjdvEDgZyGWyUm+5lXqH
+UPD9uWmwI1f8iRJgnxczlt8dbrcGjbui740Md4HzlYPsXFLxFHy/3lLe8Yf1
+/DcK9nnDaKNGSaLuSv25VurvSME+jzE2dUEi9YB9buOlfG3waBVgn7+LX94y
+GV0M2Oe7ovKrcjbeAuxzMt8xg6//JAD2+X1V85uMon6AfX5za27CeNYBwD7v
+0r8mYspdBtjn7d43zxeMFAD2ecir0QdCTemAfV793MNC2DAcsM8LDn8g/xx0
+AOzzT8ARb7UQC9jnVG+VUStdT8A+H1uO/uK9IAnY5y8k1fgCRU4C9jn/1qjP
+TuRdgH1Onnoi8zndBrDPeQ9mmR5r+kPBPrc72Xpx1FcfsM8FNV/rzSsjwD7X
+KdifbCvNC9jnznKlU8Z3tgH2OapVfDi8vpeCfW7NyvMsfHs9YJ9fudHY63yu
+CrDPRd66SXk3FQP2ua+aefpejVuAfd7NYuN74mcCYJ9f2UCesC/wA+zzM9Iy
+27oWDwL2ueiL+AvBSmWAfd6nu4Zwl/8OYJ8LxZsIzK/LAOxz/8MZCpqV4YB9
+bk1KQR0pZwD7fNppUdJLQICCfU6a2i3W5xQH2OfOSnXzoeAJ2OePx0/5uuds
+BuzzQ2Jfo5TlzQD7PFWcaSI2VRWwz9dwdGaE7LEF7HNv3fyU7qlVgH1uzifu
+t9nMALDPDfqP+exu1QDs84GADjVmHn7APr8fFWX+i2EHYJ/XR5AWt6cOUbDP
+OXv6t20OEAfsc9tvrD802eYp2Oc1BgzRnImcgH3O2JV0SntdA4X0MN/ak/Sd
+9v+jMwN/Q4nee1r//w+GRpOo
+    "], {{{}, {}, {}, 
+      {RGBColor[0.368417, 0.506779, 0.709798], Opacity[0.2], EdgeForm[None], 
+       GraphicsGroupBox[
+        PolygonBox[{{168, 3, 27, 167, 201, 194, 200, 188, 193, 199, 179, 185, 
+         192, 198, 173, 177, 183, 191, 197, 169, 171, 175, 181, 189, 195, 166,
+          165, 164, 163, 162, 161, 160, 159, 158, 157, 156, 155, 154, 153, 
+         152, 151, 150, 149, 148, 147, 146, 145, 144, 143, 172, 176, 182, 190,
+          196}}]]}, {}, 
+      {RGBColor[0.368417, 0.506779, 0.709798], Opacity[0.2], EdgeForm[None], 
+       GraphicsGroupBox[PolygonBox[CompressedData["
+1:eJwl0MVOHgEYBdDBgwQCobj9FClS3N0KhCWPwAPQbd8RlxR3Da5Bgp0/LE7u
+vd8sJjOhib/jk5FBEETwjzElvNvkphyVVQyRyjBl9PLi+YisYIBE/lBMJwee
+V8pBkimhmxv3X7I//C5yaGDZPU2W08ernSR/0sWpnSJL6eHRjpOFtLJhP5Og
+h+hg374Of5+eTT1L9gkPxNoFtLBuPxGvF9HOnn0VfP+gLOpYNI+5J8bOp5k1
+e5dLPsl0q2VBP+KOaDuPJlbtHS74IMOthnn9kFui7FwaWbG3OeedH26/mdP/
+s8UZb6S7VzOrzzDNFF8wgDiG
+         "]]]}, {}, {}, {}, {}, {}, 
+      {RGBColor[0.880722, 0.611041, 0.142051], Opacity[0.2], EdgeForm[None], 
+       GraphicsGroupBox[PolygonBox[CompressedData["
+1:eJwl02VvlFEQBtAuVKhAaWlLC3UX3F3qFAjQQCA4xa3F3d1/GG7B3YNrkGDn
+ph9OZp55N7t5750taGlrbu0QFRUV4TdnaRLOqee5wEV6maVF2j/zmvvc4BK9
+zdP5o3/DA24y0KwnHfkkP+MyfeQM/urf8pBbDDLLJprP8nOu0Ffuzj/9Ox5x
+m8FmOcTwRX7BVfrJmZH2l3qvPGa4Pp94vst3GKLPJZav8kuu0V/OIsIH+Qkj
+9AUk8EO+y1B9HnF8k8eoJSTzSh6lFpHEL3mcWkYq1+UBag86MJ5yPpqPVovp
+Qg2VPDUfqRZSRyLVVNDAT8/HqqXUk0ItVTRyz/Nh6h61kDOs5wQzOcUaDofv
+4TRtHGMaJ1nBfipYx3Gms4qD4Txp5SiNLGQbWcxiLUfCOdDMSg6EfWAGqzkU
+3pPJLGU3BeFsmcpy9lEe7p8GFrCVzLCL4a6YxBJ2kR/OlCksYy9lYd+oZz5b
+6B72POwKE1nMTvIoDftMHfPYTAZV4VxpooUd5FIS7pda5rKJdCrD3jCBRWwn
+h+Kwd9Qwh42kkU0R/ahmNhvoRiopdCWZLnQmiUQSiKcTccQSQ3T4/Uj7f/4/
+metfjw==
+         "]]]}, {}, {}, {}, {}, {}, {}}, {{}, {}, {}, 
+      TagBox[
+       {RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], 
+        Opacity[1.], LineBox[CompressedData["
+1:eJwt0VcvHgAABdCvCTVbeyXU3jS1V4zYM0YIsau09qhNa1MUv9mReDi59+0+
+3MyFrZHND4FA4JoeJUgG85EQQgkjnAgi+cRnoogmhlji+M0krXwjm1Ti2WWK
+NsrIIY1jvr/tU0MRCewxTTvl5PKFExbopZZiEtlnhg4qyCOdU37QRx0lJHHA
+LJ1Uks85PxmkkQz+sEg/9ZSSzCFzdFFFARf8YogmMvnLEgM0cMMaY3zlihVG
+aOGODSZI4Yh5uvnHJtXvn60yygPbFHLJMv8Z5p4tnmnmlnWeGOeRHV7I4oxX
+/v0n8w==
+         "]]},
+       Annotation[#, "Charting`Private`Tag$56945#1"]& ], 
+      TagBox[
+       {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
+        Opacity[1.], 
+        LineBox[{168, 196, 190, 182, 176, 172, 143, 144, 145, 146, 147, 148, 
+         149, 150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162,
+          163, 164, 165, 166, 195, 189, 181, 175, 171, 169, 197, 191, 183, 
+         177, 173, 198, 192, 185, 179, 199, 193, 188, 200, 194, 201, 167}]},
+       Annotation[#, "Charting`Private`Tag$56945#2"]& ], 
+      TagBox[
+       {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
+        Opacity[1.], Dashing[{Small, Small}], LineBox[CompressedData["
+1:eJwVz0VOBUAUBMD/0SCBQHB3d3e3sOQIHACWHBKX4O7BNUiwYlHp7jeryZ+c
+npgKBgKBGcaVWTnHPAtUuSXxpV+xxzqLVLsn861fs88GjW6ZhPJon7JEjZ3C
+j37DAZs0uWURxpN9xjK1diq/+i2HbNHslk04z/Y5K9TZaf+f4U4c0a7nEcWb
+vU2LnkMEL/YFq9Tb6QS5t4/p0POJ5t3eoVXPJZJXu0cWE8+l3SULieXT7pOl
+JLJmN8gMQuinjAf3bllEHINUcOLeKQsYJoYByhnlw3uvLGGEBIaoZIxd723y
+D7rNOIE=
+         "]]},
+       Annotation[#, "Charting`Private`Tag$56945#3"]& ]}}], {}},
+  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
+  Axes->{True, True},
+  AxesLabel->{None, None},
+  AxesOrigin->{0, 0},
+  DisplayFunction->Identity,
+  Epilog->{
+    TagBox[
+     InsetBox["Unordered", {0.5, 0.16}], "InsetString"], 
+    TagBox[
+     InsetBox["Uniform", {0.5, -0.16}], "InsetString"], 
+    TagBox[
+     InsetBox["Modulated", {-0.67, -0.06}], "InsetString"], 
+    TagBox[
+     InsetBox["(b)", {-0.9, -0.25}], "InsetString"]},
+  Frame->{{True, True}, {True, True}},
+  FrameLabel->{{
+     FormBox[
+     "\"\\!\\(\\*OverscriptBox[\\(r\\), \\(~\\)]\\) \
+\\!\\(\\*SubscriptBox[\\(D\\), \\(\[UpTee]\\)]\\)\"", TraditionalForm], 
+     None}, {
+     FormBox[
+     "\"\\!\\(\\*SubscriptBox[\\(c\\), \\(\[UpTee]\\)]\\)\"", 
+      TraditionalForm], None}},
+  FrameTicks->{{False, False}, {False, False}},
+  GridLines->{None, None},
+  GridLinesStyle->Directive[
+    GrayLevel[0.5, 0.4]],
+  ImagePadding->All,
+  ImageSize->150,
+  LabelStyle->GrayLevel[0],
+  Method->{
+   "DefaultBoundaryStyle" -> Automatic, 
+    "DefaultGraphicsInteraction" -> {
+     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
+      "Effects" -> {
+       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
+        "Droplines" -> {
+         "freeformCursorMode" -> True, 
+          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
+    AbsolutePointSize[6], "ScalingFunctions" -> None, 
+    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
+        (Identity[#]& )[
+         Part[#, 1]], 
+        (Identity[#]& )[
+         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
+        (Identity[#]& )[
+         Part[#, 1]], 
+        (Identity[#]& )[
+         Part[#, 2]]}& )}, "AxesInFront" -> True},
+  PlotRange->{{-1.1, 1.1}, {-0.3, 0.3}},
+  PlotRangeClipping->True,
+  PlotRangePadding->{{0, 0}, {0, 0}},
+  Ticks->{Automatic, Automatic}]], "Output",
+ CellChangeTimes->{{3.770583979557695*^9, 3.770584065228504*^9}, 
+   3.770584805738936*^9, {3.770584844594309*^9, 3.7705848756879*^9}, 
+   3.770585054700766*^9, 3.7705850854671583`*^9, {3.774042734852344*^9, 
+   3.774042741539318*^9}, 3.774042831455653*^9, 3.774042889186082*^9, {
+   3.7740429740628643`*^9, 3.77404298508883*^9}, 3.774043313781547*^9, {
+   3.774043658001062*^9, 3.774043668608696*^9}, {3.774043699098453*^9, 
+   3.7740438365753107`*^9}, 3.774044196466015*^9, 3.774044286339246*^9, 
+   3.774044356017893*^9},
+ CellLabel->
+  "Out[395]=",ExpressionUUID->"4b43f08e-8622-44ca-8076-254d5e591534"]
+}, Open  ]],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"pp2", "=", 
+  RowBox[{"Plot", "[", 
+   RowBox[{
+    RowBox[{"{", 
+     RowBox[{
+      RowBox[{"Piecewise", "[", 
+       RowBox[{"{", 
+        RowBox[{
+         RowBox[{"{", 
+          RowBox[{"\[ImaginaryI]", ",", 
+           RowBox[{"c", "\[GreaterEqual]", "0"}]}], "}"}], ",", 
+         RowBox[{"{", 
+          RowBox[{
+           RowBox[{
+            RowBox[{"c", "^", "2"}], "/", "4"}], ",", 
+           RowBox[{"c", "<", "0"}]}], "}"}]}], "}"}], "]"}], ",", 
+      RowBox[{"Piecewise", "[", 
+       RowBox[{"{", 
+        RowBox[{
+         RowBox[{"{", 
+          RowBox[{"0", ",", 
+           RowBox[{"c", "\[GreaterEqual]", "0"}]}], "}"}], ",", 
+         RowBox[{"{", 
+          RowBox[{"\[ImaginaryI]", ",", 
+           RowBox[{"c", "<", "0"}]}], "}"}]}], "}"}], "]"}]}], "}"}], ",", 
+    RowBox[{"{", 
+     RowBox[{"c", ",", 
+      RowBox[{"-", "1.2"}], ",", "1.2"}], "}"}], ",", 
+    RowBox[{"PlotRange", "\[Rule]", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"{", 
+        RowBox[{
+         RowBox[{"-", "1.1"}], ",", "1.1"}], "}"}], ",", 
+       RowBox[{"{", 
+        RowBox[{
+         RowBox[{"-", "0.3"}], ",", "0.3"}], "}"}]}], "}"}]}], ",", 
+    RowBox[{"Frame", "\[Rule]", "True"}], ",", 
+    RowBox[{"PlotStyle", "\[Rule]", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{
+        RowBox[{"ColorData", "[", "97", "]"}], "[", "2", "]"}], ",", 
+       RowBox[{
+        RowBox[{"ColorData", "[", "97", "]"}], "[", "1", "]"}]}], "}"}]}], 
+    ",", 
+    RowBox[{"Filling", "\[Rule]", "Bottom"}], ",", 
+    RowBox[{"PlotRange", "\[Rule]", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"-", "0.3"}], ",", "0.3"}], "}"}]}], ",", 
+    RowBox[{"FrameLabel", "\[Rule]", 
+     RowBox[{"{", 
+      RowBox[{"\"\<\!\(\*SubscriptBox[\(c\), \(\[UpTee]\)]\)\>\"", ",", 
+       RowBox[{"Style", "[", 
+        RowBox[{
+        "\"\<\!\(\*OverscriptBox[\(r\), \(~\)]\) \!\(\*SubscriptBox[\(D\), \(\
+\[UpTee]\)]\)\>\"", ",", 
+         RowBox[{"Opacity", "[", "0", "]"}]}], "]"}]}], "}"}]}], ",", 
+    RowBox[{"FrameTicks", "\[Rule]", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"{", 
+        RowBox[{"False", ",", "False"}], "}"}], ",", 
+       RowBox[{"{", 
+        RowBox[{"False", ",", "False"}], "}"}]}], "}"}]}], ",", 
+    RowBox[{"Epilog", "\[Rule]", 
+     RowBox[{"{", 
+      RowBox[{
+       RowBox[{"{", 
+        RowBox[{
+         RowBox[{"Thickness", "[", "0.012", "]"}], ",", 
+         RowBox[{
+          RowBox[{"ColorData", "[", "97", "]"}], "[", "1", "]"}], ",", 
+         RowBox[{"Line", "[", 
+          RowBox[{"{", 
+           RowBox[{
+            RowBox[{"{", 
+             RowBox[{"0", ",", "0"}], "}"}], ",", 
+            RowBox[{"{", 
+             RowBox[{"0", ",", 
+              RowBox[{"-", "2"}]}], "}"}]}], "}"}], "]"}]}], "}"}], ",", 
+       RowBox[{"{", 
+        RowBox[{
+         RowBox[{"Inset", "[", 
+          RowBox[{"\"\<Unordered\>\"", ",", 
+           RowBox[{"{", 
+            RowBox[{"0.5", ",", "0.16"}], "}"}]}], "]"}], ",", 
+         RowBox[{"Inset", "[", 
+          RowBox[{"\"\<Uniform\>\"", ",", 
+           RowBox[{"{", 
+            RowBox[{"0.5", ",", 
+             RowBox[{"-", "0.16"}]}], "}"}]}], "]"}], ",", 
+         RowBox[{"Inset", "[", 
+          RowBox[{"\"\<Modulated\>\"", ",", 
+           RowBox[{"{", 
+            RowBox[{
+             RowBox[{"-", "0.55"}], ",", 
+             RowBox[{"-", "0.14"}]}], "}"}]}], "]"}], ",", 
+         RowBox[{"Inset", "[", 
+          RowBox[{"\"\<(c)\>\"", ",", 
+           RowBox[{"{", 
+            RowBox[{
+             RowBox[{"-", "0.9"}], ",", 
+             RowBox[{"-", "0.25"}]}], "}"}]}], "]"}]}], "}"}]}], "}"}]}], ",", 
+    RowBox[{"ImageSize", "\[Rule]", "150"}], ",", 
+    RowBox[{"LabelStyle", "\[Rule]", 
+     RowBox[{"{", "Black", "}"}]}]}], "]"}]}]], "Input",
+ CellChangeTimes->CompressedData["
+1:eJwdz08sQgEAx/EYBxeL5d9s+TfMKC4O9RxkHDDazJ+9Q7UiS2mczBySZBws
+zXRpLjbPMqzMgcrKoRnTQZZNzaywQs1crJo99HuH7z7nb41iZkiZy2KxKrNB
+3Wk6fluWFD36/KWBrIHqOBeuTwvroePTyYfD1qs26DESAij28rpg+IYegFt8
+Sg7tnSWTkM0xMp4T0mVY4aRWoWnCvAddG6YDeJ9J2eHhfpUDblOjbrg0P+uB
+1ovMA1QNCkOQ/Dl5goboGaPeWBuHVL7yAxIvsi/YH6NTUJ3nyrnDR+yS0TxW
+J1Zzk6L0YjEJBVG5HH6bVYxv5ZYp+Lowp4Uym94CCy3XbqijWoPwOdHHSPw1
+hOB4U1EYStYaI/Bot5tROtL7DjcD/ASMtBzTUEVyfmGQtrE1WX3tXsawoqAZ
++rUSHtQbvB1wRbPTA/8Bk8fcvQ==
+  "],
+ CellLabel->
+  "In[407]:=",ExpressionUUID->"b3b4cd09-afae-4511-b396-8eb0b0793210"],
+
+Cell[BoxData[
+ GraphicsBox[{GraphicsComplexBox[CompressedData["
+1:eJxd1wk01GsfB/CJe1W2W7aQLCUkUpabSvNDoltJCSVLiNK1ZN/KEhq7GlmS
+pSSkGtm5th8hE6GRNZfrJmsSypJU733fc+c/7+l/zpw5n/N9nmfm+c85833+
+UjaXDO3YSCRSyz+v/75vNOWRUVH5iCr/u15BUM8yw/VAPflf4785MO2b/WdW
+xcAM/n7Tao1XEwOceXJ/1RCfQf8nq24UaLyE3P4Mys9p02gdmvva8bdWaBSy
+3n9U5D2uhAoqTXs2Q6NOzzsxz3dY+pV6Uf96E8xip3FS7wQudE+7yww1gGOJ
+lkXLqXHM46mYzAh9Cs1pduYRpaO41qJzI69ELXiTPi+Rdo2gf1Zwpvr6Kjie
+z39ic9EwennmrdPTKgfukyN3WwXeoOuEfemEczFscJGJVY4bQq7+2MXUbQXA
+Q9o/qj04gI+XBdj+HnoEByr2yXrkvMZS+dqfI1SyYcNFkUzhBz34bK26k5B/
+BhR9teVlXOjENBMtbQWO2/Cecr14xZKB9E7qOPefcdD/ziDeR/4FUiN6rKOa
+wyHiWxotvaQR+QNXtrR5+MPYp5O+JmcR31c4Jryqc4b+kC3GaFCIvbnjn798
+MgKxhErXTu9baJqzzTawZDO8AEXSpjNJQPr3+uH3wePafn+fH+7A2AgZT72j
+ifD6Stk3dY92TCuTiOqcjAXl+YW7XDN0XJEakRC9FAIBK1vm7xc+xf41oQ3V
+fZ5gXbVvLLKnHJVqypRWDp2DzRfjXIudHqB1x+16jmxtIBmPCHv3a+F6ixxb
+iWl3MsVj/8D1qVXE57Pd/bCYkvUKeRUzb/4keQto9+aTo6NfIo/Ywa53m6hA
+cckU1Y5pxs0j5svLqykQ+tWKIaXagOmTkRe+BfvC/V69fJeNVSjBqVollm8P
+juVTMcmUPJwxWMWfKK4PDRgsn88RgVfX6T7wiGOHM8Uhz0+3t6Ko3Fu3dY1R
+MOh7WLbLsAklNF06xr8Hge+XUWG5XXXY3je2l9vXDcTvKoln3SxBjbwPfw0L
+W0CviLpxu9E9XF26jyuepAbSguf1tFooxP19fUw+pJSjE/nWTwu3rEuG1njt
+nIb5l2gW4eN8aJkKO/N8Tw2YteDOfl8NHtsw4AgIhKFPDRhkmWjT6XYZno22
+xp+rrsZ3XepnszIdYM0MNy/leT6Grks/vZB6AgL9Ag4k+sehjyW5z8FSEG7f
+2vSWntGGVevF0iMNYyC8TNLWRYqOI0aRYUeUg2Hl8s9ho7N1SEoJFzsp6QGa
+eSVst9TKkIusqWvbZQX+WnbDWtZZ2HdkUDqVqgGd2NMVQfMk9hPcXNfi71uB
+kacsH+8Zs4OmG/hdxu8RnhuSFPmUpAcCgutTL0j74SkrTvoJsxHy+beCsX/c
+SMeUZykek7sUQWnrDGPwLJVYTze7iX5gdydm5wiZjkcmwwDPBM1KnoH67vZP
+/H+PA47Wt/Py9S14KqW245tgOByqifUC40ZUEwlojna9At3v1622iK1BnZGp
+pguKTjC351bBvegCFCr7OkiDk9C2Ou1MiV8CmmRGbHneJQYVkb/PkrjacauT
+barxpliw6I5iqwmn4/n7Bqfs2oNhf5VAzUHDp3h6pWKxftkDjpcseVA5ynFe
+NDG8f60NHO+6LH95LBsXnvh2v7LVBFXByjlRt3PEfm4yfiJ7H65Emv7wVRPd
+C3BbulvXikTDzy3XruYdOgyyWQdXXLlDMLH/iN0t2gKZgyLFm9h2Fw+p5vkm
+5eyCOqVMz8DGaGI969d4epp0HdVJwx2ta3nhIe1Y6pr+TIyyX5iezFSHsGU5
+EXZGEDG+GBUDu8864saFWrtuuefkuqrx77qhKahKUbt9JlcWAgp5J76axxPj
+V04kOkQadKJAVZSPHT0Z7A5rFOTpMFB0fp5bNDMOXMxFGxkLLVjArVw5eykc
+7C2F+C5da8RqL07T1LkrYOmqXNu3UIPHTvEHH1pwgtFBC/Zd0wU48+xvETZh
+I/A8pjCq8zkRRU0MPI0CJKGZff1SKrRjYHFRQk1ILNDt6hi2VXR0sQgynxIP
+AUvzVFPP4KeoRpvCJxc9oVdddAfvpXKku8c8EXpgA3skVrclPcxBCeexnKYg
+LRAzlln8WFiJm1Tc3N5+vgAhztaTV+ppOF4hmbWSdwS8lG35qWspyNPymGZ0
+5hv5Q/vjgkJqBv5ioOTTfFgFelU+lHF8v4F+m9eYTPPwQfxiH2xNv48bi++D
+D3UveC+/ctom74EDL/0fd1/sI8cf/64bxJmGu88MzHGZyYNMHU/71pWbKKgq
+4H52QRh6N8iOGP4ShHcML3oXi8+Q+fQztNrZo7Fy+pPMI8oa8Iss2nzGwQq/
+r8xFh/blk3tzvlbRzDoxyUf8AYwlQ+74SGHXCQaqirA9Csc4kJEuUcha9wIL
+fqsp3ZYWDj47Xv9FSWvEmMaFhRsq/rBZpV5MeDvi5MPk2nIzZ7Aq0FRYEi/E
+VYOaxpPXjIA++/FRdHkSil8zjHJLkAJDTdsvk6btmCg+OqNSFgv+iuZ71Lvp
+eJDseH5aNwScFaxkY24/xSdZj/ziMjwhdflh5NHMf+7/3ivSqks2sKPhhYQM
+9wMMcPB5sZ9PG5j/33KVOwzMBD7tHxIqO6k8VokOtQpHHxnZgw+/rAbX1jx8
+eJw+d8fuKGzy/ZguGR6GqUajn2f1V8FbDzd2rY8ZyCEifcl8iyoo9JT6H6uh
+oq5+uxGZxg/hqXeoRd/vo5tQ/tnd1/fBs7C50EOy3ljPec3nQe0Q2f925GRq
+axqaVyfHdtVtBz/brninyngMXvIL3zogCkGWpS5vYq7ikWbOcw7aH8kVN3e2
+WgbHIP96nUXLX7lgQ9bM8Vt25zE5Z//imCmS2951hq2KCsUgikLAPqtlspxT
+iQmnoisWz116qurSQba5p1REkfLH7RNtOmOhk+SKv6bGd78xwrSCMGVO7RTy
+j/3KO906WNRWSLjwD4E4/okawmKLjHoq+zNWbn7l5Zz7C8LXHIuLSvcyCE82
+VP1idaWTcIHc6M5ESg9hqTlSv8/V14Sbu2gJetUDhB3eVEcyrIYI84y9N/Rk
+f0PYyjXP/2zkMOHwvLGpO2tHCL9YumhrkjtKuI5+hXxScZxwR1av5DPGBOFy
+NZreded3hA2O8DHSP0wRZnOutZGmThNWjgrc7sM1Q7hNvV8vopnlJUcz+67k
+WcLLw1sqdMzmfjzPEGaeN5mOTBee2dfEGs88DzCtrSVurEll5T/2eccXvqjd
+X1k5s3+ZTufLCNP2YuU/9ueO1oTqPW9YObMvmfaSv3FaTZj1/Zn9xfTRy2/u
+1tuz5v/YP3+q9XRXdLJyZt8Q+7nJcC3lZK3P7Bemmf3BNL3uGBRvZY3PVApR
+/MmGtf54ILup34I+YZVJ73NfWln5ZT23bwvsrPkSXh9aPomz7MaTV+iu8H/7
+jXjlft+SNf+lsUYQx69yhC116Lkpzaw8O+OpThKJNf/wS+BL2Mhy2oVvS/zb
+WNYtsnB8pcQyJT/b0pUySzyf/AfZzGA+
+    "], {{{}, 
+      {RGBColor[0.880722, 0.611041, 0.142051], Opacity[0.2], EdgeForm[None], 
+       GraphicsGroupBox[
+        PolygonBox[{{1, 27, 35, 127, 34, 143, 114, 139, 88, 141, 109, 134, 64,
+          142, 112, 137, 82, 140, 105, 130, 42, 113, 138, 85, 107, 132, 55, 
+         111, 136, 73, 97, 121, 25, 89, 110, 135, 65, 83, 106, 131, 47, 86, 
+         108, 133, 60, 78, 102, 126, 33, 63, 81, 104, 129, 41, 54, 72, 96, 
+         120, 24, 46, 59, 77, 101, 125, 32, 62, 80, 103, 128, 40, 53, 71, 95, 
+         119, 23, 45, 58, 76, 100, 124, 31, 39, 52, 70, 94, 118, 22, 44, 57, 
+         75, 99, 123, 30, 38, 51, 69, 93, 117, 21, 43, 56, 74, 98, 122, 29, 
+         37, 50, 68, 92, 116, 20, 28, 36, 49, 67, 91, 115, 19, 18, 17, 16, 15,
+          14, 13, 12, 11, 10, 9, 8, 7, 6, 5, 4, 203}}]]}, {}, 
+      {RGBColor[0.368417, 0.506779, 0.709798], Opacity[0.2], EdgeForm[None], 
+       GraphicsGroupBox[
+        PolygonBox[{{169, 2, 3, 168, 202, 195, 201, 189, 194, 200, 180, 186, 
+         193, 199, 174, 178, 184, 192, 198, 170, 172, 176, 182, 190, 196, 167,
+          166, 165, 164, 163, 162, 161, 160, 159, 158, 157, 156, 155, 154, 
+         153, 152, 151, 150, 149, 148, 147, 146, 145, 144, 173, 177, 183, 191,
+          197}}]]}, {}, {}, {}, {}, {}, {}}, {{}, {}, 
+      TagBox[
+       {RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], 
+        Opacity[1.], 
+        LineBox[{203, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 
+         19, 115, 91, 67, 49, 36, 28, 20, 116, 92, 68, 50, 37, 29, 122, 98, 
+         74, 56, 43, 21, 117, 93, 69, 51, 38, 30, 123, 99, 75, 57, 44, 22, 
+         118, 94, 70, 52, 39, 31, 124, 100, 76, 58, 45, 23, 119, 95, 71, 53, 
+         40, 128, 103, 80, 62, 32, 125, 101, 77, 59, 46, 24, 120, 96, 72, 54, 
+         41, 129, 104, 81, 63, 33, 126, 102, 78, 60, 133, 108, 86, 47, 131, 
+         106, 83, 65, 135, 110, 89, 25, 121, 97, 73, 136, 111, 55, 132, 107, 
+         85, 138, 113, 42, 130, 105, 140, 82, 137, 112, 142, 64, 134, 109, 
+         141, 88, 139, 114, 143, 34, 127}]},
+       Annotation[#, "Charting`Private`Tag$57747#1"]& ], 
+      TagBox[
+       {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
+        Opacity[1.], 
+        LineBox[{169, 197, 191, 183, 177, 173, 144, 145, 146, 147, 148, 149, 
+         150, 151, 152, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163,
+          164, 165, 166, 167, 196, 190, 182, 176, 172, 170, 198, 192, 184, 
+         178, 174, 199, 193, 186, 180, 200, 194, 189, 201, 195, 202, 168}]},
+       Annotation[#, "Charting`Private`Tag$57747#2"]& ]}}], {}},
+  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
+  Axes->{True, True},
+  AxesLabel->{None, None},
+  AxesOrigin->{0, 0},
+  DisplayFunction->Identity,
+  Epilog->{{
+     Thickness[0.012], 
+     RGBColor[0.368417, 0.506779, 0.709798], 
+     LineBox[{{0, 0}, {0, -2}}]}, {
+     TagBox[
+      InsetBox["Unordered", {0.5, 0.16}], "InsetString"], 
+     TagBox[
+      InsetBox["Uniform", {0.5, -0.16}], "InsetString"], 
+     TagBox[
+      InsetBox["Modulated", {-0.55, -0.14}], "InsetString"], 
+     TagBox[
+      InsetBox["(c)", {-0.9, -0.25}], "InsetString"]}},
+  Frame->{{True, True}, {True, True}},
+  FrameLabel->{{
+     FormBox[
+      StyleBox[
+      "\"\\!\\(\\*OverscriptBox[\\(r\\), \\(~\\)]\\) \
+\\!\\(\\*SubscriptBox[\\(D\\), \\(\[UpTee]\\)]\\)\"", 
+       Opacity[0], StripOnInput -> False], TraditionalForm], None}, {
+     FormBox[
+     "\"\\!\\(\\*SubscriptBox[\\(c\\), \\(\[UpTee]\\)]\\)\"", 
+      TraditionalForm], None}},
+  FrameTicks->{{False, False}, {False, False}},
+  GridLines->{None, None},
+  GridLinesStyle->Directive[
+    GrayLevel[0.5, 0.4]],
+  ImagePadding->All,
+  ImageSize->150,
+  LabelStyle->{
+    GrayLevel[0]},
+  Method->{
+   "DefaultBoundaryStyle" -> Automatic, 
+    "DefaultGraphicsInteraction" -> {
+     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
+      "Effects" -> {
+       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
+        "Droplines" -> {
+         "freeformCursorMode" -> True, 
+          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
+    AbsolutePointSize[6], "ScalingFunctions" -> None, 
+    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
+        (Identity[#]& )[
+         Part[#, 1]], 
+        (Identity[#]& )[
+         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
+        (Identity[#]& )[
+         Part[#, 1]], 
+        (Identity[#]& )[
+         Part[#, 2]]}& )}, "AxesInFront" -> True},
+  PlotRange->{{-1.1, 1.1}, {-0.3, 0.3}},
+  PlotRangeClipping->True,
+  PlotRangePadding->{{0, 0}, {0, 0}},
+  Ticks->{Automatic, Automatic}]], "Output",
+ CellChangeTimes->{
+  3.770585024257648*^9, 3.77058510255801*^9, {3.774042750307012*^9, 
+   3.7740428248367243`*^9}, 3.774042891583301*^9, {3.774042957269361*^9, 
+   3.774042988112692*^9}, {3.774043309563458*^9, 3.774043314171821*^9}, 
+   3.7740433496116037`*^9, {3.774043421473955*^9, 3.774043554482444*^9}, {
+   3.7740438106739063`*^9, 3.7740438187563467`*^9}, {3.774043851505499*^9, 
+   3.77404385681668*^9}, {3.7740439824529047`*^9, 3.7740439937194433`*^9}, 
+   3.77404414494704*^9, 3.774044294797976*^9, {3.774044507760159*^9, 
+   3.774044523599522*^9}, {3.77404466973039*^9, 3.7740447173709927`*^9}},
+ CellLabel->
+  "Out[407]=",ExpressionUUID->"a81cbc8a-e488-404d-9dc1-7ceb38914c26"]
+}, Open  ]],
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{"Export", "[", 
+   RowBox[{
+   "\"\<~/doc/research/hidden_order/doc/overleaf/phases_scalar.pdf\>\"", ",", 
+    "pp1"}], "]"}], ";"}]], "Input",
+ CellChangeTimes->{{3.774042988914166*^9, 3.774043016730563*^9}},
+ CellLabel->
+  "In[408]:=",ExpressionUUID->"d728fcff-b9a2-4c71-8171-a8abdc435944"],
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{"Export", "[", 
+   RowBox[{
+   "\"\<~/doc/research/hidden_order/doc/overleaf/phases_vector.pdf\>\"", ",", 
+    "pp2"}], "]"}], ";"}]], "Input",
+ CellChangeTimes->{{3.774042988914166*^9, 3.774043031548027*^9}},
+ CellLabel->
+  "In[409]:=",ExpressionUUID->"8bbe13fa-dece-41bc-bbeb-237206d26e81"],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{"Indexed", "[", 
+  RowBox[{"x", ",", "2"}], "]"}]], "Input",
+ CellChangeTimes->{{3.7706319874834423`*^9, 3.7706319901051188`*^9}, {
+  3.7706320530589314`*^9, 3.770632080595777*^9}},
+ CellLabel->
+  "In[235]:=",ExpressionUUID->"8ad43851-a145-489f-ad2a-d1403083b3c5"],
+
+Cell[BoxData[
+ TemplateBox[{"x","2"},
+  "IndexedDefault"]], "Output",
+ CellChangeTimes->{
+  3.770631990380934*^9, {3.770632054776945*^9, 3.7706320809112253`*^9}},
+ CellLabel->
+  "Out[235]=",ExpressionUUID->"c8927003-4b47-4284-860c-9cf89f4a8999"]
+}, Open  ]],
+
+Cell[BoxData[
+ FractionBox["\[Lambda]", 
+  RowBox[{"1", "+", 
+   RowBox[{
+    RowBox[{"c", "/", "\[Lambda]"}], "/", 
+    RowBox[{"(", 
+     RowBox[{"1", "+", 
+      RowBox[{"d", " ", "\[CapitalDelta]t"}]}], ")"}]}]}]]], "Input",
+ CellChangeTimes->{{3.770634393156534*^9, 
+  3.7706344878251953`*^9}},ExpressionUUID->"312f4d6a-918c-4fea-a12a-\
+40cac090528c"],
+
+Cell[CellGroupData[{
+
+Cell[BoxData[
+ RowBox[{
+  RowBox[{
+   SuperscriptBox[
+    RowBox[{"(", 
+     RowBox[{
+      FractionBox[
+       RowBox[{
+        RowBox[{"b", "^", "2"}], "/", 
+        RowBox[{"(", 
+         RowBox[{"4", 
+          RowBox[{"\[Lambda]", "^", "2"}]}], ")"}]}], 
+       RowBox[{
+        RowBox[{"d", " ", 
+         RowBox[{"qs", "^", "4"}]}], "+", 
+        RowBox[{"a", " ", "\[CapitalDelta]t"}]}]], "+", 
+      RowBox[{"1", "/", "\[Lambda]"}]}], ")"}], 
+    RowBox[{"-", "1"}]], "/", "\[Lambda]"}], "//", "FullSimplify"}]], "Input",\
+
+ CellChangeTimes->{{3.770634171474904*^9, 3.770634238386848*^9}, {
+  3.7706342843863277`*^9, 3.770634298943604*^9}, {3.77063434913536*^9, 
+  3.770634349567974*^9}},
+ CellLabel->
+  "In[243]:=",ExpressionUUID->"dcfed889-f347-4169-bc6d-e488d9f5ad66"],
+
+Cell[BoxData[
+ RowBox[{"1", "-", 
+  FractionBox[
+   SuperscriptBox["b", "2"], 
+   RowBox[{
+    SuperscriptBox["b", "2"], "+", 
+    RowBox[{"4", " ", "d", " ", 
+     SuperscriptBox["qs", "4"], " ", "\[Lambda]"}], "+", 
+    RowBox[{
+    "4", " ", "a", " ", "\[CapitalDelta]t", " ", "\[Lambda]"}]}]]}]], "Output",\
+
+ CellChangeTimes->{{3.770634227655197*^9, 3.770634238785038*^9}, 
+   3.770634286362319*^9, {3.770634345136135*^9, 3.770634351688837*^9}},
+ CellLabel->
+  "Out[243]=",ExpressionUUID->"2c29099a-a822-4c55-a6db-afd274e085fa"]
+}, Open  ]],
+
+Cell[BoxData[
+ RowBox[{"Integrate", "[", " ", 
+  RowBox[{
+   RowBox[{"\[Eta]", " ", 
+    RowBox[{"Cos", "[", 
+     RowBox[{"qs", " ", "x"}], "]"}]}], ",", 
+   RowBox[{"{", 
+    RowBox[{"x", ",", 
+     RowBox[{"-", "L"}], ",", "L"}], "}"}]}]}]], "Input",
+ CellChangeTimes->{{3.770637095208066*^9, 
+  3.770637109523609*^9}},ExpressionUUID->"bec159f7-aa33-4944-8826-\
+33d1e5fac3c1"]
+},
+WindowSize->{614, 764},
+WindowMargins->{{2, Automatic}, {Automatic, 2}},
+FrontEndVersion->"12.0 for Linux x86 (64-bit) (April 8, 2019)",
+StyleDefinitions->"Default.nb"
+]
+(* End of Notebook Content *)
+
+(* Internal cache information *)
+(*CellTagsOutline
+CellTagsIndex->{}
+*)
+(*CellTagsIndex
+CellTagsIndex->{}
+*)
+(*NotebookFileOutline
+Notebook[{
+Cell[558, 20, 218, 3, 31, "Input",ExpressionUUID->"11c1c939-5ea7-4e35-820b-196b7afc7932"],
+Cell[779, 25, 2633, 77, 206, "Input",ExpressionUUID->"98997921-fc2d-40a3-b750-05550f4b5c59"],
+Cell[CellGroupData[{
+Cell[3437, 106, 1188, 31, 78, "Input",ExpressionUUID->"2ebc8031-6231-481f-b049-707f66f906b6"],
+Cell[4628, 139, 353, 9, 23, "Message",ExpressionUUID->"42656e0b-6b9d-4ade-a02a-1e1a2c506fba"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[5018, 153, 282, 5, 31, "Input",ExpressionUUID->"b1ffed5d-b41b-4c8d-9f41-b0dd65eaa34f"],
+Cell[5303, 160, 1144, 32, 62, "Output",ExpressionUUID->"c2d3dc1c-b8ea-46f2-b457-b0dbb99b7d7d"]
+}, Open  ]],
+Cell[6462, 195, 400, 9, 31, "Input",ExpressionUUID->"d2d5bc38-647e-41ae-9585-f6514cd2c6c7"],
+Cell[6865, 206, 1838, 47, 124, "Input",ExpressionUUID->"0f544c33-f435-459c-a473-d005863f95f2"],
+Cell[CellGroupData[{
+Cell[8728, 257, 251, 4, 31, "Input",ExpressionUUID->"b4dcf7cc-1bf9-49f0-a9af-a4dcf6893896"],
+Cell[8982, 263, 158, 2, 35, "Output",ExpressionUUID->"639ccfaa-595b-4074-8a0f-26ae9b6470cb"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[9177, 270, 589, 17, 31, "Input",ExpressionUUID->"5b3333e0-2528-421d-98d3-c7414514d461"],
+Cell[9769, 289, 386, 10, 59, "Output",ExpressionUUID->"81ba9537-513b-4bf7-9467-d533e9b63498"]
+}, Open  ]],
+Cell[10170, 302, 1027, 32, 90, "Input",ExpressionUUID->"19f1d159-2940-4211-a220-00ea16b951c8"],
+Cell[CellGroupData[{
+Cell[11222, 338, 365, 10, 41, "Input",ExpressionUUID->"7886c291-33ec-4ace-a89c-ce01acce2fe0"],
+Cell[11590, 350, 622, 18, 62, "Output",ExpressionUUID->"691b31d2-84e5-47dc-abab-1a896d8bf522"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[12249, 373, 1710, 49, 101, "Input",ExpressionUUID->"4da03829-80f7-4d3c-b1df-3061c9d0fa86"],
+Cell[13962, 424, 502, 14, 65, "Output",ExpressionUUID->"de2415bd-ffa4-42c0-a20c-1aaf19c5b3cd"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[14501, 443, 795, 24, 55, "Input",ExpressionUUID->"deee272c-c109-48f7-a88d-1ee72a679856"],
+Cell[15299, 469, 278, 7, 58, "Output",ExpressionUUID->"d46cf980-3684-47bc-a454-7e432e809dda"]
+}, Open  ]],
+Cell[15592, 479, 359, 9, 31, "Input",ExpressionUUID->"67413114-8b38-4eff-8a22-b1615b0ea9bb"],
+Cell[CellGroupData[{
+Cell[15976, 492, 177, 3, 31, "Input",ExpressionUUID->"415b3bc6-5bea-472c-b832-a17a204ee60e"],
+Cell[16156, 497, 600, 20, 63, "Output",ExpressionUUID->"c4009c3e-79a6-4a49-aa2e-fc55f5e134e4"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[16793, 522, 737, 18, 55, "Input",ExpressionUUID->"7dbbfd98-b060-4a55-86f6-4759247a6580"],
+Cell[17533, 542, 505, 13, 63, "Output",ExpressionUUID->"de6f3c84-c35e-4336-9744-7ae8779837e0"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[18075, 560, 488, 10, 55, "Input",ExpressionUUID->"68ae3535-4b40-4f87-b02a-90936fba2195"],
+Cell[18566, 572, 4926, 157, 381, "Output",ExpressionUUID->"aa33fa6c-e061-4c62-a7cd-54f992e9b260"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[23529, 734, 385, 10, 31, "Input",ExpressionUUID->"32ce610b-fed3-46f7-bf71-317ac1393c49"],
+Cell[23917, 746, 647, 14, 59, "Output",ExpressionUUID->"29074888-bd44-4674-a7b9-f356fb297e11"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[24601, 765, 5043, 133, 360, "Input",ExpressionUUID->"5b48adc3-c181-4ba9-8919-787aedf6881e"],
+Cell[29647, 900, 11145, 212, 118, "Output",ExpressionUUID->"4b43f08e-8622-44ca-8076-254d5e591534"]
+}, Open  ]],
+Cell[CellGroupData[{
+Cell[40829, 1117, 4385, 120, 360, "Input",ExpressionUUID->"b3b4cd09-afae-4511-b396-8eb0b0793210"],
+Cell[45217, 1239, 8970, 169, 118, "Output",ExpressionUUID->"a81cbc8a-e488-404d-9dc1-7ceb38914c26"]
+}, Open  ]],
+Cell[54202, 1411, 326, 8, 78, "Input",ExpressionUUID->"d728fcff-b9a2-4c71-8171-a8abdc435944"],
+Cell[54531, 1421, 326, 8, 78, "Input",ExpressionUUID->"8bbe13fa-dece-41bc-bbeb-237206d26e81"],
+Cell[CellGroupData[{
+Cell[54882, 1433, 287, 6, 31, "Input",ExpressionUUID->"8ad43851-a145-489f-ad2a-d1403083b3c5"],
+Cell[55172, 1441, 245, 6, 35, "Output",ExpressionUUID->"c8927003-4b47-4284-860c-9cf89f4a8999"]
+}, Open  ]],
+Cell[55432, 1450, 356, 10, 50, "Input",ExpressionUUID->"312f4d6a-918c-4fea-a12a-40cac090528c"],
+Cell[CellGroupData[{
+Cell[55813, 1464, 780, 23, 52, "Input",ExpressionUUID->"dcfed889-f347-4169-bc6d-e488d9f5ad66"],
+Cell[56596, 1489, 533, 14, 62, "Output",ExpressionUUID->"2c29099a-a822-4c55-a6db-afd274e085fa"]
+}, Open  ]],
+Cell[57144, 1506, 379, 11, 31, "Input",ExpressionUUID->"bec159f7-aa33-4944-8826-33d1e5fac3c1"]
+}
+]
+*)
+