thrustfuellim2.nb

(* Content-type: application/vnd.wolfram.mathematica *)

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

(* CreatedBy='Mathematica 12.3' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    190649,       4129]
NotebookOptionsPosition[    186062,       4063]
NotebookOutlinePosition[    186505,       4080]
CellTagsIndexPosition[    186462,       4077]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.939752165554883*^9, 3.939752168886167*^9}},
 CellLabel->
  "In[785]:=",ExpressionUUID->"46a113b6-085b-4448-9759-458b764a7e5d"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{3.9370534187414293`*^9, 3.937160323218102*^9, 
  3.937161003812353*^9},
 CellLabel->
  "In[786]:=",ExpressionUUID->"21c1531e-bfb5-4ac1-a49d-5b0746559004"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"getFinalEnergyPos", "[", 
     RowBox[{"tt_", ",", " ", "tmax_", ",", " ", "nt_", ",", " ", "delta_"}], 
     "]"}], " ", ":=", " ", 
    RowBox[{"Module", "[", 
     RowBox[{
      RowBox[{"{", "}"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{
        RowBox[{"energyform", "[", 
         RowBox[{"tff_", ",", " ", "na_"}], "]"}], " ", ":=", " ", 
        RowBox[{"Module", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
           "c", ",", " ", "r", ",", " ", "sig", ",", " ", "gamma", ",", " ", 
            "nmat"}], "}"}], ",", "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{
            RowBox[{"c", "[", "n_", "]"}], " ", "=", " ", 
            RowBox[{"{", "\[IndentingNewLine]", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{"0", ",", " ", 
                RowBox[{
                 RowBox[{"-", "2"}], "*", " ", "n"}], ",", " ", "0", ",", " ",
                 "1", ",", " ", "0", ",", " ", "0"}], "}"}], ",", 
              "\[IndentingNewLine]", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"2", " ", "*", " ", "n"}], ",", " ", "0", ",", " ", 
                "0", ",", " ", "0", ",", " ", "1", ",", " ", "0"}], "}"}], 
              ",", "\[IndentingNewLine]", 
              RowBox[{"{", 
               RowBox[{
               "0", ",", " ", "0", ",", " ", "0", ",", " ", "0", ",", "0", 
                " ", ",", " ", "1"}], "}"}], ",", "\[IndentingNewLine]", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"-", "1"}], ",", " ", "0", ",", " ", "0", ",", " ", 
                "0", ",", " ", "0", ",", " ", "0"}], "}"}], ",", 
              "\[IndentingNewLine]", 
              RowBox[{"{", 
               RowBox[{"0", ",", " ", 
                RowBox[{"-", "1"}], ",", " ", "0", ",", " ", "0", ",", "0", 
                ",", "0"}], " ", "}"}], ",", " ", "\[IndentingNewLine]", 
              RowBox[{"{", 
               RowBox[{"0", ",", " ", "0", ",", " ", 
                RowBox[{"-", "1"}], ",", " ", "0", ",", " ", "0", ",", " ", 
                "0"}], " ", "}"}]}], "\[IndentingNewLine]", "}"}]}], ";", " ", 
           RowBox[{"(*", "Good", "*)"}], "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"c", "[", "n", "]"}], " ", "//", " ", "MatrixForm"}], ";",
            "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"r", "[", 
             RowBox[{"t_", ",", " ", "h_"}], "]"}], " ", "=", " ", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{
                 RowBox[{"72", " ", "t", " ", 
                  RowBox[{"Cos", "[", 
                   RowBox[{"h", " ", "t"}], "]"}]}], "+", 
                 FractionBox[
                  RowBox[{
                   RowBox[{"154", " ", "h", " ", "t"}], "+", 
                   RowBox[{"48", " ", 
                    SuperscriptBox["h", "3"], " ", 
                    SuperscriptBox["t", "3"]}], "-", 
                   RowBox[{"384", " ", 
                    RowBox[{"Sin", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], "-", 
                   RowBox[{"27", " ", 
                    RowBox[{"Sin", "[", 
                    RowBox[{"2", " ", "h", " ", "t"}], "]"}]}]}], 
                  RowBox[{"4", " ", "h"}]]}], ",", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "3"}], " ", "h", " ", 
                  SuperscriptBox["t", "2"]}], "+", 
                 FractionBox[
                  RowBox[{"6", "-", 
                   RowBox[{"6", " ", 
                    RowBox[{"Cos", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}]}], "h"]}], ",", "0", ",", 
                FractionBox[
                 RowBox[{"8", "+", 
                  RowBox[{"12", " ", 
                   SuperscriptBox["h", "2"], " ", 
                   SuperscriptBox["t", "2"]}], "-", 
                  RowBox[{"8", " ", 
                   RowBox[{"Cos", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], "-", 
                  RowBox[{"24", " ", "h", " ", "t", " ", 
                   RowBox[{"Sin", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], "+", 
                  RowBox[{"9", " ", 
                   SuperscriptBox[
                    RowBox[{"Sin", "[", 
                    RowBox[{"h", " ", "t"}], "]"}], "2"]}]}], 
                 RowBox[{"2", " ", 
                  SuperscriptBox["h", "2"]}]], ",", 
                FractionBox[
                 RowBox[{
                  RowBox[{"23", " ", "h", " ", "t"}], "+", 
                  RowBox[{"6", " ", 
                   SuperscriptBox["h", "3"], " ", 
                   SuperscriptBox["t", "3"]}], "+", 
                  RowBox[{"42", " ", "h", " ", "t", " ", 
                   RowBox[{"Cos", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], "-", 
                  RowBox[{"56", " ", 
                   RowBox[{"Sin", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], "-", 
                  RowBox[{
                   FractionBox["9", "2"], " ", 
                   RowBox[{"Sin", "[", 
                    RowBox[{"2", " ", "h", " ", "t"}], "]"}]}]}], 
                 SuperscriptBox["h", "2"]], ",", "0"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "3"}], " ", "h", " ", 
                  SuperscriptBox["t", "2"]}], "+", 
                 FractionBox[
                  RowBox[{"6", "-", 
                   RowBox[{"6", " ", 
                    RowBox[{"Cos", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}]}], "h"]}], ",", "t", ",",
                 "0", ",", 
                FractionBox[
                 RowBox[{"2", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{
                    RowBox[{"-", "h"}], " ", "t"}], "+", 
                    RowBox[{"Sin", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], ")"}]}], 
                 SuperscriptBox["h", "2"]], ",", 
                RowBox[{
                 RowBox[{"-", 
                  FractionBox[
                   RowBox[{"3", " ", 
                    SuperscriptBox["t", "2"]}], "2"]}], "+", 
                 FractionBox[
                  RowBox[{"4", "-", 
                   RowBox[{"4", " ", 
                    RowBox[{"Cos", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}]}], 
                  SuperscriptBox["h", "2"]]}], ",", "0"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"0", ",", "0", ",", 
                FractionBox[
                 RowBox[{
                  RowBox[{"2", " ", "h", " ", "t"}], "+", 
                  RowBox[{"Sin", "[", 
                   RowBox[{"2", " ", "h", " ", "t"}], "]"}]}], 
                 RowBox[{"4", " ", "h"}]], ",", "0", ",", "0", ",", 
                FractionBox[
                 SuperscriptBox[
                  RowBox[{"Sin", "[", 
                   RowBox[{"h", " ", "t"}], "]"}], "2"], 
                 RowBox[{"2", " ", 
                  SuperscriptBox["h", "2"]}]]}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                FractionBox[
                 RowBox[{"8", "+", 
                  RowBox[{"12", " ", 
                   SuperscriptBox["h", "2"], " ", 
                   SuperscriptBox["t", "2"]}], "-", 
                  RowBox[{"8", " ", 
                   RowBox[{"Cos", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], "-", 
                  RowBox[{"24", " ", "h", " ", "t", " ", 
                   RowBox[{"Sin", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], "+", 
                  RowBox[{"9", " ", 
                   SuperscriptBox[
                    RowBox[{"Sin", "[", 
                    RowBox[{"h", " ", "t"}], "]"}], "2"]}]}], 
                 RowBox[{"2", " ", 
                  SuperscriptBox["h", "2"]}]], ",", 
                FractionBox[
                 RowBox[{"2", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{
                    RowBox[{"-", "h"}], " ", "t"}], "+", 
                    RowBox[{"Sin", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], ")"}]}], 
                 SuperscriptBox["h", "2"]], ",", "0", ",", 
                FractionBox[
                 RowBox[{
                  RowBox[{"26", " ", "h", " ", "t"}], "-", 
                  RowBox[{"32", " ", 
                   RowBox[{"Sin", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], "+", 
                  RowBox[{"3", " ", 
                   RowBox[{"Sin", "[", 
                    RowBox[{"2", " ", "h", " ", "t"}], "]"}]}]}], 
                 RowBox[{"4", " ", 
                  SuperscriptBox["h", "3"]}]], ",", 
                FractionBox[
                 RowBox[{"3", " ", 
                  SuperscriptBox[
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"-", "h"}], " ", "t"}], "+", 
                    RowBox[{"Sin", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], ")"}], "2"]}], 
                 SuperscriptBox["h", "3"]], ",", "0"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                FractionBox[
                 RowBox[{
                  RowBox[{"23", " ", "h", " ", "t"}], "+", 
                  RowBox[{"6", " ", 
                   SuperscriptBox["h", "3"], " ", 
                   SuperscriptBox["t", "3"]}], "+", 
                  RowBox[{"42", " ", "h", " ", "t", " ", 
                   RowBox[{"Cos", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], "-", 
                  RowBox[{"56", " ", 
                   RowBox[{"Sin", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], "-", 
                  RowBox[{
                   FractionBox["9", "2"], " ", 
                   RowBox[{"Sin", "[", 
                    RowBox[{"2", " ", "h", " ", "t"}], "]"}]}]}], 
                 SuperscriptBox["h", "2"]], ",", 
                RowBox[{
                 RowBox[{"-", 
                  FractionBox[
                   RowBox[{"3", " ", 
                    SuperscriptBox["t", "2"]}], "2"]}], "+", 
                 FractionBox[
                  RowBox[{"4", "-", 
                   RowBox[{"4", " ", 
                    RowBox[{"Cos", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}]}], 
                  SuperscriptBox["h", "2"]]}], ",", "0", ",", 
                FractionBox[
                 RowBox[{"3", " ", 
                  SuperscriptBox[
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"-", "h"}], " ", "t"}], "+", 
                    RowBox[{"Sin", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], ")"}], "2"]}], 
                 SuperscriptBox["h", "3"]], ",", 
                FractionBox[
                 RowBox[{
                  RowBox[{"14", " ", "h", " ", "t"}], "+", 
                  RowBox[{"3", " ", 
                   SuperscriptBox["h", "3"], " ", 
                   SuperscriptBox["t", "3"]}], "+", 
                  RowBox[{"24", " ", "h", " ", "t", " ", 
                   RowBox[{"Cos", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], "-", 
                  RowBox[{"32", " ", 
                   RowBox[{"Sin", "[", 
                    RowBox[{"h", " ", "t"}], "]"}]}], "-", 
                  RowBox[{"3", " ", 
                   RowBox[{"Sin", "[", 
                    RowBox[{"2", " ", "h", " ", "t"}], "]"}]}]}], 
                 SuperscriptBox["h", "3"]], ",", "0"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"0", ",", "0", ",", 
                FractionBox[
                 SuperscriptBox[
                  RowBox[{"Sin", "[", 
                   RowBox[{"h", " ", "t"}], "]"}], "2"], 
                 RowBox[{"2", " ", 
                  SuperscriptBox["h", "2"]}]], ",", "0", ",", "0", ",", 
                RowBox[{"-", 
                 FractionBox[
                  RowBox[{
                   RowBox[{
                    RowBox[{"-", "2"}], " ", "h", " ", "t"}], "+", 
                   RowBox[{"Sin", "[", 
                    RowBox[{"2", " ", "h", " ", "t"}], "]"}]}], 
                  RowBox[{"4", " ", 
                   SuperscriptBox["h", "3"]}]]}]}], "}"}]}], "}"}]}], ";", 
           RowBox[{"(*", "Good", "*)"}], "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"r", "[", 
             RowBox[{"t", ",", " ", "n"}], "]"}], " ", "//", " ", 
            "MatrixForm"}], ";", "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"sig", "[", 
             RowBox[{"t_", ",", " ", "n_"}], "]"}], " ", "=", " ", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"4", "-", 
                 RowBox[{"3", " ", 
                  RowBox[{"Cos", "[", 
                   RowBox[{"n", " ", "t"}], "]"}]}]}], ",", "0", ",", "0", 
                ",", 
                FractionBox[
                 RowBox[{"Sin", "[", 
                  RowBox[{"n", " ", "t"}], "]"}], "n"], ",", 
                FractionBox[
                 RowBox[{"2", "-", 
                  RowBox[{"2", " ", 
                   RowBox[{"Cos", "[", 
                    RowBox[{"n", " ", "t"}], "]"}]}]}], "n"], ",", "0"}], 
               "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"6", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"-", "n"}], " ", "t"}], "+", 
                   RowBox[{"Sin", "[", 
                    RowBox[{"n", " ", "t"}], "]"}]}], ")"}]}], ",", "1", ",", 
                "0", ",", 
                FractionBox[
                 RowBox[{"2", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"-", "1"}], "+", 
                    RowBox[{"Cos", "[", 
                    RowBox[{"n", " ", "t"}], "]"}]}], ")"}]}], "n"], ",", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "3"}], " ", "t"}], "+", 
                 FractionBox[
                  RowBox[{"4", " ", 
                   RowBox[{"Sin", "[", 
                    RowBox[{"n", " ", "t"}], "]"}]}], "n"]}], ",", "0"}], 
               "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"0", ",", "0", ",", 
                RowBox[{"Cos", "[", 
                 RowBox[{"n", " ", "t"}], "]"}], ",", "0", ",", "0", ",", 
                FractionBox[
                 RowBox[{"Sin", "[", 
                  RowBox[{"n", " ", "t"}], "]"}], "n"]}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"3", " ", "n", " ", 
                 RowBox[{"Sin", "[", 
                  RowBox[{"n", " ", "t"}], "]"}]}], ",", "0", ",", "0", ",", 
                RowBox[{"Cos", "[", 
                 RowBox[{"n", " ", "t"}], "]"}], ",", 
                RowBox[{"2", " ", 
                 RowBox[{"Sin", "[", 
                  RowBox[{"n", " ", "t"}], "]"}]}], ",", "0"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"6", " ", "n", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"-", "1"}], "+", 
                   RowBox[{"Cos", "[", 
                    RowBox[{"n", " ", "t"}], "]"}]}], ")"}]}], ",", "0", ",", 
                "0", ",", 
                RowBox[{
                 RowBox[{"-", "2"}], " ", 
                 RowBox[{"Sin", "[", 
                  RowBox[{"n", " ", "t"}], "]"}]}], ",", 
                RowBox[{
                 RowBox[{"-", "3"}], "+", 
                 RowBox[{"4", " ", 
                  RowBox[{"Cos", "[", 
                   RowBox[{"n", " ", "t"}], "]"}]}]}], ",", "0"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"0", ",", "0", ",", 
                RowBox[{
                 RowBox[{"-", "n"}], " ", 
                 RowBox[{"Sin", "[", 
                  RowBox[{"n", " ", "t"}], "]"}]}], ",", "0", ",", "0", ",", 
                RowBox[{"Cos", "[", 
                 RowBox[{"n", " ", "t"}], "]"}]}], "}"}]}], "}"}]}], ";", " ", 
           RowBox[{"(*", " ", "Good", " ", "*)"}], "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"gamma", "[", 
             RowBox[{"t_", ",", " ", "t0_", ",", " ", "n_"}], "]"}], " ", "=",
             " ", 
            RowBox[{
             RowBox[{"sig", "[", 
              RowBox[{"t", ",", " ", "n"}], "]"}], ".", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"sig", "[", 
                RowBox[{"t0", ",", " ", "n"}], "]"}], " ", "//", " ", 
               "Inverse"}], ")"}]}]}], ";", "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"nmat", "[", 
             RowBox[{"t0_", ",", " ", "tf_", ",", " ", "n_"}], "]"}], " ", 
            "=", " ", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"Join", "[", 
                RowBox[{
                 RowBox[{
                  RowBox[{"sig", "[", 
                   RowBox[{"t0", ",", " ", "n"}], "]"}], " ", "//", " ", 
                  "Inverse"}], ",", " ", 
                 RowBox[{
                  RowBox[{"-", " ", 
                   RowBox[{"sig", "[", 
                    RowBox[{"tf", ",", " ", "n"}], "]"}]}], " ", "//", " ", 
                  "Inverse"}], ",", " ", "2"}], "]"}], " ", "//", " ", 
               "Transpose"}], ")"}], ".", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"c", "[", "n", "]"}], " ", "//", " ", "Transpose"}], 
              ")"}], ".", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"r", "[", 
                RowBox[{"tf", ",", " ", "n"}], "]"}], " ", "//", " ", 
               "Inverse"}], ")"}], ".", 
             RowBox[{"c", "[", "n", "]"}], ".", 
             RowBox[{"Join", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"sig", "[", 
                 RowBox[{"t0", ",", " ", "n"}], "]"}], " ", "//", " ", 
                "Inverse"}], ",", " ", 
               RowBox[{
                RowBox[{"-", " ", 
                 RowBox[{"sig", "[", 
                  RowBox[{"tf", ",", " ", "n"}], "]"}]}], " ", "//", " ", 
                "Inverse"}], ",", " ", "2"}], "]"}]}]}], ";", 
           "\[IndentingNewLine]", 
           RowBox[{"nmat", "[", 
            RowBox[{"0", ",", " ", "tff", ",", " ", "na"}], "]"}]}]}], 
         "\[IndentingNewLine]", "]"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"getAxes", "[", 
         RowBox[{
         "terminaltime_", ",", " ", "quadform_", ",", " ", "thrustmax_"}], 
         "]"}], " ", ":=", " ", 
        RowBox[{"Module", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"testbounds", ",", " ", "teststats"}], "}"}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{
            RowBox[{"energysetaxes", "[", "quad_", "]"}], " ", ":=", 
            "\[IndentingNewLine]", 
            RowBox[{"Eigenvectors", "[", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"Join", "[", 
                RowBox[{
                 RowBox[{"Join", "[", 
                  RowBox[{
                   RowBox[{"IdentityMatrix", "[", "2", "]"}], ",", " ", 
                   RowBox[{"{", 
                    RowBox[{
                    RowBox[{"{", 
                    RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", 
                    RowBox[{"{", 
                    RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}]}], 
                    "}"}], ",", " ", "2"}], "]"}], ",", " ", 
                 RowBox[{"{", 
                  RowBox[{
                   RowBox[{"{", 
                    RowBox[{
                    "0", ",", "0", ",", "0", ",", " ", "0", ",", "0", ",", 
                    "0"}], "}"}], ",", 
                   RowBox[{"{", 
                    RowBox[{
                    "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], 
                    "}"}], ",", 
                   RowBox[{"{", 
                    RowBox[{
                    "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], 
                    "}"}], ",", 
                   RowBox[{"{", 
                    RowBox[{
                    "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], 
                    "}"}]}], "}"}]}], "]"}], ",", " ", 
               RowBox[{"quad", "[", 
                RowBox[{"[", 
                 RowBox[{
                  RowBox[{"7", ";;", " ", "12"}], ",", " ", 
                  RowBox[{"7", " ", ";;", " ", "12"}]}], "]"}], "]"}]}], 
              "}"}], "]"}]}], ";", "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"energyseteigenvalues", "[", "quad_", "]"}], " ", ":=", 
            "  ", 
            RowBox[{"Eigenvalues", "[", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"Join", "[", 
                RowBox[{
                 RowBox[{"Join", "[", 
                  RowBox[{
                   RowBox[{"IdentityMatrix", "[", "2", "]"}], ",", " ", 
                   RowBox[{"{", 
                    RowBox[{
                    RowBox[{"{", 
                    RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", 
                    RowBox[{"{", 
                    RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}]}], 
                    "}"}], ",", " ", "2"}], "]"}], ",", " ", 
                 RowBox[{"{", 
                  RowBox[{
                   RowBox[{"{", 
                    RowBox[{
                    "0", ",", "0", ",", "0", ",", " ", "0", ",", "0", ",", 
                    "0"}], "}"}], ",", 
                   RowBox[{"{", 
                    RowBox[{
                    "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], 
                    "}"}], ",", 
                   RowBox[{"{", 
                    RowBox[{
                    "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], 
                    "}"}], ",", 
                   RowBox[{"{", 
                    RowBox[{
                    "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], 
                    "}"}]}], "}"}]}], "]"}], ",", " ", 
               RowBox[{"quad", "[", 
                RowBox[{"[", 
                 RowBox[{
                  RowBox[{"7", ";;", " ", "12"}], ",", " ", 
                  RowBox[{"7", " ", ";;", " ", "12"}]}], "]"}], "]"}]}], 
              "}"}], "]"}]}], ";", "\[IndentingNewLine]", 
           RowBox[{"vecs", " ", "=", " ", 
            RowBox[{"energysetaxes", "[", "quadform", "]"}]}], ";", 
           "\[IndentingNewLine]", 
           RowBox[{"vals", " ", "=", " ", 
            RowBox[{"energyseteigenvalues", "[", "quadform", "]"}]}], ";", 
           "\[IndentingNewLine]", 
           RowBox[{"ax1", " ", "=", " ", 
            RowBox[{
             RowBox[{
              RowBox[{"vecs", "[", 
               RowBox[{"[", "1", "]"}], "]"}], "[", 
              RowBox[{"[", 
               RowBox[{"1", ";;", "2"}], "]"}], "]"}], " ", "*", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"terminaltime", " ", 
                 RowBox[{"thrustmax", "^", "2"}], "*", " ", 
                 RowBox[{"vals", "[", 
                  RowBox[{"[", "1", "]"}], "]"}]}], ")"}], " ", "//", " ", 
               "Sqrt"}], ")"}]}]}], ";", "\[IndentingNewLine]", 
           RowBox[{"ax2", " ", "=", " ", 
            RowBox[{
             RowBox[{
              RowBox[{"vecs", "[", 
               RowBox[{"[", "2", "]"}], "]"}], "[", 
              RowBox[{"[", 
               RowBox[{"1", ";;", "2"}], "]"}], "]"}], " ", "*", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"terminaltime", " ", 
                 RowBox[{"thrustmax", "^", "2"}], "*", " ", 
                 RowBox[{"vals", "[", 
                  RowBox[{"[", "2", "]"}], "]"}]}], ")"}], " ", "//", " ", 
               "Sqrt"}], ")"}]}]}], ";", "\[IndentingNewLine]", 
           RowBox[{"{", 
            RowBox[{"ax1", ",", "ax2"}], "}"}]}]}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"axes", " ", "=", " ", 
        RowBox[{"getAxes", "[", 
         RowBox[{"tt", ",", " ", 
          RowBox[{"energyform", "[", 
           RowBox[{"tt", ",", " ", "nt"}], "]"}], ",", " ", "tmax"}], "]"}]}],
        ";", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"ellipseVecFromUnitVec", "[", 
         RowBox[{"uv_", ",", "ax1_", ",", " ", "ax2_"}], "]"}], ":=", 
        RowBox[{
         RowBox[{"1", "/", 
          RowBox[{"Sqrt", "[", 
           RowBox[{
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"Normalize", "[", "ax1", "]"}], ".", "uv"}], ")"}], 
              "^", "2"}], "/", 
             RowBox[{
              RowBox[{"Norm", "[", "ax1", "]"}], "^", "2"}]}], "+", 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"Normalize", "[", "ax2", "]"}], ".", "uv"}], ")"}], 
              "^", "2"}], "/", 
             RowBox[{
              RowBox[{"Norm", "[", "ax2", "]"}], "^", "2"}]}]}], "]"}]}], "*",
          "uv"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"pos", " ", "=", " ", 
        RowBox[{"ellipseVecFromUnitVec", "[", 
         RowBox[{"delta", ",", " ", 
          RowBox[{"axes", "[", 
           RowBox[{"[", "1", "]"}], "]"}], ",", " ", 
          RowBox[{"axes", "[", 
           RowBox[{"[", "2", "]"}], "]"}]}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"pos6", "=", " ", 
        RowBox[{"Maximize", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Join", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}],
               "}"}], ",", " ", 
             RowBox[{"{", 
              RowBox[{
              "x1", ",", " ", "x2", ",", " ", "x3", ",", "x4", ",", "x5", ",",
                "x6"}], "}"}]}], "]"}], ".", 
           RowBox[{"energyform", "[", 
            RowBox[{"tt", ",", " ", "nt"}], "]"}], ".", 
           RowBox[{"Join", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}],
               "}"}], ",", 
             RowBox[{"{", 
              RowBox[{
              "x1", ",", " ", "x2", ",", " ", "x3", ",", "x4", ",", "x5", ",",
                "x6"}], "}"}]}], "]"}]}], ",", " ", 
          RowBox[{
           RowBox[{"x1", " ", "==", 
            RowBox[{"pos", "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], " ", "&&", " ", 
           RowBox[{"x2", "==", 
            RowBox[{"pos", "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], " ", "&&", " ", 
           RowBox[{
            RowBox[{
             RowBox[{"Join", "[", 
              RowBox[{
               RowBox[{"{", 
                RowBox[{
                "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], 
                "}"}], ",", " ", 
               RowBox[{"{", 
                RowBox[{
                "x1", ",", " ", "x2", ",", " ", "x3", ",", "x4", ",", "x5", 
                 ",", "x6"}], "}"}]}], "]"}], ".", 
             RowBox[{"energyform", "[", 
              RowBox[{"tt", ",", " ", "nt"}], "]"}], ".", 
             RowBox[{"Join", "[", 
              RowBox[{
               RowBox[{"{", 
                RowBox[{
                "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], 
                "}"}], ",", 
               RowBox[{"{", 
                RowBox[{
                "x1", ",", " ", "x2", ",", " ", "x3", ",", "x4", ",", "x5", 
                 ",", "x6"}], "}"}]}], "]"}]}], " ", "==", " ", 
            RowBox[{"tt", " ", 
             RowBox[{"tmax", "^", "2"}]}]}]}], ",", " ", 
          RowBox[{"{", 
           RowBox[{
           "x1", ",", "x2", ",", "x3", ",", "x4", ",", "x5", ",", "x6"}], 
           "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", "pos6"}]}], 
     "\[IndentingNewLine]", "]"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"propStateCostatev2", "[", 
    RowBox[{"ics_", ",", " ", "isp_", ",", " ", "tf_"}], "]"}], " ", ":=", 
   " ", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"vars", ",", " ", "dyn"}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{
       RowBox[{"mag", "[", 
        RowBox[{"lxv_", ",", " ", "lyv_", ",", " ", "lzv_"}], "]"}], " ", ":=", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"lxv", " ", "^", " ", "2"}], "+", " ", 
          RowBox[{"lyv", " ", "^", "2"}], " ", "+", " ", 
          RowBox[{"lzv", " ", "^", "2"}]}], ")"}], " ", "//", " ", "Sqrt"}]}],
       ";", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"switchingfunc", " ", "[", 
        RowBox[{
        "lxv_", ",", " ", "lyv_", ",", " ", "lzv_", ",", " ", "lmag_", ",", 
         " ", "mas_", ",", " ", "c_"}], "]"}], " ", ":=", " ", 
       RowBox[{
        RowBox[{"mag", "[", 
         RowBox[{"lxv", ",", " ", "lyv", ",", " ", "lzv"}], "]"}], " ", "-", 
        " ", 
        RowBox[{"lmag", " ", 
         RowBox[{"mas", " ", "/", "c"}]}]}]}], ";", " ", 
      "\[IndentingNewLine]", 
      RowBox[{"vars", "=", 
       RowBox[{"{", 
        RowBox[{
        "x", ",", "y", ",", "z", ",", "vx", ",", "vy", ",", "vz", ",", "m", 
         ",", "lxr", ",", "lyr", ",", "lzr", ",", "lxv", ",", "lyv", ",", 
         "lzv", ",", "lm"}], "}"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"eqnsOpt", " ", "=", " ", 
       RowBox[{"Join", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", " ", "A"}], " ", "//", " ", "Transpose"}], ")"}], 
          ".", 
          RowBox[{"{", 
           RowBox[{
           "lxr", ",", " ", "lyr", ",", " ", "lzr", ",", " ", "lxv", ",", " ",
             "lyv", ",", " ", "lzv"}], "}"}]}], ",", " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"(", " ", 
            RowBox[{
             RowBox[{"1", "/", "2"}], " ", "+", 
             RowBox[{
              RowBox[{"Sign", "[", 
               RowBox[{"switchingfunc", "[", 
                RowBox[{
                "lxv", ",", " ", "lyv", ",", " ", "lzv", ",", " ", "lm", ",", 
                 " ", "m", ",", " ", "isp"}], "]"}], "]"}], " ", "/", " ", 
              "2"}]}], ")"}], "tmax", " ", 
           RowBox[{
            RowBox[{"mag", "[", 
             RowBox[{"lxv", ",", " ", "lyv", ",", " ", "lzv"}], "]"}], "  ", 
            "/", " ", 
            RowBox[{"(", 
             RowBox[{"m", " ", "^", "2"}], ")"}]}]}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"eqnsStateOpt", "=", 
       RowBox[{
        RowBox[{"{", "  ", 
         RowBox[{
          RowBox[{
           RowBox[{"2", " ", "vy", " ", "n"}], " ", "+", " ", 
           RowBox[{"3", 
            RowBox[{"n", "^", "2"}], " ", "x"}]}], ",", "    ", 
          RowBox[{
           RowBox[{"-", "2"}], " ", "n", " ", "vx"}], ",", " ", 
          RowBox[{
           RowBox[{"-", " ", 
            RowBox[{"n", "^", "2"}]}], " ", "z"}]}], "  ", "}"}], " ", "+", 
        " ", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"1", "/", "m"}], ")"}], 
         RowBox[{"(", " ", 
          RowBox[{
           RowBox[{"1", "/", "2"}], " ", "+", 
           RowBox[{
            RowBox[{"Sign", "[", 
             RowBox[{"switchingfunc", "[", 
              RowBox[{
              "lxv", ",", " ", "lyv", ",", " ", "lzv", ",", " ", "lm", ",", 
               " ", "m", ",", " ", "isp"}], "]"}], "]"}], " ", "/", " ", 
            "2"}]}], ")"}], " ", "tmax", " ", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"lxv", ",", " ", "lyv", ",", " ", "lzv"}], "}"}], " ", "/",
           " ", 
          RowBox[{"mag", "[", 
           RowBox[{"lxv", ",", " ", "lyv", ",", " ", "lzv"}], "]"}]}]}]}]}], 
      ";", "\[IndentingNewLine]", 
      RowBox[{"(*", 
       RowBox[{
        RowBox[{
         RowBox[{"*", "eqnsStateOpt"}], " ", "=", " ", 
         RowBox[{"{", "  ", 
          RowBox[{
           RowBox[{
            RowBox[{"2", " ", "vy", " ", "n"}], " ", "+", " ", 
            RowBox[{"3", 
             RowBox[{"n", "^", "2"}], " ", "x"}], " ", "-", "lxv"}], ",", 
           "    ", 
           RowBox[{
            RowBox[{
             RowBox[{"-", "2"}], " ", "n", " ", "vx"}], " ", "-", "lyv"}], 
           ",", " ", 
           RowBox[{
            RowBox[{
             RowBox[{"-", " ", 
              RowBox[{"n", "^", "2"}]}], " ", "z"}], " ", "-", " ", "lzv"}]}],
           " ", "}"}]}], ";"}], " ", "**)"}], "\[IndentingNewLine]", 
      RowBox[{"dyn", "=", 
       RowBox[{"Join", "[", 
        RowBox[{
         RowBox[{"Join", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"vx", ",", "vy", ",", "vz"}], "}"}], ",", 
           RowBox[{"Join", "[", 
            RowBox[{"eqnsStateOpt", ",", " ", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"-", " ", 
                RowBox[{"(", " ", 
                 RowBox[{
                  RowBox[{"1", "/", "2"}], " ", "+", 
                  RowBox[{
                   RowBox[{"Sign", "[", 
                    RowBox[{"switchingfunc", "[", 
                    RowBox[{
                    "lxv", ",", " ", "lyv", ",", " ", "lzv", " ", ",", " ", 
                    "lm", ",", " ", "m", ",", " ", "isp"}], "]"}], "]"}], " ",
                    "/", " ", "2"}]}], ")"}]}], 
               RowBox[{"tmax", "/", "isp"}]}], "}"}]}], "]"}]}], "]"}], ",", 
         "eqnsOpt"}], "]"}]}], ";", "\[IndentingNewLine]", 
      "\[IndentingNewLine]", 
      RowBox[{"soln", " ", "=", "  ", 
       RowBox[{"NDSolveValue", "[", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"Join", "[", 
          RowBox[{
           RowBox[{"Thread", "[", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{
                RowBox[{
                 RowBox[{"#", "'"}], "[", "t", "]"}], "&"}], "/@", 
               RowBox[{"vars", "[", 
                RowBox[{"[", 
                 RowBox[{"1", ";;", "14"}], "]"}], "]"}]}], ")"}], "==", 
             RowBox[{"(", 
              RowBox[{"dyn", "/.", 
               RowBox[{"Thread", "[", 
                RowBox[{"vars", "->", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"#", "[", "t", "]"}], "&"}], "/@", "vars"}], 
                  ")"}]}], "]"}]}], ")"}]}], "]"}], ",", 
           "\[IndentingNewLine]", 
           RowBox[{"Thread", "[", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{
                RowBox[{"#", "[", "0", "]"}], "&"}], "/@", "vars"}], ")"}], 
             "\[Equal]", "ics"}], "]"}]}], "]"}], ",", "\[IndentingNewLine]", 
         "vars", ",", "  ", 
         RowBox[{"{", 
          RowBox[{"t", ",", " ", "tf", ",", " ", "tf"}], "}"}]}], 
        "\[IndentingNewLine]", "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{
        RowBox[{"#", "[", "tf", "]"}], "&"}], "/@", "soln"}]}]}], 
    "\[IndentingNewLine]", "]"}]}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"findInitialCostates", "[", 
    RowBox[{
    "icRel_", ",", " ", "fcRel_", ",", " ", "n_", ",", " ", "A_", ",", " ", 
     "td_", ",", " ", "mmax_", ",", " ", "m0_", ",", " ", "tmax_"}], "]"}], 
   " ", ":=", " ", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"vars", ",", " ", "dyn"}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"eqnsStateOpt", "=", 
       RowBox[{
        RowBox[{"{", "  ", 
         RowBox[{
          RowBox[{
           RowBox[{"2", " ", "vy", " ", "n"}], " ", "+", " ", 
           RowBox[{"3", 
            RowBox[{"n", "^", "2"}], " ", "x"}]}], ",", "    ", 
          RowBox[{
           RowBox[{"-", "2"}], " ", "n", " ", "vx"}], ",", " ", 
          RowBox[{
           RowBox[{"-", " ", 
            RowBox[{"n", "^", "2"}]}], " ", "z"}]}], "  ", "}"}], " ", "-", 
        " ", 
        RowBox[{"{", 
         RowBox[{"lxv", ",", " ", "lyv", ",", " ", "lzv"}], "}"}]}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"eqnsOpt", " ", "=", " ", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", " ", "A"}], " ", "//", " ", "Transpose"}], ")"}], ".", 
        
        RowBox[{"{", 
         RowBox[{
         "lxr", ",", " ", "lyr", ",", " ", "lzr", ",", " ", "lxv", ",", " ", 
          "lyv", ",", " ", "lzv"}], "}"}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"getInitCostates", "[", 
        RowBox[{"varMat_", ",", "icsRel_", ",", "fcsRel_"}], "]"}], ":=", 
       RowBox[{"LinearSolve", "[", 
        RowBox[{
         RowBox[{"varMat", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"1", ";;", "6"}], ",", 
            RowBox[{"7", ";;", "12"}]}], "]"}], "]"}], ",", 
         RowBox[{"fcsRel", "-", 
          RowBox[{
           RowBox[{"varMat", "[", 
            RowBox[{"[", 
             RowBox[{
              RowBox[{"1", ";;", "6"}], ",", 
              RowBox[{"1", ";;", "6"}]}], "]"}], "]"}], ".", "icsRel"}]}]}], 
        "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"topLeftMatrix", "=", "A"}], ";", "\[IndentingNewLine]", 
      RowBox[{"topRightMatrix", "=", 
       RowBox[{"ConstantArray", "[", 
        RowBox[{"0", ",", 
         RowBox[{"{", 
          RowBox[{"3", ",", "6"}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"middleLeftMatrix", "=", 
       RowBox[{"ConstantArray", "[", 
        RowBox[{"0", ",", 
         RowBox[{"{", 
          RowBox[{"3", ",", "3"}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"middleRightMatrix", "=", 
       RowBox[{"-", 
        RowBox[{"IdentityMatrix", "[", "3", "]"}]}]}], ";", "\n", "     ", 
      RowBox[{"bottomRightMatrix", "=", 
       RowBox[{
        RowBox[{"-", "A"}], " ", "//", " ", "Transpose"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"bottomLeftMatrix", "=", 
       RowBox[{"ConstantArray", "[", 
        RowBox[{"0", ",", 
         RowBox[{"{", 
          RowBox[{"6", ",", "6"}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"fullmatrix", " ", "=", " ", 
       RowBox[{"Join", "[", 
        RowBox[{
         RowBox[{"Join", "[", 
          RowBox[{"A", ",", " ", 
           RowBox[{"Join", "[", " ", 
            RowBox[{"topRightMatrix", ",", " ", 
             RowBox[{"Join", "[", 
              RowBox[{
              "middleLeftMatrix", ",", " ", "middleRightMatrix", ",", " ", 
               "2"}], "]"}]}], "]"}], ",", " ", "2"}], "]"}], ",", " ", 
         RowBox[{"Join", "[", 
          RowBox[{
          "bottomLeftMatrix", ",", " ", "bottomRightMatrix", ",", " ", "2"}], 
          " ", "]"}]}], " ", "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"stm", "[", "tf_", "]"}], " ", ":=", " ", 
       RowBox[{"MatrixExp", "[", 
        RowBox[{"fullmatrix", " ", "tf"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"initCo", " ", "=", " ", 
       RowBox[{"getInitCostates", "[", 
        RowBox[{
         RowBox[{"stm", "[", "td", "]"}], ",", " ", "icRel", ",", " ", 
         "fcRel"}], "]"}]}], ";"}]}], "\[IndentingNewLine]", 
    "\[IndentingNewLine]", "]"}]}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"numerPropAndSTM2", "[", 
    RowBox[{
    "vars_", ",", "ics_", ",", "dyn_", ",", "jac_", ",", "jac2_", ",", 
     "time_"}], "]"}], ":=", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"soln", ",", "zeros", ",", "tpzeros"}], "}"}], ",", 
     RowBox[{
      RowBox[{"zeros", "=", 
       RowBox[{"ConstantArray", "[", 
        RowBox[{"0", ",", 
         RowBox[{"Length", "[", "vars", "]"}]}], "]"}]}], ";", 
      RowBox[{"tpzeros", "=", 
       RowBox[{"TensorProduct", "[", 
        RowBox[{
         RowBox[{"TensorProduct", "[", 
          RowBox[{"zeros", ",", "zeros"}], "]"}], ",", "zeros"}], "]"}]}], 
      ";", 
      RowBox[{"soln", "=", 
       RowBox[{"NDSolve", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Thread", "[", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{
                RowBox[{
                 RowBox[{"#", "'"}], "[", "t", "]"}], "&"}], "/@", "vars"}], 
              ")"}], "==", 
             RowBox[{"(", 
              RowBox[{"dyn", "/.", 
               RowBox[{"Thread", "[", 
                RowBox[{"vars", "->", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"#", "[", "t", "]"}], "&"}], "/@", "vars"}], 
                  ")"}]}], "]"}]}], ")"}]}], "]"}], ",", 
           RowBox[{"Thread", "[", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{
                RowBox[{"#", "[", "0", "]"}], "&"}], "/@", "vars"}], ")"}], 
             "\[Equal]", "ics"}], "]"}], ",", 
           RowBox[{
            RowBox[{
             RowBox[{"Phi", "'"}], "[", "t", "]"}], "==", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"jac", "/.", 
               RowBox[{"Thread", "[", 
                RowBox[{"vars", "->", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"#", "[", "t", "]"}], "&"}], "/@", "vars"}], 
                  ")"}]}], "]"}]}], ")"}], ".", 
             RowBox[{"Phi", "[", "t", "]"}]}]}], ",", 
           RowBox[{
            RowBox[{
             RowBox[{"Phi2", "'"}], "[", "t", "]"}], "==", 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{"jac", "/.", 
                RowBox[{"Thread", "[", 
                 RowBox[{"vars", "->", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{
                    RowBox[{"#", "[", "t", "]"}], "&"}], "/@", "vars"}], 
                   ")"}]}], "]"}]}], ")"}], ".", 
              RowBox[{"Phi2", "[", "t", "]"}]}], "+", 
             RowBox[{
              RowBox[{"Transpose", "[", 
               RowBox[{
                RowBox[{"Transpose", "[", 
                 RowBox[{"Phi", "[", "t", "]"}], "]"}], ".", 
                RowBox[{"Transpose", "[", 
                 RowBox[{"(", 
                  RowBox[{"jac2", "/.", 
                   RowBox[{"Thread", "[", 
                    RowBox[{"vars", "->", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"#", "[", "t", "]"}], "&"}], "/@", "vars"}], 
                    ")"}]}], "]"}]}], ")"}], "]"}]}], "]"}], ".", 
              RowBox[{"Phi", "[", "t", "]"}]}]}]}], ",", 
           RowBox[{
            RowBox[{"Phi", "[", "0", "]"}], "==", 
            RowBox[{"IdentityMatrix", "[", 
             RowBox[{"Length", "[", "vars", "]"}], "]"}]}], ",", 
           RowBox[{
            RowBox[{"Phi2", "[", "0", "]"}], "==", "tpzeros"}]}], "}"}], ",", 
         
         RowBox[{"Join", "[", 
          RowBox[{"vars", ",", 
           RowBox[{"{", 
            RowBox[{"Phi", ",", "Phi2"}], "}"}]}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", "time", ",", "time"}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"#", "[", "time", "]"}], "&"}], "/@", "vars"}], ")"}], "/.", 
         RowBox[{"soln", "[", 
          RowBox[{"[", "1", "]"}], "]"}]}], ",", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Phi", "[", "time", "]"}], "/.", 
          RowBox[{"soln", "[", 
           RowBox[{"[", "1", "]"}], "]"}]}], ")"}], ",", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"Phi2", "[", "time", "]"}], "/.", 
          RowBox[{"soln", "[", 
           RowBox[{"[", "1", "]"}], "]"}]}], ")"}]}], "}"}]}]}], "]"}]}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"findInitialCostatesThrust", "[", 
    RowBox[{
    "icRel_", ",", " ", "deltas_", ",", " ", "n_", ",", " ", "A_", ",", " ", 
     "td_", ",", " ", "mmax_", ",", " ", "m0_", ",", " ", "tmax_", ",", " ", 
     "isp_"}], "]"}], " ", ":=", " ", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"vars", ",", " ", "dyn"}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"vars", "=", 
       RowBox[{"{", 
        RowBox[{
        "x", ",", "y", ",", "z", ",", "vx", ",", "vy", ",", "vz", ",", "lxr", 
         ",", "lyr", ",", "lzr", ",", "lxv", ",", "lyv", ",", "lzv"}], 
        "}"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"eqnsStateOpt", "=", 
       RowBox[{"{", "  ", 
        RowBox[{
         RowBox[{
          RowBox[{"2", " ", "vy", " ", "n"}], " ", "+", " ", 
          RowBox[{"3", 
           RowBox[{"n", "^", "2"}], " ", "x"}]}], ",", "    ", 
         RowBox[{
          RowBox[{"-", "2"}], " ", "n", " ", "vx"}], ",", " ", 
         RowBox[{
          RowBox[{"-", " ", 
           RowBox[{"n", "^", "2"}]}], " ", "z"}]}], "  ", "}"}]}], " ", ";", 
      "\[IndentingNewLine]", 
      RowBox[{"eqnsOpt", " ", "=", " ", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", " ", "A"}], " ", "//", " ", "Transpose"}], ")"}], ".", 
        
        RowBox[{"{", 
         RowBox[{
         "lxr", ",", " ", "lyr", ",", " ", "lzr", ",", " ", "lxv", ",", " ", 
          "lyv", ",", " ", "lzv"}], "}"}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"getInitCostates", "[", 
        RowBox[{"varMat_", ",", "d_"}], "]"}], ":=", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"varMat", "[", 
           RowBox[{"[", 
            RowBox[{
             RowBox[{"7", ";;", "12"}], ",", " ", 
             RowBox[{"7", ";;", "12"}]}], "]"}], "]"}], " ", "//", " ", 
          "Inverse"}], ")"}], ".", 
        RowBox[{"Join", "[", 
         RowBox[{"d", ",", " ", 
          RowBox[{"ConstantArray", "[", 
           RowBox[{"0", ",", "4"}], " ", "]"}]}], "]"}]}]}], ";", 
      "\[IndentingNewLine]", "\[IndentingNewLine]", 
      RowBox[{"dyn", "=", 
       RowBox[{"Join", "[", 
        RowBox[{
         RowBox[{"Join", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"vx", ",", "vy", ",", "vz"}], "}"}], ",", 
           "eqnsStateOpt"}], "]"}], ",", "eqnsOpt"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"jac", "=", 
       RowBox[{"D", "[", 
        RowBox[{"dyn", ",", 
         RowBox[{"{", "vars", "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"jac2", "=", 
       RowBox[{"D", "[", 
        RowBox[{"jac", ",", 
         RowBox[{"{", "vars", "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
      "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"stm", "[", "tf_", "]"}], " ", ":=", 
       RowBox[{
        RowBox[{"numerPropAndSTM2", "[", 
         RowBox[{"vars", ",", " ", 
          RowBox[{"ConstantArray", "[", 
           RowBox[{"0", ",", " ", "12"}], "]"}], ",", " ", "dyn", ",", " ", 
          "jac", ",", " ", "jac2", ",", " ", "tf"}], "]"}], "[", 
        RowBox[{"[", "2", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"initCo", " ", "=", " ", 
       RowBox[{"getInitCostates", "[", 
        RowBox[{
         RowBox[{"stm", "[", "td", "]"}], ",", " ", "deltas"}], "]"}]}], ";", 
      "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"lambdastm", "[", "ta_", "]"}], " ", "=", " ", 
       RowBox[{"MatrixExp", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "A"}], " ", "//", " ", "Transpose"}], ")"}], " ", 
         "ta"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"Print", "[", 
       RowBox[{"mmax", " ", 
        RowBox[{"isp", "/", "tmax"}]}], "]"}], ";", "\[IndentingNewLine]", 
      "\[IndentingNewLine]", 
      RowBox[{"lma", " ", "=", " ", 
       RowBox[{
        RowBox[{"isp", " ", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"lambdastm", "[", 
                RowBox[{"mmax", " ", 
                 RowBox[{"isp", " ", "/", " ", "tmax"}]}], "]"}], ".", 
               "initCo"}], ")"}], "[", 
             RowBox[{"[", 
              RowBox[{"4", ";;", "6"}], "]"}], "]"}], "   ", "//", " ", 
            "Norm"}], ")"}], "/", 
          RowBox[{"(", 
           RowBox[{"m0", " ", "-", " ", "mmax"}], ")"}]}]}], " ", "-", " ", 
        RowBox[{"NIntegrate", "[", "  ", 
         RowBox[{
          RowBox[{"tmax", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"lambdastm", "[", "t", "]"}], ".", "initCo"}], 
                 ")"}], "[", 
                RowBox[{"[", 
                 RowBox[{"4", ";;", "6"}], "]"}], "]"}], ")"}], " ", "//", 
              " ", "Norm"}], ")"}], " ", "/", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"m0", " ", "-", " ", 
                RowBox[{"tmax", " ", 
                 RowBox[{"t", " ", "/", " ", "isp"}]}]}], ")"}], " ", "^", 
              "2"}], ")"}]}]}], ",", " ", 
          RowBox[{"{", 
           RowBox[{"t", ",", " ", "0", ",", " ", 
            RowBox[{"mmax", " ", 
             RowBox[{"isp", " ", "/", " ", "tmax"}]}]}], "}"}]}], "]"}]}]}], 
      ";", "\[IndentingNewLine]", 
      RowBox[{"initCoM", " ", "=", " ", 
       RowBox[{"Join", "[", 
        RowBox[{"initCo", ",", " ", 
         RowBox[{"{", "lma", "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
      "\[IndentingNewLine]", "initCoM"}]}], "\[IndentingNewLine]", "]"}]}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"getInitMass", "[", 
   RowBox[{
   "mmax_", ",", " ", "isp_", ",", " ", "tmax_", ",", " ", "m0_", ",", " ", 
    "A_"}], "]"}], " ", ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"vars", ",", " ", "dyn"}], "}"}], ",", " ", 
    "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"lambdastm", "[", "ta_", "]"}], " ", "=", " ", 
      RowBox[{"MatrixExp", "[", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "A"}], " ", "//", " ", "Transpose"}], ")"}], " ", 
        "ta"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"lma", " ", "=", " ", 
      RowBox[{
       RowBox[{"isp", " ", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"lambdastm", "[", 
               RowBox[{"mmax", " ", 
                RowBox[{"isp", " ", "/", " ", "tmax"}]}], "]"}], ".", 
              "initCo"}], ")"}], "[", 
            RowBox[{"[", 
             RowBox[{"4", ";;", "6"}], "]"}], "]"}], "   ", "//", " ", 
           "Norm"}], ")"}], "/", 
         RowBox[{"(", 
          RowBox[{"m0", " ", "-", " ", "mmax"}], ")"}]}]}], " ", "-", " ", 
       RowBox[{"NIntegrate", "[", "  ", 
        RowBox[{
         RowBox[{"tmax", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 RowBox[{"lambdastm", "[", "t", "]"}], ".", "initCo"}], ")"}],
                "[", 
               RowBox[{"[", 
                RowBox[{"4", ";;", "6"}], "]"}], "]"}], ")"}], " ", "//", " ",
              "Norm"}], ")"}], " ", "/", " ", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"m0", " ", "-", " ", 
               RowBox[{"tmax", " ", 
                RowBox[{"t", " ", "/", " ", "isp"}]}]}], ")"}], " ", "^", 
             "2"}], ")"}]}]}], ",", " ", 
         RowBox[{"{", 
          RowBox[{"t", ",", " ", "0", ",", " ", 
           RowBox[{"mmax", " ", 
            RowBox[{"isp", " ", "/", " ", "tmax"}]}]}], "}"}]}], "]"}]}]}], 
     ";", "\[IndentingNewLine]", "lma"}]}], "\[IndentingNewLine]", 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"propInitialCostatesThrust", "[", 
    RowBox[{
    "icRel_", ",", " ", "n_", ",", " ", "A_", ",", " ", "td_", ",", " ", 
     "mmax_", ",", " ", "m0_", ",", " ", "tmax_"}], "]"}], " ", ":=", " ", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"vars", ",", " ", "dyn"}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"vars", "=", 
       RowBox[{"{", 
        RowBox[{
        "x", ",", "y", ",", "z", ",", "vx", ",", "vy", ",", "vz", ",", "lxr", 
         ",", "lyr", ",", "lzr", ",", "lxv", ",", "lyv", ",", "lzv"}], 
        "}"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"eqnsStateOpt", "=", 
       RowBox[{
        RowBox[{"{", "  ", 
         RowBox[{
          RowBox[{
           RowBox[{"2", " ", "vy", " ", "n"}], " ", "+", " ", 
           RowBox[{"3", 
            RowBox[{"n", "^", "2"}], " ", "x"}]}], ",", "    ", 
          RowBox[{
           RowBox[{"-", "2"}], " ", "n", " ", "vx"}], ",", " ", 
          RowBox[{
           RowBox[{"-", " ", 
            RowBox[{"n", "^", "2"}]}], " ", "z"}]}], "  ", "}"}], " ", "+", 
        " ", 
        RowBox[{"tmax", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"lxv", ",", " ", "lyv", ",", " ", "lzv"}], "}"}], " ", "/",
           " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"lxv", "^", "2"}], " ", "+", " ", 
              RowBox[{"lyv", "^", "2"}], " ", "+", 
              RowBox[{"lzv", "^", "2"}]}], " ", ")"}], " ", "//", " ", 
            "Sqrt"}], ")"}]}]}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"eqnsOpt", " ", "=", " ", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", " ", "A"}], " ", "//", " ", "Transpose"}], ")"}], ".", 
        
        RowBox[{"{", 
         RowBox[{
         "lxr", ",", " ", "lyr", ",", " ", "lzr", ",", " ", "lxv", ",", " ", 
          "lyv", ",", " ", "lzv"}], "}"}]}]}], ";", "\[IndentingNewLine]", 
      "\[IndentingNewLine]", 
      RowBox[{"dyn", "=", 
       RowBox[{"Join", "[", 
        RowBox[{
         RowBox[{"Join", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"vx", ",", "vy", ",", "vz"}], "}"}], ",", 
           "eqnsStateOpt"}], "]"}], ",", "eqnsOpt"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"jac", "=", 
       RowBox[{"D", "[", 
        RowBox[{"dyn", ",", 
         RowBox[{"{", "vars", "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"jac2", "=", 
       RowBox[{"D", "[", 
        RowBox[{"jac", ",", 
         RowBox[{"{", "vars", "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
      "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"numerPropAndSTM2", "[", 
        RowBox[{"vars", ",", " ", 
         RowBox[{"Join", "[", 
          RowBox[{
           RowBox[{"ConstantArray", "[", 
            RowBox[{"0", ",", "6"}], "]"}], ",", " ", "icRel"}], "]"}], ",", 
         " ", "dyn", ",", " ", "jac", ",", " ", "jac2", ",", " ", "td"}], 
        "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], "\[IndentingNewLine]", "]"}]}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"mu1", " ", "=", " ", 
   RowBox[{"3.986", " ", 
    RowBox[{"10", "^", "14"}]}]}], " ", ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"endTime", " ", "=", " ", "3600"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"alpha", " ", "=", " ", "7780000"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"n", " ", "=", " ", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"mu1", " ", "/", " ", 
      RowBox[{"alpha", "^", "3"}]}], ")"}], "^", 
    RowBox[{"(", 
     RowBox[{"1", "/", "2"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"A", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "0", ",", " ", "0", ",", " ", "0", ",", " ", "1", ",", " ", "0", ",", 
       " ", "0"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{
      "0", ",", " ", "0", ",", " ", "0", ",", " ", "0", ",", " ", "1", ",", 
       "0"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1"}], "}"}], 
     ",", " ", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"3", " ", 
        RowBox[{"n", "^", "2"}]}], ",", " ", "0", ",", " ", "0", ",", " ", 
       "0", ",", " ", 
       RowBox[{"2", "n"}], ",", " ", "0"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"0", ",", " ", "0", ",", " ", "0", ",", " ", 
       RowBox[{
        RowBox[{"-", "2"}], "n"}], ",", " ", "0", ",", "0"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"0", ",", " ", "0", ",", " ", 
       RowBox[{"-", 
        RowBox[{"n", "^", "2"}]}], ",", " ", "0", ",", " ", "0", ",", " ", 
       "0"}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"isp", " ", "=", " ", "10000"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"tmax", " ", "=", " ", "0.002"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"mmax", " ", "=", " ", 
   RowBox[{"0.002", " ", 
    RowBox[{"endTime", " ", "/", " ", "isp"}]}]}], 
  ";"}], "\[IndentingNewLine]", "mmax", "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"m0", " ", "=", " ", "1"}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", 
   RowBox[{
    RowBox[{"tmax", "^", "2"}], " ", "endTime"}], "]"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"delta", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"1", "/", " ", 
      RowBox[{"(", 
       RowBox[{"2", " ", "//", " ", "Sqrt"}], ")"}]}], ",", " ", 
     RowBox[{"1", "/", " ", 
      RowBox[{"(", 
       RowBox[{"2", " ", "//", " ", "Sqrt"}], ")"}]}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"energyPosF", " ", "=", " ", 
   RowBox[{"getFinalEnergyPos", "[", 
    RowBox[{"endTime", ",", " ", "tmax", ",", " ", "n", ",", " ", "delta"}], 
    "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"posf", " ", "=", " ", 
    RowBox[{"energyPosF", "[", 
     RowBox[{"[", 
      RowBox[{"2", ",", " ", "All", ",", " ", "2"}], "]"}], "]"}]}], " ", 
   ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"*", "initCostates"}], " ", "=", " ", 
     RowBox[{"findInitialCostates", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], 
        "}"}], ",", " ", "posf", ",", " ", "n", ",", " ", "A", ",", " ", 
       "endTime", ",", " ", "mmax", ",", " ", "m0", ",", " ", 
       RowBox[{"tmax", " ", "/", " ", "m0"}]}], " ", "]"}]}], ";"}], "**)"}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"initCostatesThrust", " ", "=", " ", 
   RowBox[{"findInitialCostatesThrust", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}], 
     ",", " ", "delta", ",", " ", "n", ",", " ", "A", ",", " ", "endTime", 
     ",", " ", "mmax", ",", " ", "m0", ",", " ", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"tmax", " ", 
         RowBox[{"endTime", " ", "/", "isp"}]}], " ", ">", " ", "mmax"}], ",",
        " ", 
       RowBox[{"mmax", " ", 
        RowBox[{"isp", "/", " ", "endTime"}]}], ",", " ", "tmax"}], " ", 
      "]"}], ",", " ", "isp"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", "\"\<initial costates thrust:\>\"", " ", 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", "initCostatesThrust", "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"thrustPos", " ", "=", "  ", 
    RowBox[{"propInitialCostatesThrust", "[", 
     RowBox[{
      RowBox[{"initCostatesThrust", "[", 
       RowBox[{"[", 
        RowBox[{"1", ";;", "6"}], "]"}], "]"}], ",", " ", "n", ",", " ", "A", 
      ",", " ", "endTime", ",", " ", "mmax", ",", " ", "m0", ",", " ", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"tmax", " ", 
          RowBox[{"endTime", " ", "/", " ", "isp"}]}], " ", ">", " ", 
         "mmax"}], ",", " ", 
        RowBox[{"mmax", " ", 
         RowBox[{"isp", "/", " ", "endTime"}]}], ",", " ", "tmax"}], " ", 
       "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
  
  RowBox[{"(*", 
   RowBox[{
   "*", " ", "the", " ", "actual", " ", "terminal", " ", "direction", " ", 
    "from", " ", "the", " ", "Gong", " ", "and", " ", "Gong", " ", "problem", 
    " ", "that", " ", "we", " ", "are", " ", 
    RowBox[{"targeting", " ", "--"}], " ", "to", " ", "be", " ", "used", " ", 
    "in", " ", "our", " ", "new", " ", "terminal", " ", "position", " ", 
    "constraint"}], "**)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"sf", " ", "=", " ", 
    RowBox[{"thrustPos", "[", 
     RowBox[{"[", 
      RowBox[{"1", ";;", "2"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
   "*", " ", "just", " ", "getting", " ", "final", " ", "energy", " ", 
    "position", " ", "and", " ", "first", " ", "order", " ", "approximate", 
    " ", "solution", " ", "for", " ", "comparison"}], 
   "**)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"ics", " ", "=", " ", 
    RowBox[{"Join", "[", 
     RowBox[{"Join", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
        "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "m0"}], 
        "}"}], ",", "initCostatesThrust"}], "]"}], "]"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  "\"\<final states for energy constrained trajectory:\>\"", 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"getFinalEnergyPos", "[", 
   RowBox[{"endTime", ",", " ", "tmax", ",", " ", "n", ",", " ", 
    RowBox[{"sf", " ", "/", " ", 
     RowBox[{"Norm", "[", "sf", "]"}]}]}], "]"}], "[", 
  RowBox[{"[", 
   RowBox[{"2", ",", " ", "All", ",", " ", "2"}], "]"}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"getFinalEnergyPos", "[", 
     RowBox[{"endTime", ",", " ", "tmax", ",", " ", "n", ",", " ", 
      RowBox[{"sf", " ", "/", " ", 
       RowBox[{"Norm", "[", "sf", "]"}]}]}], "]"}], "[", 
    RowBox[{"[", 
     RowBox[{"2", ",", " ", "All", ",", " ", "2"}], "]"}], "]"}], "[", 
   RowBox[{"[", 
    RowBox[{"1", ";;", "2"}], "]"}], "]"}], "  ", "/", " ", 
  RowBox[{"Norm", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"getFinalEnergyPos", "[", 
      RowBox[{"endTime", ",", " ", "tmax", ",", " ", "n", ",", " ", 
       RowBox[{"sf", " ", "/", " ", 
        RowBox[{"Norm", "[", "sf", "]"}]}]}], "]"}], "[", 
     RowBox[{"[", 
      RowBox[{"2", ",", " ", "All", ",", " ", "2"}], "]"}], "]"}], "[", 
    RowBox[{"[", 
     RowBox[{"1", ";;", "2"}], "]"}], "]"}], " ", 
   "]"}]}], "\[IndentingNewLine]", "thrustPos", "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"thrustPos", "[", 
    RowBox[{"[", 
     RowBox[{"1", ";;", "2"}], "]"}], "]"}], " ", "/", " ", 
   RowBox[{"Norm", "[", 
    RowBox[{"thrustPos", "[", 
     RowBox[{"[", 
      RowBox[{"1", ";;", "2"}], "]"}], "]"}], "]"}]}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"propStateCostatev2", "[", 
   RowBox[{
    RowBox[{"Join", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
       "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "m0"}], 
       "}"}], ",", " ", "initCostatesThrust"}], 
     RowBox[{"(*", 
      RowBox[{"*", " ", 
       RowBox[{"Join", "[", 
        RowBox[{
         RowBox[{"initCostatesThrust", "[", 
          RowBox[{"[", 
           RowBox[{"1", ";;", "6"}], "]"}], "]"}], ",", " ", 
         RowBox[{"{", 
          RowBox[{"poss", "[", 
           RowBox[{"[", 
            RowBox[{"1", ",", "2"}], "]"}], "]"}], "}"}]}], "]"}]}], "**)"}], 
     "]"}], ",", " ", "isp", ",", " ", "endTime"}], "]"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Print", "[", "\"\<hello\>\"", "]"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"*", 
    RowBox[{"findInitialCostatesMass", "[", 
     RowBox[{
     "initCostatesThrust", ",", " ", "n", ",", " ", "A", ",", " ", "endTime", 
      ",", " ", "mmax", ",", " ", "m0", ",", " ", 
      RowBox[{
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"tmax", " ", "endTime"}], " ", ">", " ", "mmax"}], ",", " ", 
         RowBox[{"mmax", "/", " ", "endTime"}], ",", " ", 
         RowBox[{"tmax", " ", "endTime"}]}], " ", "]"}], " ", "/", " ", 
       "m0"}], ",", " ", "thrustPos"}], "]"}]}], "**)"}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
   "*", " ", "constraint", " ", "and", " ", "root", " ", "finding", " ", 
    "function"}], " ", "**)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"cons", "[", 
    RowBox[{
    "dey_", ",", " ", "lf_", ",", " ", "xf_", ",", " ", "m0_", ",", " ", 
     "mmax_", ",", " ", "isp_"}], "]"}], " ", ":=", " ", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"udd", " ", "=", " ", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"2", " ", "//", " ", "IdentityMatrix"}], ")"}], " ", "-", 
        " ", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"{", "dey", "}"}], " ", "//", " ", "Transpose"}], ")"}], 
         ".", 
         RowBox[{"{", "dey", "}"}]}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"c1", " ", "=", 
       RowBox[{
        RowBox[{"{", " ", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"udd", ".", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xf", "[", 
               RowBox[{"[", "1", "]"}], "]"}], ",", " ", 
              RowBox[{"xf", "[", 
               RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], ")"}], "[", 
          RowBox[{"[", "1", "]"}], "]"}], "}"}], " ", "/", " ", 
        RowBox[{"Norm", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xf", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", " ", 
           RowBox[{"xf", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}], "]"}]}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"(*", 
       RowBox[{
        RowBox[{"*", 
         RowBox[{"Print", "[", 
          RowBox[{
           RowBox[{"lf", "[", 
            RowBox[{"[", 
             RowBox[{"1", ";;", "2"}], "]"}], "]"}], " ", ".", " ", "dey"}], 
          "]"}]}], ";"}], "**)"}], "\[IndentingNewLine]", 
      RowBox[{"c2", " ", "=", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"lf", "[", 
            RowBox[{"[", 
             RowBox[{"1", ";;", "2"}], "]"}], "]"}], " ", ".", " ", "dey"}], 
          " ", "/", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xf", "[", 
               RowBox[{"[", "1", "]"}], "]"}], ",", " ", 
              RowBox[{"xf", "[", 
               RowBox[{"[", "2", "]"}], "]"}]}], "}"}], " ", "//", " ", 
            "Norm"}], " ", ")"}]}], " ", "-", " ", "1"}], "}"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"c3", " ", "=", " ", 
       RowBox[{"lf", "[", 
        RowBox[{"[", 
         RowBox[{"3", ";;", "6"}], "]"}], "]"}]}], " ", ";", 
      "\[IndentingNewLine]", 
      RowBox[{"c4", " ", "=", "  ", 
       RowBox[{"If", " ", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"endTime", " ", 
           RowBox[{"tmax", " ", "/", " ", "isp"}]}], "<=", "mmax"}], ",", " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"m0", " ", "-", " ", 
             RowBox[{"xf", "[", 
              RowBox[{"[", "7", "]"}], "]"}]}], ")"}], " ", "-", 
           RowBox[{"endTime", " ", 
            RowBox[{"tmax", " ", "/", " ", "isp"}]}]}], "}"}], ",", " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"m0", " ", "-", " ", 
             RowBox[{"xf", "[", 
              RowBox[{"[", "7", "]"}], "]"}]}], ")"}], " ", "-", "mmax"}], 
          "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"Join", "[", 
        RowBox[{
         RowBox[{"Join", "[", 
          RowBox[{
           RowBox[{"Join", "[", 
            RowBox[{"c1", ",", " ", "c2"}], "]"}], ",", " ", "c3"}], "]"}], 
         ",", " ", "c4"}], "]"}], " ", "//", " ", "Flatten"}]}]}], 
    "\[IndentingNewLine]", 
    RowBox[{"(*", 
     RowBox[{"*", "c4"}], "**)"}], "\[IndentingNewLine]", "]"}]}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"f", "[", "joint_", "]"}], " ", ":=", " ", 
    RowBox[{"Module", "[", 
     RowBox[{
      RowBox[{"{", "out", "}"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"out", "=", " ", 
        RowBox[{"propStateCostatev2", "[", 
         RowBox[{
          RowBox[{"Join", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
             "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
              "m0"}], "}"}], ",", 
            RowBox[{"joint", "[", 
             RowBox[{"[", 
              RowBox[{"1", ";;", "7"}], "]"}], "]"}]}], " ", 
           RowBox[{"(*", 
            RowBox[{"*", 
             RowBox[{"Join", "[", 
              RowBox[{
               RowBox[{"initCostatesThrust", "[", 
                RowBox[{"[", 
                 RowBox[{"1", ";;", "6"}], "]"}], "]"}], ",", " ", 
               RowBox[{"{", "joint", "}"}]}], "]"}]}], "**)"}], "]"}], ",", 
          " ", "isp", ",", " ", "endTime"}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"cons", "[", 
        RowBox[{
         RowBox[{"sf", " ", "/", " ", 
          RowBox[{"Norm", "[", "sf", "]"}]}], ",", " ", 
         RowBox[{"out", "[", 
          RowBox[{"[", 
           RowBox[{"8", ";;", "14"}], "]"}], "]"}], ",", " ", 
         RowBox[{"out", "[", 
          RowBox[{"[", 
           RowBox[{"1", ";;", "7"}], "]"}], "]"}], " ", ",", " ", "m0", ",", 
         " ", "mmax", ",", " ", "isp"}], "]"}]}]}], "\[IndentingNewLine]", 
     "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
  "\[IndentingNewLine]", "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{
  "f", "\"\<initial costates of thrust trajectory from root finding\>\""}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"poss", " ", "=", " ", 
  RowBox[{"FindRoot", "[", 
   RowBox[{
    RowBox[{"f", "[", "l", "]"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"{", 
      RowBox[{"l", ",", " ", 
       RowBox[{"initCostatesThrust", "[", 
        RowBox[{"[", 
         RowBox[{"1", ";;", "7"}], "]"}], "]"}]}], 
      RowBox[{"(*", 
       RowBox[{"*", 
        RowBox[{"initCostatesThrust", "[", 
         RowBox[{"[", "7", "]"}], "]"}]}], "**)"}], "}"}], "}"}], ",", " ", 
    RowBox[{"Evaluated", "\[Rule]", "False"}], ",", " ", 
    RowBox[{"PrecisionGoal", "\[Rule]", "10"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  "\"\<final states of thrust-constrained trajectory:\>\"", 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"deb", " ", "=", " ", 
  RowBox[{"propStateCostatev2", "[", 
   RowBox[{
    RowBox[{"Join", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
       "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "m0"}], 
       "}"}], ",", " ", 
      RowBox[{
       RowBox[{"poss", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "2"}], "]"}], "]"}], "[", 
       RowBox[{"[", 
        RowBox[{"1", ";;", "7"}], "]"}], "]"}]}], 
     RowBox[{"(*", 
      RowBox[{"*", " ", 
       RowBox[{"Join", "[", 
        RowBox[{
         RowBox[{"initCostatesThrust", "[", 
          RowBox[{"[", 
           RowBox[{"1", ";;", "6"}], "]"}], "]"}], ",", " ", 
         RowBox[{"{", 
          RowBox[{"poss", "[", 
           RowBox[{"[", 
            RowBox[{"1", ",", "2"}], "]"}], "]"}], "}"}]}], "]"}]}], "**)"}], 
     "]"}], ",", " ", "isp", ",", " ", "endTime"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  "\"\<final constraint values of thrust-constrained trajectory:\>\"", 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"cons", "[", 
   RowBox[{
    RowBox[{"sf", " ", "/", " ", 
     RowBox[{"Norm", "[", "sf", "]"}]}], ",", " ", 
    RowBox[{"deb", "[", 
     RowBox[{"[", 
      RowBox[{"8", ";;", "14"}], "]"}], "]"}], ",", " ", 
    RowBox[{"deb", "[", 
     RowBox[{"[", 
      RowBox[{"1", ";;", "7"}], "]"}], "]"}], ",", " ", "m0", ",", " ", 
    "mmax", ",", " ", "isp"}], "]"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"mmax", " ", "=", " ", 
    RowBox[{"0.002", " ", "0.5", " ", 
     RowBox[{"endTime", " ", "/", " ", "isp"}]}]}], ";"}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
   "*", " ", "constraint", " ", "and", " ", "root", " ", "finding", " ", 
    "function"}], " ", "**)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"cons", "[", 
    RowBox[{
    "dey_", ",", " ", "lf_", ",", " ", "xf_", ",", " ", "m0_", ",", " ", 
     "mmax_", ",", " ", "isp_"}], "]"}], " ", ":=", " ", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"udd", " ", "=", " ", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"2", " ", "//", " ", "IdentityMatrix"}], ")"}], " ", "-", 
        " ", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"{", "dey", "}"}], " ", "//", " ", "Transpose"}], ")"}], 
         ".", 
         RowBox[{"{", "dey", "}"}]}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"c1", " ", "=", 
       RowBox[{
        RowBox[{"{", " ", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"udd", ".", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xf", "[", 
               RowBox[{"[", "1", "]"}], "]"}], ",", " ", 
              RowBox[{"xf", "[", 
               RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], ")"}], "[", 
          RowBox[{"[", "1", "]"}], "]"}], "}"}], " ", "/", " ", 
        RowBox[{"Norm", "[", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xf", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", " ", 
           RowBox[{"xf", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}], "]"}]}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"c2", " ", "=", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"lf", "[", 
            RowBox[{"[", 
             RowBox[{"1", ";;", "2"}], "]"}], "]"}], " ", ".", " ", "dey"}], 
          " ", "/", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xf", "[", 
               RowBox[{"[", "1", "]"}], "]"}], ",", " ", 
              RowBox[{"xf", "[", 
               RowBox[{"[", "2", "]"}], "]"}]}], "}"}], " ", "//", " ", 
            "Norm"}], " ", ")"}]}], " ", "-", " ", "1"}], "}"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"c3", " ", "=", " ", 
       RowBox[{"lf", "[", 
        RowBox[{"[", 
         RowBox[{"3", ";;", "6"}], "]"}], "]"}]}], " ", ";", 
      "\[IndentingNewLine]", 
      RowBox[{"c4", " ", "=", "  ", 
       RowBox[{"If", " ", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"endTime", " ", 
           RowBox[{"tmax", " ", "/", " ", "isp"}]}], "<=", "mmax"}], ",", " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"m0", " ", "-", " ", 
             RowBox[{"xf", "[", 
              RowBox[{"[", "7", "]"}], "]"}]}], ")"}], " ", "-", 
           RowBox[{"endTime", " ", 
            RowBox[{"tmax", " ", "/", " ", "isp"}]}]}], "}"}], ",", " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"m0", " ", "-", " ", 
             RowBox[{"xf", "[", 
              RowBox[{"[", "7", "]"}], "]"}]}], ")"}], " ", "-", "mmax"}], 
          "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"Join", "[", 
        RowBox[{
         RowBox[{"Join", "[", 
          RowBox[{
           RowBox[{"Join", "[", 
            RowBox[{"c1", ",", " ", "c2"}], "]"}], ",", " ", "c3"}], "]"}], 
         ",", " ", "c4"}], "]"}], " ", "//", " ", "Flatten"}]}]}], 
    "\[IndentingNewLine]", 
    RowBox[{"(*", 
     RowBox[{"*", "c4"}], "**)"}], "\[IndentingNewLine]", "]"}]}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"f", "[", "joint_", "]"}], " ", ":=", " ", 
    RowBox[{"Module", "[", 
     RowBox[{
      RowBox[{"{", "out", "}"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"out", "=", " ", 
        RowBox[{"propStateCostatev2", "[", 
         RowBox[{
          RowBox[{"Join", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
             "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
              "m0"}], "}"}], ",", 
            RowBox[{"joint", "[", 
             RowBox[{"[", 
              RowBox[{"1", ";;", "7"}], "]"}], "]"}]}], "]"}], ",", " ", 
          "isp", ",", " ", "endTime"}], "]"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"Join", "[", 
        RowBox[{"cons", "[", 
         RowBox[{
          RowBox[{"sf", " ", "/", " ", 
           RowBox[{"Norm", "[", "sf", "]"}]}], ",", 
          RowBox[{"out", "[", 
           RowBox[{"[", 
            RowBox[{"8", ";;", "14"}], "]"}], "]"}], ",", " ", 
          RowBox[{"out", "[", 
           RowBox[{"[", 
            RowBox[{"1", ";;", "7"}], "]"}], "]"}], " ", ",", " ", "m0", ",", 
          " ", "mmax", ",", " ", "isp"}], "]"}], "]"}]}]}], 
     "\[IndentingNewLine]", "]"}]}], ";"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"*", "initCostatesThrust"}], " ", "=", " ", 
    RowBox[{"findInitialCostatesThrust", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}],
       ",", " ", "delta", ",", " ", "n", ",", " ", "A", ",", " ", "endTime", 
      ",", " ", "mmax", ",", " ", "m0", ",", " ", 
      RowBox[{"tmax", " ", "/", " ", "2"}], ",", " ", "isp"}], "]"}]}], 
   "**)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"lm2", " ", "=", " ", 
    RowBox[{"getInitMass", "[", 
     RowBox[{
     "mmax", ",", " ", "isp", ",", " ", "tmax", ",", " ", "m0", ",", " ", 
      "A"}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  "\"\<final states and costates of fuel-constrained trajectory:\>\"", 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"poss2", " ", "=", " ", 
   RowBox[{"FindRoot", "[", 
    RowBox[{
     RowBox[{"f", "[", "l", "]"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"{", 
       RowBox[{"l", ",", " ", 
        RowBox[{"Join", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"poss", "[", 
            RowBox[{"[", 
             RowBox[{"1", ",", "2"}], "]"}], "]"}], "[", 
           RowBox[{"[", 
            RowBox[{"1", ";;", "6"}], "]"}], "]"}], ",", " ", 
          RowBox[{"{", "lm2", "}"}]}], "]"}]}], 
       RowBox[{"(*", 
        RowBox[{"*", 
         RowBox[{
          RowBox[{"poss", "[", 
           RowBox[{"[", 
            RowBox[{"1", ",", "2"}], "]"}], "]"}], "[", 
          RowBox[{"[", 
           RowBox[{"1", ";;", "7"}], "]"}], "]"}]}], "**)"}], " ", 
       RowBox[{"(*", 
        RowBox[{"*", 
         RowBox[{"Join", "[", 
          RowBox[{
           RowBox[{"initCostatesThrust", "[", 
            RowBox[{"[", 
             RowBox[{"1", ";;", "6"}], "]"}], "]"}], ",", " ", 
           RowBox[{"{", 
            RowBox[{"poss", "[", 
             RowBox[{"[", 
              RowBox[{"1", ",", "2"}], "]"}], "]"}], "}"}]}], "]"}]}], 
        "**)"}], " ", "}"}], "}"}], ",", " ", 
     RowBox[{"Evaluated", "\[Rule]", "False"}], ",", " ", 
     RowBox[{"PrecisionGoal", "\[Rule]", "10"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"deb", " ", "=", " ", 
  RowBox[{"propStateCostatev2", "[", 
   RowBox[{
    RowBox[{"Join", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
       "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "m0"}], 
       "}"}], ",", " ", 
      RowBox[{
       RowBox[{"poss2", "[", 
        RowBox[{"[", 
         RowBox[{"1", ",", "2"}], "]"}], "]"}], "[", 
       RowBox[{"[", 
        RowBox[{"1", ";;", "7"}], "]"}], "]"}]}], "]"}], ",", " ", "isp", ",",
     " ", "endTime"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"deb", "[", 
    RowBox[{"[", 
     RowBox[{"1", ";;", "2"}], "]"}], "]"}], " ", "//", " ", "Normalize"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", "\"\<final constraint satisfaction:\>\"", 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"cons", "[", 
   RowBox[{
    RowBox[{"sf", " ", "/", " ", 
     RowBox[{"Norm", "[", "sf", "]"}]}], ",", " ", 
    RowBox[{"deb", "[", 
     RowBox[{"[", 
      RowBox[{"8", ";;", "14"}], "]"}], "]"}], ",", " ", 
    RowBox[{"deb", "[", 
     RowBox[{"[", 
      RowBox[{"1", ";;", "7"}], "]"}], "]"}], ",", " ", "m0", ",", " ", 
    "mmax", ",", " ", "isp"}], "]"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"getFinalEnergyPos", "[", 
     RowBox[{"endTime", ",", " ", "tmax", ",", " ", "n", ",", " ", 
      RowBox[{"{", 
       RowBox[{"0.2395253917501401", ",", 
        RowBox[{"-", "0.9708901002208962"}]}], "}"}]}], "]"}], "[", 
    RowBox[{"[", 
     RowBox[{"2", ",", " ", "All", ",", " ", "2"}], "]"}], "]"}], "[", 
   RowBox[{"[", 
    RowBox[{"1", ";;", "2"}], "]"}], "]"}], " ", 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.93705350835292*^9, 3.937053558858734*^9}, {
   3.9370535972118607`*^9, 3.937053657572639*^9}, {3.93706183443497*^9, 
   3.937061872834649*^9}, {3.9370619084015293`*^9, 3.937061929977111*^9}, {
   3.937061989344798*^9, 3.937062032324397*^9}, {3.937062139875431*^9, 
   3.937062235594915*^9}, {3.9370627251470423`*^9, 3.937062740073971*^9}, {
   3.937062846064793*^9, 3.9370628916681757`*^9}, {3.9370635836910973`*^9, 
   3.937063584167029*^9}, {3.937063648225327*^9, 3.9370636486559353`*^9}, {
   3.937063721800565*^9, 3.9370637280857983`*^9}, {3.937063871424686*^9, 
   3.937063880156501*^9}, {3.937063961512188*^9, 3.9370639631099033`*^9}, 
   3.937063993643808*^9, {3.93706408213673*^9, 3.937064181802022*^9}, {
   3.937064255141423*^9, 3.9370643676816797`*^9}, {3.937064398192511*^9, 
   3.937064415637086*^9}, {3.9370845723303213`*^9, 3.93708457274804*^9}, 
   3.937084925365774*^9, {3.937085007680208*^9, 3.937085035298552*^9}, {
   3.937085439687063*^9, 3.937085439879554*^9}, {3.937085585182461*^9, 
   3.937085640693947*^9}, {3.937136550130637*^9, 3.9371366017443943`*^9}, {
   3.93713665139963*^9, 3.937136730161511*^9}, {3.9371368164188013`*^9, 
   3.937136917397297*^9}, {3.9371372751678143`*^9, 3.9371372850353203`*^9}, {
   3.9371375332617817`*^9, 3.937137540531474*^9}, {3.93713762552958*^9, 
   3.937137626424892*^9}, {3.937137685934916*^9, 3.9371377623629637`*^9}, {
   3.937137801206596*^9, 3.937137812920061*^9}, {3.9371382568417463`*^9, 
   3.9371382589969673`*^9}, {3.93713833710721*^9, 3.9371383459135237`*^9}, {
   3.937138443156794*^9, 3.937138537654931*^9}, {3.9371386818414927`*^9, 
   3.9371387128263693`*^9}, {3.937138857685326*^9, 3.937138864699889*^9}, {
   3.9371391215582523`*^9, 3.937139186316373*^9}, {3.937139316492498*^9, 
   3.937139345885913*^9}, {3.937139382312868*^9, 3.937139385693008*^9}, {
   3.937139794954873*^9, 3.9371398120733013`*^9}, {3.9371400443125687`*^9, 
   3.937140050048024*^9}, {3.937140105430173*^9, 3.9371401271845016`*^9}, {
   3.937140258103634*^9, 3.937140276762204*^9}, {3.937140313411147*^9, 
   3.937140318678821*^9}, {3.9371403725078697`*^9, 3.937140425905015*^9}, {
   3.9371404611096287`*^9, 3.937140591076659*^9}, {3.937140705302157*^9, 
   3.937140712692605*^9}, {3.9371409196911497`*^9, 3.937140942880443*^9}, {
   3.937140977067135*^9, 3.937140984309245*^9}, {3.937141033654162*^9, 
   3.93714114389382*^9}, {3.9371411830685577`*^9, 3.9371411883314238`*^9}, {
   3.937141253172132*^9, 3.937141256721178*^9}, {3.9371413933670177`*^9, 
   3.937141423986206*^9}, {3.937141459407975*^9, 3.937141485765575*^9}, {
   3.9371415968591557`*^9, 3.937141668210816*^9}, 3.937141707358955*^9, {
   3.93714175832408*^9, 3.9371417757611713`*^9}, {3.937141885655761*^9, 
   3.937141969608136*^9}, {3.9371420419043407`*^9, 3.937142043011734*^9}, {
   3.937142103175931*^9, 3.937142165049892*^9}, {3.937142208830531*^9, 
   3.937142246980523*^9}, {3.937142792769424*^9, 3.937142796682479*^9}, 
   3.937142829363441*^9, {3.9371428721182623`*^9, 3.937142873652424*^9}, {
   3.9371429155006638`*^9, 3.937142977224495*^9}, {3.937143008758419*^9, 
   3.937143091175146*^9}, {3.937143139949913*^9, 3.93714316970857*^9}, {
   3.937143201078314*^9, 3.93714322095926*^9}, {3.9371432865771427`*^9, 
   3.9371432934297953`*^9}, {3.9371433242837343`*^9, 3.937143330394905*^9}, {
   3.9371433676087103`*^9, 3.937143371482441*^9}, {3.937143434100151*^9, 
   3.9371434840248957`*^9}, {3.9371435576004877`*^9, 3.937143605274767*^9}, 
   3.93714370233375*^9, {3.937144186957676*^9, 3.937144187600604*^9}, {
   3.937147801889592*^9, 3.937147823688785*^9}, {3.93714793771751*^9, 
   3.937147992278372*^9}, {3.9371480853182487`*^9, 3.937148199476568*^9}, {
   3.937148234430827*^9, 3.9371482488804626`*^9}, {3.937148280698773*^9, 
   3.937148297849609*^9}, {3.937148354132831*^9, 3.937148372952359*^9}, {
   3.937148410029887*^9, 3.937148416712673*^9}, {3.9371485498626966`*^9, 
   3.937148580842741*^9}, {3.937148611955162*^9, 3.937148612005842*^9}, {
   3.937148665518873*^9, 3.937148699379434*^9}, {3.9371488992848597`*^9, 
   3.937148902373657*^9}, {3.937148941173057*^9, 3.937148955495116*^9}, {
   3.93715583205921*^9, 3.937155835941169*^9}, {3.937155883685265*^9, 
   3.937155937478775*^9}, {3.9371569125308313`*^9, 3.937156986198902*^9}, {
   3.937157041812933*^9, 3.937157072432756*^9}, {3.937157127041087*^9, 
   3.937157152413331*^9}, {3.937157187991329*^9, 3.937157205501108*^9}, {
   3.937157407355515*^9, 3.937157407382217*^9}, {3.937157464857601*^9, 
   3.937157506169724*^9}, {3.937157844453987*^9, 3.9371578452444344`*^9}, {
   3.937158430398716*^9, 3.9371584316497927`*^9}, {3.937158483461738*^9, 
   3.9371584870088663`*^9}, {3.93715866745231*^9, 3.9371586690206127`*^9}, {
   3.937159948975848*^9, 3.937160010690188*^9}, {3.93716010148298*^9, 
   3.937160140876852*^9}, 3.9371602719520082`*^9, {3.937160404899146*^9, 
   3.937160406973774*^9}, {3.937160488586479*^9, 3.937160493882042*^9}, {
   3.937160594793235*^9, 3.937160778878606*^9}, {3.9371608187145157`*^9, 
   3.9371608230635147`*^9}, {3.9371609192146797`*^9, 
   3.9371609234178333`*^9}, {3.93716096325671*^9, 3.9371609726254597`*^9}, {
   3.937161028630447*^9, 3.9371610308550673`*^9}, 3.937161084174088*^9, {
   3.937161114781172*^9, 3.937161132568988*^9}, {3.9371750383381157`*^9, 
   3.937175041485235*^9}, {3.93722127111056*^9, 3.937221303446209*^9}, {
   3.937221343416337*^9, 3.937221406149288*^9}, {3.937221436323106*^9, 
   3.937221502484054*^9}, {3.937221541819149*^9, 3.937221546750074*^9}, {
   3.937221595601055*^9, 3.937221630099284*^9}, {3.9372220117040167`*^9, 
   3.93722201615558*^9}, {3.937222050074133*^9, 3.937222075086882*^9}, {
   3.937222106544096*^9, 3.937222108691905*^9}, {3.937222145986369*^9, 
   3.93722217376891*^9}, {3.93722234942669*^9, 3.93722242651884*^9}, {
   3.937223147062192*^9, 3.937223148609968*^9}, {3.9372242665105877`*^9, 
   3.937224273116838*^9}, {3.937226269523575*^9, 3.9372263346997232`*^9}, {
   3.937226413075775*^9, 3.937226457341488*^9}, {3.937226491028414*^9, 
   3.937226515961446*^9}, {3.93722659215732*^9, 3.9372266382022457`*^9}, {
   3.937226669276916*^9, 3.937226695706571*^9}, {3.937226777189086*^9, 
   3.9372267812901993`*^9}, 3.9372268196002693`*^9, {3.937226924639823*^9, 
   3.937226930040144*^9}, {3.937226969182375*^9, 3.937226969614772*^9}, {
   3.937227025983151*^9, 3.9372270292817793`*^9}, {3.9372270886293707`*^9, 
   3.937227156793589*^9}, {3.9372273248760023`*^9, 3.937227348137269*^9}, {
   3.937227384793066*^9, 3.937227448479545*^9}, {3.93722748667642*^9, 
   3.937227534199085*^9}, {3.9372277617083483`*^9, 3.937227762970385*^9}, {
   3.9372278264429913`*^9, 3.937227942253145*^9}, {3.9372280529270983`*^9, 
   3.9372280554723587`*^9}, {3.937228125465045*^9, 3.937228131229741*^9}, 
   3.937228190275086*^9, {3.937228231307666*^9, 3.9372282773746862`*^9}, {
   3.937228311650166*^9, 3.937228312073114*^9}, {3.9372283733156977`*^9, 
   3.937228385293358*^9}, {3.9372284238741217`*^9, 3.9372284759352407`*^9}, {
   3.93722850712556*^9, 3.937228724527253*^9}, {3.9372287786566057`*^9, 
   3.937228781024295*^9}, {3.937229001568819*^9, 3.9372290596137037`*^9}, {
   3.9372290967259817`*^9, 3.9372290970218363`*^9}, {3.937229141993535*^9, 
   3.937229242853209*^9}, {3.937229275492755*^9, 3.937229276656104*^9}, {
   3.937229332939745*^9, 3.937229363832939*^9}, {3.9372294074909286`*^9, 
   3.9372294079303913`*^9}, {3.9372294918574047`*^9, 3.93722949707275*^9}, 
   3.937229540340839*^9, {3.937230019302351*^9, 3.9372300318473797`*^9}, {
   3.937230137597727*^9, 3.937230137833198*^9}, {3.937230182640107*^9, 
   3.937230183283309*^9}, {3.9372306770302362`*^9, 3.937230677324362*^9}, {
   3.937230741461863*^9, 3.937230807349839*^9}, {3.9372308738061237`*^9, 
   3.937230887938032*^9}, {3.937231001914177*^9, 3.937231071782484*^9}, {
   3.937231470950163*^9, 3.937231471566681*^9}, {3.937231521789637*^9, 
   3.937231524563328*^9}, {3.937231583179676*^9, 3.937231730078823*^9}, {
   3.937231782855981*^9, 3.937231906452795*^9}, {3.937231978349945*^9, 
   3.937232026873981*^9}, {3.9372320606619787`*^9, 3.937232157990294*^9}, 
   3.9372330720485764`*^9, {3.937233111045189*^9, 3.937233124214448*^9}, {
   3.93723318039574*^9, 3.937233187023243*^9}, 3.937233249508154*^9, {
   3.937233489027205*^9, 3.93723350532703*^9}, {3.9372336132104187`*^9, 
   3.9372336392504597`*^9}, 3.937233776460741*^9, {3.937233869629094*^9, 
   3.9372339263370323`*^9}, {3.937233964574546*^9, 3.937234137756032*^9}, {
   3.937234273566703*^9, 3.937234277546871*^9}, {3.937234312247201*^9, 
   3.9372343295561647`*^9}, {3.937234406926111*^9, 3.937234466136793*^9}, {
   3.937234497987656*^9, 3.937234544380595*^9}, {3.937234590657741*^9, 
   3.937234758038455*^9}, {3.937234843805668*^9, 3.937234852070725*^9}, 
   3.937234954030533*^9, {3.9372349858638573`*^9, 3.937235005657782*^9}, {
   3.937235048386322*^9, 3.937235375453539*^9}, {3.937235735669389*^9, 
   3.937235844095662*^9}, {3.937235874768655*^9, 3.937235959976728*^9}, {
   3.9372360391585484`*^9, 3.937236057265679*^9}, {3.937236167079422*^9, 
   3.937236175585017*^9}, {3.937236292030177*^9, 3.937236307575183*^9}, {
   3.9372363478781767`*^9, 3.93723644508202*^9}, {3.937236479758645*^9, 
   3.937236497054875*^9}, {3.937236608936*^9, 3.937236618269876*^9}, {
   3.937236664236661*^9, 3.9372367094114113`*^9}, {3.937236780412383*^9, 
   3.937236783895102*^9}, {3.937236933502548*^9, 3.937236941315362*^9}, {
   3.937237051394486*^9, 3.9372370629376783`*^9}, {3.937237574316998*^9, 
   3.937237579468853*^9}, 3.937237609919063*^9, {3.9372376712823772`*^9, 
   3.937237799490637*^9}, {3.937237858493752*^9, 3.9372378857442007`*^9}, {
   3.937237991360865*^9, 3.937238046470437*^9}, {3.9372380868822308`*^9, 
   3.937238110659472*^9}, {3.937238417883522*^9, 3.93723848310644*^9}, {
   3.937240993563719*^9, 3.937240995607195*^9}, 3.937241419528412*^9, {
   3.937241485879045*^9, 3.9372415648105307`*^9}, {3.937242083678384*^9, 
   3.9372420839743013`*^9}, {3.937242188648161*^9, 3.9372422038108587`*^9}, {
   3.937242286849629*^9, 3.937242355919815*^9}, {3.937242600753927*^9, 
   3.9372426010537577`*^9}, {3.937242643375333*^9, 3.937242648262897*^9}, {
   3.937242793295389*^9, 3.937242927326079*^9}, {3.937243094378413*^9, 
   3.937243096956916*^9}, {3.937243377350292*^9, 3.937243386491987*^9}, {
   3.937243427268278*^9, 3.937243453789185*^9}, {3.937243484573009*^9, 
   3.937243488458584*^9}, {3.93724378395986*^9, 3.9372437850626383`*^9}, {
   3.937243838183195*^9, 3.937243875217252*^9}, 3.9372444259565763`*^9, {
   3.9372447192036963`*^9, 3.9372448338364677`*^9}, {3.937244940405778*^9, 
   3.937244960445355*^9}, {3.937245007375515*^9, 3.937245020605474*^9}, {
   3.937245053135069*^9, 3.937245060958953*^9}, {3.93724514906931*^9, 
   3.937245190087422*^9}, {3.937245226226473*^9, 3.9372452272544394`*^9}, {
   3.9372452915109463`*^9, 3.9372453151314077`*^9}, {3.937245527704549*^9, 
   3.9372455496073713`*^9}, {3.9372456903209333`*^9, 3.937245708777264*^9}, {
   3.937245759192192*^9, 3.937245828102579*^9}, {3.937245906542501*^9, 
   3.937245909580248*^9}, {3.937245977347066*^9, 3.937245999742633*^9}, {
   3.9372462898878098`*^9, 3.937246342641831*^9}, {3.9372463910544167`*^9, 
   3.937246394684618*^9}, {3.937246439957123*^9, 3.937246442937044*^9}, {
   3.937246474377736*^9, 3.937246495355443*^9}, {3.937246533059181*^9, 
   3.9372465339878407`*^9}, {3.937246564675798*^9, 3.9372465693979053`*^9}, {
   3.937246605210766*^9, 3.937246611176898*^9}, {3.937246684891202*^9, 
   3.9372466958481417`*^9}, {3.937246740280349*^9, 3.937246764859291*^9}, {
   3.9372468008955727`*^9, 3.937246909946329*^9}, {3.937246940897829*^9, 
   3.937246968832357*^9}, {3.937247008889876*^9, 3.937247011579167*^9}, {
   3.937247061243184*^9, 3.937247129883552*^9}, {3.9372471621218576`*^9, 
   3.937247168344054*^9}, {3.93724725123955*^9, 3.9372472749246407`*^9}, {
   3.937247311896473*^9, 3.937247313169403*^9}, {3.9372473627423477`*^9, 
   3.9372474277384644`*^9}, {3.937247566240055*^9, 3.937247584383545*^9}, {
   3.9372476242749577`*^9, 3.9372477572574167`*^9}, {3.9372478194656477`*^9, 
   3.93724782288511*^9}, {3.9372479197321177`*^9, 3.937247926129434*^9}, {
   3.93724796037096*^9, 3.9372480083816566`*^9}, {3.937248109853971*^9, 
   3.9372481117801237`*^9}, {3.9372482832333107`*^9, 3.937248339641569*^9}, {
   3.937248434125031*^9, 3.937248434474708*^9}, {3.937248495632597*^9, 
   3.93724856834566*^9}, {3.937248662741152*^9, 3.937248674428624*^9}, {
   3.9372487365532923`*^9, 3.9372487430207033`*^9}, {3.9372488117402687`*^9, 
   3.937248861602325*^9}, {3.937392902355701*^9, 3.937392906572844*^9}, {
   3.9373945752720947`*^9, 3.937394625628478*^9}, {3.9379286582181463`*^9, 
   3.937928697282323*^9}, {3.937928731304757*^9, 3.937928735755253*^9}, {
   3.937928783433587*^9, 3.9379288001427402`*^9}, {3.937928875740234*^9, 
   3.937928878665063*^9}, {3.937928920008093*^9, 3.937928960898484*^9}, {
   3.937928997496633*^9, 3.937929097761903*^9}, {3.9379291330417643`*^9, 
   3.93792913353296*^9}, {3.937929550374154*^9, 3.937929579381473*^9}, {
   3.937929770315551*^9, 3.9379297899002*^9}, {3.937929846201851*^9, 
   3.937929847617475*^9}, {3.9379298840548143`*^9, 3.9379299039100037`*^9}, 
   3.937929959917014*^9, {3.937929991448649*^9, 3.9379300434487867`*^9}, {
   3.937930170643105*^9, 3.937930191058443*^9}, {3.937930231448262*^9, 
   3.9379302480230017`*^9}, {3.937930289252945*^9, 3.937930311574604*^9}, {
   3.937930533942113*^9, 3.937930561569619*^9}, {3.937930658757688*^9, 
   3.937930707605255*^9}, {3.937930737656764*^9, 3.9379307956561937`*^9}, {
   3.937930949954968*^9, 3.9379310367339687`*^9}, {3.9379311221127653`*^9, 
   3.937931124971847*^9}, {3.937931164112298*^9, 3.937931253070757*^9}, 
   3.937931292211236*^9, 3.9379314527045727`*^9, {3.93793158757267*^9, 
   3.9379315889710093`*^9}, {3.937931718848564*^9, 3.937931757225237*^9}, {
   3.937931804626539*^9, 3.937931845685945*^9}, {3.9379319377364683`*^9, 
   3.937931939873837*^9}, {3.937931988233344*^9, 3.9379320677733183`*^9}, {
   3.937932119023587*^9, 3.937932135539732*^9}, {3.937932234189713*^9, 
   3.9379322545578003`*^9}, {3.937932288678462*^9, 3.937932292578457*^9}, {
   3.937932381873559*^9, 3.937932382381572*^9}, 3.937932554976591*^9, 
   3.9379326636830683`*^9, {3.9379328591693077`*^9, 3.9379330210583057`*^9}, {
   3.9379331240060577`*^9, 3.937933124901989*^9}, {3.937933353970059*^9, 
   3.937933355437975*^9}, {3.937933451307913*^9, 3.937933452392623*^9}, {
   3.937933528126402*^9, 3.937933560168275*^9}, {3.9379336208110857`*^9, 
   3.937933621795575*^9}, {3.937934815092249*^9, 3.937934819802092*^9}, {
   3.937934862473784*^9, 3.937934864152856*^9}, {3.937935475956504*^9, 
   3.937935480636264*^9}, {3.9379356398871717`*^9, 3.9379356871665173`*^9}, {
   3.937935918204529*^9, 3.937936036845602*^9}, {3.937936091725575*^9, 
   3.937936113853488*^9}, {3.937936144925297*^9, 3.9379361511461267`*^9}, {
   3.937936397367529*^9, 3.937936416457609*^9}, {3.9379364741033297`*^9, 
   3.937936504086844*^9}, {3.93793658393305*^9, 3.937936596825048*^9}, {
   3.9379366343870373`*^9, 3.937936641984309*^9}, {3.937936694593205*^9, 
   3.9379366959767733`*^9}, {3.937937226073986*^9, 3.937937227109765*^9}, {
   3.9379381711226387`*^9, 3.937938194694655*^9}, {3.937938866987706*^9, 
   3.937938870459466*^9}, {3.937939018490645*^9, 3.9379390194916773`*^9}, 
   3.937939081105111*^9, {3.937939143135133*^9, 3.937939199748076*^9}, {
   3.937939241079938*^9, 3.9379392633255444`*^9}, 3.9379393601311207`*^9, 
   3.937939904066613*^9, {3.937940179829056*^9, 3.937940227013883*^9}, 
   3.937940279791051*^9, {3.9379403534222116`*^9, 3.937940355947411*^9}, {
   3.937940494744897*^9, 3.9379404955272503`*^9}, 3.93794060408525*^9, {
   3.937940652372191*^9, 3.9379406941752996`*^9}, {3.937940727597788*^9, 
   3.9379408047241592`*^9}, {3.9379408384983263`*^9, 3.937940858593943*^9}, {
   3.9379409602168913`*^9, 3.937940960756592*^9}, {3.937941122067276*^9, 
   3.93794112254879*^9}, {3.937941170825095*^9, 3.93794120341905*^9}, {
   3.937942054977233*^9, 3.937942055668906*^9}, {3.9379427593497477`*^9, 
   3.9379427717096767`*^9}, {3.937942989159569*^9, 3.937943008068905*^9}, {
   3.9379430399776373`*^9, 3.937943041615728*^9}, {3.937943080445533*^9, 
   3.937943081024728*^9}, {3.937943176148539*^9, 3.937943176548811*^9}, {
   3.9379433092780447`*^9, 3.937943314104207*^9}, {3.937943449927746*^9, 
   3.937943452633654*^9}, {3.93794348788608*^9, 3.937943490107112*^9}, {
   3.937943781311935*^9, 3.937943789970656*^9}, {3.937944158399091*^9, 
   3.937944158528042*^9}, {3.937944611322379*^9, 3.9379446344954243`*^9}, {
   3.9379447037617693`*^9, 3.937944704465185*^9}, {3.937945065380196*^9, 
   3.937945065733341*^9}, {3.937945221913816*^9, 3.937945223043703*^9}, {
   3.937945346407764*^9, 3.937945348398875*^9}, {3.937945423577631*^9, 
   3.937945424177815*^9}, {3.9379454562085743`*^9, 3.9379454589957314`*^9}, {
   3.937945755523217*^9, 3.937945765764037*^9}, {3.937945815250793*^9, 
   3.9379458161980047`*^9}, {3.9379459314880733`*^9, 3.937945935641095*^9}, {
   3.9379464527213383`*^9, 3.93794645746675*^9}, {3.937946507397752*^9, 
   3.937946549549014*^9}, {3.937946594167945*^9, 3.937946597886277*^9}, {
   3.937946695739216*^9, 3.937946702920233*^9}, {3.937946869915038*^9, 
   3.937946879438891*^9}, {3.937947011126465*^9, 3.9379470619320374`*^9}, {
   3.937947461620597*^9, 3.937947461760099*^9}, {3.9379475399209948`*^9, 
   3.937947559811679*^9}, {3.9379476803209057`*^9, 3.9379477149555407`*^9}, {
   3.937947750507724*^9, 3.937947782900285*^9}, {3.9379478165748653`*^9, 
   3.937947851238202*^9}, {3.9379479288897333`*^9, 3.937947930604753*^9}, {
   3.9379480015585737`*^9, 3.937948006456456*^9}, {3.93794804117202*^9, 
   3.937948041518631*^9}, {3.937949077732122*^9, 3.9379490780296383`*^9}, {
   3.9379491495653048`*^9, 3.9379491512677517`*^9}, {3.9379491987212143`*^9, 
   3.937949205058626*^9}, {3.937949270582224*^9, 3.937949304314286*^9}, {
   3.9379494566809196`*^9, 3.93794947519915*^9}, 3.9379495286337633`*^9, {
   3.9379495763358393`*^9, 3.937949594301819*^9}, {3.937949646066321*^9, 
   3.937949761169353*^9}, {3.9379497922389727`*^9, 3.9379497929876966`*^9}, 
   3.9379498340766897`*^9, {3.937949944406872*^9, 3.937949945709919*^9}, {
   3.937950106651745*^9, 3.937950111473778*^9}, {3.9379501460964317`*^9, 
   3.937950148928978*^9}, 3.937950210770644*^9, {3.937950245513818*^9, 
   3.937950276730234*^9}, {3.937950327597857*^9, 3.937950331922773*^9}, {
   3.9379503763565903`*^9, 3.937950379726297*^9}, {3.937950508518075*^9, 
   3.937950510458166*^9}, {3.937950542131851*^9, 3.937950544900867*^9}, {
   3.9379505981169643`*^9, 3.9379506060110693`*^9}, {3.9379506640986156`*^9, 
   3.9379506819243917`*^9}, {3.937950721161377*^9, 3.937950757703631*^9}, {
   3.937951416203043*^9, 3.937951504745661*^9}, {3.937951549568096*^9, 
   3.937951625974906*^9}, {3.937951661996087*^9, 3.937951824050748*^9}, {
   3.937951855818864*^9, 3.9379518889903307`*^9}, {3.937951921587183*^9, 
   3.9379519296190987`*^9}, {3.937951982875657*^9, 3.937951992925568*^9}, {
   3.937952044192618*^9, 3.93795208673603*^9}, {3.937952119356409*^9, 
   3.937952130335836*^9}, {3.93795233811819*^9, 3.9379523940014887`*^9}, {
   3.937952427327332*^9, 3.9379524651907997`*^9}, {3.937952604589239*^9, 
   3.93795267166609*^9}, {3.937952709894182*^9, 3.937952733326742*^9}, {
   3.937952776203478*^9, 3.9379528980156937`*^9}, {3.937953003712031*^9, 
   3.9379530215588303`*^9}, {3.9379531029486027`*^9, 3.937953122238484*^9}, {
   3.937953337958297*^9, 3.9379533638493843`*^9}, {3.9379534195996532`*^9, 
   3.937953423998402*^9}, {3.937953455951891*^9, 3.937953496327292*^9}, {
   3.937953645461586*^9, 3.937953647169009*^9}, {3.937953680186326*^9, 
   3.93795368152885*^9}, {3.937953815895062*^9, 3.93795381602892*^9}, {
   3.937954254078802*^9, 3.9379542779831753`*^9}, {3.937954330205964*^9, 
   3.937954385035489*^9}, {3.937954726239669*^9, 3.937954770885026*^9}, {
   3.937954808336512*^9, 3.937954810668487*^9}, 3.937955002198905*^9, {
   3.937955714869948*^9, 3.937955738127133*^9}, {3.9379561936836967`*^9, 
   3.937956195682014*^9}, {3.937956283421863*^9, 3.937956287023099*^9}, {
   3.9379563310171432`*^9, 3.937956334668744*^9}, {3.937956368398637*^9, 
   3.937956373206222*^9}, {3.93795645955682*^9, 3.937956477512946*^9}, {
   3.9379565093179417`*^9, 3.937956566096622*^9}, {3.937956617079165*^9, 
   3.9379566194124537`*^9}, {3.937956659918696*^9, 3.937956742844537*^9}, {
   3.937956817582711*^9, 3.937956819624219*^9}, {3.93795691218547*^9, 
   3.9379569172190228`*^9}, {3.9379569802523203`*^9, 3.937957009141124*^9}, {
   3.937957056926897*^9, 3.937957057450499*^9}, {3.937957188198333*^9, 
   3.9379572260599403`*^9}, 3.937957325628159*^9, {3.937957483320032*^9, 
   3.93795752327733*^9}, {3.9379577347959757`*^9, 3.937957776680243*^9}, {
   3.937958678704185*^9, 3.937958685482545*^9}, {3.937958847042952*^9, 
   3.937958864865059*^9}, {3.937958925224598*^9, 3.93795905604384*^9}, {
   3.9379590953650846`*^9, 3.937959141849354*^9}, {3.937959206513473*^9, 
   3.937959294831327*^9}, {3.9379593356711493`*^9, 3.937959377942114*^9}, {
   3.9379594110991507`*^9, 3.93795942643568*^9}, {3.937959490077298*^9, 
   3.937959491735262*^9}, {3.937959523664856*^9, 3.937959590359892*^9}, {
   3.937959785246048*^9, 3.9379598077042637`*^9}, {3.937959840419393*^9, 
   3.937959861426804*^9}, {3.937959941067325*^9, 3.937959941261526*^9}, {
   3.937960141741815*^9, 3.937960155286171*^9}, {3.937960293020919*^9, 
   3.937960305977549*^9}, {3.937960383556004*^9, 3.93796038568758*^9}, {
   3.937960433773074*^9, 3.937960437479886*^9}, {3.937960468068597*^9, 
   3.937960531328085*^9}, {3.93796060662391*^9, 3.9379606070772257`*^9}, {
   3.937960706437233*^9, 3.937960715293477*^9}, {3.937960769330121*^9, 
   3.93796081004469*^9}, {3.937997501882469*^9, 3.937997523734627*^9}, {
   3.937997605715549*^9, 3.937997609211095*^9}, {3.938002997070312*^9, 
   3.9380030830611467`*^9}, {3.938003276838945*^9, 3.938003278676323*^9}, {
   3.938003441789331*^9, 3.938003483623494*^9}, {3.938003516573255*^9, 
   3.938003516766095*^9}, {3.938003624089608*^9, 3.938003642570149*^9}, {
   3.938003766991377*^9, 3.9380038029311323`*^9}, {3.938004001251601*^9, 
   3.9380040346272087`*^9}, 3.938004116974658*^9, {3.938004155442512*^9, 
   3.938004289807262*^9}, {3.938004533490635*^9, 3.938004568686741*^9}, {
   3.938005433508244*^9, 3.93800546181595*^9}, {3.938005564419766*^9, 
   3.938005567772811*^9}, {3.938005618449395*^9, 3.938005626853236*^9}, {
   3.938005727313665*^9, 3.9380057596269617`*^9}, {3.938005859896553*^9, 
   3.938005887501047*^9}, {3.938005933950965*^9, 3.938005940174878*^9}, {
   3.9380060176791277`*^9, 3.938006035424382*^9}, {3.938006066785885*^9, 
   3.938006073479762*^9}, {3.938006203525161*^9, 3.938006257654694*^9}, {
   3.9380063519066467`*^9, 3.938006354209611*^9}, {3.938006662806581*^9, 
   3.93800670927528*^9}, {3.938006788985238*^9, 3.938006907299941*^9}, {
   3.93800696776521*^9, 3.938007027161718*^9}, {3.938007111444522*^9, 
   3.938007126773397*^9}, {3.938007179551375*^9, 3.93800720149144*^9}, {
   3.938007434140874*^9, 3.938007510624452*^9}, {3.938007585399374*^9, 
   3.938007592269741*^9}, {3.938007754081888*^9, 3.938007974305126*^9}, {
   3.938008042765717*^9, 3.938008052888371*^9}, {3.938008107399741*^9, 
   3.9380081555630713`*^9}, {3.938008220306551*^9, 3.9380083174635067`*^9}, {
   3.938008363930333*^9, 3.9380084150274353`*^9}, {3.938008470495253*^9, 
   3.93800847086775*^9}, {3.938008608011409*^9, 3.9380086247172403`*^9}, {
   3.9380087042006207`*^9, 3.938008748029833*^9}, 3.9380088038050423`*^9, {
   3.9380089167665577`*^9, 3.9380089410206947`*^9}, {3.9380090131391*^9, 
   3.938009060090501*^9}, {3.938009104010318*^9, 3.93800913678295*^9}, {
   3.938009172055335*^9, 3.938009188381671*^9}, {3.938009375006949*^9, 
   3.938009405263756*^9}, {3.938009439400552*^9, 3.9380094406846724`*^9}, {
   3.938009472452663*^9, 3.938009497197611*^9}, 3.938009531051103*^9, 
   3.93800957917478*^9, {3.938009630693367*^9, 3.9380096314327173`*^9}, {
   3.9380096787791986`*^9, 3.938009679168891*^9}, {3.9380097148630466`*^9, 
   3.938009728604886*^9}, {3.9380098096846113`*^9, 3.93800986781102*^9}, {
   3.938009909197123*^9, 3.938009993077168*^9}, {3.938010053770741*^9, 
   3.9380100540253983`*^9}, {3.9380101330672617`*^9, 
   3.9380101332800617`*^9}, {3.938010289225956*^9, 3.9380104215375566`*^9}, {
   3.9380104571013803`*^9, 3.938010499733975*^9}, {3.9380105509066877`*^9, 
   3.9380105515113583`*^9}, {3.938010630903558*^9, 3.938010631648871*^9}, {
   3.938010700137034*^9, 3.938010715457485*^9}, {3.938010762462555*^9, 
   3.938010762698612*^9}, {3.938010804881495*^9, 3.938010833978918*^9}, 
   3.93801089647488*^9, 3.938010979602172*^9, {3.938011015271311*^9, 
   3.938011027163248*^9}, {3.938011072733461*^9, 3.9380110775100107`*^9}, 
   3.9380111820523167`*^9, {3.938011305515286*^9, 3.938011305791918*^9}, {
   3.938012015896798*^9, 3.938012067299283*^9}, {3.9380121861880007`*^9, 
   3.938012189220265*^9}, {3.93801227968174*^9, 3.9380123567321253`*^9}, {
   3.938012429389269*^9, 3.938012460743844*^9}, {3.938012498834777*^9, 
   3.938012602681164*^9}, {3.938023515713524*^9, 3.938023540082715*^9}, {
   3.93802374544678*^9, 3.938023797209293*^9}, {3.938023848081275*^9, 
   3.9380238491066093`*^9}, {3.938023905731139*^9, 3.938023966567552*^9}, {
   3.9380240212368402`*^9, 3.938024050906294*^9}, {3.938024097730019*^9, 
   3.93802409805333*^9}, {3.9380241448481703`*^9, 3.938024179499345*^9}, {
   3.938024247828623*^9, 3.93802435848884*^9}, {3.9380244453356047`*^9, 
   3.938024561094759*^9}, 3.938024628847625*^9, {3.938024661167695*^9, 
   3.938024664508931*^9}, {3.938024954739627*^9, 3.938024967076111*^9}, 
   3.938025014443465*^9, {3.938025112873294*^9, 3.938025117045535*^9}, {
   3.938025370162318*^9, 3.938025376780686*^9}, {3.938025431572445*^9, 
   3.938025436119864*^9}, {3.938025495967811*^9, 3.938025497449992*^9}, {
   3.938025593700362*^9, 3.938025630535735*^9}, {3.938025663236372*^9, 
   3.938025675828327*^9}, {3.938025746997197*^9, 3.938025768069715*^9}, {
   3.9380258411681013`*^9, 3.9380258785025787`*^9}, {3.938025919060666*^9, 
   3.938025926261546*^9}, {3.938026031651608*^9, 3.938026034388341*^9}, {
   3.93802607695573*^9, 3.9380260796314077`*^9}, {3.938026116229885*^9, 
   3.938026126043518*^9}, {3.938026209039234*^9, 3.938026235377016*^9}, {
   3.938026287041013*^9, 3.938026315689879*^9}, {3.938026427479209*^9, 
   3.9380264523241987`*^9}, 3.938026506766323*^9, {3.938026595172036*^9, 
   3.9380266071289663`*^9}, {3.938026663534281*^9, 3.938026686143022*^9}, {
   3.938026750504992*^9, 3.9380268037848167`*^9}, {3.938026843514221*^9, 
   3.938026845214637*^9}, {3.938026877061681*^9, 3.9380268997854357`*^9}, {
   3.938026935715006*^9, 3.938026960967476*^9}, {3.938027007767294*^9, 
   3.938027074996011*^9}, {3.938027119184142*^9, 3.938027120323926*^9}, {
   3.9380271713025703`*^9, 3.938027296039206*^9}, {3.938028213649768*^9, 
   3.938028249880805*^9}, {3.938028285945223*^9, 3.938028291905273*^9}, {
   3.938035027024746*^9, 3.938035069457697*^9}, {3.938035120955669*^9, 
   3.9380351237397423`*^9}, {3.938035173226851*^9, 3.938035199696966*^9}, {
   3.938037779709429*^9, 3.9380377871313667`*^9}, {3.938037848513259*^9, 
   3.9380378539750357`*^9}, {3.93803824654399*^9, 3.938038269448393*^9}, {
   3.9380966342682133`*^9, 3.938096719917303*^9}, {3.938096782668829*^9, 
   3.938096786914235*^9}, {3.938096860591511*^9, 3.938096874268532*^9}, {
   3.938096911648013*^9, 3.9380969399489613`*^9}, {3.938097123478121*^9, 
   3.938097142985695*^9}, {3.9380972448388*^9, 3.938097253689465*^9}, {
   3.938097324816733*^9, 3.938097354793654*^9}, {3.938097385932271*^9, 
   3.938097452019113*^9}, {3.938101404383659*^9, 3.9381014050250683`*^9}, {
   3.93810192557918*^9, 3.938101934731495*^9}, {3.938101989620722*^9, 
   3.938101993350664*^9}, {3.938102860538067*^9, 3.938102943695158*^9}, {
   3.938102996452153*^9, 3.938103010005265*^9}, {3.938103048964369*^9, 
   3.9381030572653723`*^9}, {3.9381032257752953`*^9, 3.938103230436812*^9}, {
   3.938104083130435*^9, 3.9381041206546297`*^9}, {3.9381041537252893`*^9, 
   3.938104267282956*^9}, {3.938104639336255*^9, 3.9381046439667883`*^9}, {
   3.9381049770878963`*^9, 3.938104997049534*^9}, {3.938105029736331*^9, 
   3.9381050403921824`*^9}, {3.938106772429142*^9, 3.93810677289529*^9}, {
   3.93810691329412*^9, 3.938106924044578*^9}, {3.93810699300624*^9, 
   3.938107086845159*^9}, {3.9381071210300922`*^9, 3.938107175294312*^9}, {
   3.9381072143839703`*^9, 3.938107271574939*^9}, {3.938107738614127*^9, 
   3.938107763129703*^9}, {3.938114873529941*^9, 3.93811488274017*^9}, {
   3.938115241100313*^9, 3.938115241492303*^9}, {3.93811596056842*^9, 
   3.938115963942329*^9}, {3.9381161687735558`*^9, 3.938116168823407*^9}, {
   3.938116310228982*^9, 3.9381163283945627`*^9}, {3.938116431873468*^9, 
   3.938116434546831*^9}, {3.938116789608926*^9, 3.938116790081843*^9}, {
   3.938116824312313*^9, 3.938116827448729*^9}, 3.938116870003018*^9, {
   3.9381169104312887`*^9, 3.938116910970998*^9}, {3.938117049058552*^9, 
   3.938117117202886*^9}, {3.938117183427526*^9, 3.938117211341226*^9}, 
   3.938117258083767*^9, {3.938117292161469*^9, 3.938117307088067*^9}, {
   3.938118113779326*^9, 3.938118147435178*^9}, {3.938118192590496*^9, 
   3.938118232165169*^9}, {3.938118277740942*^9, 3.938118291129589*^9}, {
   3.938118507661171*^9, 3.938118514353874*^9}, {3.938118728511302*^9, 
   3.9381187580172157`*^9}, {3.938118857667012*^9, 3.938118944515983*^9}, {
   3.938118984700889*^9, 3.9381190349213676`*^9}, {3.938119069304469*^9, 
   3.938119085878209*^9}, {3.938119126283184*^9, 3.938119129274827*^9}, {
   3.938119202591611*^9, 3.938119209761936*^9}, {3.938119269557374*^9, 
   3.938119270551217*^9}, {3.938119423492381*^9, 3.9381194378374157`*^9}, {
   3.938119485181884*^9, 3.938119550681847*^9}, {3.9381195875208197`*^9, 
   3.938119590283517*^9}, {3.938119766962553*^9, 3.93811976701239*^9}, {
   3.938119885537759*^9, 3.93811988605335*^9}, {3.9381199615814466`*^9, 
   3.938119961680441*^9}, 3.938120022985626*^9, {3.938120083702619*^9, 
   3.938120104884576*^9}, {3.938120147719503*^9, 3.938120160855044*^9}, {
   3.938120206138413*^9, 3.938120211127398*^9}, {3.938120291506155*^9, 
   3.938120302031168*^9}, {3.938120339650607*^9, 3.9381203414808493`*^9}, 
   3.938120384793126*^9, {3.938120473683936*^9, 3.938120479126684*^9}, {
   3.938120611683885*^9, 3.938120637009308*^9}, {3.938120735320735*^9, 
   3.938120767599473*^9}, {3.938120841763009*^9, 3.938120931125325*^9}, {
   3.938121003992687*^9, 3.93812105991593*^9}, {3.938121122586055*^9, 
   3.938121122926982*^9}, {3.9381211667085667`*^9, 3.9381212469032183`*^9}, {
   3.938683287724983*^9, 3.9386833038688107`*^9}, 3.938683387047976*^9, {
   3.938683433336697*^9, 3.938683440321169*^9}, {3.938683509821046*^9, 
   3.938683510161796*^9}, {3.938685296856278*^9, 3.938685308097847*^9}, {
   3.938685361264468*^9, 3.938685367937632*^9}, {3.9386854371046886`*^9, 
   3.9386854440043297`*^9}, 3.93868548315968*^9, {3.938691107929818*^9, 
   3.9386911082116117`*^9}, {3.938691196897572*^9, 3.938691197757142*^9}, {
   3.938691349767489*^9, 3.9386913543344717`*^9}, {3.938691608366197*^9, 
   3.938691612107622*^9}, {3.938691681746273*^9, 3.938691697323954*^9}, {
   3.93869173799718*^9, 3.938691819127214*^9}, {3.938691873511807*^9, 
   3.938691880536426*^9}, {3.9386919690078793`*^9, 3.938691982125866*^9}, {
   3.9386922358571367`*^9, 3.938692236231312*^9}, {3.9386923568422403`*^9, 
   3.9386923799524717`*^9}, {3.938692413894853*^9, 3.938692419253274*^9}, {
   3.93869245086935*^9, 3.938692464302335*^9}, {3.938692501735408*^9, 
   3.938692505142548*^9}, {3.9386925885653276`*^9, 3.9386925925549927`*^9}, {
   3.9386929384887533`*^9, 3.938692990354384*^9}, 3.9386930388246717`*^9, {
   3.938693133743705*^9, 3.9386931426382027`*^9}, {3.938693570433908*^9, 
   3.938693574333301*^9}, {3.938693610393961*^9, 3.938693614658547*^9}, {
   3.938693662951879*^9, 3.9386936634447193`*^9}, {3.93869388943139*^9, 
   3.938693890164185*^9}, {3.9386939960883217`*^9, 3.9386939964774*^9}, {
   3.938694172657921*^9, 3.938694195213561*^9}, {3.9386942364905767`*^9, 
   3.938694269226692*^9}, {3.9386943390097027`*^9, 3.938694380439823*^9}, {
   3.938694411806199*^9, 3.938694444249743*^9}, 3.938694520204052*^9, {
   3.938694554276915*^9, 3.938694555532546*^9}, {3.93869460670665*^9, 
   3.938694619877657*^9}, {3.938694651240756*^9, 3.938694661810508*^9}, {
   3.938694773843067*^9, 3.9386947826009083`*^9}, {3.938695560917501*^9, 
   3.9386955674346867`*^9}, {3.938695645146678*^9, 3.938695663191742*^9}, {
   3.938695702448855*^9, 3.938695717900813*^9}, 3.938695753109356*^9, 
   3.9386958405824738`*^9, {3.938695976737164*^9, 3.938696004463256*^9}, {
   3.93869606378337*^9, 3.938696066669855*^9}, {3.938696304303516*^9, 
   3.93869632097484*^9}, {3.9386964086996307`*^9, 3.938696418407651*^9}, {
   3.938698872483584*^9, 3.938698939673051*^9}, {3.938699035408197*^9, 
   3.93869904704913*^9}, {3.938699086990196*^9, 3.938699159573947*^9}, {
   3.938699339950082*^9, 3.938699341358225*^9}, {3.9386994750485697`*^9, 
   3.938699494249486*^9}, {3.938699531866549*^9, 3.938699535620966*^9}, {
   3.938699667098558*^9, 3.938699672954994*^9}, {3.938699761932259*^9, 
   3.938699764253763*^9}, {3.938699814676527*^9, 3.938699815460908*^9}, {
   3.93869985436761*^9, 3.9386999020240717`*^9}, {3.938699982900897*^9, 
   3.938699982984728*^9}, {3.938700019958934*^9, 3.938700029276306*^9}, {
   3.938700088430114*^9, 3.938700097518919*^9}, {3.9387001381017637`*^9, 
   3.938700143011744*^9}, {3.938700180267355*^9, 3.938700244462339*^9}, {
   3.9387002804413767`*^9, 3.938700296615444*^9}, {3.938700330903137*^9, 
   3.938700368186722*^9}, {3.938707092233663*^9, 3.938707145344391*^9}, {
   3.938707215706867*^9, 3.938707312871168*^9}, {3.938707492232202*^9, 
   3.938707519961013*^9}, {3.9387139133425694`*^9, 3.938713920243813*^9}, {
   3.938713983060566*^9, 3.938714011277969*^9}, {3.938717138995234*^9, 
   3.9387171522095613`*^9}, {3.938724225703261*^9, 3.9387242388997183`*^9}, {
   3.938724373559597*^9, 3.9387243873708067`*^9}, {3.938724480077979*^9, 
   3.938724483174242*^9}, {3.938724522931614*^9, 3.9387245436647463`*^9}, {
   3.93872463310157*^9, 3.938724646294544*^9}, {3.9388060222723503`*^9, 
   3.938806040484552*^9}, 3.9388065797174263`*^9, 3.938806672288148*^9, {
   3.9388069524617434`*^9, 3.938806994396418*^9}, 3.93880704382791*^9, {
   3.9388074700072823`*^9, 3.938807481689423*^9}, {3.938807539811397*^9, 
   3.9388075520229397`*^9}, 3.9388076227264853`*^9, {3.938807764853909*^9, 
   3.938807767281958*^9}, {3.938807879784671*^9, 3.9388078874895277`*^9}, {
   3.9388079748167677`*^9, 3.9388079773373337`*^9}, {3.9388081183609943`*^9, 
   3.938808118822205*^9}, {3.93880826349724*^9, 3.938808265392721*^9}, {
   3.9388095893219833`*^9, 3.938809601145488*^9}, {3.938809798503036*^9, 
   3.938809857655003*^9}, {3.938810059318918*^9, 3.9388100868812733`*^9}, {
   3.938811125831841*^9, 3.938811190704961*^9}, {3.9388112387972717`*^9, 
   3.9388113445393047`*^9}, {3.938811520619246*^9, 3.938811567602542*^9}, {
   3.938811707394759*^9, 3.938811741659225*^9}, {3.9388147074682617`*^9, 
   3.93881470766185*^9}, 3.9388166476793203`*^9, {3.938817473907086*^9, 
   3.9388174918871117`*^9}, {3.938817580996648*^9, 3.938817630652711*^9}, {
   3.9388177761983137`*^9, 3.938817777852166*^9}, {3.938817851886751*^9, 
   3.9388178562611723`*^9}, {3.9388179199259253`*^9, 3.938818063725128*^9}, {
   3.938818122402212*^9, 3.938818123819234*^9}, {3.938818229096074*^9, 
   3.9388182298984423`*^9}, {3.939770497363111*^9, 3.939770526249392*^9}},
 CellLabel->
  "In[787]:=",ExpressionUUID->"c387f004-7059-4de4-92e0-ddeb48546c32"],

