Small typo and remove functional time dependence

xmds will handle time dependence internally

1 files changed, 278 insertions, 224 deletions

diff --git a/xmds2/Shahriar_system/GenerateShahriarSystem.nb b/xmds2/Shahriar_system/GenerateShahriarSystem.nb
index f5dcb38..9968348 100644
--- a/xmds2/Shahriar_system/GenerateShahriarSystem.nb
+++ b/xmds2/Shahriar_system/GenerateShahriarSystem.nb
@@ -10,10 +10,10 @@
 NotebookFileLineBreakTest

 NotebookFileLineBreakTest

 NotebookDataPosition[       145,          7]

-NotebookDataLength[     63326,       1879]

-NotebookOptionsPosition[     61101,       1799]

-NotebookOutlinePosition[     61466,       1815]

-CellTagsIndexPosition[     61423,       1812]

+NotebookDataLength[     65931,       1933]

+NotebookOptionsPosition[     63713,       1853]

+NotebookOutlinePosition[     64070,       1869]

+CellTagsIndexPosition[     64027,       1866]

 WindowFrame->Normal*)

 

 (* Beginning of Notebook Content *)

@@ -50,29 +50,22 @@ frequencies ",
   FormBox[

    StyleBox[

     SubscriptBox["\[Omega]", "2"], "InlineMath"], TraditionalForm]]],

- " (near resonant with the ",

- Cell[BoxData[

-  StyleBox[

-   RowBox[{

-    RowBox[{

-     RowBox[{

-      RowBox[{"|", "2"}], "\[RightAngleBracket]"}], "\[Rule]", 

-     RowBox[{"|", "3"}]}], "\[RightAngleBracket]"}], "InlineMath"]]],

- " transition), and  ",

+ " and ",

  Cell[BoxData[

   FormBox[

    StyleBox[

-    SubscriptBox["\[Omega]", "c"], "InlineMath"], TraditionalForm]]],

- " (near resonant with the ",

+    SubscriptBox["\[Omega]", "3"], "InlineMath"], TraditionalForm]]],

+ "(near resonant with the ",

  Cell[BoxData[

   StyleBox[

    RowBox[{

     RowBox[{

      RowBox[{

       RowBox[{"|", "2"}], "\[RightAngleBracket]"}], "\[Rule]", 

-     RowBox[{"|", "4"}]}], "\[RightAngleBracket]"}], "InlineMath"]]],

- " transition)."

+     RowBox[{"|", "3"}]}], "\[RightAngleBracket]"}], "InlineMath"]]],

+ " transition)\[AliasDelimiter]."

 }], "Text",

+ CellChangeTimes->{{3.524013673755671*^9, 3.524013697746722*^9}},

  CellID->525777075],

 

 Cell["\<\

@@ -190,7 +183,8 @@ Cell[BoxData[
       RowBox[{"2", ",", 

        RowBox[{"{", 

         RowBox[{"Dipole", ",", "1"}], "}"}], ",", "3"}], "]"}]]}], ",", "0", 

-   ",", "0"}], "}"}]], "Output"]

+   ",", "0"}], "}"}]], "Output",

+ CellChangeTimes->{3.524014012716285*^9}]

 }, Open  ]],

 

 Cell["\<\

@@ -268,7 +262,8 @@ Cell[BoxData[
          Offset[0.4]}, 

         Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],

   Function[BoxForm`e$, 

-   MatrixForm[BoxForm`e$]]]], "Output"]

+   MatrixForm[BoxForm`e$]]]], "Output",

+ CellChangeTimes->{3.524014012961052*^9}]

 }, Open  ]],

 

 Cell["The level diagram for the system.", "MathCaption",

@@ -289,21 +284,23 @@ Cell[BoxData[
        RowBox[{"Energy", "[", "2", "]"}], "\[Rule]", 

        RowBox[{"-", "1"}]}], ",", 

       RowBox[{

-       RowBox[{"Energy", "[", "4", "]"}], "\[Rule]", ".5"}]}], "}"}]}]}], 

+       RowBox[{"Energy", "[", "3", "]"}], "\[Rule]", ".5"}]}], "}"}]}]}], 

   "]"}]], "Input",

+ CellChangeTimes->{{3.524013707327201*^9, 3.524013712245067*^9}},

  CellID->167259034],

 

 Cell[BoxData[

  GraphicsBox[{{{{}, 

      LineBox[{{-0.9, -1.5}, {-0.09999999999999998, -1.5}}]}, {{}, 

-     LineBox[{{-0.9, -1}, {-0.09999999999999998, -1}}]}, {{}, 

-     LineBox[{{0.09999999999999998, 0}, {0.9, 0}}]}}, {{}, {}, {}}, 

+     LineBox[{{-0.9, -1.}, {-0.09999999999999998, -1.}}]}, {{}, 

+     LineBox[{{0.09999999999999998, 0.}, {0.9, 0.}}]}}, {{}, {}, {}}, 

    {Arrowheads[{-0.07659574468085106, 0.07659574468085106}], 

     ArrowBox[{{-0.45999999999999996`, -1.5}, {0.45999999999999996`, 0.}}], 

-    ArrowBox[{{-0.45999999999999996`, -1}, {0.45999999999999996`, 0.}}]}, 

+    ArrowBox[{{-0.45999999999999996`, -1.}, {0.45999999999999996`, 0.}}]}, 

    {PointSize[0.0225]}},

-  ImagePadding->{{2, 2}, {2, 2}},

-  ImageSize->94.]], "Output"]

+  ImagePadding->{{2., 2}, {2., 2.}},

+  ImageSize->94.]], "Output",

+ CellChangeTimes->{3.524014013119625*^9}]

 }, Open  ]],

 

 Cell["\<\

@@ -424,7 +421,8 @@ Cell[BoxData[
          Offset[0.4]}, 

         Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],

   Function[BoxForm`e$, 

-   MatrixForm[BoxForm`e$]]]], "Output"]

+   MatrixForm[BoxForm`e$]]]], "Output",

