path: root/mathemathica_fwm/fwm_degenerate.nb

(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 7.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[    102785,       2526]
NotebookOptionsPosition[     99915,       2423]
NotebookOutlinePosition[    100292,       2440]
CellTagsIndexPosition[    100249,       2437]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["General setup", "Section"],

Cell["This loads the package.", "MathCaption",
 CellID->836781195],

Cell[BoxData[
 RowBox[{"<<", "AtomicDensityMatrix`"}]], "Input",
 CellChangeTimes->{{3.522532598595615*^9, 3.522532603186735*^9}},
 CellID->2058623809],

Cell[TextData[{
 "We define an atomic system consisting of two even-parity lower states and \
two odd-parity upper states. We apply a light field with components at \
frequencies ",
 Cell[BoxData[
  FormBox[
   StyleBox[
    SubscriptBox["\[Omega]", "1"], "InlineMath"], TraditionalForm]]],
 " (near resonant with the ",
 Cell[BoxData[
  StyleBox[
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"|", "1"}], "\[RightAngleBracket]"}], "\[Rule]", 
     RowBox[{"|", "3"}]}], "\[RightAngleBracket]"}], "InlineMath"]]],
 " transition), ",
 Cell[BoxData[
  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),  ",
 Cell[BoxData[
  FormBox[
   StyleBox[
    SubscriptBox["\[Omega]", "3"], "InlineMath"], TraditionalForm]]],
 " (near resonant with the ",
 Cell[BoxData[
  StyleBox[
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"|", "2"}], "\[RightAngleBracket]"}], "\[Rule]", 
     RowBox[{"|", "4"}]}], "\[RightAngleBracket]"}], "InlineMath"]]],
 " transition), and ",
 Cell[BoxData[
  FormBox[
   StyleBox[
    SubscriptBox["\[Omega]", "4"], "InlineMath"], TraditionalForm]]],
 " (near resonant with the ",
 Cell[BoxData[
  StyleBox[
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"|", "1"}], "\[RightAngleBracket]"}], "\[Rule]", 
     RowBox[{"|", "4"}]}], "\[RightAngleBracket]"}], "InlineMath"]]],
 " transition)"
}], "Text",
 CellChangeTimes->{{3.522532429821333*^9, 3.522532492302448*^9}, {
  3.522540911043191*^9, 3.522540911147507*^9}},
 CellID->525777075],

Cell["Define the atomic system.", "MathCaption",
 CellID->429217524],

Cell[BoxData[
 RowBox[{
  RowBox[{"system", "=", 
   RowBox[{"{", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"AtomicState", "[", 
      RowBox[{"1", ",", 
       RowBox[{"NaturalWidth", "\[Rule]", "0"}], ",", 
       RowBox[{"Parity", "\[Rule]", "Even"}]}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"AtomicState", "[", 
      RowBox[{"2", ",", 
       RowBox[{"NaturalWidth", "\[Rule]", "0"}], ",", 
       RowBox[{"Parity", "\[Rule]", "Even"}]}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"AtomicState", "[", 
      RowBox[{"3", ",", 
       RowBox[{"Energy", "\[Rule]", "0"}], ",", 
       RowBox[{"NaturalWidth", "\[Rule]", 
        SubscriptBox["\[CapitalGamma]", "3"]}], ",", 
       RowBox[{"Parity", "\[Rule]", "Odd"}]}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"AtomicState", "[", 
      RowBox[{"4", ",", 
       RowBox[{"NaturalWidth", "\[Rule]", 
        SubscriptBox["\[CapitalGamma]", "4"]}], ",", 
       RowBox[{"Parity", "\[Rule]", "Odd"}]}], "]"}]}], "\[IndentingNewLine]",
     "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.522536311483895*^9, 3.52253631268178*^9}},
 CellID->433132487],

Cell[TextData[{
 "Define the optical field with three frequencies, ",
 Cell[BoxData[
  FormBox[
   StyleBox[
    SubscriptBox["\[Omega]", "1"], "InlineMath"], TraditionalForm]]],
 ", ",
 Cell[BoxData[
  FormBox[
   StyleBox[
    SubscriptBox["\[Omega]", "2"], "InlineMath"], TraditionalForm]]],
 ", ",
 Cell[BoxData[
  FormBox[
   StyleBox[
    SubscriptBox["\[Omega]", "3"], "InlineMath"], TraditionalForm]]],
 ", and ",
 Cell[BoxData[
  FormBox[
   StyleBox[
    SubscriptBox["\[Omega]", "4"], "InlineMath"], TraditionalForm]]],
 "."
}], "MathCaption",
 CellChangeTimes->{{3.522540939999547*^9, 3.522540967294499*^9}, {
  3.522541679413973*^9, 3.522541681294852*^9}},
 CellID->133602844],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"field", "=", 
  RowBox[{
   RowBox[{"OpticalField", "[", 
    RowBox[{"\[Omega]1", ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"\[CapitalOmega]1", "/", 
        RowBox[{"ReducedME", "[", 
         RowBox[{"1", ",", 
          RowBox[{"{", 
           RowBox[{"Dipole", ",", "1"}], "}"}], ",", "3"}], "]"}]}], ",", 
       "\[Phi]1"}], "}"}]}], "]"}], "+", 
   RowBox[{"OpticalField", "[", 
    RowBox[{"\[Omega]2", ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"\[CapitalOmega]2", "/", 
        RowBox[{"ReducedME", "[", 
         RowBox[{"2", ",", 
          RowBox[{"{", 
           RowBox[{"Dipole", ",", "1"}], "}"}], ",", "3"}], "]"}]}], ",", 
       "\[Phi]2"}], "}"}]}], "]"}], "+", 
   RowBox[{"OpticalField", "[", 
    RowBox[{"\[Omega]3", ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"\[CapitalOmega]3", "/", 
        RowBox[{"ReducedME", "[", 
         RowBox[{"2", ",", 
          RowBox[{"{", 
           RowBox[{"Dipole", ",", "1"}], "}"}], ",", "4"}], "]"}]}], ",", 
       "\[Phi]3"}], "}"}]}], "]"}], "+", "\[IndentingNewLine]", 
   RowBox[{"OpticalField", "[", 
    RowBox[{"\[Omega]4", ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"\[CapitalOmega]4", "/", 
        RowBox[{"ReducedME", "[", 
         RowBox[{"1", ",", 
          RowBox[{"{", 
           RowBox[{"Dipole", ",", "1"}], "}"}], ",", "4"}], "]"}]}], ",", 
       "\[Phi]4"}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.522534194542496*^9, 3.522534219533982*^9}, {
  3.522540920436486*^9, 3.522540924714358*^9}},
 CellID->534530029],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    FractionBox[
     RowBox[{
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{"\[ImaginaryI]", " ", 
        RowBox[{"(", 
         RowBox[{"\[Phi]1", "-", 
          RowBox[{"t", " ", "\[Omega]1"}]}], ")"}]}]], " ", 
      "\[CapitalOmega]1"}], 
     RowBox[{"ReducedME", "[", 
      RowBox[{"1", ",", 
       RowBox[{"{", 
        RowBox[{"Dipole", ",", "1"}], "}"}], ",", "3"}], "]"}]], "+", 
    FractionBox[
     RowBox[{
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{"\[ImaginaryI]", " ", 
        RowBox[{"(", 
         RowBox[{"\[Phi]4", "-", 
          RowBox[{"t", " ", "\[Omega]4"}]}], ")"}]}]], " ", 
      "\[CapitalOmega]4"}], 
     RowBox[{"ReducedME", "[", 
      RowBox[{"1", ",", 
       RowBox[{"{", 
        RowBox[{"Dipole", ",", "1"}], "}"}], ",", "4"}], "]"}]], "+", 
    FractionBox[
     RowBox[{
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{"\[ImaginaryI]", " ", 
        RowBox[{"(", 
         RowBox[{"\[Phi]2", "-", 
          RowBox[{"t", " ", "\[Omega]2"}]}], ")"}]}]], " ", 
      "\[CapitalOmega]2"}], 
     RowBox[{"ReducedME", "[", 
      RowBox[{"2", ",", 
       RowBox[{"{", 
        RowBox[{"Dipole", ",", "1"}], "}"}], ",", "3"}], "]"}]], "+", 
    FractionBox[
     RowBox[{
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{"\[ImaginaryI]", " ", 
        RowBox[{"(", 
         RowBox[{"\[Phi]3", "-", 
          RowBox[{"t", " ", "\[Omega]3"}]}], ")"}]}]], " ", 
      "\[CapitalOmega]3"}], 
     RowBox[{"ReducedME", "[", 
      RowBox[{"2", ",", 
       RowBox[{"{", 
        RowBox[{"Dipole", ",", "1"}], "}"}], ",", "4"}], "]"}]]}], ",", "0", 
   ",", "0"}], "}"}]], "Output",
 CellChangeTimes->{3.522533257349164*^9, 3.52253422045873*^9, 
  3.522536315439293*^9, 3.522540739248025*^9, 3.522540970145529*^9, 
  3.522541564921328*^9, 3.522542172549766*^9, 3.522542936599973*^9, 
  3.522544786198204*^9, 3.522544854264734*^9, 3.52254522556669*^9, 
  3.522545943915771*^9}]
}, Open  ]],

Cell["\<\
The Hamiltonian for the system subject to the optical field. Each field is \
assumed to interact with only one transition\[LongDash]the other terms are \
set to zero.\
\>", "MathCaption",
 CellID->462076121],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"MatrixForm", "[", 
  RowBox[{"H", "=", 
   RowBox[{
    RowBox[{"Expand", "@", 
     RowBox[{"Hamiltonian", "[", 
      RowBox[{"system", ",", 
       RowBox[{"ElectricField", "\[Rule]", "field"}]}], "]"}]}], "/.", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"Cos", "[", "_", "]"}], " ", 
      RowBox[{"ReducedME", "[", 
       RowBox[{"_", ",", 
        RowBox[{"{", 
         RowBox[{"Dipole", ",", "1"}], "}"}], ",", "_"}], "]"}]}], "\[Rule]", 
     "0"}]}]}], "]"}]], "Input",
 CellID->494599775],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"Energy", "[", "1", "]"}], "0", 
      RowBox[{
       RowBox[{"-", "\[CapitalOmega]1"}], " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Phi]1", "-", 
         RowBox[{"t", " ", "\[Omega]1"}]}], "]"}]}], 
      RowBox[{
       RowBox[{"-", "\[CapitalOmega]4"}], " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Phi]4", "-", 
         RowBox[{"t", " ", "\[Omega]4"}]}], "]"}]}]},
     {"0", 
      RowBox[{"Energy", "[", "2", "]"}], 
      RowBox[{
       RowBox[{"-", "\[CapitalOmega]2"}], " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Phi]2", "-", 
         RowBox[{"t", " ", "\[Omega]2"}]}], "]"}]}], 
      RowBox[{
       RowBox[{"-", "\[CapitalOmega]3"}], " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Phi]3", "-", 
         RowBox[{"t", " ", "\[Omega]3"}]}], "]"}]}]},
     {
      RowBox[{
       RowBox[{"-", "\[CapitalOmega]1"}], " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Phi]1", "-", 
         RowBox[{"t", " ", "\[Omega]1"}]}], "]"}]}], 
      RowBox[{
       RowBox[{"-", "\[CapitalOmega]2"}], " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Phi]2", "-", 
         RowBox[{"t", " ", "\[Omega]2"}]}], "]"}]}], "0", "0"},
     {
      RowBox[{
       RowBox[{"-", "\[CapitalOmega]4"}], " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Phi]4", "-", 
         RowBox[{"t", " ", "\[Omega]4"}]}], "]"}]}], 
      RowBox[{
       RowBox[{"-", "\[CapitalOmega]3"}], " ", 
       RowBox[{"Cos", "[", 
        RowBox[{"\[Phi]3", "-", 
         RowBox[{"t", " ", "\[Omega]3"}]}], "]"}]}], "0", 
      RowBox[{"Energy", "[", "4", "]"}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.522533260351258*^9, 3.522534225599567*^9, 
  3.522536320272253*^9, 3.522540741905562*^9, 3.522540981472922*^9, 
  3.522541568259588*^9, 3.522542173099244*^9, 3.522542936745729*^9, 
  3.522544786457194*^9, 3.522544854514756*^9, 3.522545225737479*^9, 
  3.522545944334553*^9}]
}, Open  ]],

Cell["The level diagram for the system.", "MathCaption",
 CellID->358620443],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"LevelDiagram", "[", 
  RowBox[{"system", ",", 
   RowBox[{"H", "/.", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"Energy", "[", "1", "]"}], "\[Rule]", 
       RowBox[{"-", "1.5"}]}], ",", 
      RowBox[{
       RowBox[{"Energy", "[", "2", "]"}], "\[Rule]", 
       RowBox[{"-", "1"}]}], ",", 
      RowBox[{
       RowBox[{"Energy", "[", "4", "]"}], "\[Rule]", ".5"}]}], "}"}]}]}], 
  "]"}]], "Input",
 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.09999999999999998, 0.5}, {0.9, 0.5}}]}}, {{}, {}, {}}, 
   {Arrowheads[{-0.07659574468085106, 0.07659574468085106}], 
    ArrowBox[{{-0.45999999999999996`, -1.5}, {0.45999999999999996`, 0.}}], 
    ArrowBox[{{-0.45999999999999996`, -1.5}, {0.45999999999999996`, 0.5}}], 
    ArrowBox[{{-0.45999999999999996`, -1.}, {0.45999999999999996`, 0.}}], 
    ArrowBox[{{-0.45999999999999996`, -1.}, {0.45999999999999996`, 0.5}}]}, 
   {PointSize[0.0225]}},
  ImagePadding->{{2., 2}, {2., 2.}},
  ImageSize->94.]], "Output",
 CellChangeTimes->{3.522533274929177*^9, 3.522534231598043*^9, 
  3.522536287987677*^9, 3.522536326037291*^9, 3.522540747140555*^9, 
  3.522540990177467*^9, 3.522541571244472*^9, 3.522542173445689*^9, 
  3.522542937562187*^9, 3.522544786610477*^9, 3.522544854751983*^9, 
  3.522545225868822*^9, 3.522545945164071*^9}]
}, Open  ]],