Cell[BoxData["0.0007199999999999999`"], "Output",
 CellChangeTimes->{
  3.9386993435075912`*^9, 3.938699496785405*^9, 3.938699680589446*^9, 
   3.938699768659811*^9, 3.9386998180243273`*^9, 3.938699864049407*^9, 
   3.9386999043412647`*^9, 3.938699987126308*^9, {3.938700022290429*^9, 
   3.938700037164534*^9}, 3.938700099555936*^9, 3.9387001452733088`*^9, 
   3.938700248885579*^9, 3.93870033813731*^9, 3.938700370290875*^9, 
   3.938707638223526*^9, 3.938713922998406*^9, 3.9387140148682747`*^9, 
   3.938724269330748*^9, 3.938724396103724*^9, 3.938724488440359*^9, 
   3.938724546678544*^9, 3.9387246482973013`*^9, 3.938772452339512*^9, 
   3.938806593423192*^9, 3.938806695697021*^9, 3.938807631176777*^9, {
   3.938807770268972*^9, 3.938807797517008*^9}, 3.938807891147723*^9, 
   3.93880806479662*^9, 3.938808122366763*^9, 3.9388082813046923`*^9, 
   3.9388098303431997`*^9, {3.938809863725388*^9, 3.93880987656607*^9}, {
   3.938810078448085*^9, 3.938810104391127*^9}, 3.938810179289878*^9, 
   3.938811974316352*^9, 3.938814710090949*^9, 3.938816649996491*^9, 
   3.938817633049756*^9, 3.9388177819587717`*^9, 3.9388178577747927`*^9, 
   3.938818077210712*^9, 3.938818127214962*^9, 3.938865457490511*^9, 
   3.939770533295917*^9},
 CellLabel->
  "Out[802]=",ExpressionUUID->"c78c004f-444d-4c3c-bed2-bc3fc65dccbd"],