+ CellChangeTimes->{3.52401401316455*^9}]

 }, Open  ]],

 

 Cell[TextData[{

@@ -491,7 +489,8 @@ Cell[BoxData[
          Offset[0.4]}, 

         Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],

   Function[BoxForm`e$, 

-   MatrixForm[BoxForm`e$]]]], "Output"]

+   MatrixForm[BoxForm`e$]]]], "Output",

+ CellChangeTimes->{3.524014013274888*^9}]

 }, Open  ]],

 

 Cell[TextData[{

@@ -577,7 +576,8 @@ Cell[BoxData[
          Offset[0.4]}, 

         Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],

   Function[BoxForm`e$, 

-   MatrixForm[BoxForm`e$]]]], "Output"]

+   MatrixForm[BoxForm`e$]]]], "Output",

+ CellChangeTimes->{3.524014013584079*^9}]

 }, Open  ]],

 

 Cell["Here are the evolution equations.", "MathCaption",

@@ -601,9 +601,8 @@ Cell[BoxData[
       RowBox[{

        RowBox[{

         SuperscriptBox[

-         RowBox[{"(", 

-          SubscriptBox["\[Rho]", 

-           RowBox[{"1", ",", "1"}]], ")"}], "\[Prime]",

+         SubscriptBox["\[Rho]", 

+          RowBox[{"1", ",", "1"}]], "\[Prime]",

          MultilineFunction->None], "[", "t", "]"}], "\[Equal]", 

        RowBox[{

         FractionBox["\[Gamma]t", "2"], "-", 

@@ -641,9 +640,8 @@ Cell[BoxData[
       RowBox[{

        RowBox[{

         SuperscriptBox[

-         RowBox[{"(", 

-          SubscriptBox["\[Rho]", 

-           RowBox[{"1", ",", "2"}]], ")"}], "\[Prime]",

+         SubscriptBox["\[Rho]", 

+          RowBox[{"1", ",", "2"}]], "\[Prime]",

          MultilineFunction->None], "[", "t", "]"}], "\[Equal]", 

        RowBox[{

         RowBox[{

@@ -701,9 +699,8 @@ Cell[BoxData[
       RowBox[{

        RowBox[{

         SuperscriptBox[

-         RowBox[{"(", 

-          SubscriptBox["\[Rho]", 

-           RowBox[{"1", ",", "3"}]], ")"}], "\[Prime]",

+         SubscriptBox["\[Rho]", 

+          RowBox[{"1", ",", "3"}]], "\[Prime]",

          MultilineFunction->None], "[", "t", "]"}], "\[Equal]", 

        RowBox[{

         RowBox[{

@@ -771,9 +768,8 @@ Cell[BoxData[
       RowBox[{

        RowBox[{

         SuperscriptBox[

-         RowBox[{"(", 

-          SubscriptBox["\[Rho]", 

-           RowBox[{"2", ",", "1"}]], ")"}], "\[Prime]",

+         SubscriptBox["\[Rho]", 

+          RowBox[{"2", ",", "1"}]], "\[Prime]",

          MultilineFunction->None], "[", "t", "]"}], "\[Equal]", 

        RowBox[{

         RowBox[{

@@ -830,9 +826,8 @@ Cell[BoxData[
       RowBox[{

        RowBox[{

         SuperscriptBox[

-         RowBox[{"(", 

-          SubscriptBox["\[Rho]", 

-           RowBox[{"2", ",", "2"}]], ")"}], "\[Prime]",

+         SubscriptBox["\[Rho]", 

+          RowBox[{"2", ",", "2"}]], "\[Prime]",

          MultilineFunction->None], "[", "t", "]"}], "\[Equal]", 

        RowBox[{

         FractionBox["\[Gamma]t", "2"], "-", 

@@ -897,9 +892,8 @@ Cell[BoxData[
       RowBox[{

        RowBox[{

         SuperscriptBox[

-         RowBox[{"(", 

-          SubscriptBox["\[Rho]", 

-           RowBox[{"2", ",", "3"}]], ")"}], "\[Prime]",

+         SubscriptBox["\[Rho]", 

+          RowBox[{"2", ",", "3"}]], "\[Prime]",

          MultilineFunction->None], "[", "t", "]"}], "\[Equal]", 

        RowBox[{

         RowBox[{

@@ -977,9 +971,8 @@ Cell[BoxData[
       RowBox[{

        RowBox[{

         SuperscriptBox[

-         RowBox[{"(", 

-          SubscriptBox["\[Rho]", 

-           RowBox[{"3", ",", "1"}]], ")"}], "\[Prime]",

+         SubscriptBox["\[Rho]", 

+          RowBox[{"3", ",", "1"}]], "\[Prime]",

          MultilineFunction->None], "[", "t", "]"}], "\[Equal]", 

        RowBox[{

         RowBox[{

@@ -1044,9 +1037,8 @@ Cell[BoxData[
       RowBox[{

        RowBox[{

         SuperscriptBox[

-         RowBox[{"(", 

-          SubscriptBox["\[Rho]", 

-           RowBox[{"3", ",", "2"}]], ")"}], "\[Prime]",

+         SubscriptBox["\[Rho]", 

+          RowBox[{"3", ",", "2"}]], "\[Prime]",

          MultilineFunction->None], "[", "t", "]"}], "\[Equal]", 

        RowBox[{

         RowBox[{

@@ -1122,9 +1114,8 @@ Cell[BoxData[
       RowBox[{

        RowBox[{

         SuperscriptBox[

-         RowBox[{"(", 

-          SubscriptBox["\[Rho]", 

-           RowBox[{"3", ",", "3"}]], ")"}], "\[Prime]",

+         SubscriptBox["\[Rho]", 

+          RowBox[{"3", ",", "3"}]], "\[Prime]",

          MultilineFunction->None], "[", "t", "]"}], "\[Equal]", 

        RowBox[{

         RowBox[{"\[Gamma]p", " ", 

@@ -1217,7 +1208,8 @@ Cell[BoxData[
         Offset[0.2]}, "RowsIndexed" -> {}}],

    Column],

   Function[BoxForm`e$, 

-   TableForm[BoxForm`e$]]]], "Output"]

+   TableForm[BoxForm`e$]]]], "Output",

+ CellChangeTimes->{3.524014013741555*^9}]

 }, Open  ]],

 

 Cell["Convert to c form.", "Text"],

@@ -1252,19 +1244,19 @@ Cell[BoxData[{
            FractionBox[

             RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Delta]"}], "2"], "+", 

            RowBox[{"\[ImaginaryI]", " ", "\[Phi]2"}]}]], " ", "\[Rule]", 

-         RowBox[{"2", 

-          RowBox[{"d", "[", "t", "]"}], 

+         RowBox[{"2", "d", " ", 

           RowBox[{"E2", "/", "\[CapitalOmega]2"}]}]}], ",", 

         RowBox[{

          SuperscriptBox["\[ExponentialE]", 

           RowBox[{

            RowBox[{"-", 

             FractionBox[

-             RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Delta]"}], "2"]}], 

-           "-", 

+             RowBox[{GridBox[{

+                {"\[ImaginaryI]"},

+                {"\[Placeholder]"}

+               }], " ", "t", " ", "\[Delta]"}], "2"]}], "-", 

            RowBox[{"\[ImaginaryI]", " ", "\[Phi]2"}]}]], " ", "\[Rule]", 

-         RowBox[{"2", 

-          RowBox[{"dc", "[", "t", "]"}], 

+         RowBox[{"2", "dc", " ", 

           RowBox[{"E2c", "/", "\[CapitalOmega]2"}]}]}], ",", 

         RowBox[{

          SuperscriptBox["\[ExponentialE]", 

@@ -1274,8 +1266,7 @@ Cell[BoxData[{
              RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Delta]"}], "2"]}], 

            "+", 

            RowBox[{"\[ImaginaryI]", " ", "\[Phi]3"}]}]], " ", "\[Rule]", 

-         RowBox[{"2", 

-          RowBox[{"dc", "[", "t", "]"}], 

+         RowBox[{"2", "d", " ", 

           RowBox[{"E3", "/", "\[CapitalOmega]3"}]}]}], ",", 

         RowBox[{

          SuperscriptBox["\[ExponentialE]", 

@@ -1283,8 +1274,7 @@ Cell[BoxData[{
            FractionBox[

             RowBox[{"\[ImaginaryI]", " ", "t", " ", "\[Delta]"}], "2"], "-", 

            RowBox[{"\[ImaginaryI]", " ", "\[Phi]3"}]}]], " ", "\[Rule]", 

-         RowBox[{"2", 

-          RowBox[{"d", "[", "t", "]"}], 

+         RowBox[{"2", "d", " ", 

           RowBox[{"E3c", "/", "\[CapitalOmega]3"}]}]}], ",", 

         RowBox[{"\[CapitalGamma]", "\[Rule]", 

          RowBox[{"2", "G"}]}], ",", 

@@ -1349,7 +1339,8 @@ Cell[BoxData[{
   RowBox[{

    RowBox[{

     RowBox[{"NotebookDirectory", "[", "]"}], "<>", "\"\<\\\\code.txt\>\""}], 

-   ",", "%"}], "]"}]}], "Input"],

+   ",", "%"}], "]"}]}], "Input",

+ CellChangeTimes->{{3.524014168318809*^9, 3.524014208362378*^9}}],

 

 Cell[BoxData[

  RowBox[{"{", 

@@ -1385,13 +1376,9 @@ Cell[BoxData[
         RowBox[{"da", " ", "i"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"1", ",", "2"}]]}], "-", 

-     RowBox[{"i", " ", 

+     RowBox[{"d", " ", 

       RowBox[{"(", 

-       RowBox[{

-        RowBox[{"E2", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E3", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+       RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"1", ",", "3"}]]}], "+", 

      RowBox[{"E1c", " ", "i", " ", 

@@ -1405,13 +1392,21 @@ Cell[BoxData[
       RowBox[{"-", "E1c"}], " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"1", ",", "1"}]]}], "-", 

-     RowBox[{"i", " ", 

+     RowBox[{

+      FractionBox["1", "2"], " ", "i", " ", 

       RowBox[{"(", 

        RowBox[{

-        RowBox[{"E3c", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E2c", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+        RowBox[{"2", " ", "d", " ", "E3c"}], "+", 

+        RowBox[{

+         SuperscriptBox["\[ExponentialE]", 

+          RowBox[{

+           RowBox[{"-", 

+            FractionBox["1", "2"]}], " ", "i", " ", 

+           RowBox[{"(", 

+            RowBox[{

+             RowBox[{"t", " ", "\[Delta]"}], "+", 

+             RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ", 

+         "\[CapitalOmega]2"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"1", ",", "2"}]]}], "+", 

      RowBox[{

@@ -1441,13 +1436,21 @@ Cell[BoxData[
      RowBox[{"E1", " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"2", ",", "3"}]]}], "+", 

-     RowBox[{"i", " ", 

+     RowBox[{

+      FractionBox["1", "2"], " ", "i", " ", 

       RowBox[{"(", 

        RowBox[{

-        RowBox[{"E3c", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E2c", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+        RowBox[{"2", " ", "d", " ", "E3c"}], "+", 

+        RowBox[{

+         SuperscriptBox["\[ExponentialE]", 

+          RowBox[{

+           RowBox[{"-", 

+            FractionBox["1", "2"]}], " ", "i", " ", 

+           RowBox[{"(", 

+            RowBox[{

+             RowBox[{"t", " ", "\[Delta]"}], "+", 

+             RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ", 

+         "\[CapitalOmega]2"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"3", ",", "1"}]]}]}]}], ",", 

    RowBox[{

@@ -1457,22 +1460,26 @@ Cell[BoxData[
      RowBox[{"2", " ", "gt", " ", 

       SubscriptBox["r", 

        RowBox[{"2", ",", "2"}]]}], "-", 

-     RowBox[{"i", " ", 

+     RowBox[{"d", " ", 

       RowBox[{"(", 

-       RowBox[{

-        RowBox[{"E2", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E3", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+       RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"2", ",", "3"}]]}], "+", 

-     RowBox[{"i", " ", 

+     RowBox[{

+      FractionBox["1", "2"], " ", "i", " ", 

       RowBox[{"(", 

        RowBox[{

-        RowBox[{"E3c", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E2c", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+        RowBox[{"2", " ", "d", " ", "E3c"}], "+", 

+        RowBox[{

+         SuperscriptBox["\[ExponentialE]", 

+          RowBox[{

+           RowBox[{"-", 

+            FractionBox["1", "2"]}], " ", "i", " ", 

+           RowBox[{"(", 

+            RowBox[{

+             RowBox[{"t", " ", "\[Delta]"}], "+", 

+             RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ", 

+         "\[CapitalOmega]2"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"3", ",", "2"}]]}], "+", 

      RowBox[{"G", " ", 

@@ -1486,13 +1493,21 @@ Cell[BoxData[
       RowBox[{"-", "E1c"}], " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"2", ",", "1"}]]}], "-", 

-     RowBox[{"i", " ", 

+     RowBox[{

+      FractionBox["1", "2"], " ", "i", " ", 

       RowBox[{"(", 

        RowBox[{

-        RowBox[{"E3c", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E2c", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+        RowBox[{"2", " ", "d", " ", "E3c"}], "+", 

+        RowBox[{

+         SuperscriptBox["\[ExponentialE]", 

+          RowBox[{

+           RowBox[{"-", 

+            FractionBox["1", "2"]}], " ", "i", " ", 

+           RowBox[{"(", 

+            RowBox[{

+             RowBox[{"t", " ", "\[Delta]"}], "+", 

+             RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ", 

+         "\[CapitalOmega]2"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"2", ",", "2"}]]}], "+", 

      RowBox[{

@@ -1503,13 +1518,21 @@ Cell[BoxData[
         RowBox[{"da", " ", "i"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"2", ",", "3"}]]}], "+", 

-     RowBox[{"i", " ", 

+     RowBox[{

+      FractionBox["1", "2"], " ", "i", " ", 

       RowBox[{"(", 

        RowBox[{

-        RowBox[{"E3c", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E2c", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+        RowBox[{"2", " ", "d", " ", "E3c"}], "+", 

+        RowBox[{

+         SuperscriptBox["\[ExponentialE]", 

+          RowBox[{

+           RowBox[{"-", 

+            FractionBox["1", "2"]}], " ", "i", " ", 

+           RowBox[{"(", 

+            RowBox[{

+             RowBox[{"t", " ", "\[Delta]"}], "+", 

+             RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ", 

+         "\[CapitalOmega]2"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"3", ",", "3"}]]}]}]}], ",", 

    RowBox[{

@@ -1519,13 +1542,9 @@ Cell[BoxData[
      RowBox[{"E1", " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"1", ",", "1"}]]}], "+", 

-     RowBox[{"i", " ", 

+     RowBox[{"d", " ", 

       RowBox[{"(", 

-       RowBox[{

-        RowBox[{"E2", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E3", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+       RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"2", ",", "1"}]]}], "+", 

      RowBox[{

@@ -1546,13 +1565,9 @@ Cell[BoxData[
      RowBox[{"E1", " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"1", ",", "2"}]]}], "+", 

-     RowBox[{"i", " ", 

+     RowBox[{"d", " ", 

       RowBox[{"(", 

-       RowBox[{

-        RowBox[{"E2", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E3", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+       RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"2", ",", "2"}]]}], "+", 

      RowBox[{

@@ -1563,13 +1578,9 @@ Cell[BoxData[
         RowBox[{"da", " ", "i"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"3", ",", "2"}]]}], "-", 

-     RowBox[{"i", " ", 

+     RowBox[{"d", " ", 

       RowBox[{"(", 

-       RowBox[{

-        RowBox[{"E2", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E3", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+       RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"3", ",", "3"}]]}]}]}], ",", 

    RowBox[{

@@ -1582,32 +1593,38 @@ Cell[BoxData[
      RowBox[{"E1", " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"1", ",", "3"}]]}], "+", 

-     RowBox[{"i", " ", 

+     RowBox[{"d", " ", 

       RowBox[{"(", 

-       RowBox[{

-        RowBox[{"E2", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E3", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+       RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"2", ",", "3"}]]}], "-", 

      RowBox[{"E1c", " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"3", ",", "1"}]]}], "-", 

-     RowBox[{"i", " ", 

+     RowBox[{

+      FractionBox["1", "2"], " ", "i", " ", 

       RowBox[{"(", 

        RowBox[{

-        RowBox[{"E3c", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E2c", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+        RowBox[{"2", " ", "d", " ", "E3c"}], "+", 

+        RowBox[{

+         SuperscriptBox["\[ExponentialE]", 

+          RowBox[{

+           RowBox[{"-", 

+            FractionBox["1", "2"]}], " ", "i", " ", 

+           RowBox[{"(", 

+            RowBox[{

+             RowBox[{"t", " ", "\[Delta]"}], "+", 

+             RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ", 

+         "\[CapitalOmega]2"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"3", ",", "2"}]]}], "-", 

      RowBox[{"2", " ", 

       RowBox[{"(", 

        RowBox[{"G", "+", "gt"}], ")"}], " ", 

       SubscriptBox["r", 

-       RowBox[{"3", ",", "3"}]]}]}]}]}], "}"}]], "Output"],

+       RowBox[{"3", ",", "3"}]]}]}]}]}], "}"}]], "Output",

+ CellChangeTimes->{3.524014014687035*^9, 3.524014182827821*^9, 

+  3.524014214903372*^9}],

 

 Cell[BoxData[

  RowBox[{"{", 

@@ -1643,13 +1660,9 @@ Cell[BoxData[
         RowBox[{"da", " ", "i"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"1", ",", "2"}]]}], "-", 

-     RowBox[{"i", " ", 

+     RowBox[{"d", " ", 

       RowBox[{"(", 

-       RowBox[{

-        RowBox[{"E2", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E3", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+       RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"1", ",", "3"}]]}], "+", 

      RowBox[{"E1c", " ", "i", " ", 

@@ -1663,13 +1676,21 @@ Cell[BoxData[
       RowBox[{"-", "E1c"}], " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"1", ",", "1"}]]}], "-", 

-     RowBox[{"i", " ", 

+     RowBox[{

+      FractionBox["1", "2"], " ", "i", " ", 

       RowBox[{"(", 

        RowBox[{

-        RowBox[{"E3c", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E2c", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+        RowBox[{"2", " ", "d", " ", "E3c"}], "+", 

+        RowBox[{

+         SuperscriptBox["\[ExponentialE]", 

+          RowBox[{

+           RowBox[{"-", 

+            FractionBox["1", "2"]}], " ", "i", " ", 

+           RowBox[{"(", 

+            RowBox[{

+             RowBox[{"t", " ", "\[Delta]"}], "+", 

+             RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ", 

+         "\[CapitalOmega]2"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"1", ",", "2"}]]}], "+", 

      RowBox[{

@@ -1690,22 +1711,26 @@ Cell[BoxData[
      RowBox[{"2", " ", "gt", " ", 

       SubscriptBox["r", 

        RowBox[{"2", ",", "2"}]]}], "-", 

-     RowBox[{"i", " ", 

+     RowBox[{"d", " ", 

       RowBox[{"(", 

-       RowBox[{

-        RowBox[{"E2", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E3", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+       RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"2", ",", "3"}]]}], "+", 

-     RowBox[{"i", " ", 

+     RowBox[{

+      FractionBox["1", "2"], " ", "i", " ", 

       RowBox[{"(", 

        RowBox[{

-        RowBox[{"E3c", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E2c", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+        RowBox[{"2", " ", "d", " ", "E3c"}], "+", 

+        RowBox[{

+         SuperscriptBox["\[ExponentialE]", 

+          RowBox[{

+           RowBox[{"-", 

+            FractionBox["1", "2"]}], " ", "i", " ", 

+           RowBox[{"(", 

+            RowBox[{

+             RowBox[{"t", " ", "\[Delta]"}], "+", 

+             RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ", 

+         "\[CapitalOmega]2"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"3", ",", "2"}]]}], "+", 

      RowBox[{"G", " ", 

@@ -1719,13 +1744,21 @@ Cell[BoxData[
       RowBox[{"-", "E1c"}], " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"2", ",", "1"}]]}], "-", 

-     RowBox[{"i", " ", 

+     RowBox[{

+      FractionBox["1", "2"], " ", "i", " ", 

       RowBox[{"(", 

        RowBox[{

-        RowBox[{"E3c", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E2c", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+        RowBox[{"2", " ", "d", " ", "E3c"}], "+", 

+        RowBox[{

+         SuperscriptBox["\[ExponentialE]", 

+          RowBox[{

+           RowBox[{"-", 

+            FractionBox["1", "2"]}], " ", "i", " ", 

+           RowBox[{"(", 

+            RowBox[{

+             RowBox[{"t", " ", "\[Delta]"}], "+", 

+             RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ", 

+         "\[CapitalOmega]2"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"2", ",", "2"}]]}], "+", 

      RowBox[{

@@ -1736,13 +1769,21 @@ Cell[BoxData[
         RowBox[{"da", " ", "i"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"2", ",", "3"}]]}], "+", 

-     RowBox[{"i", " ", 

+     RowBox[{

+      FractionBox["1", "2"], " ", "i", " ", 

       RowBox[{"(", 

        RowBox[{

-        RowBox[{"E3c", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E2c", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+        RowBox[{"2", " ", "d", " ", "E3c"}], "+", 

+        RowBox[{

+         SuperscriptBox["\[ExponentialE]", 

+          RowBox[{

+           RowBox[{"-", 

+            FractionBox["1", "2"]}], " ", "i", " ", 

+           RowBox[{"(", 

+            RowBox[{

+             RowBox[{"t", " ", "\[Delta]"}], "+", 

+             RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ", 

+         "\[CapitalOmega]2"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"3", ",", "3"}]]}]}]}], ",", 

    RowBox[{

@@ -1755,52 +1796,65 @@ Cell[BoxData[
      RowBox[{"E1", " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"1", ",", "3"}]]}], "+", 

-     RowBox[{"i", " ", 

+     RowBox[{"d", " ", 

       RowBox[{"(", 

-       RowBox[{

-        RowBox[{"E2", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E3", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+       RowBox[{"E2", "+", "E3"}], ")"}], " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"2", ",", "3"}]]}], "-", 

      RowBox[{"E1c", " ", "i", " ", 

       SubscriptBox["r", 

        RowBox[{"3", ",", "1"}]]}], "-", 

-     RowBox[{"i", " ", 

+     RowBox[{

+      FractionBox["1", "2"], " ", "i", " ", 

       RowBox[{"(", 

        RowBox[{

-        RowBox[{"E3c", " ", 

-         RowBox[{"d", "[", "t", "]"}]}], "+", 

-        RowBox[{"E2c", " ", 

-         RowBox[{"dc", "[", "t", "]"}]}]}], ")"}], " ", 

+        RowBox[{"2", " ", "d", " ", "E3c"}], "+", 

+        RowBox[{

+         SuperscriptBox["\[ExponentialE]", 

+          RowBox[{

+           RowBox[{"-", 

+            FractionBox["1", "2"]}], " ", "i", " ", 

+           RowBox[{"(", 

+            RowBox[{

+             RowBox[{"t", " ", "\[Delta]"}], "+", 

+             RowBox[{"2", " ", "\[Phi]2"}]}], ")"}]}]], " ", 

+         "\[CapitalOmega]2"}]}], ")"}], " ", 

       SubscriptBox["r", 

        RowBox[{"3", ",", "2"}]]}], "-", 

      RowBox[{"2", " ", 

       RowBox[{"(", 

        RowBox[{"G", "+", "gt"}], ")"}], " ", 

       SubscriptBox["r", 

-       RowBox[{"3", ",", "3"}]]}]}]}]}], "}"}]], "Output"],

+       RowBox[{"3", ",", "3"}]]}]}]}]}], "}"}]], "Output",

+ CellChangeTimes->{3.524014014687035*^9, 3.524014182827821*^9, 

+  3.524014214922399*^9}],

 

 Cell[BoxData["\<\"dr11_dt = gt - 2*(gp + gt)*r11 - E1*i*r13 + E1c*i*r31 + \

-G*r33;\\ndr12_dt = (-gp - 2*gt - d1*i + da*i)*r12 - i*(E2*d(t) + \

-E3*dc(t))*r13 + E1c*i*r32;\\ndr13_dt = -(E1c*i*r11) - i*(E3c*d(t) + \

-E2c*dc(t))*r12 + (-G - gp - 2*gt - d1*i)*r13 + E1c*i*r33;\\ndr22_dt = gt - \

-2*gt*r22 - i*(E2*d(t) + E3*dc(t))*r23 + i*(E3c*d(t) + E2c*dc(t))*r32 + G*r33;\

-\\ndr23_dt = -(E1c*i*r21) - i*(E3c*d(t) + E2c*dc(t))*r22 + (-G - 2*gt - \

-da*i)*r23 + i*(E3c*d(t) + E2c*dc(t))*r33;\\ndr33_dt = 2*gp*r11 + E1*i*r13 + \

-i*(E2*d(t) + E3*dc(t))*r23 - E1c*i*r31 - i*(E3c*d(t) + E2c*dc(t))*r32 - 2*(G \

-+ gt)*r33;\\n\"\>"], "Output"],

+G*r33;\\ndr12_dt = (-gp - 2*gt - d1*i + da*i)*r12 - d*(E2 + E3)*i*r13 + \

+E1c*i*r32;\\ndr13_dt = -(E1c*i*r11) - (i*(2*d*E3c + \

+\[CapitalOmega]2/Power(E,(i*(t*\[Delta] + 2*\[Phi]2))/2.))*r12)/2. + (-G - gp \

+- 2*gt - d1*i)*r13 + E1c*i*r33;\\ndr22_dt = gt - 2*gt*r22 - d*(E2 + E3)*i*r23 \

++ (i*(2*d*E3c + \[CapitalOmega]2/Power(E,(i*(t*\[Delta] + \

+2*\[Phi]2))/2.))*r32)/2. + G*r33;\\ndr23_dt = -(E1c*i*r21) - (i*(2*d*E3c + \

+\[CapitalOmega]2/Power(E,(i*(t*\[Delta] + 2*\[Phi]2))/2.))*r22)/2. + (-G - \

+2*gt - da*i)*r23 + (i*(2*d*E3c + \[CapitalOmega]2/Power(E,(i*(t*\[Delta] + 2*\

+\[Phi]2))/2.))*r33)/2.;\\ndr33_dt = 2*gp*r11 + E1*i*r13 + d*(E2 + E3)*i*r23 - \

+E1c*i*r31 - (i*(2*d*E3c + \[CapitalOmega]2/Power(E,(i*(t*\[Delta] + \

+2*\[Phi]2))/2.))*r32)/2. - 2*(G + gt)*r33;\\n\"\>"], "Output",

+ CellChangeTimes->{3.524014014687035*^9, 3.524014182827821*^9, 

+  3.524014214936878*^9}],

 

-Cell[BoxData["\<\"C:\\\\cygwin\\\\home\\\\Simon\\\\Nresonances\\\\xmds2\\\\\

-Shahriar_system\\\\\\\\code.txt\"\>"], "Output"]

+Cell[BoxData["\<\"/mnt/light_huge_archive/home/evmik/src/my_src/Nresonances/\

+xmds2/Shahriar_system/\\\\code.txt\"\>"], "Output",

+ CellChangeTimes->{3.524014014687035*^9, 3.524014182827821*^9, 

+  3.52401421495504*^9}]

 }, Open  ]]

 }, Open  ]]

 },

-WindowSize->{1010, 875},

-WindowMargins->{{0, Automatic}, {Automatic, 0}},

+WindowSize->{960, 501},

+WindowMargins->{{Automatic, 0}, {24, Automatic}},

 ShowSelection->True,

-FrontEndVersion->"7.0 for Microsoft Windows (64-bit) (February 18, 2009)",

+FrontEndVersion->"8.0 for Linux x86 (32-bit) (February 23, 2011)",

 StyleDefinitions->"Default.nb"

 ]

 (* End of Notebook Content *)

@@ -1815,70 +1869,70 @@ CellTagsIndex->{}
 (*NotebookFileOutline

 Notebook[{

 Cell[CellGroupData[{

-Cell[567, 22, 36, 0, 71, "Section"],

-Cell[606, 24, 66, 1, 43, "MathCaption",

+Cell[567, 22, 36, 0, 74, "Section"],

+Cell[606, 24, 66, 1, 44, "MathCaption",

  CellID->836781195],

-Cell[675, 27, 85, 2, 31, "Input",

+Cell[675, 27, 85, 2, 30, "Input",

  CellID->2058623809],

-Cell[763, 31, 1290, 44, 47, "Text",

+Cell[763, 31, 1128, 37, 52, "Text",

  CellID->525777075],

-Cell[2056, 77, 112, 3, 43, "MathCaption",

+Cell[1894, 70, 112, 3, 44, "MathCaption",

  CellID->429217524],

-Cell[2171, 82, 1057, 27, 112, "Input",

+Cell[2009, 75, 1057, 27, 126, "Input",

  CellID->433132487],

-Cell[3231, 111, 125, 3, 43, "MathCaption",

+Cell[3069, 104, 125, 3, 44, "MathCaption",

  CellID->133602844],

 Cell[CellGroupData[{

-Cell[3381, 118, 1085, 32, 52, "Input"],

-Cell[4469, 152, 1306, 40, 53, "Output"]

+Cell[3219, 111, 1085, 32, 69, "Input"],

+Cell[4307, 145, 1348, 41, 52, "Output"]

 }, Open  ]],

-Cell[5790, 195, 166, 4, 43, "MathCaption",

+Cell[5670, 189, 166, 4, 44, "MathCaption",

  CellID->462076121],

 Cell[CellGroupData[{

-Cell[5981, 203, 508, 15, 31, "Input"],

-Cell[6492, 220, 1755, 50, 71, "Output"]

+Cell[5861, 197, 508, 15, 50, "Input"],

+Cell[6372, 214, 1797, 51, 73, "Output"]

 }, Open  ]],

-Cell[8262, 273, 76, 1, 43, "MathCaption",

+Cell[8184, 268, 76, 1, 44, "MathCaption",

  CellID->358620443],

 Cell[CellGroupData[{

-Cell[8363, 278, 464, 15, 31, "Input",

+Cell[8285, 273, 530, 16, 30, "Input",

  CellID->167259034],

-Cell[8830, 295, 535, 10, 94, "Output"]

+Cell[8818, 291, 585, 11, 94, "Output"]

 }, Open  ]],

-Cell[9380, 308, 182, 4, 43, "MathCaption",

+Cell[9418, 305, 182, 4, 44, "MathCaption",

  CellID->577766068],

 Cell[CellGroupData[{

-Cell[9587, 316, 1379, 38, 72, "Input"],

-Cell[10969, 356, 2462, 70, 117, "Output"]

+Cell[9625, 313, 1379, 38, 69, "Input"],

+Cell[11007, 353, 2503, 71, 112, "Output"]

 }, Open  ]],

-Cell[13446, 429, 456, 13, 43, "MathCaption",

+Cell[13525, 427, 456, 13, 44, "MathCaption",

  CellID->610306692],

 Cell[CellGroupData[{

-Cell[13927, 446, 933, 25, 57, "Input",

+Cell[14006, 444, 933, 25, 58, "Input",

  CellID->645617687],

-Cell[14863, 473, 717, 20, 71, "Output"]

+Cell[14942, 471, 759, 21, 72, "Output"]

 }, Open  ]],

-Cell[15595, 496, 390, 12, 43, "MathCaption",

+Cell[15716, 495, 390, 12, 44, "MathCaption",

  CellID->854192725],

 Cell[CellGroupData[{

-Cell[16010, 512, 1070, 29, 59, "Input",

+Cell[16131, 511, 1070, 29, 58, "Input",

  CellID->465762594],

-Cell[17083, 543, 1226, 36, 92, "Output"]

+Cell[17204, 542, 1268, 37, 91, "Output"]

 }, Open  ]],

-Cell[18324, 582, 76, 1, 43, "MathCaption",

+Cell[18487, 582, 76, 1, 44, "MathCaption",

  CellID->314466782],

 Cell[CellGroupData[{

-Cell[18425, 587, 250, 7, 31, "Input",

+Cell[18588, 587, 250, 7, 30, "Input",

  CellID->298399236],

-Cell[18678, 596, 23733, 623, 315, "Output"]

+Cell[18841, 596, 23487, 615, 298, "Output"]

 }, Open  ]],

-Cell[42426, 1222, 34, 0, 29, "Text"],

+Cell[42343, 1214, 34, 0, 30, "Text"],

 Cell[CellGroupData[{

-Cell[42485, 1226, 4646, 125, 233, "Input"],

-Cell[47134, 1353, 7898, 256, 164, "Output"],

-Cell[55035, 1611, 5277, 171, 107, "Output"],

-Cell[60315, 1784, 630, 8, 145, "Output"],

-Cell[60948, 1794, 125, 1, 30, "Output"]

+Cell[42402, 1218, 4657, 124, 255, "Input"],

+Cell[47062, 1344, 8995, 282, 278, "Output"],

+Cell[56060, 1628, 6444, 201, 205, "Output"],

+Cell[62507, 1831, 958, 13, 210, "Output"],

+Cell[63468, 1846, 217, 3, 29, "Output"]

 }, Open  ]]

 }, Open  ]]

 }