Cell["Apply the rotating-wave approximation to the Hamiltonian.", \
"MathCaption",
 CellID->577766068],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{"Hrwa", "=", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"RotatingWaveApproximation", "[", 
        RowBox[{"system", ",", 
         RowBox[{"H", "/.", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"\[Omega]3", "\[Rule]", 
             RowBox[{"\[Omega]2", "+", "\[Omega]43"}]}], ",", 
            RowBox[{"\[Omega]4", "\[Rule]", 
             RowBox[{"\[Omega]1", "+", "\[Omega]41"}]}]}], "}"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"\[Omega]1", ",", "\[Omega]2", ",", "\[Omega]43"}], "}"}], 
         ",", "\[IndentingNewLine]", 
         RowBox[{"TransformMatrix", "\[Rule]", 
          RowBox[{"MatrixExp", "[", 
           RowBox[{
            RowBox[{"-", "\[ImaginaryI]"}], " ", "t", " ", 
            RowBox[{"DiagonalMatrix", "[", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"-", "\[Omega]1"}], ",", 
               RowBox[{"-", "\[Omega]2"}], ",", "0", ",", "\[Omega]43"}], 
              "}"}], "]"}]}], "]"}]}]}], "]"}], "/.", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"\[Omega]1", "\[Rule]", 
          RowBox[{
           RowBox[{"-", 
            RowBox[{"Energy", "[", "1", "]"}]}], "+", "\[Delta]1"}]}], ",", 
         RowBox[{"\[Omega]2", "\[Rule]", 
          RowBox[{
           RowBox[{"-", 
            RowBox[{"Energy", "[", "2", "]"}]}], "+", "\[Delta]2"}]}], ",", 
         RowBox[{"\[Omega]43", "\[Rule]", 
          RowBox[{
           RowBox[{"Energy", "[", "4", "]"}], "+", "\[Delta]3", "-", 
           "\[Delta]2"}]}], ",", 
         RowBox[{"\[Omega]41", "\[Rule]", 
          RowBox[{
           RowBox[{"Energy", "[", "4", "]"}], "+", "\[Delta]4", "-", 
           "\[Delta]1"}]}]}], "}"}]}]}], ")"}], "//", "MatrixForm"}], " ", "//",
    "Simplify"}], " ", "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.522537910365833*^9, 3.522537961926054*^9}, 
   3.522538015680879*^9, {3.522538077788326*^9, 3.522538108632129*^9}, {
   3.522538220981776*^9, 3.522538294286786*^9}, {3.522538364562097*^9, 
   3.522538365489798*^9}, 3.522538408544395*^9, {3.522538642590504*^9, 
   3.522538692268066*^9}, {3.522538753179414*^9, 3.522538789067889*^9}, 
   3.522538841634358*^9, {3.522539728181287*^9, 3.52253981610592*^9}, {
   3.522539877916903*^9, 3.522539880715579*^9}, {3.522540013440722*^9, 
   3.52254005705363*^9}, {3.522540998041653*^9, 3.522541007397957*^9}, {
   3.522541371243734*^9, 3.522541473321752*^9}, {3.522541586733216*^9, 
   3.52254160084184*^9}, {3.522542443165734*^9, 3.522542553763046*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"\[Delta]1", "0", 
      RowBox[{
       RowBox[{"-", 
        FractionBox["1", "2"]}], " ", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{
         RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]1"}]], " ", 
       "\[CapitalOmega]1"}], 
      RowBox[{
       RowBox[{"-", 
        FractionBox["1", "2"]}], " ", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{
         RowBox[{"-", "\[ImaginaryI]"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"t", " ", 
            RowBox[{"(", 
             RowBox[{
             "\[Delta]1", "-", "\[Delta]2", "+", "\[Delta]3", "-", 
              "\[Delta]4"}], ")"}]}], "+", "\[Phi]4"}], ")"}]}]], " ", 
       "\[CapitalOmega]4"}]},
     {"0", "\[Delta]2", 
      RowBox[{
       RowBox[{"-", 
        FractionBox["1", "2"]}], " ", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{
         RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]2"}]], " ", 
       "\[CapitalOmega]2"}], 
      RowBox[{
       RowBox[{"-", 
        FractionBox["1", "2"]}], " ", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{
         RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]3"}]], " ", 
       "\[CapitalOmega]3"}]},
     {
      RowBox[{
       RowBox[{"-", 
        FractionBox["1", "2"]}], " ", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{"\[ImaginaryI]", " ", "\[Phi]1"}]], " ", "\[CapitalOmega]1"}], 
      RowBox[{
       RowBox[{"-", 
        FractionBox["1", "2"]}], " ", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{"\[ImaginaryI]", " ", "\[Phi]2"}]], " ", "\[CapitalOmega]2"}],
       "0", "0"},
     {
      RowBox[{
       RowBox[{"-", 
        FractionBox["1", "2"]}], " ", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{"\[ImaginaryI]", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"t", " ", 
            RowBox[{"(", 
             RowBox[{
             "\[Delta]1", "-", "\[Delta]2", "+", "\[Delta]3", "-", 
              "\[Delta]4"}], ")"}]}], "+", "\[Phi]4"}], ")"}]}]], " ", 
       "\[CapitalOmega]4"}], 
      RowBox[{
       RowBox[{"-", 
        FractionBox["1", "2"]}], " ", 
       SuperscriptBox["\[ExponentialE]", 
        RowBox[{"\[ImaginaryI]", " ", "\[Phi]3"}]], " ", "\[CapitalOmega]3"}],
       "0", 
      RowBox[{"\[Delta]2", "-", "\[Delta]3"}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.522541474095441*^9, 3.522541485568013*^9}, {
   3.522541574936735*^9, 3.522541601520124*^9}, 3.522542174052798*^9, {
   3.522542462309133*^9, 3.522542468302471*^9}, {3.522542523892669*^9, 
   3.522542555317527*^9}, 3.522542937866057*^9, 3.522544786796598*^9, 
   3.522544855009029*^9, 3.522545226303203*^9, 3.522545945473287*^9}]
}, Open  ]],

Cell[TextData[{
 StyleBox["Assume that it is degenerate four-wave mixing ",
  FontWeight->"Bold",
  FontColor->RGBColor[1, 0, 0]],
 " ",
 Cell[BoxData[
  FormBox[
   StyleBox[
    SubscriptBox["\[Omega]", "1"], "InlineMath"], TraditionalForm]]],
 "-",
 Cell[BoxData[
  FormBox[
   StyleBox[
    SubscriptBox["\[Omega]", "2"], "InlineMath"], TraditionalForm]]],
 " =  ",
 Cell[BoxData[
  FormBox[
   StyleBox[
    SubscriptBox["\[Omega]", "4"], "InlineMath"], TraditionalForm]]],
 "-",
 Cell[BoxData[
  FormBox[
   StyleBox[
    SubscriptBox["\[Omega]", "3"], "InlineMath"], TraditionalForm]]],
 "."
}], "Text",
 CellChangeTimes->{{3.522541626868227*^9, 3.522541723219086*^9}, {
  3.522541835758232*^9, 3.52254185271491*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Hrwa", "=", 
  RowBox[{"Hrwa", "/.", " ", 
   RowBox[{"\[Delta]4", "\[Rule]", 
    RowBox[{"\[Delta]1", "-", "\[Delta]2", "+", "\[Delta]3"}]}]}]}]], "Input",
 CellChangeTimes->{{3.522541725461055*^9, 3.522541787947193*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"\[Delta]1", ",", "0", ",", 
     RowBox[{
      RowBox[{"-", 
       FractionBox["1", "2"]}], " ", 
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{
        RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]1"}]], " ", 
      "\[CapitalOmega]1"}], ",", 
     RowBox[{
      RowBox[{"-", 
       FractionBox["1", "2"]}], " ", 
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{
        RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]4"}]], " ", 
      "\[CapitalOmega]4"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", "\[Delta]2", ",", 
     RowBox[{
      RowBox[{"-", 
       FractionBox["1", "2"]}], " ", 
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{
        RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]2"}]], " ", 
      "\[CapitalOmega]2"}], ",", 
     RowBox[{
      RowBox[{"-", 
       FractionBox["1", "2"]}], " ", 
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{
        RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]3"}]], " ", 
      "\[CapitalOmega]3"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       FractionBox["1", "2"]}], " ", 
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{"\[ImaginaryI]", " ", "\[Phi]1"}]], " ", "\[CapitalOmega]1"}], 
     ",", 
     RowBox[{
      RowBox[{"-", 
       FractionBox["1", "2"]}], " ", 
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{"\[ImaginaryI]", " ", "\[Phi]2"}]], " ", "\[CapitalOmega]2"}], 
     ",", "0", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"-", 
       FractionBox["1", "2"]}], " ", 
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{"\[ImaginaryI]", " ", "\[Phi]4"}]], " ", "\[CapitalOmega]4"}], 
     ",", 
     RowBox[{
      RowBox[{"-", 
       FractionBox["1", "2"]}], " ", 
      SuperscriptBox["\[ExponentialE]", 
       RowBox[{"\[ImaginaryI]", " ", "\[Phi]3"}]], " ", "\[CapitalOmega]3"}], 
     ",", "0", ",", 
     RowBox[{"\[Delta]2", "-", "\[Delta]3"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.522541788900106*^9, 3.522542174067565*^9, 
  3.52254256103613*^9, 3.522542937980884*^9, 3.522544786810239*^9, 
  3.522544855387083*^9, 3.522545226372991*^9, 3.522545945925881*^9}]
}, Open  ]],

Cell[TextData[{
 Cell[BoxData[
  ButtonBox["IntrinsicRelaxation",
   BaseStyle->"Link",
   ButtonData->"paclet:AtomicDensityMatrix/ref/IntrinsicRelaxation"]]],
 " and ",
 Cell[BoxData[
  ButtonBox["TransitRelaxation",
   BaseStyle->"Link",
   ButtonData->"paclet:AtomicDensityMatrix/ref/TransitRelaxation"]]],
 " supply the relaxation matrices."
}], "MathCaption",
 CellID->610306692],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"MatrixForm", "[", 
  RowBox[{"relax", "=", 
   RowBox[{
    RowBox[{"IntrinsicRelaxation", "[", "system", "]"}], "+", 
    RowBox[{"TransitRelaxation", "[", 
     RowBox[{"system", ",", "\[Gamma]t"}], "]"}]}]}], "]"}]], "Input",
 CellID->645617687],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"\[Gamma]t", "0", "0", "0"},
     {"0", "\[Gamma]t", "0", "0"},
     {"0", "0", 
      RowBox[{"\[Gamma]t", "+", 
       SubscriptBox["\[CapitalGamma]", "3"]}], "0"},
     {"0", "0", "0", 
      RowBox[{"\[Gamma]t", "+", 
       SubscriptBox["\[CapitalGamma]", "4"]}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.52253328148618*^9, 3.522541927263808*^9, 
  3.522542174082849*^9, 3.522542571871882*^9, 3.522542938071843*^9, 
  3.522544786827352*^9, 3.522544855406043*^9, 3.522545226441606*^9, 
  3.522545945953746*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.522542008327447*^9, 3.522542028419855*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Version using complex DM variables", "Section"],

Cell["Remove explict time dependence from the density matrix.", "MathCaption",
 CellID->690131918],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"SetOptions", "[", 
  RowBox[{"DensityMatrix", ",", 
   RowBox[{"TimeDependence", "\[Rule]", "False"}], ",", 
   RowBox[{"ComplexExpandVariables", "\[Rule]", "False"}]}], "]"}]], "Input",
 CellID->718931880],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"TimeDependence", "\[Rule]", "False"}], ",", 
   RowBox[{"Representation", "\[Rule]", "Zeeman"}], ",", 
   RowBox[{"DMSymbol", "\[Rule]", "\[Rho]"}], ",", 
   RowBox[{"Label", "\[Rule]", "None"}], ",", 
   RowBox[{"ComplexExpandVariables", "\[Rule]", "False"}], ",", 
   RowBox[{"TimeVariable", "\[Rule]", "t"}]}], "}"}]], "Output",
 CellChangeTimes->{3.522542122980978*^9, 3.522542174200228*^9, 
  3.52254257611812*^9, 3.522542938179186*^9, 3.522544788879686*^9, 
  3.522544855426596*^9, 3.522545229935206*^9, 3.522545945983281*^9},
 ImageSize->{432, 33},
 ImageMargins->{{0, 0}, {0, 0}},
 ImageRegion->{{0, 1}, {0, 1}}]
}, Open  ]],

Cell[TextData[{
 Cell[BoxData[
  ButtonBox["OpticalRepopulation",
   BaseStyle->"Link",
   ButtonData->"paclet:AtomicDensityMatrix/ref/OpticalRepopulation"]]],
 " and ",
 Cell[BoxData[
  ButtonBox["TransitRepopulation",
   BaseStyle->"Link",
   ButtonData->"paclet:AtomicDensityMatrix/ref/TransitRepopulation"]]],
 " supply the repopulation matrices."
}], "MathCaption",
 CellID->854192725],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"MatrixForm", "[", 
  RowBox[{"repop", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"OpticalRepopulation", "[", "system", "]"}], "+", 
     RowBox[{"TransitRepopulation", "[", 
      RowBox[{"system", ",", "\[Gamma]t"}], "]"}]}], "/.", 
    RowBox[{
     RowBox[{"BranchingRatio", "[", 
      RowBox[{"a_", ",", "b_"}], "]"}], "\[Rule]", 
     SubscriptBox["R", 
      RowBox[{"a", ",", "b"}]]}]}]}], "]"}]], "Input",
 CellID->465762594],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{
       FractionBox["\[Gamma]t", "2"], "+", 
       RowBox[{
        SubscriptBox["\[CapitalGamma]", "3"], " ", 
        SubscriptBox["R", 
         RowBox[{"3", ",", "1"}]], " ", 
        SubscriptBox["\[Rho]", 
         RowBox[{"3", ",", "3"}]]}], "+", 
       RowBox[{
        SubscriptBox["\[CapitalGamma]", "4"], " ", 
        SubscriptBox["R", 
         RowBox[{"4", ",", "1"}]], " ", 
        SubscriptBox["\[Rho]", 
         RowBox[{"4", ",", "4"}]]}]}], "0", "0", "0"},
     {"0", 
      RowBox[{
       FractionBox["\[Gamma]t", "2"], "+", 
       RowBox[{
        SubscriptBox["\[CapitalGamma]", "3"], " ", 
        SubscriptBox["R", 
         RowBox[{"3", ",", "2"}]], " ", 
        SubscriptBox["\[Rho]", 
         RowBox[{"3", ",", "3"}]]}], "+", 
       RowBox[{
        SubscriptBox["\[CapitalGamma]", "4"], " ", 
        SubscriptBox["R", 
         RowBox[{"4", ",", "2"}]], " ", 
        SubscriptBox["\[Rho]", 
         RowBox[{"4", ",", "4"}]]}]}], "0", "0"},
     {"0", "0", "0", "0"},
     {"0", "0", "0", "0"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{
  3.522542126889324*^9, 3.522542174305952*^9, 3.522542579319489*^9, 
   3.522542938528408*^9, 3.522544791355964*^9, 3.522544855628949*^9, {
   3.522545232720883*^9, 3.522545241737016*^9}, 3.522545946211666*^9}]
}, Open  ]],

Cell["Here are the evolution equations.", "MathCaption",
 CellID->314466782],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"TableForm", "[", 
  RowBox[{
   RowBox[{"eqs", "=", 
    RowBox[{
     RowBox[{"LiouvilleEquation", "[", 
      RowBox[{"system", ",", "Hrwa", ",", "relax", ",", "repop"}], "]"}], "//",
      "Expand"}]}], ",", 
   RowBox[{"TableHeadings", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"DMVariables", "[", "system", "]"}], ",", "None"}], "}"}]}]}], 
  "]"}]], "Input",
 CellID->298399236],