Cell[CellGroupData[{

Cell[BoxData["0.0144`"], "Print",
 CellChangeTimes->{
  3.938699346884555*^9, 3.938699500191321*^9, 3.938699683985469*^9, 
   3.9386997719937553`*^9, 3.938699821379321*^9, 3.9386998673420267`*^9, 
   3.938699907690009*^9, 3.938699990438398*^9, {3.938700025571706*^9, 
   3.938700040706686*^9}, 3.938700102802505*^9, 3.9387001485746593`*^9, 
   3.938700252183187*^9, 3.938700341370242*^9, 3.938700373733715*^9, 
   3.938707641584076*^9, 3.938713926276226*^9, 3.93871401816787*^9, 
   3.938724272787839*^9, 3.938724399443666*^9, 3.93872449171428*^9, 
   3.9387245499196043`*^9, 3.9387246515740833`*^9, 3.938772455733667*^9, 
   3.938806596778872*^9, 3.938806698948567*^9, 3.938807634656785*^9, {
   3.93880777361756*^9, 3.938807800862137*^9}, 3.938807894481145*^9, 
   3.938808068068837*^9, 3.938808125648101*^9, 3.938808284623829*^9, 
   3.938809833823061*^9, {3.938809867155486*^9, 3.938809879839347*^9}, {
   3.938810081729137*^9, 3.938810107729225*^9}, 3.93881018283461*^9, 
   3.938811977714941*^9, 3.9388147134279623`*^9, 3.93881665328976*^9, 
   3.938817637292341*^9, 3.938817786809588*^9, 3.938817857776849*^9, 
   3.938818077212473*^9, 3.938818127216411*^9, 3.938865457491803*^9, 
   3.939770533297275*^9},
 CellLabel->
  "During evaluation of \
In[787]:=",ExpressionUUID->"5f9e968b-475f-4f11-a53d-45b1dfb70a03"],

Cell[BoxData["3599.9999999999995`"], "Print",
 CellChangeTimes->{
  3.938699346884555*^9, 3.938699500191321*^9, 3.938699683985469*^9, 
   3.9386997719937553`*^9, 3.938699821379321*^9, 3.9386998673420267`*^9, 
   3.938699907690009*^9, 3.938699990438398*^9, {3.938700025571706*^9, 
   3.938700040706686*^9}, 3.938700102802505*^9, 3.9387001485746593`*^9, 
   3.938700252183187*^9, 3.938700341370242*^9, 3.938700373733715*^9, 
   3.938707641584076*^9, 3.938713926276226*^9, 3.93871401816787*^9, 
   3.938724272787839*^9, 3.938724399443666*^9, 3.93872449171428*^9, 
   3.9387245499196043`*^9, 3.9387246515740833`*^9, 3.938772455733667*^9, 
   3.938806596778872*^9, 3.938806698948567*^9, 3.938807634656785*^9, {
   3.93880777361756*^9, 3.938807800862137*^9}, 3.938807894481145*^9, 
   3.938808068068837*^9, 3.938808125648101*^9, 3.938808284623829*^9, 
   3.938809833823061*^9, {3.938809867155486*^9, 3.938809879839347*^9}, {
   3.938810081729137*^9, 3.938810107729225*^9}, 3.93881018283461*^9, 
   3.938811977714941*^9, 3.9388147134279623`*^9, 3.93881665328976*^9, 
   3.938817637292341*^9, 3.938817786809588*^9, 3.938817857776849*^9, 
   3.938818077212473*^9, 3.938818127216411*^9, 3.938865457491803*^9, 
   3.9397705366515903`*^9},
 CellLabel->
  "During evaluation of \
In[787]:=",ExpressionUUID->"244b1158-e79b-4be1-92a5-de83b1ac5953"],

Cell[BoxData["\<\"initial costates thrust:\"\>"], "Print",
 CellChangeTimes->{
  3.938699346884555*^9, 3.938699500191321*^9, 3.938699683985469*^9, 
   3.9386997719937553`*^9, 3.938699821379321*^9, 3.9386998673420267`*^9, 
   3.938699907690009*^9, 3.938699990438398*^9, {3.938700025571706*^9, 
   3.938700040706686*^9}, 3.938700102802505*^9, 3.9387001485746593`*^9, 
   3.938700252183187*^9, 3.938700341370242*^9, 3.938700373733715*^9, 
   3.938707641584076*^9, 3.938713926276226*^9, 3.93871401816787*^9, 
   3.938724272787839*^9, 3.938724399443666*^9, 3.93872449171428*^9, 
   3.9387245499196043`*^9, 3.9387246515740833`*^9, 3.938772455733667*^9, 
   3.938806596778872*^9, 3.938806698948567*^9, 3.938807634656785*^9, {
   3.93880777361756*^9, 3.938807800862137*^9}, 3.938807894481145*^9, 
   3.938808068068837*^9, 3.938808125648101*^9, 3.938808284623829*^9, 
   3.938809833823061*^9, {3.938809867155486*^9, 3.938809879839347*^9}, {
   3.938810081729137*^9, 3.938810107729225*^9}, 3.93881018283461*^9, 
   3.938811977714941*^9, 3.9388147134279623`*^9, 3.93881665328976*^9, 
   3.938817637292341*^9, 3.938817786809588*^9, 3.938817857776849*^9, 
   3.938818077212473*^9, 3.938818127216411*^9, 3.938865457491803*^9, 
   3.9397705368710537`*^9},
 CellLabel->
  "During evaluation of \
In[787]:=",ExpressionUUID->"082f52d2-1de5-44fa-838a-e5d18fd52213"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "9.852851753752176`"}], ",", "0.7071067811865475`", ",", 
   "0.`", ",", 
   RowBox[{"-", "3182.413052702611`"}], ",", 
   RowBox[{"-", "5106.358162474497`"}], ",", "0.`", ",", 
   RowBox[{"-", "13324.351357950318`"}]}], "}"}]], "Print",
 CellChangeTimes->{
  3.938699346884555*^9, 3.938699500191321*^9, 3.938699683985469*^9, 
   3.9386997719937553`*^9, 3.938699821379321*^9, 3.9386998673420267`*^9, 
   3.938699907690009*^9, 3.938699990438398*^9, {3.938700025571706*^9, 
   3.938700040706686*^9}, 3.938700102802505*^9, 3.9387001485746593`*^9, 
   3.938700252183187*^9, 3.938700341370242*^9, 3.938700373733715*^9, 
   3.938707641584076*^9, 3.938713926276226*^9, 3.93871401816787*^9, 
   3.938724272787839*^9, 3.938724399443666*^9, 3.93872449171428*^9, 
   3.9387245499196043`*^9, 3.9387246515740833`*^9, 3.938772455733667*^9, 
   3.938806596778872*^9, 3.938806698948567*^9, 3.938807634656785*^9, {
   3.93880777361756*^9, 3.938807800862137*^9}, 3.938807894481145*^9, 
   3.938808068068837*^9, 3.938808125648101*^9, 3.938808284623829*^9, 
   3.938809833823061*^9, {3.938809867155486*^9, 3.938809879839347*^9}, {
   3.938810081729137*^9, 3.938810107729225*^9}, 3.93881018283461*^9, 
   3.938811977714941*^9, 3.9388147134279623`*^9, 3.93881665328976*^9, 
   3.938817637292341*^9, 3.938817786809588*^9, 3.938817857776849*^9, 
   3.938818077212473*^9, 3.938818127216411*^9, 3.938865457491803*^9, 
   3.9397705368719673`*^9},
 CellLabel->
  "During evaluation of \
In[787]:=",ExpressionUUID->"9359be76-0a4b-40d2-8b6a-99517fa60fc3"],

Cell[BoxData["\<\"final states for energy constrained trajectory:\"\>"], \
"Print",
 CellChangeTimes->{
  3.938699346884555*^9, 3.938699500191321*^9, 3.938699683985469*^9, 
   3.9386997719937553`*^9, 3.938699821379321*^9, 3.9386998673420267`*^9, 
   3.938699907690009*^9, 3.938699990438398*^9, {3.938700025571706*^9, 
   3.938700040706686*^9}, 3.938700102802505*^9, 3.9387001485746593`*^9, 
   3.938700252183187*^9, 3.938700341370242*^9, 3.938700373733715*^9, 
   3.938707641584076*^9, 3.938713926276226*^9, 3.93871401816787*^9, 
   3.938724272787839*^9, 3.938724399443666*^9, 3.93872449171428*^9, 
   3.9387245499196043`*^9, 3.9387246515740833`*^9, 3.938772455733667*^9, 
   3.938806596778872*^9, 3.938806698948567*^9, 3.938807634656785*^9, {
   3.93880777361756*^9, 3.938807800862137*^9}, 3.938807894481145*^9, 
   3.938808068068837*^9, 3.938808125648101*^9, 3.938808284623829*^9, 
   3.938809833823061*^9, {3.938809867155486*^9, 3.938809879839347*^9}, {
   3.938810081729137*^9, 3.938810107729225*^9}, 3.93881018283461*^9, 
   3.938811977714941*^9, 3.9388147134279623`*^9, 3.93881665328976*^9, 
   3.938817637292341*^9, 3.938817786809588*^9, 3.938817857776849*^9, 
   3.938818077212473*^9, 3.938818127216411*^9, 3.938865457491803*^9, 
   3.939770537108231*^9},
 CellLabel->
  "During evaluation of \
In[787]:=",ExpressionUUID->"ba843332-dd83-41ca-8de1-89a37b08414f"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "5249.489249956984`"}], ",", "28751.73450234287`", ",", 
   "146.69570883699035`", ",", "6.307576970560026`", ",", 
   "10.094547795031076`", ",", "0.001235440334316833`"}], "}"}]], "Output",
 CellChangeTimes->{
  3.9386993435075912`*^9, 3.938699496785405*^9, 3.938699680589446*^9, 
   3.938699768659811*^9, 3.9386998180243273`*^9, 3.938699864049407*^9, 
   3.9386999043412647`*^9, 3.938699987126308*^9, {3.938700022290429*^9, 
   3.938700037164534*^9}, 3.938700099555936*^9, 3.9387001452733088`*^9, 
   3.938700248885579*^9, 3.93870033813731*^9, 3.938700370290875*^9, 
   3.938707638223526*^9, 3.938713922998406*^9, 3.9387140148682747`*^9, 
   3.938724269330748*^9, 3.938724396103724*^9, 3.938724488440359*^9, 
   3.938724546678544*^9, 3.9387246482973013`*^9, 3.938772452339512*^9, 
   3.938806593423192*^9, 3.938806695697021*^9, 3.938807631176777*^9, {
   3.938807770268972*^9, 3.938807797517008*^9}, 3.938807891147723*^9, 
   3.93880806479662*^9, 3.938808122366763*^9, 3.9388082813046923`*^9, 
   3.9388098303431997`*^9, {3.938809863725388*^9, 3.93880987656607*^9}, {
   3.938810078448085*^9, 3.938810104391127*^9}, 3.938810179289878*^9, 
   3.938811974316352*^9, 3.938814710090949*^9, 3.938816649996491*^9, 
   3.938817633049756*^9, 3.9388177819587717`*^9, 3.9388178577747927`*^9, 
   3.938818077210712*^9, 3.938818127214962*^9, 3.938865457490511*^9, 
   3.939770540407753*^9},
 CellLabel->
  "Out[815]=",ExpressionUUID->"6acaa1b0-ba4b-4c8e-8805-5359733492af"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "0.1796107564922578`"}], ",", "0.9837377578157092`"}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.9386993435075912`*^9, 3.938699496785405*^9, 3.938699680589446*^9, 
   3.938699768659811*^9, 3.9386998180243273`*^9, 3.938699864049407*^9, 
   3.9386999043412647`*^9, 3.938699987126308*^9, {3.938700022290429*^9, 
   3.938700037164534*^9}, 3.938700099555936*^9, 3.9387001452733088`*^9, 
   3.938700248885579*^9, 3.93870033813731*^9, 3.938700370290875*^9, 
   3.938707638223526*^9, 3.938713922998406*^9, 3.9387140148682747`*^9, 
   3.938724269330748*^9, 3.938724396103724*^9, 3.938724488440359*^9, 
   3.938724546678544*^9, 3.9387246482973013`*^9, 3.938772452339512*^9, 
   3.938806593423192*^9, 3.938806695697021*^9, 3.938807631176777*^9, {
   3.938807770268972*^9, 3.938807797517008*^9}, 3.938807891147723*^9, 
   3.93880806479662*^9, 3.938808122366763*^9, 3.9388082813046923`*^9, 
   3.9388098303431997`*^9, {3.938809863725388*^9, 3.93880987656607*^9}, {
   3.938810078448085*^9, 3.938810104391127*^9}, 3.938810179289878*^9, 
   3.938811974316352*^9, 3.938814710090949*^9, 3.938816649996491*^9, 
   3.938817633049756*^9, 3.9388177819587717`*^9, 3.9388178577747927`*^9, 
   3.938818077210712*^9, 3.938818127214962*^9, 3.938865457490511*^9, 
   3.939770546895418*^9},
 CellLabel->
  "Out[816]=",ExpressionUUID->"0c7184e8-70d2-4f47-b93e-2a0285cf4161"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "4209.164974380354`"}], ",", "23053.822583013472`", ",", 
   "0.`", ",", "7.002512208981952`", ",", "10.015743461266341`", ",", "0.`", 
   ",", "0.7071059113621312`", ",", "0.7071067811865475`", ",", "0.`", ",", 
   RowBox[{"-", "0.000029307517870378246`"}], ",", 
   RowBox[{"-", "0.0006302909488393995`"}], ",", "0.`"}], "}"}]], "Output",
 CellChangeTimes->{
  3.9386993435075912`*^9, 3.938699496785405*^9, 3.938699680589446*^9, 
   3.938699768659811*^9, 3.9386998180243273`*^9, 3.938699864049407*^9, 
   3.9386999043412647`*^9, 3.938699987126308*^9, {3.938700022290429*^9, 
   3.938700037164534*^9}, 3.938700099555936*^9, 3.9387001452733088`*^9, 
   3.938700248885579*^9, 3.93870033813731*^9, 3.938700370290875*^9, 
   3.938707638223526*^9, 3.938713922998406*^9, 3.9387140148682747`*^9, 
   3.938724269330748*^9, 3.938724396103724*^9, 3.938724488440359*^9, 
   3.938724546678544*^9, 3.9387246482973013`*^9, 3.938772452339512*^9, 
   3.938806593423192*^9, 3.938806695697021*^9, 3.938807631176777*^9, {
   3.938807770268972*^9, 3.938807797517008*^9}, 3.938807891147723*^9, 
   3.93880806479662*^9, 3.938808122366763*^9, 3.9388082813046923`*^9, 
   3.9388098303431997`*^9, {3.938809863725388*^9, 3.93880987656607*^9}, {
   3.938810078448085*^9, 3.938810104391127*^9}, 3.938810179289878*^9, 
   3.938811974316352*^9, 3.938814710090949*^9, 3.938816649996491*^9, 
   3.938817633049756*^9, 3.9388177819587717`*^9, 3.9388178577747927`*^9, 
   3.938818077210712*^9, 3.938818127214962*^9, 3.938865457490511*^9, 
   3.939770546896467*^9},
 CellLabel->
  "Out[817]=",ExpressionUUID->"de398d5d-f758-4747-ad5d-ab409793fc30"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "0.17961075649225777`"}], ",", "0.983737757815709`"}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.9386993435075912`*^9, 3.938699496785405*^9, 3.938699680589446*^9, 
   3.938699768659811*^9, 3.9386998180243273`*^9, 3.938699864049407*^9, 
   3.9386999043412647`*^9, 3.938699987126308*^9, {3.938700022290429*^9, 
   3.938700037164534*^9}, 3.938700099555936*^9, 3.9387001452733088`*^9, 
   3.938700248885579*^9, 3.93870033813731*^9, 3.938700370290875*^9, 
   3.938707638223526*^9, 3.938713922998406*^9, 3.9387140148682747`*^9, 
   3.938724269330748*^9, 3.938724396103724*^9, 3.938724488440359*^9, 
   3.938724546678544*^9, 3.9387246482973013`*^9, 3.938772452339512*^9, 
   3.938806593423192*^9, 3.938806695697021*^9, 3.938807631176777*^9, {
   3.938807770268972*^9, 3.938807797517008*^9}, 3.938807891147723*^9, 
   3.93880806479662*^9, 3.938808122366763*^9, 3.9388082813046923`*^9, 
   3.9388098303431997`*^9, {3.938809863725388*^9, 3.93880987656607*^9}, {
   3.938810078448085*^9, 3.938810104391127*^9}, 3.938810179289878*^9, 
   3.938811974316352*^9, 3.938814710090949*^9, 3.938816649996491*^9, 
   3.938817633049756*^9, 3.9388177819587717`*^9, 3.9388178577747927`*^9, 
   3.938818077210712*^9, 3.938818127214962*^9, 3.938865457490511*^9, 
   3.9397705468973494`*^9},
 CellLabel->
  "Out[818]=",ExpressionUUID->"c4ca6d6d-984f-467e-8bee-ba69b59e69a4"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "4208.617318802811`"}], ",", "23057.1527117401`", ",", "0.`", 
   ",", "7.0055300406346666`", ",", "10.016496627224155`", ",", "0.`", ",", 
   "0.9992800001862455`", ",", "0.7071059183810939`", ",", 
   "0.7071067811865475`", ",", "0.`", ",", 
   RowBox[{"-", "0.00002827787565999185`"}], ",", 
   RowBox[{"-", "0.0006252048896529178`"}], ",", "0.`", ",", 
   "6.315291827962259`"}], "}"}]], "Output",
 CellChangeTimes->{
  3.9386993435075912`*^9, 3.938699496785405*^9, 3.938699680589446*^9, 
   3.938699768659811*^9, 3.9386998180243273`*^9, 3.938699864049407*^9, 
   3.9386999043412647`*^9, 3.938699987126308*^9, {3.938700022290429*^9, 
   3.938700037164534*^9}, 3.938700099555936*^9, 3.9387001452733088`*^9, 
   3.938700248885579*^9, 3.93870033813731*^9, 3.938700370290875*^9, 
   3.938707638223526*^9, 3.938713922998406*^9, 3.9387140148682747`*^9, 
   3.938724269330748*^9, 3.938724396103724*^9, 3.938724488440359*^9, 
   3.938724546678544*^9, 3.9387246482973013`*^9, 3.938772452339512*^9, 
   3.938806593423192*^9, 3.938806695697021*^9, 3.938807631176777*^9, {
   3.938807770268972*^9, 3.938807797517008*^9}, 3.938807891147723*^9, 
   3.93880806479662*^9, 3.938808122366763*^9, 3.9388082813046923`*^9, 
   3.9388098303431997`*^9, {3.938809863725388*^9, 3.93880987656607*^9}, {
   3.938810078448085*^9, 3.938810104391127*^9}, 3.938810179289878*^9, 
   3.938811974316352*^9, 3.938814710090949*^9, 3.938816649996491*^9, 
   3.938817633049756*^9, 3.9388177819587717`*^9, 3.9388178577747927`*^9, 
   3.938818077210712*^9, 3.938818127214962*^9, 3.938865457490511*^9, 
   3.939770546899796*^9},
 CellLabel->
  "Out[819]=",ExpressionUUID->"fd40a038-ce73-4160-8fce-f06c91c19535"],

Cell[CellGroupData[{

Cell[BoxData["\<\"hello\"\>"], "Print",
 CellChangeTimes->{
  3.938699346884555*^9, 3.938699500191321*^9, 3.938699683985469*^9, 
   3.9386997719937553`*^9, 3.938699821379321*^9, 3.9386998673420267`*^9, 
   3.938699907690009*^9, 3.938699990438398*^9, {3.938700025571706*^9, 
   3.938700040706686*^9}, 3.938700102802505*^9, 3.9387001485746593`*^9, 
   3.938700252183187*^9, 3.938700341370242*^9, 3.938700373733715*^9, 
   3.938707641584076*^9, 3.938713926276226*^9, 3.93871401816787*^9, 
   3.938724272787839*^9, 3.938724399443666*^9, 3.93872449171428*^9, 
   3.9387245499196043`*^9, 3.9387246515740833`*^9, 3.938772455733667*^9, 
   3.938806596778872*^9, 3.938806698948567*^9, 3.938807634656785*^9, {
   3.93880777361756*^9, 3.938807800862137*^9}, 3.938807894481145*^9, 
   3.938808068068837*^9, 3.938808125648101*^9, 3.938808284623829*^9, 
   3.938809833823061*^9, {3.938809867155486*^9, 3.938809879839347*^9}, {
   3.938810081729137*^9, 3.938810107729225*^9}, 3.93881018283461*^9, 
   3.938811977714941*^9, 3.9388147134279623`*^9, 3.93881665328976*^9, 
   3.938817637292341*^9, 3.938817786809588*^9, 3.938817857776849*^9, 
   3.938818077212473*^9, 3.938818127216411*^9, 3.938865457491803*^9, 
   3.939770546900807*^9},
 CellLabel->
  "During evaluation of \
In[787]:=",ExpressionUUID->"8070c84e-345e-48d3-8ad3-6fdfd47b706a"],

Cell[BoxData[
 RowBox[{"\<\"initial costates of thrust trajectory from root finding\"\>", 
  " ", "f"}]], "Print",
 CellChangeTimes->{
  3.938699346884555*^9, 3.938699500191321*^9, 3.938699683985469*^9, 
   3.9386997719937553`*^9, 3.938699821379321*^9, 3.9386998673420267`*^9, 
   3.938699907690009*^9, 3.938699990438398*^9, {3.938700025571706*^9, 
   3.938700040706686*^9}, 3.938700102802505*^9, 3.9387001485746593`*^9, 
   3.938700252183187*^9, 3.938700341370242*^9, 3.938700373733715*^9, 
   3.938707641584076*^9, 3.938713926276226*^9, 3.93871401816787*^9, 
   3.938724272787839*^9, 3.938724399443666*^9, 3.93872449171428*^9, 
   3.9387245499196043`*^9, 3.9387246515740833`*^9, 3.938772455733667*^9, 
   3.938806596778872*^9, 3.938806698948567*^9, 3.938807634656785*^9, {
   3.93880777361756*^9, 3.938807800862137*^9}, 3.938807894481145*^9, 
   3.938808068068837*^9, 3.938808125648101*^9, 3.938808284623829*^9, 
   3.938809833823061*^9, {3.938809867155486*^9, 3.938809879839347*^9}, {
   3.938810081729137*^9, 3.938810107729225*^9}, 3.93881018283461*^9, 
   3.938811977714941*^9, 3.9388147134279623`*^9, 3.93881665328976*^9, 
   3.938817637292341*^9, 3.938817786809588*^9, 3.938817857776849*^9, 
   3.938818077212473*^9, 3.938818127216411*^9, 3.938865457491803*^9, 
   3.939770546901627*^9},
 CellLabel->
  "During evaluation of \
In[787]:=",ExpressionUUID->"9cb42d38-ffa1-4020-9edd-2446261d20c5"]
}, Open  ]],