Cell[BoxData[
 TagBox[
  TagBox[GridBox[{
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"1", ",", "1"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        FractionBox["\[Gamma]t", "2"], "-", 
        RowBox[{"\[Gamma]t", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "1"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]1"}]], " ", "\[CapitalOmega]1",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "3"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]4"}]], " ", "\[CapitalOmega]4",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "4"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]1"}]], " ", 
         "\[CapitalOmega]1", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "1"}]]}], "+", 
        RowBox[{
         SubscriptBox["\[CapitalGamma]", "3"], " ", 
         SubscriptBox["R", 
          RowBox[{"3", ",", "1"}]], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "3"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]4"}]], " ", 
         "\[CapitalOmega]4", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "1"}]]}], "+", 
        RowBox[{
         SubscriptBox["\[CapitalGamma]", "4"], " ", 
         SubscriptBox["R", 
          RowBox[{"4", ",", "1"}]], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "4"}]]}]}]}]},
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"1", ",", "2"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        RowBox[{
         RowBox[{"-", "\[Gamma]t"}], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "2"}]]}], "-", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]1", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "2"}]]}], "+", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "2"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]2"}]], " ", "\[CapitalOmega]2",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "3"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]3"}]], " ", "\[CapitalOmega]3",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "4"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]1"}]], " ", 
         "\[CapitalOmega]1", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "2"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]4"}]], " ", 
         "\[CapitalOmega]4", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "2"}]]}]}]}]},
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"1", ",", "3"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        RowBox[{
         RowBox[{"-", 
          FractionBox["1", "2"]}], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]1"}]], " ", 
         "\[CapitalOmega]1", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "1"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]2"}]], " ", 
         "\[CapitalOmega]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "2"}]]}], "-", 
        RowBox[{"\[Gamma]t", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "3"}]]}], "-", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]1", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "3"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         SubscriptBox["\[CapitalGamma]", "3"], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "3"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]1"}]], " ", 
         "\[CapitalOmega]1", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "3"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]4"}]], " ", 
         "\[CapitalOmega]4", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "3"}]]}]}]}]},
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"1", ",", "4"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        RowBox[{
         RowBox[{"-", 
          FractionBox["1", "2"]}], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]4"}]], " ", 
         "\[CapitalOmega]4", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "1"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]3"}]], " ", 
         "\[CapitalOmega]3", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "2"}]]}], "-", 
        RowBox[{"\[Gamma]t", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "4"}]]}], "-", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]1", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "4"}]]}], "+", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "4"}]]}], "-", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]3", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "4"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         SubscriptBox["\[CapitalGamma]", "4"], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "4"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]1"}]], " ", 
         "\[CapitalOmega]1", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "4"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]4"}]], " ", 
         "\[CapitalOmega]4", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "4"}]]}]}]}]},
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"2", ",", "1"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        RowBox[{
         RowBox[{"-", "\[Gamma]t"}], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "1"}]]}], "+", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]1", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "1"}]]}], "-", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "1"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]1"}]], " ", "\[CapitalOmega]1",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "3"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]4"}]], " ", "\[CapitalOmega]4",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "4"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]2"}]], " ", 
         "\[CapitalOmega]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "1"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]3"}]], " ", 
         "\[CapitalOmega]3", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "1"}]]}]}]}]},
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"2", ",", "2"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        FractionBox["\[Gamma]t", "2"], "-", 
        RowBox[{"\[Gamma]t", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "2"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]2"}]], " ", "\[CapitalOmega]2",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "3"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]3"}]], " ", "\[CapitalOmega]3",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "4"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]2"}]], " ", 
         "\[CapitalOmega]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "2"}]]}], "+", 
        RowBox[{
         SubscriptBox["\[CapitalGamma]", "3"], " ", 
         SubscriptBox["R", 
          RowBox[{"3", ",", "2"}]], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "3"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]3"}]], " ", 
         "\[CapitalOmega]3", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "2"}]]}], "+", 
        RowBox[{
         SubscriptBox["\[CapitalGamma]", "4"], " ", 
         SubscriptBox["R", 
          RowBox[{"4", ",", "2"}]], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "4"}]]}]}]}]},
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"2", ",", "3"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        RowBox[{
         RowBox[{"-", 
          FractionBox["1", "2"]}], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]1"}]], " ", 
         "\[CapitalOmega]1", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "1"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]2"}]], " ", 
         "\[CapitalOmega]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "2"}]]}], "-", 
        RowBox[{"\[Gamma]t", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "3"}]]}], "-", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "3"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         SubscriptBox["\[CapitalGamma]", "3"], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "3"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]2"}]], " ", 
         "\[CapitalOmega]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "3"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]3"}]], " ", 
         "\[CapitalOmega]3", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "3"}]]}]}]}]},
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"2", ",", "4"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        RowBox[{
         RowBox[{"-", 
          FractionBox["1", "2"]}], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]4"}]], " ", 
         "\[CapitalOmega]4", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "1"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]3"}]], " ", 
         "\[CapitalOmega]3", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "2"}]]}], "-", 
        RowBox[{"\[Gamma]t", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "4"}]]}], "-", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]3", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "4"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         SubscriptBox["\[CapitalGamma]", "4"], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "4"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]2"}]], " ", 
         "\[CapitalOmega]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "4"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]3"}]], " ", 
         "\[CapitalOmega]3", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "4"}]]}]}]}]},
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"3", ",", "1"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]1"}]], " ", "\[CapitalOmega]1",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "1"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]2"}]], " ", "\[CapitalOmega]2",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "1"}]]}], "-", 
        RowBox[{"\[Gamma]t", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "1"}]]}], "+", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]1", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "1"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         SubscriptBox["\[CapitalGamma]", "3"], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "1"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]1"}]], " ", "\[CapitalOmega]1",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "3"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]4"}]], " ", "\[CapitalOmega]4",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "4"}]]}]}]}]},
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"3", ",", "2"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]1"}]], " ", "\[CapitalOmega]1",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "2"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]2"}]], " ", "\[CapitalOmega]2",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "2"}]]}], "-", 
        RowBox[{"\[Gamma]t", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "2"}]]}], "+", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "2"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         SubscriptBox["\[CapitalGamma]", "3"], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "2"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]2"}]], " ", "\[CapitalOmega]2",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "3"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]3"}]], " ", "\[CapitalOmega]3",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "4"}]]}]}]}]},
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"3", ",", "3"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]1"}]], " ", "\[CapitalOmega]1",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "3"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]2"}]], " ", "\[CapitalOmega]2",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "3"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]1"}]], " ", 
         "\[CapitalOmega]1", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "1"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]2"}]], " ", 
         "\[CapitalOmega]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "2"}]]}], "-", 
        RowBox[{"\[Gamma]t", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "3"}]]}], "-", 
        RowBox[{
         SubscriptBox["\[CapitalGamma]", "3"], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "3"}]]}]}]}]},
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"3", ",", "4"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]1"}]], " ", "\[CapitalOmega]1",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "4"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]2"}]], " ", "\[CapitalOmega]2",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "4"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]4"}]], " ", 
         "\[CapitalOmega]4", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "1"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]3"}]], " ", 
         "\[CapitalOmega]3", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "2"}]]}], "-", 
        RowBox[{"\[Gamma]t", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "4"}]]}], "+", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "4"}]]}], "-", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]3", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "4"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         SubscriptBox["\[CapitalGamma]", "3"], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "4"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         SubscriptBox["\[CapitalGamma]", "4"], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"3", ",", "4"}]]}]}]}]},
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"4", ",", "1"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]4"}]], " ", "\[CapitalOmega]4",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "1"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]3"}]], " ", "\[CapitalOmega]3",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "1"}]]}], "-", 
        RowBox[{"\[Gamma]t", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "1"}]]}], "+", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]1", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "1"}]]}], "-", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "1"}]]}], "+", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]3", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "1"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         SubscriptBox["\[CapitalGamma]", "4"], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "1"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]1"}]], " ", "\[CapitalOmega]1",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "3"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]4"}]], " ", "\[CapitalOmega]4",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "4"}]]}]}]}]},
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"4", ",", "2"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]4"}]], " ", "\[CapitalOmega]4",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "2"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]3"}]], " ", "\[CapitalOmega]3",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "2"}]]}], "-", 
        RowBox[{"\[Gamma]t", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "2"}]]}], "+", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]3", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "2"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         SubscriptBox["\[CapitalGamma]", "4"], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "2"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]2"}]], " ", "\[CapitalOmega]2",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "3"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]3"}]], " ", "\[CapitalOmega]3",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "4"}]]}]}]}]},
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"4", ",", "3"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]4"}]], " ", "\[CapitalOmega]4",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "3"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]3"}]], " ", "\[CapitalOmega]3",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "3"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]1"}]], " ", 
         "\[CapitalOmega]1", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "1"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]2"}]], " ", 
         "\[CapitalOmega]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "2"}]]}], "-", 
        RowBox[{"\[Gamma]t", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "3"}]]}], "-", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]2", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "3"}]]}], "+", 
        RowBox[{"\[ImaginaryI]", " ", "\[Delta]3", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "3"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         SubscriptBox["\[CapitalGamma]", "3"], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "3"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", 
         SubscriptBox["\[CapitalGamma]", "4"], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "3"}]]}]}]}]},
     {
      TagBox[
       SubscriptBox["\[Rho]", 
        RowBox[{"4", ",", "4"}]],
       HoldForm], 
      RowBox[{"0", "\[Equal]", 
       RowBox[{
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]4"}]], " ", "\[CapitalOmega]4",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"1", ",", "4"}]]}], "+", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"\[ImaginaryI]", " ", "\[Phi]3"}]], " ", "\[CapitalOmega]3",
          " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"2", ",", "4"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]4"}]], " ", 
         "\[CapitalOmega]4", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "1"}]]}], "-", 
        RowBox[{
         FractionBox["1", "2"], " ", "\[ImaginaryI]", " ", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "\[ImaginaryI]"}], " ", "\[Phi]3"}]], " ", 
         "\[CapitalOmega]3", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "2"}]]}], "-", 
        RowBox[{"\[Gamma]t", " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "4"}]]}], "-", 
        RowBox[{
         SubscriptBox["\[CapitalGamma]", "4"], " ", 
         SubscriptBox["\[Rho]", 
          RowBox[{"4", ",", "4"}]]}]}]}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxDividers->{
     "Columns" -> {False, {True}, False}, "ColumnsIndexed" -> {}, 
      "Rows" -> {{False}}, "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.5599999999999999]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}],
   OutputFormsDump`HeadedColumn],
  Function[BoxForm`e$, 
   TableForm[BoxForm`e$, TableHeadings -> {{
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 1, 1], 
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 1, 2], 
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 1, 3], 
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 1, 4], 
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 2, 1], 
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 2, 2], 
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 2, 3], 
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 2, 4], 
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 3, 1], 
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 3, 2], 
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 3, 3], 
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 3, 4], 
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 4, 1], 
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 4, 2], 
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 4, 3], 
       Subscript[AtomicDensityMatrix`AtomicDensityMatrix`\[Rho], 4, 4]}, 
      None}]]]], "Output",
 CellChangeTimes->{{3.522542134400776*^9, 3.522542138220477*^9}, 
   3.522542174324492*^9, 3.522542207187112*^9, 3.522542582215765*^9, 
   3.522542938655801*^9, 3.522544793068121*^9, 3.522544855825677*^9, 
   3.522545244285186*^9, 3.522545946433595*^9}]
}, Open  ]],

Cell["Make plots from paper.", "Text"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.522545604432086*^9, 3.52254560497677*^9}, 
   3.522545837580253*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"params", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Phi]1", "\[Rule]", "0"}], ",", 
     RowBox[{"\[Phi]2", "\[Rule]", "0"}], ",", 
     RowBox[{"\[Phi]3", "\[Rule]", "0"}], ",", 
     RowBox[{"\[Phi]4", "\[Rule]", "0"}], " ", ",", 
     RowBox[{"\[Gamma]t", "\[Rule]", 
      RowBox[{"2", "\[Pi]", " ", ".01"}]}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalGamma]", "3"], "\[Rule]", 
      RowBox[{"2", " ", "2", "\[Pi]", " ", "2.7"}]}], ",", 
     RowBox[{
      SubscriptBox["\[CapitalGamma]", "4"], "\[Rule]", 
      RowBox[{"2", " ", "2", "\[Pi]", " ", "3"}]}], ",", 
     RowBox[{
      SubscriptBox["R", 
       RowBox[{"4", ",", "1"}]], "\[Rule]", "0.5"}], ",", 
     RowBox[{
      SubscriptBox["R", 
       RowBox[{"4", ",", "2"}]], "\[Rule]", "0.5"}], ",", 
     RowBox[{
      SubscriptBox["R", 
       RowBox[{"3", ",", "_"}]], "\[Rule]", ".5"}], ",", 
     RowBox[{"\[CapitalOmega]1", "\[Rule]", 
      RowBox[{"2", "\[Pi]", " ", ".005"}]}], ",", 
     RowBox[{"\[CapitalOmega]2", "\[Rule]", 
      RowBox[{"2", "\[Pi]", " ", "15."}]}], ",", 
     RowBox[{"\[CapitalOmega]3", "\[Rule]", 
      RowBox[{"2", "\[Pi]", " ", "10."}]}], ",", 
     RowBox[{"\[CapitalOmega]4", "\[Rule]", 
      RowBox[{"2", "\[Pi]", " ", "0.001"}]}], ",", 
     RowBox[{"\[Delta]1", "\[Rule]", 
      RowBox[{"2", " ", "d"}]}], ",", 
     RowBox[{"\[Delta]2", "\[Rule]", "0"}], ",", 
     RowBox[{"\[Delta]3", "\[Rule]", "0"}]}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.522542639308531*^9, 3.522542653876242*^9}, 
   3.522542961645403*^9, 3.522543008355341*^9, {3.522543075408626*^9, 
   3.522543098002856*^9}, 3.522543145146279*^9, {3.522543184507017*^9, 
   3.522543189475381*^9}, {3.522543261197481*^9, 3.52254328797459*^9}, 
   3.522543355566535*^9, {3.522543426678711*^9, 3.522543426918095*^9}, {
   3.522543480813404*^9, 3.522543483724296*^9}, {3.522543596096022*^9, 
   3.522543597336738*^9}, 3.522543637720401*^9, 3.522543700412568*^9, {
   3.522543819190993*^9, 3.522543969803677*^9}, 3.522544019666515*^9, {
   3.52254405332497*^9, 3.522544053804354*^9}, 3.522544121248746*^9, {
   3.522544624131693*^9, 3.522544625322849*^9}, {3.522544811206698*^9, 
   3.5225448127211*^9}, 3.522545891693146*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Re", "[", 
      RowBox[{
       SubscriptBox["\[Rho]", 
        RowBox[{"1", ",", "3"}]], "/.", 
       RowBox[{
        RowBox[{"NSolve", "[", 
         RowBox[{
          RowBox[{"eqs", "/.", "params"}], ",", 
          RowBox[{"DMVariables", "[", "system", "]"}]}], "]"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}]}], "]"}], ",", 
     RowBox[{"Im", "[", 
      RowBox[{
       SubscriptBox["\[Rho]", 
        RowBox[{"1", ",", "3"}]], "/.", 
       RowBox[{
        RowBox[{"NSolve", "[", 
         RowBox[{
          RowBox[{"eqs", "/.", "params"}], ",", 
          RowBox[{"DMVariables", "[", "system", "]"}]}], "]"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}]}], "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"d", ",", 
     RowBox[{"-", "60"}], ",", "60"}], "}"}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotPoints", "\[Rule]", " ", "50"}], ",", 
   RowBox[{"MaxRecursion", "\[Rule]", "1"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "60"}], ",", "60"}], "}"}], ",", "Full"}], "}"}]}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.522542076105163*^9, 3.522542106948172*^9}, {
  3.522542247272649*^9, 3.52254227331632*^9}, {3.522542379859035*^9, 
  3.522542404138859*^9}, {3.52254260442503*^9, 3.522542738601174*^9}, {
  3.522544168650641*^9, 3.522544235160111*^9}, {3.522544324681123*^9, 
  3.522544357750749*^9}, {3.522544495634713*^9, 3.522544514917687*^9}, {
  3.522544546544331*^9, 3.522544579418693*^9}, {3.52254589412908*^9, 
  3.522545902384872*^9}, {3.522545984364466*^9, 3.52254603852176*^9}, {
  3.522546191788578*^9, 3.522546295569195*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwtlPs7E2oAgDGFjuRSWZc5dsYsc9nF5jL2fcPmMk6TkEWJmuSkTrpQdEIk
SkkiTiUKlSOVU55k36c5apVVCockSkrkUrlU5HSepx/e5/0L3pccvmXFBg01
NbXQ7/xv3rPtAzMzvjjgEWnYUGgGuir1f4dRvjikvaFA1WIGkpIrxvc1++Bv
b5jkG1vMwUHj6PvnM8V492ZSW14PFZg00Epipd44qDi0MGErDUybbVzAt/XC
e2WEJZFcS/Crj5g9dsgT+7FOkjg+dHA61sav4p0HDu9SGDknWIHr1OwRe38P
XEq4Zjhy1Ro0tX88oqgV4TLzEKlYyxY0EhbPP2sswitG28j05Qwgt4b5e2OF
2BNsiljHZIJOZlmUtNkdj2TzqEdpLPCZo8fjMt1x4KXYwSQeGzz/59PhvgI3
nOVb1dkeaAfqA551n9B0w0FLpIenMjkgLSSVeHKdK/50y+BO6wcuEEfYSgrl
ApzPkFEd1zqAtrFVOtuoAmxd4HF58JAjiEhPvu2ZBfFx34YPZxycAO1dezbv
A8B9p8cNzlc5goOKytOv7/Dxm+d+B9IfcEFFxpwd7P0umFo9Jzaz3gY88pOJ
k6AzNr2taWMUqguSGPVmcYY8/FQh4Ao3/oxY85Z+i+l1xF15XlVHVTbo/LPp
Ip1TDnhkt25udRgLEctfuJ0LsMf5stmbaP1M1EJzS+LN52LT/cETSIOBjl0o
lTc322Fbv+WtGTI66hTozBiWsfGFVMh+RrJAFh3RwH8PC2dotFyUm1CQXViA
lz2XiU9YVRzSX2uC9j1ILOW+tMXg1DyVp8QYVWstx5IUa3wz6ciTqXAD1O9q
2rHJjI6biBy1eD1tVBw8k/bgGA0HJn6bS3H6Kl+9tYtto07Fx9qMXUyDe+Tl
w3jBwT0U3Dduce4i7uR/iimZeDVmit3mFF6ZiP7Er3HSRmR7Eo57KnoQNjzN
P5CnzNqTTMTe74d6a+ZqgDWXSvev32iEE2SJ46UmBBC9V2us4E9dTI2Puu0e
pw4OhN+crF5AwGcG3txKbf3CFx3umDxF+og6sw0ex4je83Fa38J3fYUo3lJn
dGgql68fTCmX8T+Aqs+a4681BuV6NYsYCwUEeNYiqzxl1ld5joKs+zl8LlQW
rYltvKKOAg4S3yqjjWAT8q97ZEJAPQxTb2MqEVKN3J9YcTWQjKyb1iIkwfDj
s9J3tU7Lx7SnupUsMvyNbyXd4PNRXpI/pLx8nQJ19Z/UeEjb5H60nqu5TlR4
PCxCK0TWy+9Nb8ldlkODQae3Brx3/srf2a+MqxuwhNqkO2TSUi2QGlN1V2eX
NWxdU7SkXWgA7sWLdenqDFhET9aePGcMGs4BJo3JhKva6EmRAhMgf8gONF/H
giuUi6NLlL+AGFFpt3kmG05Zkr1Gb37vtY4YTf3bDq4e2j4wMmEJdEWaVimT
HLjGqu89pcQWfFEND77g2UPGjVTC43tMIIWdF3OiHOD8hOiBhjgWqL1yd6NH
niNMqRMolBkMIOXYhAZ1OcGiRiH7mxoFTNbk+EWaO0Psqbl0RyUBaZt/cVUF
ukDTPW+Nl72go0XZYXacdD7UhEt6O2o5qNtTmHnmGoArJFbFT/0dkb/agtH9
WyCkDOYu8nRxQgGvFfZvDQQwo1HidKLSEQXd27ZXXC2ACpP1eUd2OKBjlGK3
DrErlBjG1Kso9kiV8Fg76pUr7OyUWK97wkG0tS795w+4wYjV6yU/JdqhFMEF
5Uu6O7xikdlKFLJRjdfTEHKTOzSp3e0/m8lCQ34zw2tjhHBp7+ijFgcmClFd
VL+nKYIVTJudg5EMdN870MiuUATlpcGWNStt0SoJH/rwPWAjS79+Ubc16gug
bl7f4wF1dt51S8yyQn89zCKWx3vC6rLEq4MxdLTda0IxYOAFu6837VDJLFFN
umjfTmdvGNx/d9dIAQ05dejFhu4TQ2leWuWqYSqymf9veGGlDxT9MbPyca45
qt5FLfMN8oXZoUWmxW/M0I9/wx//Rv8B4k6yBw==
     "]]}, 
   {Hue[0.9060679774997897, 0.6, 0.6], LineBox[CompressedData["
1:eJw1lX04E3gcwDdTkZfkZRuhxjYrZ142htX8jG1mXCH0KqbkvS5KLyolWhQP
e6507rru5PUqruT0wvdXJC9FiV6koRd56VKW0Cl39zx3f3yez/+ffz40+bbA
LRoEAiHyH/41/1ny6OysP5Y/uaLUV87Uqy4afIdi/HGWDt+rS/y1/tDh85/S
Ov2w+E3yWq4RAebc04kwsvfDXcMFMSUcIhyjxLUVZ8vwha07Z47v0oCuP6V/
6Kp9se2WfEVKMwksG1lFSet88ShxS3SV1xzISfbpfdIoxeEihqpqcC58oUeb
COyl2KXp+gKzm1rwrZ+MM3HcB7etfnh6SF8HziSxA86PSLDAPmiKcEMPaph5
73lBEqy3POlUYaAB3Huqzm24LsaTwTvaFXMNoYlkZvwLRYx3E6wU912Nod4O
FRxIEuHZ+NVhf70gQ69jacy6Tm9cdLluf7ipGUw76/NdHL1xP4n+9DnVEp7f
/nhi8Acv3HipSRwYSYObwc/6T2p6Yc3bho+11tIhc0MG9XSEEHP148YnM1gg
i7RfVVjviUc0q6XfmrPh8cQa7R1MT9yb5XNCkuMEkYrDt3xyEE7vacx05fGA
NfI0jz/ugUM2RBSXWK2AYw0Xz7y+I8Dbtd+Th0yFcD5r/k7OkRX4yt31Aw+b
RHA/IEp2CC3HWYvdN6piJXDI4SZ9tyEfM7WXTRp2iMBpgfnXxFduWN22t/2L
hhCKn305q/2TK363oPZTjxEfqGV9XueCeTh5piMs5iUHulleh/jGLjjs6tj1
B9gO8stL6js7uVhTriiYGGVCr6f2rGEpB3+aW2cXedIKbHriPIL2OWHqkCLE
1mMxcMODpTwXR/yBRl7zjrAI0u7uL3F5YY/71AW5wVupUD1vJV6Vboc7t0Xu
30qlwLBwSU8s3RZXnlYGckrJ8Ova2cy7+SxcuV6lc4BIgfXbVRw2kYlZ/HMz
KVZUKBvDJsf2WWPs/FJ3eNwUPiYWTb6cWII7E3qIkgRzqHXXAhrPAlPGkqPM
fBfD0VMtOfsOU7FU6TrSKKNB2G8lRzZHG+E3jw+mTPOtIe7AvIkfftTF2p35
HbcYDDgqvzZVbULCDs/T51MqmSA+0TP1k4Ua8vSn69pu2ADOHCSPDBaCidaC
iIMWLDBYa10WJRj3IDeabPa5bQP6taYOZE8S6npSdWW2jAnKBprutFwPPbv0
5aC9BgOCj1GHWuKM0Kmz0vEKd2sYcFjiS2FS0er7lY3TCTSIoulmdossED3V
LS8WLYYJrZn+Ficaenjm0XDqMnMoKnjXUlljjTI/pN6h3zSFANbApe/dmYjY
dIdRTaTCK0X390uVLESQmXTK/+m3a7hld93oMpTsKmVWlJAhI7GqWTvFDvUj
5UJ9Ywq07pHp2hIdED16PsMqjQqN5zwcWY6OqCr0s6pKbQb1HZwQRoQTwuLC
4vIuS0gUl/QzsjmoptVqZ53QCizrqHHMK1yEajN2VccyQVes+U36lDO6ejKm
IU/PDj63j73t4/OQf4FwWXwsB9ah3gpljCuyeHTZrfWeG1z/vTlacsoNGUyn
7BG9RbDOmb0xVOWOZrLkYpJQBFO1yoCtjOUowlPSnRslAS3GZ2F7yAp0y5gX
ymGLwTQvnOusECDugxcfLgwIod9HlP3zZQ/EltDyMrIFEEQw+XBkG0LxTkvl
3YGuEPy6gTe00BM9Xf3cR2TOgdDWHQdk1Z6oYt7QNzaZbMi3/tWrRyZEx6VO
zRZ2S6E99YFWzEshao71fqyziAGsTSuGi496ofwi/5x912iQ7lne8sLWG2mu
od5XCSyhVtq1gXbPG1nX0fONeWbwLmB2bFOiCNmzywtzP5NhQ3sFsVVTjHKq
TvyYkGYMbb4hRtxCMcq+pC5TJxjCmlUC5CeQoJ9rzkrrswxgMJiZsHlAgrTa
48vTBvXgQkcOtWyPD3Jr4gdprNeBZOlkw+hCKaK6bwxLZmgDedjSyrFUiogP
9vrOQfOgViFO27XcF9WYhL29oJwDe36vmI5R+aJbfaRc+2kSuPfoJ21Mk6H4
SVKRMEcDuCnq2FSKH9q8suAgIZwIbOMn8sKLfihacbyubyUBqlOYpf6h/mih
4SubI/u/1v/3B/T/H/4GZmzqsA==
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 0},
  Frame->True,
  PlotRange->{{-60, 60}, {-0.0007146440279521698, 0.00040010361659988916`}},
  PlotRangeClipping->True,
  PlotRangePadding->{Automatic, 
    Scaled[0.02]}]], "Output",
 CellChangeTimes->{{3.522542156195219*^9, 3.522542182530973*^9}, 
   3.522542280444924*^9, 3.522542411408412*^9, {3.522542593479043*^9, 
   3.522542626285554*^9}, 3.522542681997767*^9, {3.522542723280386*^9, 
   3.522542748242407*^9}, {3.522542949188668*^9, 3.522542977352989*^9}, 
   3.522543024498731*^9, {3.522543088219265*^9, 3.52254311067984*^9}, 
   3.522543163269546*^9, 3.522543203476909*^9, {3.522543258287421*^9, 
   3.522543307465341*^9}, 3.522543371309052*^9, 3.522543444314203*^9, 
   3.522543500408285*^9, 3.522543615581949*^9, 3.522543654012719*^9, 
   3.522543714809282*^9, {3.52254383655232*^9, 3.522543954209074*^9}, 
   3.522543987854934*^9, 3.522544040177055*^9, 3.522544138645664*^9, 
   3.52254423961056*^9, {3.522544332311615*^9, 3.522544359410752*^9}, {
   3.522544499019484*^9, 3.522544519115584*^9}, 3.522544551529608*^9, 
   3.522544583202106*^9, 3.522544798152048*^9, 3.522544856018762*^9, 
   3.522545918022345*^9, 3.522545953245106*^9, {3.522545984075476*^9, 
   3.522546042552423*^9}, {3.522546198674039*^9, 3.522546306798113*^9}},
 ImageCache->GraphicsData["CompressedBitmap", "\<\
eJztnQe4VdWVx5/wqCZmxpRJMkbNxMTJmEkyTkwmKkWwBTWKLYpdxN4QsRsT
RUSNCnbABgiISG+CIiJNuvQqHenSRJpwZv/W2f/NuXjeue/exzeQmckXn957
z9l77bXXXn2tfdZVd17f5Jar7ryh8VVHnnbHVbdff0PjZkeeetsd7qvKB5WU
VDqspOSgWkeW8N9RSYn+RO7/8Z/D+cPXBX9xvP2rarRjx47o5ptvjho3bhzd
eOON0eWXXx6tW7euUuLnW2+9Nbryyivtn9tuuy3auXNnffu5etSxY8eoQYMG
gBU1adIkWrVqlU3kX6vEOyWV9Za9VMf+Vo7mzp0baUWV7W+V6JBDDinxYGqI
krQxapcHuopNUa+iKywWtZVy96dFixZR8+bNA8wvvPBCdO6552qQxx57LLri
iivCz8zRqlWrkz0GZsyYYZDv3r27JM8rJ/lXzjjjjOjhhx+OJk+e7BDD3yrR
q6++Gh7X65X88/71nH3JA1VFpqhf+MrqVgiVOftRM/rxj38czZ49Ozy1evXq
qFKlStGmTZvsM78PGzZMoH/wwQfRT37ykxzQZ82a5c9xFED3r4Vh/XvalBUr
VkR169aN1q5d61+pEu3atSs8rtfTpq2dPcXJ+2CKIlZWt+L4PCjemd/7n9av
X28U4c5VmHz79u323bx588Lva9as0WjuMNtXGzduDHypTZs2gWmU+r+c4Dyv
a6cWLlwYXXLJJcYnKu8FmXu1JO3dWuV87uQKzFHB5dXeN0j2W/Yr/1O1atUC
MFX8XzeCbfz8+fPtt8WLF4eZJkyYYN8xC7Px34FlVuY/7Sv3aslxZU5RPUyx
bdu2qGHDhtGWLVvssydfDZeYoaSwKSprisRJqAxOogsvvDBme+Wco1Y5nzul
AnNIXlUQHSKRCm5czqk+OFqwYIH9dO2110bLli2Lli5dGp1wwglR/frGdPJC
VdezskmTJon5BzZV3pXF7K0arDE65phjorFjx5Zr7vJu3qlFjr8v1lZ7H+F5
r23bvHlzdNBBB0VTpkwx9alGjRrRiSeeGC1atMiG27Bhg73q2H0qu4rZaWnU
smVL2P/etJjv9ZM93bVr1y467rjj9p4yiwnnG/q0woeu4GJq7yOc+i063v90
7LHHBsFcpQxhPXTo0Oioo46q5xfAhjdq1CgaNWoUo1bOfusk/zPwLlmyJPzs
eFNUpUqVshQE/3ad7J9PrtDgxS2obuGY20ub79KlS3TWWWdJsYqGDx8eXXfd
dUkV1KmdYZDLLrssevzxx6VeuXVyHmMRk/2KcF+nTp3otNNOCxP27ds3Ovro
o/NMWCf755MrNHgRq6lbIfT5PTjBc8u2bdtGl156aXTLLbdExx9/fDRgwIAg
p7DJbr/99mCTNW3aFBqR4tW7d++v6bdlvKINgJdiSZxzzjnRDTfcENWrVw9q
y/N2nfLAU+zgRa+mKNx57O8Ln8W++OLE+Av3vzmmDCa+dVQodPDPP6R8948V
eO7/2rsHOnxV/3/r/3/rv/bTGSnfnfk/8FyBIB0TfxetmLoginbvtk+HV2D6
H/nxNi5bs/ccFQH9Xj9sn5uei1ZMWWCfbk557ib/3Np5y6NBd7W3T/cXNtU9
foj1i1dFbWs3jYa3eqvkloypJr4x2J5bPWtJcVOVRmOe7xO1rdU0eu20e6Kd
W7dnLazHNU/ZoyyuwNnutr9Vol07voo6nPlA1PWiR22ouYMnpOLyRv/8ljUb
bYFd/9TCnl81c3HJfXmmtl804fyhk+3FpWNnuUlbRP1uezFzwhm9R9nzq2ct
jjqc9aCt9Z4waoakjSesGvW//cWoywWPRLt37Y4mdhgStatzR7Rpxbqbypyz
ajSoeTtg41P0SeehWui9haH4LvtbLdq4fK2hbFKHITbi5hWfGxCOVkrSoLjB
v7fjy23RK/WaO4Lobe85Yog6/vHBaGDzdiVFgVI5Gvty/6h93WbRF6s3BGPO
bYBtBAf+xjLBqRwtHD7VELF80rzwbgI5OirlBam5H2LXV7uiTmf/OXr3nleC
29ERoQ372SfzM0Ea9miX6I0z7ndbu2uP+yCBpSJBWjB8ik2/eNSMABLDciI/
fKxrKkjX++ccKEamQx9+s6TUfzel6wcxlqYv8lRZMEhQ5Jvn/tVWqmEdLxKT
EEZSQXIs0qb/9INPwrvC0uD7XysUpDv9ECLj8a8O8vbqnqnmDBqfCdLHL/Y1
0t7+xdbwLstxICXptVCQJrz2roHkQAvDOrKO3rq4ZdTv1hcyQep2yWPRgDte
jrmO/256jxG2nC1rNxYLUufz/mqkuOfxPULiy3WbNH3Ou9f55zYsWW3TOzDC
uyumfGrfLfl4VqEgNfNDuMNvQ8zoObKkmv/OHepwuLNAmtJ1mD23eeXnASS3
ibacyW8OLRYktxobFmkukBx2YjB7jcoECS2gx9V/s096F+n0/l86Jo9QoSB9
0iU+tJBodf/dGw3ui0a17pmcPufda/1z2zZucXTYLBr/yqDwrjtpUbdLW5UL
JPtFgHzQonPU+fyH7ZMG631Dm6h/05fyAuLOoS1izeyl4d2Jrw824HZu2xHP
kU+QCpB3rvpbUKM02Ifu0Dq+lBcQJ8Oi10+/14SM3l04YloQHgVu0h1+CIQH
fMTJtZIagUT3bFwWSLBUp0DlLMfpBfburL5jigXp84UrbYh5QyYGkJwgiVc6
Y1EmSChhCA4+6d3X/3BfNPLpHkmmUihIUrnWffpZGHbjsrWeU49LBamJMOz0
QzQGRzUlNf13fW953p26ZysC0rj2A6L2J90Z7dr5VQAJ0flKfbeZL/XLBClg
+L2JAaRRbXqZLIS+igTJKSBR9yuesE8atvuVT5piktikVJAWfhTrRmvnLgvv
zh4w1r5zmC4WpC4XPhIN/WunHJDe+3MH5FlekBw3tumdFhneXTNnqX3nFJws
kL4Rf8jOJaka/fznP7e4vn87N3mlak7ySqE/10z9OSXtpWrKKpr61e7Yss3L
pfdLDvbfxRpCs+irHTszsZfka8Kee8eOAgpPObCXmflR1UJybgHhZ/+4fk5m
p1T8Z5+kIh3Z56jYpyoZGFw5fWGsCI+eETA4//1JgZlkYbDPjc+abccnveuO
VzhOzfJjsIwkj6rZP+dJESnXz9UyM0iyUDazz+igIGnZ6+Yv9wr45EyUYYeP
+Fv3HJSZ9L/g4Xwo0+P5EjZKy3yuavK5r2dcFPWzO7QCIpHOkYm/kc/0QFHI
wYHTVEytd4c3FX/X+OdQu8Dz1G4fehrao7Bu27QlC3+xnaekC7cc/u5JG8kX
pxZav/5c1eRzX89SKOjnRA6EgEikQGSiNSmyhZqkktwkA62BC4yZGd5dNn5O
sBfKQZb7C381c37OyEYweKuwgts8KPgVRjzVPQdjTvtFJciLsTkDx3nRvya8
u3X9ZvtuWvfh8Rw2W/VcPOULtpeW+VzuAfxa6kCen2vm/JyREpAmb2/3QAW7
0hls3/Tf4cxEp8NvkYWxsW37R+3r3YlCH95989y/mLuFT3cULi18rLwsaZET
Sq/Yz7mx+xoZKErauVrmrH5jgp6YhSIMyLcvi2W63jUbpvFT5UVRGaHsqtk/
lysQnvfnMqLsWQSVNMG1YkIT8s1lYav75XtUDb07rt0AsxYhsnJgq4zQc9Xs
n8sVuC7z59jdmB4Sz0JV0kmg5TqRZ6j6pPPQVFQ19s+Z3VUvtrsO8d85VSVo
eRmoSjNeBFLSXaBhcTnqSF+TAdKmz7wd3v/j8K6FYIhQvDu+UJDE05Pugm/5
7/re/FzU6/rWeUFaMmamP7qfBpDA3Ksn32WYKxKkpLtAIDkdMOhAWSA5UWLv
OtES3u1x9R6sNy0OpKS74Ft7TbVlzcZMkJzExGNhn/Tu0IffxJKtCEhJd4GG
JYBkKsjEuZkgEd/ode0zOSBN7vS+aZS4zosEKeku0LBO6Q+e1SyQcAq7E5vz
7qKR04OrrkiQku4CDYv7+7VT7zEXUhpIV/vnCDYxPUEqH/BGfQknLh9I9sut
Gixhef+jCPOapwjFJfGQCsiaucuCsSRAwjlz5/a2MFuGy1KArJq5OBiwAsRh
iA3IC8i89ybGx2D+8vCu8bS749N1e2GbJJBm9//YhnUsLgw7/tVBhi/IMQsk
XLd6Lrkcp8lXBCScs9AIzjMNO2/IhLD6LJAMm17w6N3k6SoSpAHNXg6HVsNq
Mz8d9kkmSFDaAO8c17vyNa+Zs7RYkJLqqIbdvvlLG9atOBMksDvi6Xdy3t2w
dHWIuBcI0i2ansiLG2LKW8NKvu2/I/j4Yau3ktPnvHuVf07a+rR3Pio5NHHO
zAPrzlmRICkCSHxKIOEK6nPTc3lBUhAq+S4nTnrcbcWBJJfMF6vXh2GHtexC
1DkvSMo+4LTqXTtxcXZAsSCNat0zSHgNy6mRtpsF0ujnegeTSu9O6vhefOK2
7SgWJMSk26gckNxJCzIpCyTMEKdw57wrl7jj58WC1PGPf44+evLtnGHlS8P5
nwUSoRSC28l3FU91vK1QkHyGSsKK/6jkO/47YulwZUIkaSBd6Z9zyov5jzld
ejd54vKBZL8IkM8mzw/az3f9d6gVCh9dlQHI+kWxzjln4Ljw7jtX7Yl03Bxm
yxCxAmR6Tx8jX7cpDJZUubIAkc7jjLrwbny6HrVPtxa3SSihHc58wD5pWNOp
vWKaBZJszq0bvwjvTrKsoThiWiRIBKn73vJ8DkiOawf1PQuk4Y93wxuV824y
4FQkSObg8jJJwypQyYZkgcRyel/fJufdPadrYqEgKZ9ty9qNwYP0PRGrN/sg
1iyQyCEa9miXHJCCZdu2f7EgLR03J6QiCaRwkN2waSBd4Z+TeEY7+Cf/HRHF
wffGVclpiYTlAEm22LaNW8KwsLv3HnwjiZFUkJI6lN7lxCmiWCRIkKjEqYY1
oeCdpVkgBU3TmZt6N5w4p0EWCZIjz5AioWGJEr96yt2m6GaBNP6VQSZOv9q+
I7wrXkUMv0iQSIBwIOSA5EzFOF60an0mSO8/1CHEgfSuU2jCcSkQJJ/CFixW
Eh6+vxfJI1qyQEI1U8aU3kX7HP5EtySjyZZnAkTulZXTF4bByDCz9fUelQrI
5f65YOM+06PkB/47B1jU85qn9yw2nzwTIEmnowAhdRQpTSZmFiBuE0PClgD5
+KV+wemYDyNlgGQ2lt97DYt/U5boFRkguU1Uum14V869tfOWFwvSW41aRkPu
fz0HJMtwuuSxvCBJpyZLT+9KeMDlCwRJyaBKQ3O6cMkP/XdjHfLNb+SQnwXS
6Gd7eZ16dwCJTSelZFSbnsWCJINm6bjZAaRZ3v4nNpQFUlKntncrIXWLBEQ5
UVs3fBEAUSYhYb0sQJKiRe86nBhuwFdaCnM5QIKhySuoYbf4lD9iZGkgXeaf
U1ISAvif/XeKuTlKygJp/+e25LZ1kcM/rSNMWkBX6COldcgD8QEUCpAv0ugu
z0CffM7kfuldaiTkdisH+vZfcktu4xYFRJONW9LiJMryNYeN90we5r8zf5GX
Y1lokwaAiBLaQkThpX7lQdv+yGhRMnOyjUoWiiQxlo2fE1CEBOpywSN5UaR6
AFRTveukcEBvWub8AZLBUprzc5s2beIckMpknAVYK2fgTYpDcu1JZSILb+S+
dzgjtlj1rlOjUKfy4a3kwMxcqZHzcxlNQJSvnYpPEh3l2BZOFD/HGsvCJ/ap
nFmHeQ7rtGR7t0v7UVG/vrOirVt3puaVHyCpKzkINKSlNeQwBFomyZ4c4hYO
dW/koE1yFb00C2043Ia17JJ4t1q0YFzMDi6p9VRUp1bb6MwGb0QvvfhxtGLF
ptjbdSDlseTOooO9V4eMNKJTMYHFzGvHMfPDA1YexMBNYirn3UuTssWb+zH2
akYzpq+KLj6vU/RyrTuiznd3iqZNWxn99S/vR/VPah+dVKdddP99g93+rkuL
zB0A6S2VPQL3amiRhUAlhi2bMDcgcEDTl4I7MAuB8mk6QetVwqqYdtHb3aYa
wi65+K3ozYta4j+p5A/12rVbotdeHR+d80enPp3+ejR1yorUQMwBkAgjmk92
qMjCpCVw144zF4VJ0qlwYOD7z8KkXER4FYTJnj2n2wF+8P4h0RdfbI8+eORN
Inelfo/B5BWXvx2ddsqr0diPl6Q5sg+A/BgpM8kuE2nKjJDoOGHwJguJc9+N
i2LWL16ViURiRBLKP/Tc8KYb+0TXXtMTsrS9ndrtw2BYCZMbN261Z06u1z76
cNinaR7CDPtIkJOd6URgDuRJJT0Lcjwz8n3JnF29+ouobu22UdcunwSWLPsU
Va+Khx7SuPXmvsad+veflebuyIA+WQFlkcmuH5Qc4b8zt6kzxse1H5AK/SV6
bkf83JgX+gbLsPvbUw36zz7bFKAXv4VXV/XQO6Ee3dV8oD3bs8f0NCOyHNAH
D9ikeSVHBpw+GRJrLs2AXvkLBJkEPVRzXZOe/q3E/nqzWtDv2LEr+vOD7xn0
Q4bMSxPV5YDe2TbGOFA6Bb3lh/mkhyzoVbSyftHKVMoRE4N54OvDk1otAX3z
Owca3Y8auShti8sBvTNqgzda0IckbqdaZ0GPW0LsMTjKEpQj6OWD+XzBCq+5
xJRz/XW9jP85KZ02TTmgT6bgCPokNWVBb5lo17XOObVJyhH0Kt2a3nOETzaF
52yLrrqiu8nATz9dd0lh0KsAQdFNh++SH+/1HT6ZNOgb+edgiTggCNimUY6q
RaliFl8T9Dz3pwu6ROed+2a0bNnGRnmgt18Es7LVyOz4F/8dsSDlQKShotFe
mJzZZ3QqvQhmqB0tEGqv6WF2cEZnn9UhanRR12jdui1/CoBl+JgFs8ORr1vf
HWAeeGfbUI6YBfOCD+PidKom06hEMAcaX7jSq+Wl0by5a41CoBTH5y8qjkqw
MOR6EvQUdJOhhaM1C/rx7QeGUKq860kqUQGv9UzwoY1veOgnT1pup/PG63tz
Wv9UGPRKDA5BhZ4jS37iv5s3ZGLAVRb0g+5uH719+eP2SdAn6UVlpE6tCjQu
6EeOXBTVq9vepJPjlBcWB33gz053EfTrFnxm381/f1Im9CQwOUrOgT5JOYJ+
mJMVzhK0T9/00DtpZKvEdtm1a/cFxUFPJqiq6gR9Mic3DfqLU3iToN+0aZtB
1fnNyaGKVPl1SLBDPPTsEapvyxbDgP784qDvdM5DIWvzKP8d6bbu5CbxnAq9
ZABmkaCfPn2lQfXxmCV7ynITdP8tD33HDpPsucdaGvTnFQa9ck6drPNO3/El
P9VOOw3XUUVe6JOuLUE/cMBsgwo/gCokoXtVFAn6N16faM893upDRMS5xUGv
TgzkcAn6oHuv35wJfTJ8JujbvjQ2OvXkV8FngB66VxD7Hzz0r782waB/8onh
QN+wOOgH37cnMibo5WYlypoFfbdLW4WEAUF/3z3vRldfFTv6Bb3iH3CGf/TQ
v9J+vEH/1JMfAf05hUGv5EWyKKhCpxrqZzqNG78ImYdp0HvRYrlR8HtK9ySt
Lm30FnwkB3oV3iOJD/XQt2s71qB/5K9Dnem36+zioKdvECM7uVlydIITys91
cQb0q2YssncXjpgWcO8gMT7uKCPUdjnqCnQv6N/sNNl40913DUJa/bE46NHl
yUQgHUnQU4WmThJZ0IdS1BWfB+gXLvzccDr0/fkBeqP7hg/Zp2976Pv2mWkW
4C039Yk2b95+Vh7o7RfBjLtHRUlHJ9ahLLEsmNGMZHML5uHDFxjM8+etDTAn
qf07HuZhH3xqtoc7GWhkDQJgGRqZ0gWp7UZ/cXpMyb9q932uDAlNWTBbZru3
tgUzuw/+tm3bGWBO0vh3Pczjxy01nQbXltMo0xqrZVDJ1xJnx8wM0Nupdfo4
nss06L365HuV3GFeS53QR1t8gD5un1RNBI07PmafBP3MGavM83bWGW9EEycu
b1Ac9FbcUe9OS/b5uf+O6go5/9PUVEEvbdid8gC902fQsnKgxzJUw6Hveeid
5DJNGNdin94z/lAc9Oy+s9zsk3QaKHHI4LnRl1/uCPb+1Ldjl5yjgswVkahI
n5mkDYvO/vxzo3NWpEI40sfjhKAq0ZYtO6J7nWzgtLR+emT0ldPGTytsWd7f
E9LuSA+Pl1WT4aILz+9sw0O0DzvGPHrUYmOJck05BTVrdRa29nkjSf2fIfv0
mRkqqxBXpLOS4vJ9vzonqy3ywLNNb+8ffb7uS9/kvODVqQ0E7Eir+8BtGkM7
ZTJq/cxII2w+n3VGh6hr58lWne/U57TVeWU+x7bQ3k2YsMyGmTJlRSjKoeDC
qHHzl/650ujdQXMMqQ3P7hRNmrT8lOIWZl2sPC/Swm64vpdZqeg67hvbRqdz
Bj9Z91tfQiJmLUxKNckVWlivnjNsYZ9//mVYmFgmmQU/9Nu2YMG66LJLuhk7
7PDGRAOjfmGr82680L/EgRPOGhEWwOjxzjTvIakSbd/+lfmy7zijTZASaeai
ludkd/COaXmcnzMbxP65ZDEVHp5Bzdv5DIAqnHHjl8Bw5x0DovXrt/prfgpe
XuiK5shCy/vLn9+LzvjD68ZK5KKdPWt1dHad52UopK3sggTEVhHyWNfgjbyj
aX/sPvukCislShFsOMxPM2hgTJEcBGe9nlTYopLJedgFJONrUStXbjYd6sUX
xgSLnyhO61rNo243vZi6V1qRGA1SWSu64LzO2Bc5K0rSyo8SpIh2YLbgo8Pw
gegKk0KXpTqX+UMnh2W55diy3PKCUxTe2LJBy+j52s2jTRu3ppn9Wpl5OX1b
OFEhHB1o8Y2oSgceOqBZnFChleF7bd9unE2PF23sx0t0dUa20qb1jHj6ndBM
UE4dSBvaczQYDNQRr8cZHq3u6l2StZiBiexfLQa6ZTFOZoTF4BCC8FBwj/CL
4bkrLnvbmNMTrYYTUIl73efT5kJ+mONw+OvwxGpzHH+wyR2/CH5HaWTX1moV
fTR8YZpX43wNueMrcy2Pfq53IDucIgy5dOmGsJ51n8bOmFn9xngHcFVjR6+0
H2fK6R/P7EBYhVNxfGGUp8Wp7pWaRu0UohEW71h9cPbBqB7/wyN2OJwamub0
0OI++2R+SAXU4qAmAIZdq4AG6SJN98d+s5YsXh81va2fIeKmG3qblv774lbm
A25W5auVjRix0IZ20jmsjGyZV065J6pb62WcI1krE3nRBUsre/CBIdGVl8ed
PbQyFAIyCvEl/IvfNkLPaHHwP87VM0+PiDZs2Pq7whanmBwWhvzjokkkMAoV
etrBHp/yRT5yW+zNR8tJc/loheS5q7hOK3SrY5U5K1TRD9bHUX6FzrKzVaEB
c9A7dZwEqRxX3Aqli382eX5Y4bnndIpaPPKBfdIKVeUx/MUBRmF/e/Kj1BX6
7yyxVz4ErRAdhnfdsQrVHs7mDI0gf+pXiJf+gfsHGwmB6fffmwdF/2dxK6Tt
InG7pN8eHHq10evmVWjTYW2Hn3t2tFHOwgWfl6Q5vrTCkFk4bWFYoWMqNq5j
MmGFilyTT3y0XyFecSwensWa+fDDBazwPwpboWJQSgrG/6g9JGeF0d1RDCuk
YhmrceWCVWb9YW9lrdC3dbL4gFaI6cC4c2avCSuEc6MtIN9iS7aancMRHy2M
rr6y+96r1J1facvc/5m+pcmfuTdS93glrotUVC91N8gBJyeLaIl2Y0D/2O+6
ZMkGvxtV43YZceq7Jcvw+/hxS9O8lecmEE2YjDCFklm7dP7EXkXXVK2TZY7X
jmsFfp7YEDYANqMNQRbCxv49/4YcSI3x/KWRle1qs8rJ+yJTN0ahFFNCfa2B
QgFof7AksHCIn67frS+YurbDaQDIZ0jYkW2aI1Ybo0p1zrg2Bt+5E+P2SRtD
jwW8ssNadPY99auZGjB06HxTnozjXdA5eqf7NEt1TOu7fwBkJyuZ4+sXRVZO
XhSZuRs4+NFdsZu0G7hdJeS1G+qGi+sE+gVF7jhl7QYRP988PuzGrbf0xf9g
n1SKNeaFvubDIj7xi8QxGT16sWXTyJHQ9uWxFjH5edm7sb8TodNyKHNvh6ye
ArzCFQSwUKRwpCrUcfGfukYPORtiz15UC/3UfbFYdPNNfUxZgPekefo9L7O2
BdaP9o3BIb0V5+ljLYflbAi1cWwISUS/8sueOXOVpdfgieAfMuGc+Dy67L0o
OTCTq3NnkYcxeUtjWmKhtojyaFBDV3MdF6f1GUpefWW8t6CrmpqNzGarZjmT
DHPMWelZuzPJ369A3qF255T6r5iaxie5/Ag5wLwcEyv5tacI3LEvv/SxHROO
C7wSYYZ/6aj8x+XA2AvNkrx0cU9vQfmbaG1OtAHi1yFxdpMt22mofgeqBScD
Pl7ozWnS5pFZvfqLkrSAkudkVromP5W2gaxExu/Vc0YwyCXGyNI71r6rbpSA
gdrkmh72PLbBffcONkNs585dPwqLOSBTvpUfmrhHMbWMSrEb+g5CsBuWrg47
gY2geJN2ItwiEHe0NQ8W8eOWLYZl7oTSoem4rJ2YP3+djT9y5KKgC9PvBr0X
/fg3fk5nNVic+ryGnex57NDnnxtNfpdyKguR5v9z/RDLvD8xLVNXkQGsPFgC
OZGKMhN9hytt3/6Vj2lXC61OiKMwNI4VnsGKSAtTaieoN0SXxrmgnRj78RLD
7OxZq/fYXeNmBwvot/5MoFeNGbMYIeZUvFfsHc4HeSXr1289LP9m7M/Oi2Vc
qJi1GXZDRa34GiNtBp7oiy9K9j+qZlLWrrJ5fbANjceQAIZTkfwlUVXs6LBZ
TEOqEkX1XoaQNm2i6OOX+oXLovr3mxXHFdZ9GVQuvIkQB+6Z//KbQsojOe2k
g/E8RHDbrf3M0ecY5A/yb8r+SGDXz2Vcspi2KZIatGOjSJn4meLqmF733D3I
Ph3qNwVeAvf3N3gZNwdBHw1f6ONl37SWymQIojqTw4Rr0224TcvGIxWIv+HF
ICwBf/TB6TjvpnYcVdGmcVKRFOTFa9NocEU4gzuzTvCbRv4BWSkyHdEp8F4Q
dHMW7ffK3rQ0r4UwQ8YBmgrXGolcObOcVadSeB9Y9dC7HA8VmOEZ1AwUSCfb
4qjUoUbVCAYifjxOyir9BOjaBFsAAtx0OEFAYK9rn7GFchKgVHoZOBF8uocE
b5A8dMLXE48PxzqzT4Ebvdzf8Mo1QbU8xHhXySUAR+AKnIE7p6dFc+eu+XZh
+JJzXi03yCAWJcnfNHDA7OAzdEfVYtH40UPH6fFxPPOtrlNKdB0pBcXwVzyt
eOocD8ennGazVfJEigrEbVyk4gAKwpUKPQwtXTseQtzvjg8lELhM3YlPIK7U
FFSEJ6H43bt2xRGdGiaj4dDYazAHoCbL+Km/jTDzj+jPtwpD4Pl+RpImCDD0
ur41EAuHo0YuCjEGdXii/wOEQX6ycIgH+vRTX4vWOD2utschJ1KZEOCc9zav
+DyNIZT6VyiooJdDn5ueNTz1bPI0/bF0Eys7Yvr7th0BfVC7TDHR4qcffBIk
Xhy5rInSbUWSQIoXmlXhycSWhR4J0bjTc3Ax6ItLW6EVWr4oXwUOwCzsinCH
TsAhpG+kcMeuogOTQVbLixi0LAiI04NpQTVOmpFazT9PXzOcX76HgGHoJI0/
baGNM779wFC5QtDmuWfj/Athsr/jOQQrmOtkjzW80nhLOZ2cWFEdeeKPPDw0
GvzuXERTGmgZWJMP1awBB9m4dgNCngxRzPPOjXN65UNVriTRnEp7YXfSpOX+
2ugqhl2wxbPkdFMFVDNl/up6fvOXtnN2/d9Fj1rmYh1/IGD8kF+nhg8RaxGS
iHyJWQiZZJaw+STon+5fd0qUqd2kuDr9oRIvua9xHPAVuRHOPEzjJxl4O9cP
TugOHgzMjqMJdRR73NG0fwJ1NexEEeR6+/LHLQEHIbHDvU7K4uWXdkNIeP2j
CrqRoS1JRGnnoaZ/Htc6Y4M+RCqKqK6HnzNoXMibEJ4IXShV78gEM1EBWAM/
Lm4mTAUcfYQ6nP1cEiOwioU8SJ8kV2PBAtPBsuKn8nrRKBswyVb+N/8dXKBN
61H2ST5IxS0IdorSxo1bGrb8tx5EYMSRzLbjTaODwjdT5leECIZulObG7nVd
a+u+eILfTRrVIjRoK+DGFbY+mfxZSLM5MrHxhMfoDeRE1Fl+dKQ+1jYIQ6M/
95xOHmF5YuXychjDdOghIVboURJT714zgpcD7y9J48TU01j/bzxApLZT/sFy
EZ34MQ9JQY+S8ZkfbRjm7v7bx5Et1glYNG8TwbznQ9TLl230EdtSU28QXhCj
4wENw/oyKOMc/y4eA3RF7pNw6NXqpRRMduxFqycBQ7d27M28EUH/6YckERpm
zLAz+4xOu233EI8o4hyIXBBFpiimhA/HWicrtGPycFW71LVLHNBwR9jH60qt
yBgdxWnR0E+B+ehn+1HYMlQOsiLdf+viXNXYYrvIyQKLZW2k4VbaC7APh33q
/V9VLCKjxmN0d3K2Uuo1w99KPm9huGb2POr4b7Xtf+5gDiZuoREyCEY6QWSf
FExSYR1X9xWIiT8mDhntFliiUz5+4UcmrKz0qR947kp7LGiW8lTOBNyVw3j9
tb3MubFhw9aSX/rV8SzDWqD48W72/KEpYCjVnRRYUdHsAWN/o6M6LD6q3MSn
iC1HsLHPnv9JQn5x9ECm43oF5vl7D5DRJLORa+4DYdHtziRVtlXsY6huBgIz
oaHJSoOFY684dv3vHijYKyoGOJv33sTvpMx8qEcXHR5gN6CLKhscq//hhyED
HrTASb/asVNoIC3w7rti0/EoMdct26zggJSx9YtWFlgwcGYCmbrF2FHoL/2E
uPed3WOfZCOqsgQC1OEgC4HD4bRpz15jiYDfGuzC+bFHv5sCwrf982Q/K6cD
W5iGIb9K8D9MOs6vO7s6C+TJkIrAp5/6RzG5eRTHkZMgBZYgnJHAhwHvFjrm
hb7ChzIH3u421ZsyNexRTEqogzbIUAZaJd4m/ICbNm3zGnMVu5cHrZdnyd/j
TrM0A/q7/nmqGeD5MBjiVXDPX3gQwS/mE6TmtPJQnOVOp4odjg44WWuuAWjE
PVpgYUMDPwobRIABnDjD18fxQ64BNVE/8qfFYGscw+ZvEscHasFWpx3+zC+P
BJuxZmw1M72MO1a+nwLB9xIkRVoCJIV6iK/r3xI7hoBgLHcKFSWhR4IPGvpd
KLWAF5QEOtbOXZavUmKvn/6QQAf7bkfmlUG/9hOi+DLhyhWbveFR3ciXySBr
ChZBh8qUHNqO8sujVY5sIvgSHCHNJ6bsairUIHVODNnyKCFHJ9BBzrKvv5Oj
XPl8Y8cu9QZJqYlohCJ8a83cZQXe7C57AFMDTwoTTuowxCezWCUcNhUh9iM8
OtbMXmpWB8FEYudkQuIVww3pxLDS1Ej1o0EhI7LbGDv/nAJBjKKq1ggG+uSw
YIKSAe0sfTFKbKfB975qwzkiktOJNBUwgq4tWmJmMILG5mipwBKP0/wojrQs
9QgV2hl4x/oJn35qhOmvfDrSY4Qr2nFJsfPwwIXurKB0kZ3sNvJwv0Lon5/R
EzhatCN2OkWaM1u5zvQ2IAfK9HI3A54q6VWYdXAX6Gda9+FyiUgjgnEd44ch
TR6noxmz7QeYpVVg9cWpflYMJxaAReD4oc/rsgI5lbXGFFDNnFkUcgM8W9ev
59SQuHVYYtfxSpnnz0FH24QtazcengKBXgGlXBUBD+AVhnaQHOEBRIaDXRRu
p3rK20Ftm9sT+/QLMTAHoFgiHk2c3acXhpaTE7yQQhKGohk3Q0klQpTcf5/5
6736U9WMM3gdDOPJ6ztG9R3PnzF9lQ4DzjpYKwvRGlfNXHxkCgSHe2xjptHR
nnUDBc0NUHtQd1QZAOsk55MhUU5w3wbHRothWN326ZcelxQ/o+KDajd9vsoW
+6W+n8usb6fZMBfWOjTn8xxN92z99Ej7pNRIXEpcVwnsrevdF9125nPRxg1b
vyticgimBzeVA4AE8uDfKO3/kgLXkQkixM3JUWZopgD5BEi+56eGhXL/KCcJ
2jZZvWOnbC1KX5xSaZ9+5VcHd0AewB3ki6gfcJBh153kB8ChiZoNehiI7CMZ
E3BUfM9vdprsvetxPHTB8ClR50aP2TpeOPUBqzck/+UfPB9irUiVDx/rasyP
5yAHhK5TItJi+iGf2ak8HCouZqoEUO47KAR7BTtKYVl4GFRl6HJrZysQXlxT
JEsE341SUGVrcc4BK3ZltbDlnlzYUavrh8KXzqrMDeTOCKiXTerEtrEYJ8aD
nxPMdX2if9TqxHtZWq2mbgz4JsEFmAV5uiTDx1Hrb9jzgEtiOVWPKJJowCwX
jCL3YFPIe0fUP0sBVfIaSURXKWd/+amrmuikPJmz6gwfMyR27fwqnr2meaV5
Dd8jj2Bj4KFmZmwK9tIdW6m0eDbavjTWPh3rEYTBbX4vb5pQ3lBggUxtPxSx
WRw6DMW/uZzG55oHzXHChGWeicQZ3Fh7Deu9EL3bup9VFTp8VfZoB4WsH30A
nsKNHKX+N1zw+GSQkIsdhbGtaUlbUsZpYEzXLbbCe/lKNA9sE7cXCjZOyNdO
vSf8BgXSEx6PIDwHpVzKMK2YVNr0G78kzG1gIt6EXwunj9OW6haG0BP9aIZQ
L9DR6Wf1HYOWI4/SMF8St2Txeq/OHGpAoV/h016+fKMdVSwRzgE1jrAfjhad
t2DutFYGYH7jGfQEtzFpeaLxif2mKVfgAkZJuB9bkhOOgc6QMESYJseBGwuo
4pPkV4Hz8OELPAOraoFIwAHXMCLGdPIpX33MXj/9PsH6YGvcaqPoIF3c3c7o
vZCESzbbMV564c8lkADytmzZkVb9oasZxnsf6qxZq727uapdfAL9sgjUf+6C
cgtPgzRjEb/ziwDF07p/FA4mi+AmD2xIj5foxefH4PVNbE0pVb5mEKJgOTyk
FXpoEWoQgafzeL8IXBGkVKJVQ8GcFqfAF1gJc5xfhFM/zbThHLAI/EWEeNkJ
fyCivzz0Pp57+ySXjEqB2rcbl1qqohVQrYi5QaL5iX5OPAJcI2CqEGrw9a1N
8NEao8DKkGP9kJwIHH8IIuPoXBbllDe4my+UtDp8BSvEbIkhI5fdKfVUkr4M
4p2KEdX2sln7LxWHzUAMINMdDfymsJXEEMXeC+wxzjsWOyOjGaGCkm2ATang
7kUXdglJm1Lw0DSoIcKyd5SWhk8tijLn66+N71I4yS8K5zRmKjydqWmxBIUR
I3ULLrAk4hi/KDRZuAyUpmg1i+L8w/FJo1SoHweech3lbtqrNddvMxbVpHEP
klLtU30/PZ4B2PSwR7uEhAQEFzglE4MeZr8sbGWxbVbTtCf8psgQVoMKI5UC
Xku+9NYNX/iyfosJK4Mw+BZXfxGdf+6bZm07BpdWOaTFoWU/2yYOIMl8Y3FE
+Ul9cKQYKyWl5shFm8bodurevxW2vNioOdj0aJRINCtEr5aH/kLIHCpFOYFd
eLeNdchSrcuxHkg8s+r0tX791rSzoRs+pCpjU8r7h5UDh8An7xi3dD5ICaGE
XMQzBxQ/zbNQ++VHnpdy2PDF4eLGih4YF4tofDaRGDUhAkoGcaArgKXKG/Le
FIahlwNeAsxDJ+PTKqS0ys8+2xSXJg2e613fcYgPEYUkxMqmtb90TMgVdkBI
j8sfWK6CxPkMlH/yTNJnkdlSkP2oGAScnRKaMw+xBXQF4IAfOApT1gs9TQB7
3botXp0qtfx8uA39Ojdv3n5sxqqnTV0RQujn+VWTjU7ZGwFTtEX6dcfQVAum
FEIb+p47eILJPxxxPyyMomVRkddLYASqxRhm2yErelM6S7GyxwOkRb4SvJAj
RCTOkZfKJchOIFcLsXbcXnhg9zdt2pbGKhV3gp9RdD937hpkOzQf8gbVQ4GK
2Is8nQIwhg0tYND/2SEnbip7NKF5QKcUILOxSCDMWTp1FJjeVN2zTIQbKUmw
TdDEQUdvJJkND4ebvUTTI/hQgvANYr7xPOoshfBX+lUly29+6xFGAwW4/TVX
9yC1+tcp4CjwT2dJaunJZCDoQoC+ig/Q000b4fCC07WmTFlhe6IovJRIqncu
88cAbZtOXmAJhHIMUOCdBepNlxpmuvhrE80YRkrAJFCX8A3u3tPzqmzMHmy/
EOYEkxAcmhXMCm2cE44pCOY4gZoavY4Nhb0ieZEsYJT3oQOmVqtM3JHUK/Pp
dx6rY8cutbTqxle9QzQwrbov3uTvWC1/m9ajQnYSUogh6dvQs8d02yCoGjW0
Y4eJ9iy/Q+XKpqarF885rCsjlQTrs86I45bX+GNHVAySBInccIAJTmkeWSJu
9fHiv2HsB0cvlikespFP9zAtHaOKAwD+iVJjSu5heif4ibkioSTOvvTorGkW
NgQLOpEVKJ64QLCjcGjB/t3B1g2ltLxQ2e9/eYyO8w2GIDLH+MrWFKoZJugH
CHYgbXr9UV1KTmcNjwqsv6Hvz7e8iT+c9po9S7UZVWcOixKj/Hz1ld0T8r+G
uZhwYiwePcNQwqnE5UQ7KxgYLQer2sJTyfIEDyW9VfgVQYIHDvUBAw/5do+f
nlrm57yq8XuPCCwX3BNkAC5atD7N/lReHWEesp1ZHGoq0Xqn/HzLQ0BSE4wO
lCofE0uZs6scWZKS724eB1eb6+w6053gDRoyTivUK7oZJohh/9esFnR5TVpl
xX0eA3A9da5Q7g4BIzRFKoYdt/tFyutaIpHHPr1n2FZSQUSEy0mabyd+BufJ
ElQnmKWqJmO595fN6g7EitSMy26y0L1kSZy++96QeZ6Dl5qedukl3Qx/Dlne
oMl5nRnRb3XucS4jvP8pgWfONpl4KLZ0EnQ8wLvec+LDD+TH8/4rMM0zdzpO
YRkg5ZPJn/lMvdiaA1n4BZyAjk2TasnXc/rmDH53boiecSKQ87xKYh49Wnz/
KlOkFFnOQOT+rg0tqHQ0i1zJhJWiKO8ufBVvrXqdovf8q1+9KSMOf3yvTktO
nB3uX8V4QfrpVZKIlQbNb3zvJGF5ULu/6ggLKjPMQq1SfSmaquPxwzmVNwtW
CdqP9jyFIWCjyf5cR/r3sGcxFzjjFElQeiaXlCqoFy78/H8XXq3Id05CKu/v
EsdiKiBTSYP0YLf99qme32I4DwwN7o7yC9koaSJx4mw4epsp3cadPvNKWn/Z
v42wHEs5XZOe4yLEwv9kpWLm3OlYJLO/+Z1xf0mFl2nzBfbABiKXe2iO8L/t
1QDPxlTYiNwusKU6KHwTdfw8Sdd1ObC4P0sMM+dOxyLZbOrZdorHBvWY6Odg
iswW2HgsPePeFqSxYx76Hokl/qyGVBs6xUrHX7J4vQIZkLec5+VA5P4sC8yc
Oxfqe/3qICEQhg52WuLgku/L9xjGlBr/IEGOUveaNTUxKm6WM74yuubMXuOV
9g60LFfJAjhu9diH+ZCa5qcR5MqZ79tnps//qWKNBtTohwoNZ8L7BMaqVpqu
tmReucI4TbuUTxQhFRiqGTFioYxtAsOMQruDB4sDXmoFADbwNMpWUL6Ktqwa
Yfr+Y0GrqwJbQ0XmxRfFGrUzU/C4pbFtRQ5N57yht+1mv76zFJ39fN2XAX1F
LkIyfN7ctT57N17EoEFzTN7zG3UvvtwwHBGsDdgd1pkWOm3qirSuE4ocsnnc
dsaae7wzTWF75RCOGbO42EVQx6FGpLE7+rvGLah9QE/E16nav+fajLIdwjOE
AxMpxrMUJhB1JTRI6jrtSIDUMrIfHsri1AIxBwalaoM0+DXvvNV1ig/HhXzA
+fPXFbs6Ot46KOxTQz8b9x6wRURhnXFX3ZdL4ZDBkY0G6ER1WjGSIm2smfa9
PgFOCSs4zviKTPwi4QVbF10Y946O6/1qWAkVGCZH46PhCw3LVAHDzykdoZ9F
WiXmrxLQqsdqxw6TlJhEcwS3fPtUJLQUjLu9tk9pdyv/OgGBbA+Hr+MSi3Xn
pCIQcMwc5ZQJwX8kIGA6KMwxAWXzUdyqbLB8ENgvmhengLNGy5xX0S9MADrG
wnwcL1OeBezLbYl9uj8MnhHIkW8OZqas97RL6xWP0n1McG5niCvhlQPZru3Y
YjAuCNR/g7h5GgSKDVGUSMzPB26Vk60waM+e00v+XBwE9JJgCCdOUyFQVGbV
qs0mvok+rlyxWXlBc+asiYOUHy0sFgLFwtwMqRAozEGePeuna67bEaUejPSF
zgQQi4QAkaJc7jQIFBLAK81zNG5yclQ5HKH79Lovi4XAaS8hLyINAjFxx89t
qjdenxhSSMjfRYFAshc5PYdKOcpp00tCwvg49GQ0KaeR7Fgx2SKnxyJRrkHZ
05eaV0b3ojlqUKUbjkIxziIhQFNXQkAaBMd7CNas2WKnAPntzqTyYMnrFOMs
EIK7/RCqwMFHfWvKc/Jb0zgMxkPzAbfjylEGJS96FlokBEpXo+tTGgRKKFMa
PZAqk1SctPvbU0seyjO9/aJJVaJLrWzapFLx6Y8GkbvlhzYRYp5Oc4z5bz6+
r0nVt5J7ssqetNT2hJ0m8OLsP1VWiGU6tTvvYsuAgBACB4mAdRoEMhfj3W5v
xKXpxS9JZityetK1nNVin9Kml9kvdoNar8zdJLMscnpuZnTKYZ7pS800uufu
QUbcS5duUOlTklkWCQFORhUepEFQ10OAgCVREiPVaT2qRYNfOkZUEQgwV1o/
M7JMCNQDgggnW9DtrSmhVNb4pQ+uFjj9XX4IuiMxLIn0t6U857OizU2H2oOX
VhnxMEu1oChy+sWL19v0xCTKnr7UbCIidNhADv9H+tfFLCns/0txEKhMHN9N
GgRyLdN12zhb92mhkh23hdohFDm9MlLBa9r0ihiQWyj8H+G/U6sweGaR0+se
NzysZU8fyzvyQui1Euq19vTFgoEWCQF+Wney7FMaBMrNFK+GUjQ9kSt4AuGB
jOklMomZJf4nCOAsqr0uG4KYC+HywPHkCE+tamg0Lx5gEjcpfZQhQFooQwIt
9h3ME5HHXUzN/EAkmKhvVhocisuRDwjjpfOFgMARpuByXkwcbCKDX0lAwU+i
pjCIGFgcqiRKLuEuD4qxP7dZZYKn+KtTn2w8/DcqqV6T6ETx13x0ckjyJ7tt
Ai8ItjmnX7ky6D6gHscCXhO3Ozf62dThxRF3ye1lQhpfK81mkhjhXg89wEfv
sYTyApv7UxxtqWmGIoyF3qFkY6lNFh5CDGYcJY6l6V5t9blyX6cCrFRedc4g
IitoVX9KWn+B0Fb30CLc6JWNBg8+mAKIicrRQRpCdQqKbg50nDKuahm/LAPa
+LzQaQBLcePGreo9QR6O/BMFAnywB5iRsQJRS/Cm4TcAImQp/UhxIbotuNLP
qMCuEzapAPs6T0s6IKr4lDtIyYsBnZmXD9r9nyRTdGP32MlVw5yh2NC0b1K/
W9gVkh/h51heuBWDdluqJUnDqHJeuLNKlw2oCzJnA4LldD6cH6MHUh5MgZ3Z
D/WYxRELNyB8pMapsDFIFyUeJ7ESYXH0OjK2T1mY3bBhq5E7QT9hVpe30qSi
HJj9e2utHhpi7dptXlu80IpqhI7AHy20qK+Sf/Huq+V6GjqVGycDh7YB6iRF
qpf8URno3N/5L8X1Ro/5WzVTJpBSokuUbE4trWuVpI7/XsHlspFYai51fLAY
VW6T1PsNDU/pq+XA44HRNLu8DcxjDWdPO14CGOjp4BLzDPXEKc/Kn0YBUL+Z
NFT63hQhCIbBLTyipjiN5+8Lj9nNxw/z2EMTgxMqxIV8wLGFKuYwq15GKmUh
rNw0BQNK9cVKAdPE6tU3lA1CnHFn5UMBgAMkm6a4huHqIrF82UYzKHBIQjQw
RRga/S1EUKqPQXXMQt1s31SfOkShTnmcqKePFC5XDtAm3z/x2EM6k6SBqsrK
IRJueKPAWNKWMgjl66VhTwSKBs3p57Y93b6mgmXsoXJg7++kK/fP/Cs4i2Hw
rJpTS/h32tQVsljUt56EvCzEEVViCMJ3Qhxd1JHrNLgsB+L+Tjpn6wYpFDju
odeVrngW0F6I3sjnoZZ88LE07J2dwJT6zfv2fHbZjLwaGdhLM7N8ArHZ0XBl
bCDAIOBGqBwQ5RbD2j/91PgK5iwQ4fAslbiwQORSSC4rKAJE3WcDV8c4IZ/G
QLyyux1mQJTnlJCBE8Vlgqhu/NL3SEzS7arskryuLQoD8ZeeNvYGkXIrd0Lq
+hnIsFMiWNnQxdkcVLWgVLj/1qXEygd3/y4UwF8nAMSRRy88hsLgW7wnU42v
lVSVBuDeBZyU4go6hzb7jtsyC4TuWA8d5iI+DiK+mJK4rJwpqqwpZZq90n58
yR0pw0jxIjaII4KjpOuq4dbn+/hBgdCp9QX0QXUXkEHeNEfb04MghNKogMiC
jug5q2j78tgA3dUJd2CB0KnHLcq90upQrUj5UQRZYon06TTQ1OGZIAtcmOek
JqGH8R1l5Y8WBtrvPOJw4eAFw98E8sjc2rZtp9JouHlXvL9s6EoN/+AdJ6Yj
BV9CGEqm6WWaD0D75fceYwQUMH/YDNieszNDS6ZkpnIaTOok7ia1yTt1nKSa
RutbznfckPdImDUjRHmCB4gewoQdAQh3Nw5j5buQCOLsjTIBkv1GmgLuUvxf
1T3iOFaENKEyxwAK3MUTPXSMglcDvQl04TJVY7ykn6dZmdCVWsogVifmjYNE
dy9jlPI6gd6WhUGn4g14I2uEvrDDEArKVyGgKmafBp2ucVDYBwHgxbk1amZM
+nM9VhhodTxoeLVxBKtuhCx2JbLobgyCKmmgybdBSAhllT7KpR6Z1rr2ul5m
zjtWVyB0J/lREPr4VBHQTvVVggmJiuTzgceyASs1bz/sGnE1edJymdLKF+BQ
tioMMBUGcK5JToMwaMkKJMr+IG9YtzCkQRfXfle3YAInB1BQkuDRgth9frww
yJRsj37zkK/PpLs5ATVlhcDjFJxJgyy+BK2asVXsNvx6qm5hg7mcr0Cw1HQP
NouzFoSheJBYrVwNcsksaazH9JI7U8a4JMFjcUsCktNInygMErXV4K4uAsaU
QbptVO6G9BfCr2lAyOPMpSNkGvCpQFwo3EOaCLdPAoXTEZVG0d2n6nIE0yBY
Vdhsin4QVEHmgnlH+cqawLJ09GCfipnNflGHVhR0DC7Af/21CcqLoHe1Ow1l
zrEiDJQhgNTKH46CDsAcaDrEWtW+G+oCm+RsNq846nQ7E8MSd2JKCr5JjjnS
Tyk/xqCBc/bFlGd5TGJpIcQQlY7MjvCzKacQw2ofzKb+5NwAw7B+NoWzFcpD
xu+D2c5On00hY1LDG/m7s/bBbOekz6aLhJAFToHaV7M13Gs2bDxHJgqK4omG
XTk1UmH+isymWzvI6odKkNzuHChGmNQ298Fs0rxVNOkUrBCOVAQY59Q+mEoq
K1fAMixev+8nvsMlw4HfZ1PFUWEqtn3ajmJXSGKnxNunfTCblEu4F20T2LWB
A2bLlQ4Q8j7cXcRs9suFfg7E9p3NBigbXRcFKUzOZelpc5SLC0ul2ivzXJ5Z
NA6sWjhkMetIVZG+1h5M/hp6UGCUYaXeU/HZLvazUXKFmo4dt2zZRn/ZRqgN
ps/4PpitkZ+NlofqnOVMD3kq8COgmmBx31vx2aQ7YdiDRUJiTiOQbY8XTRVG
+2C2S/1sdEZBYaQ8ym3WN/xsaMbYkJzkfTDbZX62qVNW2Now/Fat2ixTl2QW
9Ru7r7DZ9n8WRW4vkdhyqRJ17NgxatCggQVomjRpYkHVtIDM5R4z1OGhxIN0
p7z7Ii2z34n7U8lUJGb2XzZEbtsPBcJmzJhhWNm9e7fM4xzwr/AYIUOLBDJ6
TDgCkZFPzgm2P72Q03qilAMj+yOLQQG6WbNmRQItbfVXevohTxQjEO5N1gKp
NLLYcX8QQIEdTJu2skAs7O/kg9Kcn9u0aVPVZx01btw4wJpWDHuV/T3YmD+9
zOD0EAj2LwKV9C5sJfxW/MbVemklveXAzf4KhNfI+Xnbtm1Rw4YNoy1bthgp
iXbsKFgsOkZJdVs6mfFQBChRUBJH25IlG9IqzJaHQaofgGu39U6aNClq3rz5
19e+11KaeLpCfcJpjxbs9JGHitv7AyOgr7PSsmVL2InWn7b8tKs40i4s2dfP
HQAx/MoeS1BOo0aNCKAe4FjaH7F6UeySJUui+vXrl3mO0u6sS7u6rbzP/V97
90CH7wDIfpAa1Lt370w16H8Lyv+fFDNJ0WxFI4Rko6qYefHn8PJ+UXLQfwMR
m8xc\
\>"]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"Re", "[", 
      RowBox[{
       SubscriptBox["\[Rho]", 
        RowBox[{"1", ",", "4"}]], "/.", 
       RowBox[{
        RowBox[{"NSolve", "[", 
         RowBox[{
          RowBox[{"eqs", "/.", "params"}], ",", 
          RowBox[{"DMVariables", "[", "system", "]"}]}], "]"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}]}], "]"}], ",", 
     RowBox[{"Im", "[", 
      RowBox[{
       SubscriptBox["\[Rho]", 
        RowBox[{"1", ",", "4"}]], "/.", 
       RowBox[{
        RowBox[{"NSolve", "[", 
         RowBox[{
          RowBox[{"eqs", "/.", "params"}], ",", 
          RowBox[{"DMVariables", "[", "system", "]"}]}], "]"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}]}], "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"d", ",", 
     RowBox[{"-", "60"}], ",", "60"}], "}"}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotPoints", "\[Rule]", " ", "50"}], ",", 
   RowBox[{"MaxRecursion", "\[Rule]", "1"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "60"}], ",", "60"}], "}"}], ",", "All"}], "}"}]}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.522542762323259*^9, 3.522542764447953*^9}, 
   3.522546267445561*^9, 3.522546298673073*^9}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwtlXtU02UDgEFHAg7lNvhxGTE253AyNrYxA937wmSwDQo+QpSUAJG4JF5A