Cell[BoxData[
 TemplateBox[{
  "FindRoot", "jsing", 
   "\"Encountered a singular Jacobian at the point \
\\!\\(\\*RowBox[{\\\"{\\\", \\\"l\\\", \\\"}\\\"}]\\) = \
\\!\\(\\*RowBox[{\\\"{\\\", RowBox[{RowBox[{\\\"-\\\", \\\"9.852851753752176`\
\\\"}], \\\",\\\", \\\"0.7071067811865475`\\\", \\\",\\\", \\\"0.`\\\", \\\",\
\\\", RowBox[{\\\"-\\\", \\\"3182.413052702611`\\\"}], \\\",\\\", \
RowBox[{\\\"-\\\", \\\"5106.358162474497`\\\"}], \\\",\\\", \\\"0.`\\\", \
\\\",\\\", RowBox[{\\\"-\\\", \\\"13324.351357950318`\\\"}]}], \
\\\"}\\\"}]\\). Try perturbing the initial point(s).\"", 2, 824, 750, 
   31053567148927707368, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.938810118065565*^9, 3.938810193175585*^9, 
  3.938811988244615*^9, 3.9388147239995213`*^9, 3.93881666373302*^9, 
  3.938817637284403*^9, 3.9388178004299793`*^9, 3.938817875563455*^9, 
  3.9388180825645933`*^9, 3.938818149086329*^9, 3.938865471426333*^9, 
  3.939770546931555*^9},
 CellLabel->
  "During evaluation of \
In[787]:=",ExpressionUUID->"e239a939-f47a-464e-bb9c-dd59710e2979"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"l", "\[Rule]", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "9.852851753752176`"}], ",", "0.7071067811865475`", ",", 
     "0.`", ",", 
     RowBox[{"-", "3182.413052702611`"}], ",", 
     RowBox[{"-", "5106.358162474497`"}], ",", "0.`", ",", 
     RowBox[{"-", "13324.351357950318`"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.9386993435075912`*^9, 3.938699496785405*^9, 3.938699680589446*^9, 
   3.938699768659811*^9, 3.9386998180243273`*^9, 3.938699864049407*^9, 
   3.9386999043412647`*^9, 3.938699987126308*^9, {3.938700022290429*^9, 
   3.938700037164534*^9}, 3.938700099555936*^9, 3.9387001452733088`*^9, 
   3.938700248885579*^9, 3.93870033813731*^9, 3.938700370290875*^9, 
   3.938707638223526*^9, 3.938713922998406*^9, 3.9387140148682747`*^9, 
   3.938724269330748*^9, 3.938724396103724*^9, 3.938724488440359*^9, 
   3.938724546678544*^9, 3.9387246482973013`*^9, 3.938772452339512*^9, 
   3.938806593423192*^9, 3.938806695697021*^9, 3.938807631176777*^9, {
   3.938807770268972*^9, 3.938807797517008*^9}, 3.938807891147723*^9, 
   3.93880806479662*^9, 3.938808122366763*^9, 3.9388082813046923`*^9, 
   3.9388098303431997`*^9, {3.938809863725388*^9, 3.93880987656607*^9}, {
   3.938810078448085*^9, 3.938810104391127*^9}, 3.938810179289878*^9, 
   3.938811974316352*^9, 3.938814710090949*^9, 3.938816649996491*^9, 
   3.938817633049756*^9, 3.9388177819587717`*^9, 3.9388178577747927`*^9, 
   3.938818077210712*^9, 3.938818127214962*^9, 3.938865457490511*^9, 
   3.9397705469387074`*^9},
 CellLabel->
  "Out[824]=",ExpressionUUID->"0d8a9c4a-62f2-443d-8cec-7fb76c722e24"],

Cell[BoxData["\<\"final states of thrust-constrained trajectory:\"\>"], \
"Print",
 CellChangeTimes->{
  3.938699346884555*^9, 3.938699500191321*^9, 3.938699683985469*^9, 
   3.9386997719937553`*^9, 3.938699821379321*^9, 3.9386998673420267`*^9, 
   3.938699907690009*^9, 3.938699990438398*^9, {3.938700025571706*^9, 
   3.938700040706686*^9}, 3.938700102802505*^9, 3.9387001485746593`*^9, 
   3.938700252183187*^9, 3.938700341370242*^9, 3.938700373733715*^9, 
   3.938707641584076*^9, 3.938713926276226*^9, 3.93871401816787*^9, 
   3.938724272787839*^9, 3.938724399443666*^9, 3.93872449171428*^9, 
   3.9387245499196043`*^9, 3.9387246515740833`*^9, 3.938772455733667*^9, 
   3.938806596778872*^9, 3.938806698948567*^9, 3.938807634656785*^9, {
   3.93880777361756*^9, 3.938807800862137*^9}, 3.938807894481145*^9, 
   3.938808068068837*^9, 3.938808125648101*^9, 3.938808284623829*^9, 
   3.938809833823061*^9, {3.938809867155486*^9, 3.938809879839347*^9}, {
   3.938810081729137*^9, 3.938810107729225*^9}, 3.93881018283461*^9, 
   3.938811977714941*^9, 3.9388147134279623`*^9, 3.93881665328976*^9, 
   3.938817637292341*^9, 3.938817786809588*^9, 3.938817857776849*^9, 
   3.938818077212473*^9, 3.938818127216411*^9, 3.938865457491803*^9, 
   3.939770546939466*^9},
 CellLabel->
  "During evaluation of \
In[787]:=",ExpressionUUID->"08417164-06f0-4485-8e4b-4b3ad0ca7bee"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "4208.617318802811`"}], ",", "23057.1527117401`", ",", "0.`", 
   ",", "7.0055300406346666`", ",", "10.016496627224155`", ",", "0.`", ",", 
   "0.9992800001862455`", ",", "0.7071059183810939`", ",", 
   "0.7071067811865475`", ",", "0.`", ",", 
   RowBox[{"-", "0.00002827787565999185`"}], ",", 
   RowBox[{"-", "0.0006252048896529178`"}], ",", "0.`", ",", 
   "6.315291827962259`"}], "}"}]], "Output",
 CellChangeTimes->{
  3.9386993435075912`*^9, 3.938699496785405*^9, 3.938699680589446*^9, 
   3.938699768659811*^9, 3.9386998180243273`*^9, 3.938699864049407*^9, 
   3.9386999043412647`*^9, 3.938699987126308*^9, {3.938700022290429*^9, 
   3.938700037164534*^9}, 3.938700099555936*^9, 3.9387001452733088`*^9, 
   3.938700248885579*^9, 3.93870033813731*^9, 3.938700370290875*^9, 
   3.938707638223526*^9, 3.938713922998406*^9, 3.9387140148682747`*^9, 
   3.938724269330748*^9, 3.938724396103724*^9, 3.938724488440359*^9, 
   3.938724546678544*^9, 3.9387246482973013`*^9, 3.938772452339512*^9, 
   3.938806593423192*^9, 3.938806695697021*^9, 3.938807631176777*^9, {
   3.938807770268972*^9, 3.938807797517008*^9}, 3.938807891147723*^9, 
   3.93880806479662*^9, 3.938808122366763*^9, 3.9388082813046923`*^9, 
   3.9388098303431997`*^9, {3.938809863725388*^9, 3.93880987656607*^9}, {
   3.938810078448085*^9, 3.938810104391127*^9}, 3.938810179289878*^9, 
   3.938811974316352*^9, 3.938814710090949*^9, 3.938816649996491*^9, 
   3.938817633049756*^9, 3.9388177819587717`*^9, 3.9388178577747927`*^9, 
   3.938818077210712*^9, 3.938818127214962*^9, 3.938865457490511*^9, 
   3.939770546940291*^9},
 CellLabel->
  "Out[826]=",ExpressionUUID->"858755d7-6d44-4408-86c4-99a10d441055"],

Cell[BoxData["\<\"final constraint values of thrust-constrained trajectory:\"\
\>"], "Print",
 CellChangeTimes->{
  3.938699346884555*^9, 3.938699500191321*^9, 3.938699683985469*^9, 
   3.9386997719937553`*^9, 3.938699821379321*^9, 3.9386998673420267`*^9, 
   3.938699907690009*^9, 3.938699990438398*^9, {3.938700025571706*^9, 
   3.938700040706686*^9}, 3.938700102802505*^9, 3.9387001485746593`*^9, 
   3.938700252183187*^9, 3.938700341370242*^9, 3.938700373733715*^9, 
   3.938707641584076*^9, 3.938713926276226*^9, 3.93871401816787*^9, 
   3.938724272787839*^9, 3.938724399443666*^9, 3.93872449171428*^9, 
   3.9387245499196043`*^9, 3.9387246515740833`*^9, 3.938772455733667*^9, 
   3.938806596778872*^9, 3.938806698948567*^9, 3.938807634656785*^9, {
   3.93880777361756*^9, 3.938807800862137*^9}, 3.938807894481145*^9, 
   3.938808068068837*^9, 3.938808125648101*^9, 3.938808284623829*^9, 
   3.938809833823061*^9, {3.938809867155486*^9, 3.938809879839347*^9}, {
   3.938810081729137*^9, 3.938810107729225*^9}, 3.93881018283461*^9, 
   3.938811977714941*^9, 3.9388147134279623`*^9, 3.93881665328976*^9, 
   3.938817637292341*^9, 3.938817786809588*^9, 3.938817857776849*^9, 
   3.938818077212473*^9, 3.938818127216411*^9, 3.938865457491803*^9, 
   3.939770546941045*^9},
 CellLabel->
  "During evaluation of \
In[787]:=",ExpressionUUID->"7b821ad8-ac4a-43d2-89a1-c9c66111c615"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.00004771666760418324`", ",", 
   RowBox[{"-", "0.9999757401966436`"}], ",", "0.`", ",", 
   RowBox[{"-", "0.00002827787565999185`"}], ",", 
   RowBox[{"-", "0.0006252048896529178`"}], ",", "0.`", ",", 
   RowBox[{"-", "1.8624552063344763`*^-10"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.9386993435075912`*^9, 3.938699496785405*^9, 3.938699680589446*^9, 
   3.938699768659811*^9, 3.9386998180243273`*^9, 3.938699864049407*^9, 
   3.9386999043412647`*^9, 3.938699987126308*^9, {3.938700022290429*^9, 
   3.938700037164534*^9}, 3.938700099555936*^9, 3.9387001452733088`*^9, 
   3.938700248885579*^9, 3.93870033813731*^9, 3.938700370290875*^9, 
   3.938707638223526*^9, 3.938713922998406*^9, 3.9387140148682747`*^9, 
   3.938724269330748*^9, 3.938724396103724*^9, 3.938724488440359*^9, 
   3.938724546678544*^9, 3.9387246482973013`*^9, 3.938772452339512*^9, 
   3.938806593423192*^9, 3.938806695697021*^9, 3.938807631176777*^9, {
   3.938807770268972*^9, 3.938807797517008*^9}, 3.938807891147723*^9, 
   3.93880806479662*^9, 3.938808122366763*^9, 3.9388082813046923`*^9, 
   3.9388098303431997`*^9, {3.938809863725388*^9, 3.93880987656607*^9}, {
   3.938810078448085*^9, 3.938810104391127*^9}, 3.938810179289878*^9, 
   3.938811974316352*^9, 3.938814710090949*^9, 3.938816649996491*^9, 
   3.938817633049756*^9, 3.9388177819587717`*^9, 3.9388178577747927`*^9, 
   3.938818077210712*^9, 3.938818127214962*^9, 3.938865457490511*^9, 
   3.939770546941827*^9},
 CellLabel->
  "Out[828]=",ExpressionUUID->"3d2677be-5a9d-47e4-94bf-e29a662b8695"],

Cell[BoxData["\<\"final states and costates of fuel-constrained trajectory:\"\
\>"], "Print",
 CellChangeTimes->{
  3.938699346884555*^9, 3.938699500191321*^9, 3.938699683985469*^9, 
   3.9386997719937553`*^9, 3.938699821379321*^9, 3.9386998673420267`*^9, 
   3.938699907690009*^9, 3.938699990438398*^9, {3.938700025571706*^9, 
   3.938700040706686*^9}, 3.938700102802505*^9, 3.9387001485746593`*^9, 
   3.938700252183187*^9, 3.938700341370242*^9, 3.938700373733715*^9, 
   3.938707641584076*^9, 3.938713926276226*^9, 3.93871401816787*^9, 
   3.938724272787839*^9, 3.938724399443666*^9, 3.93872449171428*^9, 
   3.9387245499196043`*^9, 3.9387246515740833`*^9, 3.938772455733667*^9, 
   3.938806596778872*^9, 3.938806698948567*^9, 3.938807634656785*^9, {
   3.93880777361756*^9, 3.938807800862137*^9}, 3.938807894481145*^9, 
   3.938808068068837*^9, 3.938808125648101*^9, 3.938808284623829*^9, 
   3.938809833823061*^9, {3.938809867155486*^9, 3.938809879839347*^9}, {
   3.938810081729137*^9, 3.938810107729225*^9}, 3.93881018283461*^9, 
   3.938811977714941*^9, 3.9388147134279623`*^9, 3.93881665328976*^9, 
   3.938817637292341*^9, 3.938817786809588*^9, 3.938817857776849*^9, 
   3.938818077212473*^9, 3.938818127216411*^9, 3.938865457491803*^9, 
   3.939770547148499*^9},
 CellLabel->
  "During evaluation of \
In[787]:=",ExpressionUUID->"507d2b23-3b46-40b8-925d-d6e4f8cddfab"],

Cell[BoxData[
 TemplateBox[{
  "FindRoot", "jsing", 
   "\"Encountered a singular Jacobian at the point \
\\!\\(\\*RowBox[{\\\"{\\\", \\\"l\\\", \\\"}\\\"}]\\) = \
\\!\\(\\*RowBox[{\\\"{\\\", RowBox[{RowBox[{\\\"-\\\", \\\"9.852851753752176`\
\\\"}], \\\",\\\", \\\"0.7071067811865475`\\\", \\\",\\\", \\\"0.`\\\", \\\",\
\\\", RowBox[{\\\"-\\\", \\\"3182.413052702611`\\\"}], \\\",\\\", \
RowBox[{\\\"-\\\", \\\"5106.358162474497`\\\"}], \\\",\\\", \\\"0.`\\\", \
\\\",\\\", \\\"1.2828481196279421`*^7\\\"}], \\\"}\\\"}]\\). Try perturbing \
the initial point(s).\"", 2, 834, 751, 31053567148927707368, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.938810118065565*^9, 3.938810193175585*^9, 
  3.938811988244615*^9, 3.9388147239995213`*^9, 3.93881666373302*^9, 
  3.938817637284403*^9, 3.9388178004299793`*^9, 3.938817875563455*^9, 
  3.9388180825645933`*^9, 3.938818149086329*^9, 3.938865471426333*^9, 
  3.939770547178808*^9},
 CellLabel->
  "During evaluation of \
In[787]:=",ExpressionUUID->"c0ab2dee-e1fb-4a45-a823-22975c223d10"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "6718.442384663182`"}], ",", "21616.15458197448`", ",", "0.`",
    ",", "2.035662959205049`", ",", "11.3541540696739`", ",", "0.`", ",", 
   "0.9996399999999442`", ",", "0.707105919446765`", ",", 
   "0.7071067811865475`", ",", "0.`", ",", 
   RowBox[{"-", "0.0000415238673103824`"}], ",", 
   RowBox[{"-", "0.0006244327106093388`"}], ",", "0.`", ",", 
   "1.2839016844045015`*^7"}], "}"}]], "Output",
 CellChangeTimes->{
  3.9386993435075912`*^9, 3.938699496785405*^9, 3.938699680589446*^9, 
   3.938699768659811*^9, 3.9386998180243273`*^9, 3.938699864049407*^9, 
   3.9386999043412647`*^9, 3.938699987126308*^9, {3.938700022290429*^9, 
   3.938700037164534*^9}, 3.938700099555936*^9, 3.9387001452733088`*^9, 
   3.938700248885579*^9, 3.93870033813731*^9, 3.938700370290875*^9, 
   3.938707638223526*^9, 3.938713922998406*^9, 3.9387140148682747`*^9, 
   3.938724269330748*^9, 3.938724396103724*^9, 3.938724488440359*^9, 
   3.938724546678544*^9, 3.9387246482973013`*^9, 3.938772452339512*^9, 
   3.938806593423192*^9, 3.938806695697021*^9, 3.938807631176777*^9, {
   3.938807770268972*^9, 3.938807797517008*^9}, 3.938807891147723*^9, 
   3.93880806479662*^9, 3.938808122366763*^9, 3.9388082813046923`*^9, 
   3.9388098303431997`*^9, {3.938809863725388*^9, 3.93880987656607*^9}, {
   3.938810078448085*^9, 3.938810104391127*^9}, 3.938810179289878*^9, 
   3.938811974316352*^9, 3.938814710090949*^9, 3.938816649996491*^9, 
   3.938817633049756*^9, 3.9388177819587717`*^9, 3.9388178577747927`*^9, 
   3.938818077210712*^9, 3.938818127214962*^9, 3.938865457490511*^9, 
   3.9397705471839333`*^9},
 CellLabel->
  "Out[835]=",ExpressionUUID->"dce4d94b-eb41-484a-9a62-258f53067104"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "0.29680136797496637`"}], ",", "0.9549392378409156`"}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.9386993435075912`*^9, 3.938699496785405*^9, 3.938699680589446*^9, 
   3.938699768659811*^9, 3.9386998180243273`*^9, 3.938699864049407*^9, 
   3.9386999043412647`*^9, 3.938699987126308*^9, {3.938700022290429*^9, 
   3.938700037164534*^9}, 3.938700099555936*^9, 3.9387001452733088`*^9, 
   3.938700248885579*^9, 3.93870033813731*^9, 3.938700370290875*^9, 
   3.938707638223526*^9, 3.938713922998406*^9, 3.9387140148682747`*^9, 
   3.938724269330748*^9, 3.938724396103724*^9, 3.938724488440359*^9, 
   3.938724546678544*^9, 3.9387246482973013`*^9, 3.938772452339512*^9, 
   3.938806593423192*^9, 3.938806695697021*^9, 3.938807631176777*^9, {
   3.938807770268972*^9, 3.938807797517008*^9}, 3.938807891147723*^9, 
   3.93880806479662*^9, 3.938808122366763*^9, 3.9388082813046923`*^9, 
   3.9388098303431997`*^9, {3.938809863725388*^9, 3.93880987656607*^9}, {
   3.938810078448085*^9, 3.938810104391127*^9}, 3.938810179289878*^9, 
   3.938811974316352*^9, 3.938814710090949*^9, 3.938816649996491*^9, 
   3.938817633049756*^9, 3.9388177819587717`*^9, 3.9388178577747927`*^9, 
   3.938818077210712*^9, 3.938818127214962*^9, 3.938865457490511*^9, 
   3.9397705471848373`*^9},
 CellLabel->
  "Out[836]=",ExpressionUUID->"6497c571-4391-4ba2-a2a9-e7498571bcb5"],