kWOhXCYjChMEDC8YwhwLEQ35KJP3VawcyhKUo+JYoOAngspFFPi4VOf0x3Oe
/59/HlrCrv/sWGJiYpLwN//Y/3HmyOJiGK7/Zd5dW/seMDZa74EpYZjEz8+7
VrQM5OY1vD3cHYp/bOEWmuaaA7PO5fF23qE4m9UabV1rAYoc027XFSvwSzEr
eeUTS3D/pey/5Ek5bnLVjxxgkoHbTda5jBg5zq9gjjxXWYGSzBDDw5syfK/5
+m9bzFaCeUYyRewtw37vT25PiLQGH4Yq+FNfh2BiIHHuVJgNOJPBiWh4EYwX
kjbclifaghbmsTFRZDDWP2i88VRpBzofTR5tvyrFWQNXPuHesge/LXW2/95R
ilm/j4/1fOQA2rzgiZyMIGzsrJSYWhLAwDufEtO9EdckUZbMTTqBGeEKf1/e
RpxztWsizsUV9P365ptnVRKcO0bG7pvdwPWox/0VJAkmdhk+k+xwB0e2Konv
4gNxl0A7R+hpQLHdO/xkWwDeX/u0jeDRwYOpzRZ7mQF4bneFWZofA2xX5d0I
KYGYNbotn6VkANaLR8f8JwB+tlBIZDLpoKi98czQ72I8ba84XsqigoavLPfx
Czbg6betjrGpi+K7EUmKXLgeSwX3vC72EiiXe51xwNYf+5BGUckLNvJZ6bqQ
PvgB/slzIGHBlo/qHs+ftTi9DmcqeKHr2UJEaP6U1EaJcBXFUKCyFKIeliTX
394XV44fstvxngCV1qvbursFOG7w/6OaDh9kCLBYtD3Px6JiecXRJh5a3ZsG
Ig/6YKVxvLa+hYsEcVEykS8Pv2pwydOe9UaH73yp9n3ijUt273m73IODmpd9
hMPzvXB8/EW3xndr0XCge28qg421nCOXfppmo5oti0fulLIwublaG0xio092
G/kcUyaOFYbpU13WIM1rTCk6SMctTtOj91Z7ojfp5949nXLHgxSX9DoSC7X6
mSOaiIrrd0YMOaqZqLBSV3Iwj8Cp1KbVs90MFPuDuiAx2Q6XpZeXmjE8UFrO
sqmqU2S84viaGpULFRUm/DzdTFmKqbqXl/b+YI+k3/ROn6ZOIv1YtkPiUTNk
vYWuSRJPAPNVGcuFVSSwotWJ6xCwFMrVSlp6jz0oa6eRZxKs4Lnwn6V7Wqkg
qoh4rkuzg17Hd7U5Mj3AANdd7sgk4O2lDVbu9QyQRCMf6QmiwqxYJXdeywRT
5nP9Oh8a3KxS6Lh0Fjh34pXuYgsdPgwKN5QCTxDBGrhc7seE9Zf7yl96rwGD
qp5yzzIW3EfEG5NJbLB/WHfg2sgamDobx61+wgbK9KZbFlleMKsi9zv6xFrQ
ka0gs025sJjPf/9DPgfcrAU8Fo8H1RaSXxhV3qDtD/6mVfE+0KRPiEpquCBd
qu5fVcyHZ8TjlrpKHnC7RqQxrwjgg2TvEM0FH0CWktbmTwuhqyzKxDjCB7P6
16N/+otgtfO86P6sAMRAg7YsZR3UKKoWtgQIwdVLt5KDKz+AOo3RrrOXD2KE
nG3RRj/YUGAbZs3kgOnWsojPVq2H/XstVXe7XP7uPBuo37QB7u4c9bFSvWpz
OhYnEKrEMC5DPLN2iTPqDwkqrv4RwMLjDX1WSg8UaUIZL9gF4cK3RaeHUhgo
aqhd9NwmAF65MWxdEslA0R17cxTNAfCEaH6CIqWjUnqNpFcRCDXLkytvLNCQ
/osu85SngTDWNNew54I7Yn26YbiuUAJlXQ99lXVuKD+gXveEvRF2fT5Db+a4
olbZ/a20zo2Q6oQ+7iM5o1cRi68/TQ+CldX9qosEgbbqtaYdJCm8rLb9avCQ
A7ot32QnOCmF0TmMGq4tBW0OF8NQcTB8FBY7trXJDj2LYu5MHAiGpzzFEq99
tujCHyWEJjsEGofMrbg7bVCm7F37iI0M9smKYrQ51shh2M2Dd14G37DOHjKN
WolaVdLD+9fLYcHQtZrPu61Q9iXtTIpRDg3/u37n24/JyK93Rca2wwpYZTlT
TsxYIkHWZOoXjqGwyTThwbq7Fohj/zDhZGMovKW2+bWjzRw1ZzHPh0WHwZb7
Ts5c9TL07w/gvz9AfwG21tEN
     "]]}, 
   {Hue[0.9060679774997897, 0.6, 0.6], LineBox[CompressedData["
1:eJw1lGk41HkAgAk5k4wZ4zY004gwY5wzzf83jmHMaB07EkWOtUlRsaUlOVqy
SqVLh47JzVakIvX//VrUDqGsDkeiQkOhpKiw7fPsfnif9/v74aVExPv/tEhO
Tm7jd/41uzdxbGHBB3Xk7hu5Xz2E91/W2g5ifFBF6Ejm8cARPD2j6lNapwi1
BZm7NljJcKU29XCCjQj567trh8aO4Tm6sa3FuUJUe+JsfsbsO7zrneCmxpQ3
suttra4ImsSNm+iXEoK9kT2z1ilk6D2el+jV96xJgFYH1ydekEzhc8s3Ebk2
AmQyr5dYL5nG14iEdtMHvJDbp1DNJxYz+LkEa7+qUU8UGfbex4n2Fb9BOzLp
GOCJAqvweGLCPN7WPXWosYGPVhw7vynARh7eU9DXuajLRztuj0dLkALEV4GC
1AQPtFCkfiBrdDHsY5TGBHe6I991t65xpWpw1l6T7cBwR3k1kni1yqXwefPH
g8On3VABRxxdWa8N74p7B04ouqEI/4en38aSYNb638inwl1R3jXRmrgrBlAY
aeN7BuchCelJrHcUBT6dDlLdQeOhN8wI24t9VBi5P+NPrzyATpB78zl/W0L6
aPcR9gcMFR+JDD9QzYA5jZfPDd3nolt7gfw5SwdY9bvaL3b7ViOu1m1Dk+su
8KFftDAdcFD2nqMT0XEcmG57d3mSNhtFeLYviux2hsylhvNxr53RB6PCJUu3
smBx79wF1UIn9MCSRQG+VpBc9sKtSOyI/Cer0jeuMoOP6W7pbB0HFO+pl0Ls
IsL88hK8s5OF/MupXNWG7314qgvapXbotW52mn7IJL6iJxYLSGaitPhvdwbJ
PVzWRrHA0YGBQKVZ85bJL9y0B3tKHF7aoCSR2DEoZDEmczXt2bzcErXISPWa
3dqYZN1C1oN8OlJV0au9qK+LhWzrt7OWpyHxr6JgYqUBVjaBiDnJ5qh5m0Wz
jtAU+xh36fOraVO0t/eY5aCCOVbnogIpjkZIN9PwJamNimWflOYlZ5DRYUua
+l4pHQutLNkXtYmACr/sHEUmVlhsqvL06bMaqO/x8VZTXxssO+LWTC1RAekZ
OZQXyGwx/sGemUKjKciz1le/0cjAUNYwaXT4DBwUUxSylzExrXXmZdHcD9jL
S+NBB1oYmGadni2JpwAmTbfvLnltix1tpGjMRiwBOxmSuWSiDSbOIb+RxhLA
TL3BjVwjK2zQ1tRbl0YGxTqzWwuULLBoikbWYw8jUP5MclCtiYpNq3wbkDIp
gDQfY2wCzbBLBePSKzfMgZ/A436DpinmRx+sOe5CA1uuOy2ZBgbYTpk06c7Y
SkA837FnnEXAWnYLNSzlbYE1dUmxbEIJayrCGHQGA6TkX5/1ZX/l4h12gdRw
JqDLEVomCge4cfySAWquHZg8xMmJlH+LG98hx9Kus0Dqqcz9KoZKUIOvaJU5
Yw98nheLt0h14Jf2ibcv2I5Ak/BjUlMaBQaDvoqjMU7g5vPuu5TslbCh+q9N
niedgbLyhpi4ISYMtrfesLbfBUT1Tj2tVXeGM3VH/X6mckBLalsIJ5QDVahf
XNsDV4PhJDK2x4YN9Y5sZNnv5wL2ibBDZjRHOODlkXv+GgbkQm7P8U4yYYAc
8f2+eACKCqOBDncVFA81Or5ZxgOdxm01vHwaXNuyI1VYywMmRvfab5dSYL65
xK1H6AqinM6aGWQawvaURyoxr1xBJxi+dzhMF9LDVsuKs93AIEmnZDGNADN5
5dKXlu5A2T/uHN69FNYJutZT2txB6ByvjP9ODY77LUyExXkAbZ/wq7cWK8P1
7RXyLYp8oNT5yl+LqAhbvQMJrDN8sGzlD05bS+VhkC8XiLieYFaxxihxah4f
FtO2Rg16AlG+fJV64Ff8j448ctluL5DiU26nnzqDJwo+N44tEwB0UKlLXesT
TpIZmzFKBUA6Yuodv+gjXrefn7aT4w3keIS0R84f8N3VFbMx/d5gpl2/yiN3
Enfp0UzYkCYE5JTmls8W4zhr19TmFF0RIF/d0p/XNIZb6zyLOHNZBDqY3MLW
Ihleu4tW6rPWB2wzUJvqyh3B//s9+P/3/wBZ5dLn
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 0},
  Frame->True,
  PlotRange->{{-60, 60}, All},
  PlotRangeClipping->True,
  PlotRangePadding->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.522542775492355*^9, 3.522542960547408*^9, 3.522542995770412*^9, 
   3.522543035134631*^9, {3.522543296264867*^9, 3.522543319510829*^9}, 
   3.522543384500983*^9, 3.522543455525264*^9, 3.522543514229852*^9, 
   3.522543626701997*^9, 3.522543667475726*^9, 3.522543726018032*^9, 
   3.522544871173848*^9, 3.522545966092067*^9, 3.522546271481462*^9, 
   3.522546301931841*^9}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{960, 1029},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
ShowSelection->True,
Magnification->1.5,
FrontEndVersion->"8.0 for Linux x86 (32-bit) (February 23, 2011)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[567, 22, 32, 0, 109, "Section"],
Cell[602, 24, 66, 1, 63, "MathCaption",
 CellID->836781195],
Cell[671, 27, 151, 3, 43, "Input",
 CellID->2058623809],
Cell[825, 32, 1762, 59, 121, "Text",
 CellID->525777075],
Cell[2590, 93, 68, 1, 63, "MathCaption",
 CellID->429217524],
Cell[2661, 96, 1155, 29, 179, "Input",
 CellID->433132487],
Cell[3819, 127, 689, 25, 63, "MathCaption",
 CellID->133602844],
Cell[CellGroupData[{
Cell[4533, 156, 1586, 45, 125, "Input",
 CellID->534530029],
Cell[6122, 203, 2005, 57, 124, "Output"]
}, Open  ]],
Cell[8142, 263, 217, 5, 86, "MathCaption",
 CellID->462076121],
Cell[CellGroupData[{
Cell[8384, 272, 528, 16, 98, "Input",
 CellID->494599775],
Cell[8915, 290, 2450, 67, 124, "Output"]
}, Open  ]],
Cell[11380, 360, 76, 1, 63, "MathCaption",
 CellID->358620443],
Cell[CellGroupData[{
Cell[11481, 365, 464, 15, 71, "Input",
 CellID->167259034],
Cell[11948, 382, 1077, 18, 177, "Output"]
}, Open  ]],
Cell[13040, 403, 102, 2, 63, "MathCaption",
 CellID->577766068],
Cell[CellGroupData[{
Cell[13167, 409, 2629, 57, 233, "Input"],
Cell[15799, 468, 3242, 90, 172, "Output"]
}, Open  ]],
Cell[19056, 561, 724, 27, 42, "Text"],
Cell[CellGroupData[{
Cell[19805, 592, 247, 5, 43, "Input"],
Cell[20055, 599, 2272, 66, 106, "Output"]
}, Open  ]],
Cell[22342, 668, 384, 12, 64, "MathCaption",
 CellID->610306692],
Cell[CellGroupData[{
Cell[22751, 684, 272, 7, 71, "Input",
 CellID->645617687],
Cell[23026, 693, 1053, 27, 120, "Output"]
}, Open  ]],
Cell[24094, 723, 92, 1, 43, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[24223, 729, 53, 0, 109, "Section"],
Cell[24279, 731, 98, 1, 63, "MathCaption",
 CellID->690131918],
Cell[CellGroupData[{
Cell[24402, 736, 230, 5, 71, "Input",
 CellID->718931880],
Cell[24635, 743, 672, 14, 71, "Output"]
}, Open  ]],
Cell[25322, 760, 390, 12, 64, "MathCaption",
 CellID->854192725],
Cell[CellGroupData[{
Cell[25737, 776, 459, 13, 98, "Input",
 CellID->465762594],
Cell[26199, 791, 1836, 51, 146, "Output"]
}, Open  ]],
Cell[28050, 845, 76, 1, 63, "MathCaption",
 CellID->314466782],
Cell[CellGroupData[{
Cell[28151, 850, 424, 13, 71, "Input",
 CellID->298399236],
Cell[28578, 865, 32297, 840, 580, "Output"]
}, Open  ]],
Cell[60890, 1708, 38, 0, 42, "Text"],
Cell[60931, 1710, 117, 2, 43, "Input"],
Cell[61051, 1714, 2257, 48, 98, "Input"],
Cell[CellGroupData[{
Cell[63333, 1766, 1777, 45, 125, "Input"],
Cell[65113, 1813, 28311, 470, 339, 4676, 83, "CachedBoxData", "BoxData", \
"Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[93461, 2288, 1347, 39, 125, "Input"],
Cell[94811, 2329, 5076, 90, 336, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)