Cell[BoxData["\<\"final constraint satisfaction:\"\>"], "Print",
 CellChangeTimes->{
  3.938699346884555*^9, 3.938699500191321*^9, 3.938699683985469*^9, 
   3.9386997719937553`*^9, 3.938699821379321*^9, 3.9386998673420267`*^9, 
   3.938699907690009*^9, 3.938699990438398*^9, {3.938700025571706*^9, 
   3.938700040706686*^9}, 3.938700102802505*^9, 3.9387001485746593`*^9, 
   3.938700252183187*^9, 3.938700341370242*^9, 3.938700373733715*^9, 
   3.938707641584076*^9, 3.938713926276226*^9, 3.93871401816787*^9, 
   3.938724272787839*^9, 3.938724399443666*^9, 3.93872449171428*^9, 
   3.9387245499196043`*^9, 3.9387246515740833`*^9, 3.938772455733667*^9, 
   3.938806596778872*^9, 3.938806698948567*^9, 3.938807634656785*^9, {
   3.93880777361756*^9, 3.938807800862137*^9}, 3.938807894481145*^9, 
   3.938808068068837*^9, 3.938808125648101*^9, 3.938808284623829*^9, 
   3.938809833823061*^9, {3.938809867155486*^9, 3.938809879839347*^9}, {
   3.938810081729137*^9, 3.938810107729225*^9}, 3.93881018283461*^9, 
   3.938811977714941*^9, 3.9388147134279623`*^9, 3.93881665328976*^9, 
   3.938817637292341*^9, 3.938817786809588*^9, 3.938817857776849*^9, 
   3.938818077212473*^9, 3.938818127216411*^9, 3.938865457491803*^9, 
   3.939770547185609*^9},
 CellLabel->
  "During evaluation of \
In[787]:=",ExpressionUUID->"c19d1af2-38cf-4a8e-ad8b-690b86a9db53"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "0.11849844668275983`"}], ",", 
   RowBox[{"-", "0.9999748807269476`"}], ",", "0.`", ",", 
   RowBox[{"-", "0.0000415238673103824`"}], ",", 
   RowBox[{"-", "0.0006244327106093388`"}], ",", "0.`", ",", 
   "5.576019247027775`*^-14"}], "}"}]], "Output",
 CellChangeTimes->{
  3.9386993435075912`*^9, 3.938699496785405*^9, 3.938699680589446*^9, 
   3.938699768659811*^9, 3.9386998180243273`*^9, 3.938699864049407*^9, 
   3.9386999043412647`*^9, 3.938699987126308*^9, {3.938700022290429*^9, 
   3.938700037164534*^9}, 3.938700099555936*^9, 3.9387001452733088`*^9, 
   3.938700248885579*^9, 3.93870033813731*^9, 3.938700370290875*^9, 
   3.938707638223526*^9, 3.938713922998406*^9, 3.9387140148682747`*^9, 
   3.938724269330748*^9, 3.938724396103724*^9, 3.938724488440359*^9, 
   3.938724546678544*^9, 3.9387246482973013`*^9, 3.938772452339512*^9, 
   3.938806593423192*^9, 3.938806695697021*^9, 3.938807631176777*^9, {
   3.938807770268972*^9, 3.938807797517008*^9}, 3.938807891147723*^9, 
   3.93880806479662*^9, 3.938808122366763*^9, 3.9388082813046923`*^9, 
   3.9388098303431997`*^9, {3.938809863725388*^9, 3.93880987656607*^9}, {
   3.938810078448085*^9, 3.938810104391127*^9}, 3.938810179289878*^9, 
   3.938811974316352*^9, 3.938814710090949*^9, 3.938816649996491*^9, 
   3.938817633049756*^9, 3.9388177819587717`*^9, 3.9388178577747927`*^9, 
   3.938818077210712*^9, 3.938818127214962*^9, 3.938865457490511*^9, 
   3.9397705471863747`*^9},
 CellLabel->
  "Out[838]=",ExpressionUUID->"d9dc851d-afd2-4556-a8a7-87655164e8df"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"7894.271299659034`", ",", 
   RowBox[{"-", "31998.56932617843`"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.9386993435075912`*^9, 3.938699496785405*^9, 3.938699680589446*^9, 
   3.938699768659811*^9, 3.9386998180243273`*^9, 3.938699864049407*^9, 
   3.9386999043412647`*^9, 3.938699987126308*^9, {3.938700022290429*^9, 
   3.938700037164534*^9}, 3.938700099555936*^9, 3.9387001452733088`*^9, 
   3.938700248885579*^9, 3.93870033813731*^9, 3.938700370290875*^9, 
   3.938707638223526*^9, 3.938713922998406*^9, 3.9387140148682747`*^9, 
   3.938724269330748*^9, 3.938724396103724*^9, 3.938724488440359*^9, 
   3.938724546678544*^9, 3.9387246482973013`*^9, 3.938772452339512*^9, 
   3.938806593423192*^9, 3.938806695697021*^9, 3.938807631176777*^9, {
   3.938807770268972*^9, 3.938807797517008*^9}, 3.938807891147723*^9, 
   3.93880806479662*^9, 3.938808122366763*^9, 3.9388082813046923`*^9, 
   3.9388098303431997`*^9, {3.938809863725388*^9, 3.93880987656607*^9}, {
   3.938810078448085*^9, 3.938810104391127*^9}, 3.938810179289878*^9, 
   3.938811974316352*^9, 3.938814710090949*^9, 3.938816649996491*^9, 
   3.938817633049756*^9, 3.9388177819587717`*^9, 3.9388178577747927`*^9, 
   3.938818077210712*^9, 3.938818127214962*^9, 3.938865457490511*^9, 
   3.939770550414667*^9},
 CellLabel->
  "Out[839]=",ExpressionUUID->"c8f65ab7-92b0-4fde-a273-7996c68e8cc7"]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.93868329772152*^9, 3.938683307630372*^9}},
 CellLabel->
  "In[840]:=",ExpressionUUID->"6427b661-08ba-4d9a-a811-91ed11cd1e49"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.938011021411188*^9, 3.938011031176703*^9}},
 CellLabel->
  "In[841]:=",ExpressionUUID->"4f4db268-f463-4a62-ae8f-862e4787f150"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{
  3.937148257956636*^9, {3.937955183255477*^9, 3.937955184558403*^9}, 
   3.938009568725326*^9, 3.938009746047697*^9},
 NumberMarks->False,
 CellLabel->
  "In[842]:=",ExpressionUUID->"39d767f6-6c65-4254-8d89-109a1544d8f7"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.937239054848371*^9, 3.9372390584635153`*^9}},
 CellLabel->
  "In[843]:=",ExpressionUUID->"eb8aeb14-6a69-4d67-bc73-2139e0fdc365"],

Cell[CellGroupData[{

Cell[BoxData["\[AliasDelimiter]"], "Input",
 CellChangeTimes->{3.938003504306855*^9},
 CellLabel->
  "In[844]:=",ExpressionUUID->"77176642-87f7-4d52-a6cf-187e77d311b5"],

Cell[BoxData["\[AliasDelimiter]"], "Output",
 CellChangeTimes->{
  3.9386993589742393`*^9, 3.938699511350532*^9, 3.938699695950616*^9, 
   3.938699784034501*^9, 3.938699833386594*^9, 3.9386998794544067`*^9, 
   3.938699919818776*^9, 3.9387000019526777`*^9, 3.938700052690093*^9, 
   3.9387001144839497`*^9, 3.938700160217535*^9, 3.938700268039221*^9, {
   3.938700357551981*^9, 3.938700384611143*^9}, 3.938707652433707*^9, 
   3.938713937115394*^9, 3.93871403234087*^9, 3.938724287764496*^9, 
   3.938724418972959*^9, 3.938724510851812*^9, 3.9387245637687607`*^9, 
   3.938724666641243*^9, 3.938806612000188*^9, 3.938806712541069*^9, 
   3.93880764828864*^9, {3.9388077875622263`*^9, 3.93880781472069*^9}, 
   3.938807908169424*^9, 3.93880808181244*^9, 3.938808139414295*^9, 
   3.938808298269663*^9, 3.9388098474228687`*^9, 3.9388098937608624`*^9, 
   3.938810121322361*^9, 3.938810196547358*^9, 3.9388119917676086`*^9, 
   3.9388147275788183`*^9, 3.93881666718377*^9, 3.938817804997229*^9, 
   3.938818105209723*^9, 3.938818154774643*^9, 3.9388654749272623`*^9, 
   3.939770550469075*^9},
 CellLabel->
  "Out[844]=",ExpressionUUID->"04259096-5f92-43b0-9b9b-f8b73a6ac26e"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"A", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "0", ",", " ", "0", ",", " ", "0", ",", " ", "1", ",", " ", "0", ",", 
       "0"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "1", ",", "0"}], "}"}], 
     ",", " ", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "1"}], "}"}], 
     ",", " ", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"3", 
        RowBox[{"k", "^", "2"}]}], ",", " ", "0", ",", "0", ",", "0", ",", 
       RowBox[{"2", "k"}], ",", "0"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", 
       RowBox[{
        RowBox[{"-", "2"}], "k"}], ",", "0", ",", "0"}], "}"}], ",", " ", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", 
       RowBox[{"-", 
        RowBox[{"k", "^", "2"}]}], ",", "0", ",", "0", ",", "0"}], "}"}]}], 
    "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"A", " ", "//", " ", "MatrixForm"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"A", " ", "//", " ", "Transpose"}], " ", "//", " ", "MatrixForm"}],
   ";"}], "\[IndentingNewLine]", 
 RowBox[{"s", " ", "=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"MatrixExp", "[", 
      RowBox[{
       RowBox[{"-", 
        RowBox[{"(", 
         RowBox[{"A", " ", "//", "Transpose"}], ")"}]}], " ", "t"}], "]"}], 
     " ", "/.", 
     RowBox[{"{", 
      RowBox[{"k", "->", "10"}], "}"}]}], "//", " ", "MatrixForm"}], "  ", "//",
    " ", "N"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"s", " ", "=", " ", 
   RowBox[{
    RowBox[{"MatrixExp", "[", 
     RowBox[{
      RowBox[{"-", 
       RowBox[{"(", 
        RowBox[{"A", " ", "//", "Transpose"}], ")"}]}], " ", "t"}], "]"}], 
    " ", "/.", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"k", "->", "10"}], ",", " ", 
      RowBox[{"t", "->", " ", "3600"}]}], "}"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"s", "[", 
    RowBox[{"[", 
     RowBox[{"4", " ", ";;", "5"}], "]"}], "]"}], " ", "//", " ", 
   "MatrixForm"}], " ", "//", "N"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Solve", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"s", "[", 
        RowBox[{"[", 
         RowBox[{"4", ";;", "5"}], "]"}], "]"}], ".", 
       RowBox[{"{", 
        RowBox[{
        "v1", ",", " ", "v2", ",", " ", "v3", ",", " ", "v4", ",", " ", "v5", 
         ",", " ", "v6"}], "}"}]}], " ", "==", 
      RowBox[{"{", 
       RowBox[{"0", ",", "0"}], "}"}]}], ",", " ", 
     RowBox[{"{", 
      RowBox[{
      "v1", ",", " ", "v2", ",", " ", "v3", ",", " ", "v4", ",", " ", "v5", 
       ",", " ", "v6"}], "}"}]}], "]"}], "[", 
   RowBox[{"[", "1", "]"}], "]"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.938772344034968*^9, 3.9387724978756*^9}, {
  3.9387725454968243`*^9, 3.9387725504652357`*^9}, {3.938812595327209*^9, 
  3.938812654901479*^9}, {3.938813536287328*^9, 3.9388136906806927`*^9}, {
  3.938813835639288*^9, 3.938813860622498*^9}, {3.938813904963852*^9, 
  3.9388139058151417`*^9}, {3.9388139609809523`*^9, 3.938813997586751*^9}, {
  3.938857620430049*^9, 3.938857632478833*^9}, {3.938857692941237*^9, 
  3.9388576929912786`*^9}, {3.938857735890628*^9, 3.9388577602407207`*^9}, {
  3.938857804583836*^9, 3.938857817429688*^9}, {3.9397580985592833`*^9, 
  3.939758098998438*^9}},
 CellLabel->
  "In[845]:=",ExpressionUUID->"18ba1bb9-5359-4203-9b17-149badf437d3"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", "0", "0", "1", "0", "0"},
     {"0", "0", "0", "0", "1", "0"},
     {"0", "0", "0", "0", "0", "1"},
     {
      RowBox[{"3", " ", 
       SuperscriptBox["k", "2"]}], "0", "0", "0", 
      RowBox[{"2", " ", "k"}], "0"},
     {"0", "0", "0", 
      RowBox[{
       RowBox[{"-", "2"}], " ", "k"}], "0", "0"},
     {"0", "0", 
      RowBox[{"-", 
       SuperscriptBox["k", "2"]}], "0", "0", "0"}
    },
    GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.938772441343375*^9, 3.9387724763229313`*^9}, 
   3.9387725099305162`*^9, 3.938772554329268*^9, 3.938806612014324*^9, 
   3.93880671255451*^9, 3.938807648307616*^9, {3.938807787581027*^9, 
   3.938807814739203*^9}, 3.9388079081837263`*^9, 3.9388080818256617`*^9, 
   3.938808139434301*^9, 3.9388082982902193`*^9, 3.938809847438273*^9, 
   3.938809893780755*^9, 3.938810121337162*^9, 3.938810196562296*^9, 
   3.9388119917704163`*^9, {3.938812606435754*^9, 3.938812655778474*^9}, {
   3.938813540301264*^9, 3.9388135482794847`*^9}, {3.938813579051676*^9, 
   3.93881365953008*^9}, 3.938813691947966*^9, {3.938813843113999*^9, 
   3.938813860965724*^9}, 3.938814002065447*^9, 3.9388147275977097`*^9, 
   3.938816667200032*^9, 3.938817805014983*^9, 3.938818105222166*^9, 
   3.9388181547875347`*^9, 3.938857634390337*^9, 3.93885769398458*^9, {
   3.938857736458275*^9, 3.93885776053476*^9}, {3.938857806458063*^9, 
   3.938857817680295*^9}, 3.938865474941268*^9, 3.939758101259748*^9, 
   3.939770550487207*^9},
 CellLabel->
  "Out[846]//MatrixForm=",ExpressionUUID->"b5820aeb-56f1-4930-b0f0-\
f9b6c8683e75"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{
       RowBox[{"-", "0.5`"}], " ", 
       SuperscriptBox["2.718281828459045`", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"0.`", "\[VeryThinSpace]", "-", 
           RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]], " ", 
       
       RowBox[{"(", 
        RowBox[{"3.`", "\[VeryThinSpace]", "-", 
         RowBox[{"8.`", " ", 
          SuperscriptBox["2.718281828459045`", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"0.`", "\[VeryThinSpace]", "+", 
              RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}], 
         "+", 
         RowBox[{"3.`", " ", 
          SuperscriptBox["2.718281828459045`", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"0.`", "\[VeryThinSpace]", "+", 
              RowBox[{"20.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
            "t"}]]}]}], ")"}]}], 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"0.`", "\[VeryThinSpace]", "+", 
          RowBox[{"3.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
        SuperscriptBox["2.718281828459045`", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"0.`", "\[VeryThinSpace]", "-", 
            RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]], " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1.`"}], "+", 
          SuperscriptBox["2.718281828459045`", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"0.`", "\[VeryThinSpace]", "+", 
              RowBox[{"20.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}], 
         ")"}]}], "+", 
       RowBox[{"60.`", " ", "t"}]}], "0.`", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"0.`", "\[VeryThinSpace]", "+", 
         RowBox[{"15.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
       SuperscriptBox["2.718281828459045`", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"0.`", "\[VeryThinSpace]", "-", 
           RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]], " ", 
       
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1.`"}], "+", 
         SuperscriptBox["2.718281828459045`", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"0.`", "\[VeryThinSpace]", "+", 
             RowBox[{"20.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}], 
        ")"}]}], 
      RowBox[{"30.`", " ", 
       SuperscriptBox["2.718281828459045`", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"0.`", "\[VeryThinSpace]", "-", 
           RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]], " ", 
       
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1.`"}], "+", 
          SuperscriptBox["2.718281828459045`", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"0.`", "\[VeryThinSpace]", "+", 
              RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}], 
         ")"}], "2"]}], "0.`"},
     {"0.`", "1.`", "0.`", "0.`", "0.`", "0.`"},
     {"0.`", "0.`", 
      RowBox[{
       RowBox[{"0.5`", " ", 
        SuperscriptBox["2.718281828459045`", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"0.`", "\[VeryThinSpace]", "-", 
            RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}], 
       "+", 
       RowBox[{"0.5`", " ", 
        SuperscriptBox["2.718281828459045`", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"0.`", "\[VeryThinSpace]", "+", 
            RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}]}], 
      "0.`", "0.`", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"0.`", "\[VeryThinSpace]", "+", 
          RowBox[{"5.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
        SuperscriptBox["2.718281828459045`", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"0.`", "\[VeryThinSpace]", "-", 
            RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}], 
       "-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"0.`", "\[VeryThinSpace]", "+", 
          RowBox[{"5.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
        SuperscriptBox["2.718281828459045`", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"0.`", "\[VeryThinSpace]", "+", 
            RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
          "t"}]]}]}]},
     {
      RowBox[{
       RowBox[{"(", 
        RowBox[{"0.`", "\[VeryThinSpace]", "+", 
         RowBox[{"0.05`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
       SuperscriptBox["2.718281828459045`", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"0.`", "\[VeryThinSpace]", "-", 
           RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]], " ", 
       
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1.`"}], "+", 
         SuperscriptBox["2.718281828459045`", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"0.`", "\[VeryThinSpace]", "+", 
             RowBox[{"20.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}], 
        ")"}]}], 
      RowBox[{"0.1`", " ", 
       SuperscriptBox["2.718281828459045`", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"0.`", "\[VeryThinSpace]", "-", 
           RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]], " ", 
       
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1.`"}], "+", 
          SuperscriptBox["2.718281828459045`", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"0.`", "\[VeryThinSpace]", "+", 
              RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}], 
         ")"}], "2"]}], "0.`", 
      RowBox[{"0.5`", " ", 
       SuperscriptBox["2.718281828459045`", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"0.`", "\[VeryThinSpace]", "-", 
           RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]], " ", 
       
       RowBox[{"(", 
        RowBox[{"1.`", "\[VeryThinSpace]", "+", 
         SuperscriptBox["2.718281828459045`", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"0.`", "\[VeryThinSpace]", "+", 
             RowBox[{"20.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}], 
        ")"}]}], 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"0.`", "\[VeryThinSpace]", "-", 
         RowBox[{"1.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
       SuperscriptBox["2.718281828459045`", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"0.`", "\[VeryThinSpace]", "-", 
           RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]], " ", 
       
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1.`"}], "+", 
         SuperscriptBox["2.718281828459045`", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"0.`", "\[VeryThinSpace]", "+", 
             RowBox[{"20.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}], 
        ")"}]}], "0.`"},
     {
      RowBox[{
       RowBox[{"-", "0.1`"}], " ", 
       SuperscriptBox["2.718281828459045`", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"0.`", "\[VeryThinSpace]", "-", 
           RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]], " ", 
       
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1.`"}], "+", 
          SuperscriptBox["2.718281828459045`", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"0.`", "\[VeryThinSpace]", "+", 
              RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}], 
         ")"}], "2"]}], 
      RowBox[{"0.1`", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{"0.`", "\[VeryThinSpace]", "+", 
            RowBox[{"2.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
          SuperscriptBox["2.718281828459045`", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"0.`", "\[VeryThinSpace]", "-", 
              RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]], 
          " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "1.`"}], "+", 
            SuperscriptBox["2.718281828459045`", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"0.`", "\[VeryThinSpace]", "+", 
                RowBox[{"20.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
              "t"}]]}], ")"}]}], "+", 
         RowBox[{"30.`", " ", "t"}]}], ")"}]}], "0.`", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"0.`", "\[VeryThinSpace]", "+", 
         RowBox[{"1.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
       SuperscriptBox["2.718281828459045`", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"0.`", "\[VeryThinSpace]", "-", 
           RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]], " ", 
       
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1.`"}], "+", 
         SuperscriptBox["2.718281828459045`", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"0.`", "\[VeryThinSpace]", "+", 
             RowBox[{"20.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}], 
        ")"}]}], 
      RowBox[{
       SuperscriptBox["2.718281828459045`", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"0.`", "\[VeryThinSpace]", "-", 
           RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]], " ", 
       
       RowBox[{"(", 
        RowBox[{"2.`", "\[VeryThinSpace]", "-", 
         RowBox[{"3.`", " ", 
          SuperscriptBox["2.718281828459045`", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"0.`", "\[VeryThinSpace]", "+", 
              RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}], 
         "+", 
         RowBox[{"2.`", " ", 
          SuperscriptBox["2.718281828459045`", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"0.`", "\[VeryThinSpace]", "+", 
              RowBox[{"20.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
            "t"}]]}]}], ")"}]}], "0.`"},
     {"0.`", "0.`", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"0.`", "\[VeryThinSpace]", "-", 
          RowBox[{"0.05`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
        SuperscriptBox["2.718281828459045`", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"0.`", "\[VeryThinSpace]", "-", 
            RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}], 
       "+", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"0.`", "\[VeryThinSpace]", "+", 
          RowBox[{"0.05`", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
        SuperscriptBox["2.718281828459045`", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"0.`", "\[VeryThinSpace]", "+", 
            RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}]}], 
      "0.`", "0.`", 
      RowBox[{
       RowBox[{"0.5`", " ", 
        SuperscriptBox["2.718281828459045`", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"0.`", "\[VeryThinSpace]", "-", 
            RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}], 
       "+", 
       RowBox[{"0.5`", " ", 
        SuperscriptBox["2.718281828459045`", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"0.`", "\[VeryThinSpace]", "+", 
            RowBox[{"10.`", " ", "\[ImaginaryI]"}]}], ")"}], " ", "t"}]]}]}]}
    },
    GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.938772441343375*^9, 3.9387724763229313`*^9}, 
   3.9387725099305162`*^9, 3.938772554329268*^9, 3.938806612014324*^9, 
   3.93880671255451*^9, 3.938807648307616*^9, {3.938807787581027*^9, 
   3.938807814739203*^9}, 3.9388079081837263`*^9, 3.9388080818256617`*^9, 
   3.938808139434301*^9, 3.9388082982902193`*^9, 3.938809847438273*^9, 
   3.938809893780755*^9, 3.938810121337162*^9, 3.938810196562296*^9, 
   3.9388119917704163`*^9, {3.938812606435754*^9, 3.938812655778474*^9}, {
   3.938813540301264*^9, 3.9388135482794847`*^9}, {3.938813579051676*^9, 
   3.93881365953008*^9}, 3.938813691947966*^9, {3.938813843113999*^9, 
   3.938813860965724*^9}, 3.938814002065447*^9, 3.9388147275977097`*^9, 
   3.938816667200032*^9, 3.938817805014983*^9, 3.938818105222166*^9, 
   3.9388181547875347`*^9, 3.938857634390337*^9, 3.93885769398458*^9, {
   3.938857736458275*^9, 3.93885776053476*^9}, {3.938857806458063*^9, 
   3.938857817680295*^9}, 3.938865474941268*^9, 3.939758101259748*^9, 
   3.939770550489135*^9},
 CellLabel->
  "Out[848]//MatrixForm=",ExpressionUUID->"0d25dbab-13db-4813-b4ef-\
04d2c4c516bd"],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"0.0470433982491772`", "\[VeryThinSpace]", "+", 
       RowBox[{"3.469446951953614`*^-18", " ", "\[ImaginaryI]"}]}], 
      RowBox[{
       RowBox[{"-", "0.3764870384041764`"}], "-", 
       RowBox[{"2.7755575615628914`*^-17", " ", "\[ImaginaryI]"}]}], "0.`", 
      RowBox[{
       RowBox[{"-", "0.8824351920208822`"}], "+", 
       RowBox[{"0.`", " ", "\[ImaginaryI]"}]}], 
      RowBox[{
       RowBox[{"-", "0.9408679649835441`"}], "-", 
       RowBox[{"5.551115123125783`*^-17", " ", "\[ImaginaryI]"}]}], "0.`"},
     {
      RowBox[{"0.3764870384041764`", "\[VeryThinSpace]", "+", 
       RowBox[{"2.7755575615628914`*^-17", " ", "\[ImaginaryI]"}]}], 
      RowBox[{"10800.188173592998`", "\[VeryThinSpace]", "+", 
       RowBox[{"1.1102230246251566`*^-17", " ", "\[ImaginaryI]"}]}], "0.`", 
      RowBox[{"0.9408679649835441`", "\[VeryThinSpace]", "+", 
       RowBox[{"5.551115123125783`*^-17", " ", "\[ImaginaryI]"}]}], 
      RowBox[{
       RowBox[{"-", "6.5297407680835295`"}], "+", 
       RowBox[{"4.440892098500626`*^-16", " ", "\[ImaginaryI]"}]}], "0.`"}
    },
    GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.938772441343375*^9, 3.9387724763229313`*^9}, 
   3.9387725099305162`*^9, 3.938772554329268*^9, 3.938806612014324*^9, 
   3.93880671255451*^9, 3.938807648307616*^9, {3.938807787581027*^9, 
   3.938807814739203*^9}, 3.9388079081837263`*^9, 3.9388080818256617`*^9, 
   3.938808139434301*^9, 3.9388082982902193`*^9, 3.938809847438273*^9, 
   3.938809893780755*^9, 3.938810121337162*^9, 3.938810196562296*^9, 
   3.9388119917704163`*^9, {3.938812606435754*^9, 3.938812655778474*^9}, {
   3.938813540301264*^9, 3.9388135482794847`*^9}, {3.938813579051676*^9, 
   3.93881365953008*^9}, 3.938813691947966*^9, {3.938813843113999*^9, 
   3.938813860965724*^9}, 3.938814002065447*^9, 3.9388147275977097`*^9, 
   3.938816667200032*^9, 3.938817805014983*^9, 3.938818105222166*^9, 
   3.9388181547875347`*^9, 3.938857634390337*^9, 3.93885769398458*^9, {
   3.938857736458275*^9, 3.93885776053476*^9}, {3.938857806458063*^9, 
   3.938857817680295*^9}, 3.938865474941268*^9, 3.939758101259748*^9, 
   3.9397705504916277`*^9},
 CellLabel->
  "Out[850]//MatrixForm=",ExpressionUUID->"159fe6a8-2fa8-4c3c-9ff1-\
e04830f20ba1"],

Cell[BoxData[
 TemplateBox[{
  "Solve", "svars", 
   "\"Equations may not give solutions for all \\\"solve\\\" variables.\"", 2,
    851, 752, 31053567148927707368, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{
  3.938814002122614*^9, 3.9388147276157417`*^9, 3.9388166672161694`*^9, 
   3.9388178050570583`*^9, 3.938818105243607*^9, 3.938818154804042*^9, 
   3.9388576344130363`*^9, 3.938857694000311*^9, {3.938857736478403*^9, 
   3.9388577605666027`*^9}, {3.938857806479974*^9, 3.938857817703627*^9}, 
   3.938865474949017*^9, 3.9397581013821297`*^9, 3.9397705504927073`*^9},
 CellLabel->
  "During evaluation of \
In[845]:=",ExpressionUUID->"bde679bc-2817-4f21-ac9a-b3170e3036a9"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"v4", "\[Rule]", 
    RowBox[{
     FractionBox[
      RowBox[{"\[ImaginaryI]", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "+", 
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{"72000", " ", "\[ImaginaryI]"}]]}], ")"}], " ", "v1"}], 
      RowBox[{"10", " ", 
       RowBox[{"(", 
        RowBox[{"3", "-", 
         RowBox[{"8", " ", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"36000", " ", "\[ImaginaryI]"}]]}], "+", 
         RowBox[{"3", " ", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"72000", " ", "\[ImaginaryI]"}]]}]}], ")"}]}]], "-", 
     FractionBox[
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"7", "+", 
           RowBox[{"108000", " ", "\[ImaginaryI]"}]}], ")"}], "-", 
         RowBox[{"14", " ", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"36000", " ", "\[ImaginaryI]"}]]}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"7", "-", 
            RowBox[{"108000", " ", "\[ImaginaryI]"}]}], ")"}], " ", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"72000", " ", "\[ImaginaryI]"}]]}]}], ")"}], " ", "v2"}], 
      RowBox[{"5", " ", 
       RowBox[{"(", 
        RowBox[{"3", "-", 
         RowBox[{"8", " ", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"36000", " ", "\[ImaginaryI]"}]]}], "+", 
         RowBox[{"3", " ", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"72000", " ", "\[ImaginaryI]"}]]}]}], ")"}]}]]}]}], ",", 
   RowBox[{"v5", "\[Rule]", 
    RowBox[{
     FractionBox[
      RowBox[{
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1"}], "+", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"36000", " ", "\[ImaginaryI]"}]]}], ")"}], "2"], " ", 
       "v1"}], 
      RowBox[{"5", " ", 
       RowBox[{"(", 
        RowBox[{"3", "-", 
         RowBox[{"8", " ", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"36000", " ", "\[ImaginaryI]"}]]}], "+", 
         RowBox[{"3", " ", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"72000", " ", "\[ImaginaryI]"}]]}]}], ")"}]}]], "+", 
     FractionBox[
      RowBox[{"2", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"27000", "-", "\[ImaginaryI]"}], ")"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"27000", "+", "\[ImaginaryI]"}], ")"}], " ", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"72000", " ", "\[ImaginaryI]"}]]}]}], ")"}], " ", "v2"}], 
      RowBox[{"5", " ", 
       RowBox[{"(", 
        RowBox[{"3", "-", 
         RowBox[{"8", " ", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"36000", " ", "\[ImaginaryI]"}]]}], "+", 
         RowBox[{"3", " ", 
          SuperscriptBox["\[ExponentialE]", 
           RowBox[{"72000", " ", "\[ImaginaryI]"}]]}]}], ")"}]}]]}]}]}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.938772441343375*^9, 3.9387724763229313`*^9}, 
   3.9387725099305162`*^9, 3.938772554329268*^9, 3.938806612014324*^9, 
   3.93880671255451*^9, 3.938807648307616*^9, {3.938807787581027*^9, 
   3.938807814739203*^9}, 3.9388079081837263`*^9, 3.9388080818256617`*^9, 
   3.938808139434301*^9, 3.9388082982902193`*^9, 3.938809847438273*^9, 
   3.938809893780755*^9, 3.938810121337162*^9, 3.938810196562296*^9, 
   3.9388119917704163`*^9, {3.938812606435754*^9, 3.938812655778474*^9}, {
   3.938813540301264*^9, 3.9388135482794847`*^9}, {3.938813579051676*^9, 
   3.93881365953008*^9}, 3.938813691947966*^9, {3.938813843113999*^9, 
   3.938813860965724*^9}, 3.938814002065447*^9, 3.9388147275977097`*^9, 
   3.938816667200032*^9, 3.938817805014983*^9, 3.938818105222166*^9, 
   3.9388181547875347`*^9, 3.938857634390337*^9, 3.93885769398458*^9, {
   3.938857736458275*^9, 3.93885776053476*^9}, {3.938857806458063*^9, 
   3.938857817680295*^9}, 3.938865474941268*^9, 3.939758101259748*^9, 
   3.939770550495954*^9},
 CellLabel->
  "Out[851]=",ExpressionUUID->"19a2d4db-f35e-497c-8f3e-793eaeac6dce"]
}, Open  ]],

Cell[BoxData["|"], "Input",
 CellChangeTimes->{3.938813625910653*^9},
 EmphasizeSyntaxErrors->True,
 CellLabel->
  "In[852]:=",ExpressionUUID->"0bb5789d-3499-4e74-b76f-ab38b17aa87a"]
},
WindowSize->{Full, Full},
WindowMargins->{{Automatic, 0}, {Automatic, 0}},
TaggingRules->Association["TryRealOnly" -> False],
FrontEndVersion->"12.3 for Mac OS X ARM (64-bit) (July 9, 2021)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"f90b9c7c-7060-47ee-b491-40eaabf8b6a7"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 175, 3, 30, "Input",ExpressionUUID->"46a113b6-085b-4448-9759-458b764a7e5d"],
Cell[736, 25, 200, 4, 30, "Input",ExpressionUUID->"21c1531e-bfb5-4ac1-a49d-5b0746559004"],
Cell[CellGroupData[{
Cell[961, 33, 119130, 2625, 5066, "Input",ExpressionUUID->"c387f004-7059-4de4-92e0-ddeb48546c32"],
Cell[120094, 2660, 1322, 20, 34, "Output",ExpressionUUID->"c78c004f-444d-4c3c-bed2-bc3fc65dccbd"],
Cell[CellGroupData[{
Cell[121441, 2684, 1319, 21, 24, "Print",ExpressionUUID->"5f9e968b-475f-4f11-a53d-45b1dfb70a03"],
Cell[122763, 2707, 1333, 21, 24, "Print",ExpressionUUID->"244b1158-e79b-4be1-92a5-de83b1ac5953"],
Cell[124099, 2730, 1346, 21, 24, "Print",ExpressionUUID->"082f52d2-1de5-44fa-838a-e5d18fd52213"],
Cell[125448, 2753, 1581, 28, 24, "Print",ExpressionUUID->"9359be76-0a4b-40d2-8b6a-99517fa60fc3"],
Cell[127032, 2783, 1369, 22, 24, "Print",ExpressionUUID->"ba843332-dd83-41ca-8de1-89a37b08414f"]
}, Open  ]],
Cell[128416, 2808, 1520, 25, 34, "Output",ExpressionUUID->"6acaa1b0-ba4b-4c8e-8805-5359733492af"],
Cell[129939, 2835, 1404, 24, 34, "Output",ExpressionUUID->"0c7184e8-70d2-4f47-b93e-2a0285cf4161"],
Cell[131346, 2861, 1669, 27, 34, "Output",ExpressionUUID->"de398d5d-f758-4747-ad5d-ab409793fc30"],
Cell[133018, 2890, 1406, 24, 34, "Output",ExpressionUUID->"c4ca6d6d-984f-467e-8bee-ba69b59e69a4"],
Cell[134427, 2916, 1730, 29, 34, "Output",ExpressionUUID->"fd40a038-ce73-4160-8fce-f06c91c19535"],
Cell[CellGroupData[{
Cell[136182, 2949, 1325, 21, 24, "Print",ExpressionUUID->"8070c84e-345e-48d3-8ad3-6fdfd47b706a"],
Cell[137510, 2972, 1400, 23, 24, "Print",ExpressionUUID->"9cb42d38-ffa1-4020-9edd-2446261d20c5"]
}, Open  ]],
Cell[138925, 2998, 1083, 20, 28, "Message",ExpressionUUID->"e239a939-f47a-464e-bb9c-dd59710e2979"],
Cell[140011, 3020, 1634, 29, 34, "Output",ExpressionUUID->"0d8a9c4a-62f2-443d-8cec-7fb76c722e24"],
Cell[141648, 3051, 1368, 22, 24, "Print",ExpressionUUID->"08417164-06f0-4485-8e4b-4b3ad0ca7bee"],
Cell[143019, 3075, 1730, 29, 34, "Output",ExpressionUUID->"858755d7-6d44-4408-86c4-99a10d441055"],
Cell[144752, 3106, 1379, 22, 24, "Print",ExpressionUUID->"7b821ad8-ac4a-43d2-89a1-c9c66111c615"],
Cell[146134, 3130, 1582, 26, 37, "Output",ExpressionUUID->"3d2677be-5a9d-47e4-94bf-e29a662b8695"],
Cell[147719, 3158, 1379, 22, 24, "Print",ExpressionUUID->"507d2b23-3b46-40b8-925d-d6e4f8cddfab"],
Cell[149101, 3182, 1061, 19, 31, "Message",ExpressionUUID->"c0ab2dee-e1fb-4a45-a823-22975c223d10"],
Cell[150165, 3203, 1732, 29, 37, "Output",ExpressionUUID->"dce4d94b-eb41-484a-9a62-258f53067104"],
Cell[151900, 3234, 1407, 24, 34, "Output",ExpressionUUID->"6497c571-4391-4ba2-a2a9-e7498571bcb5"],
Cell[153310, 3260, 1350, 21, 24, "Print",ExpressionUUID->"c19d1af2-38cf-4a8e-ad8b-690b86a9db53"],
Cell[154663, 3283, 1583, 27, 37, "Output",ExpressionUUID->"d9dc851d-afd2-4556-a8a7-87655164e8df"],
Cell[156249, 3312, 1399, 23, 34, "Output",ExpressionUUID->"c8f65ab7-92b0-4fde-a273-7996c68e8cc7"]
}, Open  ]],
Cell[157663, 3338, 174, 3, 30, "Input",ExpressionUUID->"6427b661-08ba-4d9a-a811-91ed11cd1e49"],
Cell[157840, 3343, 175, 3, 30, "Input",ExpressionUUID->"4f4db268-f463-4a62-ae8f-862e4787f150"],
Cell[158018, 3348, 269, 6, 30, "Input",ExpressionUUID->"39d767f6-6c65-4254-8d89-109a1544d8f7"],
Cell[158290, 3356, 177, 3, 30, "Input",ExpressionUUID->"eb8aeb14-6a69-4d67-bc73-2139e0fdc365"],
Cell[CellGroupData[{
Cell[158492, 3363, 168, 3, 30, "Input",ExpressionUUID->"77176642-87f7-4d52-a6cf-187e77d311b5"],
Cell[158663, 3368, 1173, 18, 34, "Output",ExpressionUUID->"04259096-5f92-43b0-9b9b-f8b73a6ac26e"]
}, Open  ]],
Cell[CellGroupData[{
Cell[159873, 3391, 3551, 98, 199, "Input",ExpressionUUID->"18ba1bb9-5359-4203-9b17-149badf437d3"],
Cell[163427, 3491, 1981, 44, 134, "Output",ExpressionUUID->"b5820aeb-56f1-4930-b0f0-f9b6c8683e75"],
Cell[165411, 3537, 12887, 341, 152, "Output",ExpressionUUID->"0d25dbab-13db-4813-b4ef-04d2c4c516bd"],
Cell[178301, 3880, 2670, 53, 66, "Output",ExpressionUUID->"159fe6a8-2fa8-4c3c-9ff1-e04830f20ba1"],
Cell[180974, 3935, 708, 14, 28, "Message",ExpressionUUID->"bde679bc-2817-4f21-ac9a-b3170e3036a9"],
Cell[181685, 3951, 4176, 103, 62, "Output",ExpressionUUID->"19a2d4db-f35e-497c-8f3e-793eaeac6dce"]
}, Open  ]],
Cell[185876, 4057, 182, 4, 30, "Input",ExpressionUUID->"0bb5789d-3499-4e74-b76f-ab38b17aa87a"]
}
]
*)