energythrustcomparisondebugging.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[    115983,       2611]
NotebookOptionsPosition[    111612,       2545]
NotebookOutlinePosition[    112024,       2562]
CellTagsIndexPosition[    111981,       2559]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "*", " ", "This", " ", "is", " ", "code", " ", "modified", " ", 
    "specifically", " ", "for", " ", "the", " ", "2", "d", " ", "case"}], " ",
    "**)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
   "*", " ", "Set", " ", "the", " ", "directory", " ", "that", " ", "you", 
    " ", "want", " ", "your", " ", "plots", " ", "to", " ", "go", " ", "to", 
    " ", "here"}], " ", "**)"}], "\[IndentingNewLine]", 
  RowBox[{
  "SetDirectory", "[", "\"\</users/maxzweig/energythrust\>\"", 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.9257345618408537`*^9, 3.9257345767008457`*^9}, {
  3.926591746591373*^9, 3.926591767031913*^9}},
 CellLabel->
  "In[440]:=",ExpressionUUID->"344c9343-693b-47e9-b882-096618076cd6"],

Cell[BoxData["\<\"/Users/maxzweig/energythrust\"\>"], "Output",
 CellChangeTimes->{3.926595031261936*^9, 3.926595151448574*^9, 
  3.926595447385315*^9, 3.926684006417115*^9},
 CellLabel->
  "Out[440]=",ExpressionUUID->"5dc626c2-b9aa-4372-87c6-912baed454f4"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "*", " ", "Various", " ", "helper", " ", "functions", " ", "for", " ", 
    "standard", " ", "parts", " ", "of", " ", "the", " ", 
    RowBox[{"paper", " ", "/", " ", "reachable"}], " ", "set", " ", 
    "results"}], " ", "**)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"CW", "[", "n_", "]"}], " ", ":=", " ", 
     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[{
     RowBox[{"xrf", "[", 
      RowBox[{
      "tf_", ",", " ", "delta_", ",", " ", "h_", ",", " ", "deltaT_", ",", 
       " ", "n_", ",", " ", "Rr_", ",", " ", "emax_"}], "]"}], " ", ":=", " ", 
     RowBox[{"NIntegrate", "[", "  ", 
      RowBox[{
       RowBox[{"emax", " ", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"-", 
                 FractionBox[
                  RowBox[{"Sin", "[", 
                   RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}], "n"]}], ",", 
                FractionBox[
                 RowBox[{"2", " ", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"-", "1"}], "+", 
                    RowBox[{"Cos", "[", 
                    RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}]}], ")"}]}], 
                 "n"], ",", "0"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"-", 
                 FractionBox[
                  RowBox[{"2", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "1"}], "+", 
                    RowBox[{"Cos", "[", 
                    RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}]}], ")"}]}], 
                  "n"]}], ",", 
                RowBox[{
                 RowBox[{"3", " ", "tau"}], "-", 
                 RowBox[{"3", " ", "tf"}], "-", 
                 FractionBox[
                  RowBox[{"4", " ", 
                   RowBox[{"Sin", "[", 
                    RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}]}], "n"]}], ",",
                 "0"}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{"0", ",", "0", ",", 
                RowBox[{"-", 
                 FractionBox[
                  RowBox[{"Sin", "[", 
                   RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}], "n"]}]}], 
               "}"}]}], "}"}], " ", "//", " ", "Transpose"}], ")"}], ".", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"-", 
               FractionBox[
                RowBox[{"Sin", "[", 
                 RowBox[{"n", " ", 
                  RowBox[{"(", 
                   RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}], "n"]}], ",", 
              FractionBox[
               RowBox[{"2", " ", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"-", "1"}], "+", 
                  RowBox[{"Cos", "[", 
                   RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}]}], ")"}]}], 
               "n"], ",", "0"}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"-", 
               FractionBox[
                RowBox[{"2", " ", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"-", "1"}], "+", 
                   RowBox[{"Cos", "[", 
                    RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}]}], ")"}]}], 
                "n"]}], ",", 
              RowBox[{
               RowBox[{"3", " ", "tau"}], "-", 
               RowBox[{"3", " ", "tf"}], "-", 
               FractionBox[
                RowBox[{"4", " ", 
                 RowBox[{"Sin", "[", 
                  RowBox[{"n", " ", 
                   RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}]}], "n"]}], ",",
               "0"}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{"0", ",", "0", ",", 
              RowBox[{"-", 
               FractionBox[
                RowBox[{"Sin", "[", 
                 RowBox[{"n", " ", 
                  RowBox[{"(", 
                   RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}], "n"]}]}], 
             "}"}]}], "}"}], ".", "delta"}], " ", "/", " ", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"deltaT", ".", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"{", 
                RowBox[{
                 RowBox[{"{", 
                  RowBox[{
                   RowBox[{"-", 
                    FractionBox[
                    RowBox[{"Sin", "[", 
                    RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}], "n"]}], ",", 
                   FractionBox[
                    RowBox[{"2", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "1"}], "+", 
                    RowBox[{"Cos", "[", 
                    RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}]}], ")"}]}], 
                    "n"], ",", "0"}], "}"}], ",", 
                 RowBox[{"{", 
                  RowBox[{
                   RowBox[{"-", 
                    FractionBox[
                    RowBox[{"2", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "1"}], "+", 
                    RowBox[{"Cos", "[", 
                    RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}]}], ")"}]}], 
                    "n"]}], ",", 
                   RowBox[{
                    RowBox[{"3", " ", "tau"}], "-", 
                    RowBox[{"3", " ", "tf"}], "-", 
                    FractionBox[
                    RowBox[{"4", " ", 
                    RowBox[{"Sin", "[", 
                    RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}]}], "n"]}], ",",
                    "0"}], "}"}], ",", 
                 RowBox[{"{", 
                  RowBox[{"0", ",", "0", ",", 
                   RowBox[{"-", 
                    FractionBox[
                    RowBox[{"Sin", "[", 
                    RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}], "n"]}]}], 
                  "}"}]}], "}"}], " ", "//", " ", "Transpose"}], ")"}], ".", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"{", 
                RowBox[{
                 RowBox[{"-", 
                  FractionBox[
                   RowBox[{"Sin", "[", 
                    RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}], "n"]}], ",", 
                 FractionBox[
                  RowBox[{"2", " ", 
                   RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "1"}], "+", 
                    RowBox[{"Cos", "[", 
                    RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}]}], ")"}]}], 
                  "n"], ",", "0"}], "}"}], ",", 
               RowBox[{"{", 
                RowBox[{
                 RowBox[{"-", 
                  FractionBox[
                   RowBox[{"2", " ", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"-", "1"}], "+", 
                    RowBox[{"Cos", "[", 
                    RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}]}], ")"}]}], 
                   "n"]}], ",", 
                 RowBox[{
                  RowBox[{"3", " ", "tau"}], "-", 
                  RowBox[{"3", " ", "tf"}], "-", 
                  FractionBox[
                   RowBox[{"4", " ", 
                    RowBox[{"Sin", "[", 
                    RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}]}], "n"]}], ",",
                  "0"}], "}"}], ",", 
               RowBox[{"{", 
                RowBox[{"0", ",", "0", ",", 
                 RowBox[{"-", 
                  FractionBox[
                   RowBox[{"Sin", "[", 
                    RowBox[{"n", " ", 
                    RowBox[{"(", 
                    RowBox[{"tau", "-", "tf"}], ")"}]}], "]"}], "n"]}]}], 
                "}"}]}], "}"}], ".", "delta"}], " ", "//", " ", "Sqrt"}], 
           ")"}], "[", 
          RowBox[{"[", "1", "]"}], "]"}]}]}], "  ", ",", 
       RowBox[{"{", 
        RowBox[{"tau", ",", " ", "0", ",", " ", "tf"}], "}"}]}], "]"}]}], 
    ";"}]}]}]], "Input",
 CellChangeTimes->{{3.923729826947426*^9, 3.923729830569713*^9}, {
  3.9265877766574917`*^9, 3.926587780260474*^9}, {3.926587830785193*^9, 
  3.926587849312529*^9}, {3.926588296592351*^9, 3.9265883476794243`*^9}, {
  3.926588415244143*^9, 3.926588478228994*^9}, {3.926588616827889*^9, 
  3.92658861953994*^9}, {3.926588711400141*^9, 3.926588714934411*^9}, {
  3.926588855759663*^9, 3.926588858769762*^9}, {3.926588913526075*^9, 
  3.926588931437234*^9}, {3.926588970299714*^9, 3.9265890380452433`*^9}, {
  3.9265906008374147`*^9, 3.926590601099334*^9}, {3.926590655925519*^9, 
  3.926590667696477*^9}, {3.926591550358121*^9, 3.926591554915353*^9}, {
  3.926591590248148*^9, 3.92659165468366*^9}, {3.926592138336666*^9, 
  3.9265921542392*^9}},
 CellLabel->
  "In[441]:=",ExpressionUUID->"b895b316-708f-4092-a477-0f3a66f14cfd"],

Cell[BoxData[
 TemplateBox[{
  "SetDelayed", "write", 
   "\"Tag \\!\\(\\*RowBox[{\\\"List\\\"}]\\) in \
\\!\\(\\*RowBox[{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\\\", \\\" \\\", SuperscriptBox[\\\"x\\\", \\\"2\\\"]}], \
\\\",\\\", \\\"0\\\", \\\",\\\", \\\"0\\\", \\\",\\\", \\\"0\\\", \\\",\\\", \
RowBox[{\\\"2\\\", \\\" \\\", \\\"x\\\"}], \\\",\\\", \\\"0\\\"}], \
\\\"}\\\"}], \\\",\\\", RowBox[{\\\"{\\\", RowBox[{\\\"0\\\", \\\",\\\", \
\\\"0\\\", \\\",\\\", \\\"0\\\", \\\",\\\", RowBox[{RowBox[{\\\"-\\\", \
\\\"2\\\"}], \\\" \\\", \\\"x\\\"}], \\\",\\\", \\\"0\\\", \\\",\\\", \\\"0\\\
\"}], \\\"}\\\"}], \\\",\\\", RowBox[{\\\"{\\\", RowBox[{\\\"0\\\", \
\\\",\\\", \\\"0\\\", \\\",\\\", RowBox[{\\\"-\\\", SuperscriptBox[\\\"x\\\", \
\\\"2\\\"]}], \\\",\\\", \\\"0\\\", \\\",\\\", \\\"0\\\", \\\",\\\", \
\\\"0\\\"}], \\\"}\\\"}]}], \\\"}\\\"}], \\\"[\\\", \\\"n_\\\", \
\\\"]\\\"}]\\) is Protected.\"", 2, 441, 76, 30967495211780212173, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.9266840064805193`*^9},
 CellLabel->
  "During evaluation of \
In[441]:=",ExpressionUUID->"a6b82077-2ddc-44eb-9426-d2e8355a1a8a"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"*", " ", "Computes", " ", "the", " ", "thrust"}], "-", 
    RowBox[{
    "limited", " ", "reachable", " ", "set", " ", "for", " ", "a", " ", 
     "given", " ", "final", " ", "time"}]}], " ", "**)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"thrustreachable", "[", 
    RowBox[{"xrf_", ",", "tf_", ",", " ", "na_", ",", " ", "emax_"}], "]"}], 
   " ", ":=", " ", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "points", ",", " ", "tfs", ",", " ", "hs", ",", " ", "As", ",", " ", 
       "ns", ",", " ", "Rrs", ",", " ", "pointsT", ",", " ", "rd", ",", " ", 
       "rdpoints", ",", " ", "rdboundary", ",", " ", "first2"}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"points", " ", "=", " ", 
       RowBox[{"CirclePoints", "[", "1000", "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"points", " ", "=", " ", 
       RowBox[{
        RowBox[{
         RowBox[{"Append", "[", 
          RowBox[{"#", ",", " ", "0"}], "]"}], " ", "&"}], " ", "/@", " ", 
        "points"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"points", " ", "//", "  ", "N"}], ";", "\[IndentingNewLine]", 
      RowBox[{"points", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"tfs", " ", "=", " ", 
       RowBox[{"ConstantArray", "[", 
        RowBox[{"tf", ",", " ", "1000"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"tfs", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"hs", " ", "=", " ", 
       RowBox[{"ConstantArray", "[", 
        RowBox[{"1", ",", " ", "1000"}], "]"}]}], ";", "\[IndentingNewLine]", 
      
      RowBox[{"hs", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"ns", " ", "=", " ", 
       RowBox[{"ConstantArray", "[", 
        RowBox[{"na", ",", " ", "1000"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"ns", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"Rrs", " ", "=", " ", 
       RowBox[{"ConstantArray", "[", 
        RowBox[{"1", ",", " ", "1000"}], "]"}]}], ";", "\[IndentingNewLine]", 
      
      RowBox[{"Rrs", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"emaxs", " ", "=", " ", 
       RowBox[{"ConstantArray", "[", 
        RowBox[{"emax", ",", " ", "1000"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"columntorow", "[", "x_", "]"}], " ", ":=", " ", 
       RowBox[{"{", "x", "}"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"pointsT", " ", "=", " ", 
       RowBox[{"columntorow", " ", "/@", " ", "points"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"pointsT", "[", 
        RowBox[{"[", "1", "]"}], "]"}], " ", "//", " ", "MatrixForm"}], ";", 
      "\[IndentingNewLine]", "\[IndentingNewLine]", 
      RowBox[{"rd", " ", "=", " ", 
       RowBox[{"MapThread", "[", 
        RowBox[{"xrf", ",", " ", 
         RowBox[{"{", 
          RowBox[{
          "tfs", ",", " ", "points", ",", " ", "hs", ",", " ", "pointsT", ",",
            " ", "ns", ",", " ", "Rrs", ",", " ", "emaxs"}], "}"}]}], "]"}]}],
       ";", "\[IndentingNewLine]", 
      RowBox[{"rdpoints", " ", "=", " ", 
       RowBox[{"Map", "[", 
        RowBox[{"Re", ",", " ", "rd", ",", " ", 
         RowBox[{"{", "2", "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"(*", " ", 
       RowBox[{
        RowBox[{
         RowBox[{"first2", "[", "x_", "]"}], " ", ":=", " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"x", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", " ", 
           RowBox[{"x", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], ";"}], " ", "*)"}], 
      "\[IndentingNewLine]", 
      RowBox[{"(*", " ", 
       RowBox[{
        RowBox[{"rboundary", " ", "=", " ", 
         RowBox[{"Map", "[", 
          RowBox[{"first2", ",", " ", "rdboundary"}], "]"}]}], ";"}], "*)"}], 
      "\[IndentingNewLine]", "rdpoints"}]}], "\[IndentingNewLine]", 
    "]"}]}]}]], "Input",
 CellChangeTimes->{{3.922125709580538*^9, 3.9221257398332644`*^9}, {
   3.922130104240691*^9, 3.9221301062117233`*^9}, {3.922130190198127*^9, 
   3.922130197137838*^9}, {3.922805236499282*^9, 3.922805236655064*^9}, {
   3.923729856169129*^9, 3.923729888079521*^9}, 3.925732706211501*^9, {
   3.92658918522209*^9, 3.92658918648283*^9}, 3.9265893379964743`*^9, {
   3.9265895024492807`*^9, 3.926589505424794*^9}, {3.926589580250574*^9, 
   3.926589580610242*^9}, {3.926589669672036*^9, 3.9265896761799917`*^9}, {
   3.926589766435689*^9, 3.926589781174829*^9}, {3.926589887658388*^9, 
   3.926589919496035*^9}, {3.926591669139657*^9, 3.926591688870654*^9}, {
   3.926592182618226*^9, 3.926592194707149*^9}},
 CellLabel->
  "In[443]:=",ExpressionUUID->"38e6d68a-01ae-469f-8372-72472f6edb77"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "*", " ", "Computes", " ", "the", " ", "energy", " ", "reachable", " ", 
    "quadratic", " ", "form", " ", "for", " ", "an", " ", "associated", " ", 
    "n", " ", "and", " ", "final", " ", "time"}], " ", "**)"}], 
  "\[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_", ",", " ", "n_"}], "]"}], " ", "=", " ", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{"72", "  ", "*", " ", "t", "  ", "*", " ", 
              RowBox[{"Cos", "[", 
               RowBox[{"n", " ", "*", " ", "t"}], "]"}]}], " ", "+", " ", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"154", " ", "*", " ", "n", " ", "*", " ", "t"}], " ", 
                "+", " ", 
                RowBox[{"48", " ", "*", " ", 
                 RowBox[{"n", " ", "^", " ", "3"}], " ", "*", " ", 
                 RowBox[{"t", " ", "^", " ", "3"}]}], " ", "-", " ", 
                RowBox[{"384", " ", "*", " ", 
                 RowBox[{"Sin", "[", 
                  RowBox[{"n", " ", "*", " ", "t"}], "]"}]}], " ", "-", " ", 
                RowBox[{"27", " ", "*", " ", 
                 RowBox[{"Sin", "[", 
                  RowBox[{"2", " ", "*", " ", "n", " ", "*", " ", "t"}], 
                  "]"}]}]}], ")"}], " ", "/", " ", 
              RowBox[{"(", 
               RowBox[{"4", " ", "*", " ", "n"}], ")"}]}]}], ",", " ", 
            RowBox[{
             RowBox[{
              RowBox[{"-", "3"}], " ", "*", " ", "n", " ", "*", " ", 
              RowBox[{"t", " ", "^", " ", "2"}]}], " ", "+", " ", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"6", " ", "-", " ", 
                RowBox[{"6", " ", "*", " ", 
                 RowBox[{"Cos", "[", 
                  RowBox[{"n", " ", "*", " ", "t"}], "]"}]}]}], " ", ")"}], 
              " ", "/", " ", "n"}]}], ",", " ", "0", ",", " ", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"8", "+", 
               RowBox[{"12", " ", "*", " ", 
                RowBox[{"n", "^", "2"}], "*", " ", 
                RowBox[{"t", "^", "2"}]}], "\[Minus]", 
               RowBox[{"8", "*", " ", 
                RowBox[{"Cos", "[", 
                 RowBox[{"n", " ", "*", " ", "t"}], "]"}]}], "\[Minus]", 
               RowBox[{"24", "*", "n", "*", " ", "t", "*", " ", 
                RowBox[{"Sin", "[", 
                 RowBox[{"n", " ", "*", " ", "t"}], "]"}]}], "+", 
               RowBox[{"9", " ", "*", 
                RowBox[{
                 RowBox[{"Sin", "[", " ", 
                  RowBox[{"n", " ", "*", " ", "t"}], "]"}], "^", "2"}]}]}], 
              " ", ")"}], "/", " ", 
             RowBox[{"(", 
              RowBox[{"2", " ", "*", " ", 
               RowBox[{"n", " ", "^", " ", "2"}]}], ")"}]}], ",", " ", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"23", " ", "*", "n", " ", "*", " ", "t"}], "+", 
               RowBox[{"6", " ", "*", 
                RowBox[{"n", "^", "3"}], " ", "*", 
                RowBox[{"t", "^", "3"}]}], "+", 
               RowBox[{"42", "*", "n", "*", "t", "*", " ", 
                RowBox[{"Cos", "[", " ", 
                 RowBox[{"n", " ", "*", " ", "t"}], "]"}]}], "\[Minus]", 
               RowBox[{"56", " ", "*", " ", 
                RowBox[{"Sin", "[", " ", 
                 RowBox[{"n", " ", "*", " ", "t"}], "]"}]}], "\[Minus]", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"9", "/", "2"}], ")"}], " ", "*", " ", 
                RowBox[{"Sin", " ", "[", 
                 RowBox[{"2", "*", "n", "*", "t"}], "]"}]}]}], ")"}], "/", 
             RowBox[{"n", "^", "2"}]}], ",", " ", "0"}], 
           "\[IndentingNewLine]", "}"}], ",", "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{
              RowBox[{"-", "3"}], " ", "*", " ", "n", " ", "*", " ", 
              RowBox[{"t", " ", "^", " ", "2"}]}], " ", "+", " ", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"6", " ", "-", " ", 
                RowBox[{"6", " ", "*", " ", 
                 RowBox[{"Cos", "[", 
                  RowBox[{"n", " ", "*", " ", "t"}], "]"}]}]}], " ", ")"}], 
              " ", "/", " ", "n"}]}], ",", " ", "t", ",", " ", "0", ",", "  ", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{
                RowBox[{"-", "2"}], "*", "n", " ", "*", " ", "t"}], " ", "+", 
               " ", 
               RowBox[{"2", " ", "*", " ", 
                RowBox[{"Sin", "[", 
                 RowBox[{"n", " ", "*", " ", "t"}], "]"}]}]}], ")"}], " ", 
             "/", " ", 
             RowBox[{"n", "^", "2"}]}], ",", " ", 
            RowBox[{
             RowBox[{
              RowBox[{"-", " ", "3"}], " ", "*", " ", 
              RowBox[{
               RowBox[{"t", "^", "2"}], " ", "/", " ", "2"}]}], " ", "+", " ", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"4", " ", "-", " ", 
                RowBox[{"4", " ", 
                 RowBox[{"Cos", "[", 
                  RowBox[{"n", " ", "*", " ", "t"}], "]"}]}]}], ")"}], "/", 
              RowBox[{"n", " ", "^", " ", "2"}]}]}], ",", " ", "0"}], " ", 
           "\[IndentingNewLine]", "}"}], ",", "\[IndentingNewLine]", 
          RowBox[{"{", "\[IndentingNewLine]", 
           RowBox[{"0", ",", " ", "0", ",", " ", 
            RowBox[{
             RowBox[{"(", " ", 
              RowBox[{
               RowBox[{"2", " ", "*", " ", "n", " ", "*", " ", "t"}], " ", 
               "+", " ", 
               RowBox[{"Sin", "[", 
                RowBox[{"2", " ", "*", " ", "n", " ", "*", " ", "t"}], 
                "]"}]}], ")"}], " ", "/", " ", 
             RowBox[{"(", 
              RowBox[{"4", " ", "*", "n"}], ")"}]}], " ", ",", " ", "0", ",", 
            " ", "0", ",", " ", 
            RowBox[{
             RowBox[{
              RowBox[{"Sin", "[", 
               RowBox[{"n", " ", "*", " ", "t"}], "]"}], " ", "^", " ", "2"}],
              " ", "/", " ", 
             RowBox[{"(", 
              RowBox[{"2", " ", "*", " ", 
               RowBox[{"n", "^", "2"}]}], ")"}]}]}], "\[IndentingNewLine]", 
           "}"}], ",", "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{"8", "+", 
               RowBox[{"12", " ", "*", " ", 
                RowBox[{"n", "^", "2"}], "*", " ", 
                RowBox[{"t", "^", "2"}]}], "\[Minus]", 
               RowBox[{"8", "*", " ", 
                RowBox[{"Cos", "[", 
                 RowBox[{"n", " ", "*", " ", "t"}], "]"}]}], "\[Minus]", 
               RowBox[{"24", "*", "n", "*", " ", "t", "*", " ", 
                RowBox[{"Sin", "[", 
                 RowBox[{"n", " ", "*", " ", "t"}], "]"}]}], "+", 
               RowBox[{"9", " ", "*", 
                RowBox[{
                 RowBox[{"Sin", "[", " ", 
                  RowBox[{"n", " ", "*", " ", "t"}], "]"}], "^", "2"}]}]}], 
              " ", ")"}], "/", " ", 
             RowBox[{"(", 
              RowBox[{"2", " ", "*", " ", 
               RowBox[{"n", " ", "^", " ", "2"}]}], ")"}]}], ",", " ", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{
                RowBox[{"-", "2"}], "*", "n", " ", "*", " ", "t"}], " ", "+", 
               " ", 
               RowBox[{"2", " ", "*", " ", 
                RowBox[{"Sin", "[", 
                 RowBox[{"n", " ", "*", " ", "t"}], "]"}]}]}], ")"}], " ", 
             "/", " ", 
             RowBox[{"n", "^", "2"}]}], ",", " ", "0", ",", " ", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"26", " ", "*", " ", "n", " ", "*", " ", "t"}], " ", 
               "-", " ", 
               RowBox[{"32", " ", "*", " ", 
                RowBox[{"Sin", "[", 
                 RowBox[{"n", " ", "*", " ", "t"}], "]"}]}], " ", "+", " ", 
               RowBox[{"3", " ", "*", " ", 
                RowBox[{"Sin", "[", 
                 RowBox[{"2", " ", "*", " ", "n", " ", "*", " ", "t"}], 
                 "]"}]}]}], ")"}], " ", "/", " ", 
             RowBox[{"(", 
              RowBox[{"4", " ", "*", " ", 
               RowBox[{"n", "^", "3"}]}], ")"}]}], ",", " ", 
            RowBox[{"3", " ", "*", " ", 
             RowBox[{
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", " ", "n"}], " ", "*", " ", "t"}], " ", "+", 
                 " ", 
                 RowBox[{"Sin", "[", 
                  RowBox[{"n", " ", "*", " ", "t"}], "]"}]}], ")"}], "^", 
               "2"}], "/", " ", 
              RowBox[{"n", " ", "^", " ", "3"}]}]}], ",", " ", "0"}], 
           "\[IndentingNewLine]", "}"}], ",", "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"23", " ", "*", "n", " ", "*", " ", "t"}], "+", 
               RowBox[{"6", " ", "*", 
                RowBox[{"n", "^", "3"}], " ", "*", 
                RowBox[{"t", "^", "3"}]}], "+", 
               RowBox[{"42", "*", "n", "*", "t", "*", " ", 
                RowBox[{"Cos", "[", " ", 
                 RowBox[{"n", " ", "*", " ", "t"}], "]"}]}], "\[Minus]", 
               RowBox[{"56", " ", "*", " ", 
                RowBox[{"Sin", "[", " ", 
                 RowBox[{"n", " ", "*", " ", "t"}], "]"}]}], "\[Minus]", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"9", "/", "2"}], ")"}], " ", "*", " ", 
                RowBox[{"Sin", " ", "[", 
                 RowBox[{"2", "*", "n", "*", "t"}], "]"}]}]}], ")"}], "/", 
             RowBox[{"n", "^", "2"}]}], ",", " ", 
            RowBox[{
             RowBox[{
              RowBox[{"-", " ", "3"}], " ", "*", " ", 
              RowBox[{
               RowBox[{"t", "^", "2"}], " ", "/", " ", "2"}]}], " ", "+", " ", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"4", " ", "-", " ", 
                RowBox[{"4", " ", 
                 RowBox[{"Cos", "[", 
                  RowBox[{"n", " ", "*", " ", "t"}], "]"}]}]}], ")"}], "/", 
              RowBox[{"n", " ", "^", " ", "2"}]}]}], ",", " ", "0", ",", " ", 
            
            RowBox[{"3", "*", " ", 
             RowBox[{
              RowBox[{
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "n"}], " ", "*", "t"}], " ", "+", " ", 
                 RowBox[{"Sin", "[", 
                  RowBox[{"n", " ", "*", " ", "t"}], "]"}]}], ")"}], "^", 
               "2"}], " ", "/", " ", 
              RowBox[{"n", "^", "3"}]}]}], ",", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"14", " ", "*", " ", "n", " ", "*", " ", "t"}], " ", 
               "+", " ", 
               RowBox[{"3", " ", "*", " ", 
                RowBox[{"n", "^", "3"}], " ", "*", " ", 
                RowBox[{"t", "^", "3"}]}], " ", "+", " ", 
               RowBox[{
               "24", " ", "*", " ", "n", " ", "*", " ", "t", " ", "*", " ", 
                RowBox[{"Cos", "[", 
                 RowBox[{"n", " ", "*", " ", "t"}], "]"}]}], " ", "-", " ", 
               RowBox[{"32", " ", "*", " ", 
                RowBox[{"Sin", "[", 
                 RowBox[{"n", " ", "*", " ", "t"}], "]"}]}], " ", "-", " ", 
               RowBox[{"3", " ", "*", " ", 
                RowBox[{"Sin", "[", " ", 
                 RowBox[{"2", " ", "*", "n", " ", "*", "t"}], "]"}]}]}], 
              ")"}], " ", "/", " ", 
             RowBox[{"n", " ", "^", "3"}]}], " ", ",", " ", "0"}], " ", 
           "\[IndentingNewLine]", "}"}], ",", " ", "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{"0", ",", " ", "0", ",", "  ", 
            RowBox[{
             RowBox[{
              RowBox[{"Sin", "[", 
               RowBox[{"n", " ", "*", " ", "t"}], "]"}], " ", "^", " ", "2"}],
              " ", "/", " ", 
             RowBox[{"(", 
              RowBox[{"2", " ", "*", " ", 
               RowBox[{"n", "^", "2"}]}], ")"}]}], ",", " ", "0", ",", " ", 
            "0", ",", " ", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"2", " ", "*", " ", "n", " ", "*", " ", "t"}], " ", 
               "-", " ", 
               RowBox[{"Sin", "[", 
                RowBox[{"2", " ", "*", " ", "n", " ", "*", "t"}], " ", 
                "]"}]}], ")"}], " ", "/", " ", 
             RowBox[{"(", 
              RowBox[{"4", " ", "*", " ", 
               RowBox[{"n", "^", "3"}]}], ")"}]}]}], "\[IndentingNewLine]", 
           "}"}]}], "\[IndentingNewLine]", "}"}]}], ";", 
       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]", "\[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]"}]}]], "Input",
 CellChangeTimes->{{3.922125494138694*^9, 3.922125528969516*^9}, {
   3.922125718901556*^9, 3.922125751610413*^9}, {3.922126978187389*^9, 
   3.9221269916052513`*^9}, {3.92280286749896*^9, 3.922802893835882*^9}, {
   3.922802933355482*^9, 3.92280298334046*^9}, {3.922803101441881*^9, 
   3.922803141051291*^9}, {3.922805107573048*^9, 3.922805161495988*^9}, {
   3.922805240214939*^9, 3.922805240325514*^9}, {3.923729935517523*^9, 
   3.923729936238494*^9}, {3.923731435264474*^9, 3.9237314544649553`*^9}, {
   3.9237316096822863`*^9, 3.923731610560288*^9}, {3.923735874979231*^9, 
   3.923735878507802*^9}, {3.92373609001552*^9, 3.9237360921921463`*^9}, {
   3.923736169308229*^9, 3.92373617964412*^9}, {3.923737104093314*^9, 
   3.9237371045879517`*^9}, {3.923737267148921*^9, 3.923737267952726*^9}, {
   3.925224670660235*^9, 3.92522467839744*^9}, {3.925390776789486*^9, 
   3.925390777437336*^9}, 3.925392605723798*^9, {3.925396851597014*^9, 
   3.925396856478215*^9}, {3.925398796765225*^9, 3.925398817962708*^9}, {
   3.925399177132155*^9, 3.9253992517095327`*^9}, 3.925399425013669*^9, {
   3.925400412071677*^9, 3.925400425115653*^9}, {3.925401383881632*^9, 
   3.9254013862167807`*^9}, 3.925402083479985*^9, 3.925402122722498*^9, 
   3.9257324022950687`*^9},
 CellLabel->
  "In[444]:=",ExpressionUUID->"e8b154ed-e19c-47e3-a98d-7cc40f6a965d"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "*", " ", "Given", " ", "points", " ", "on", " ", "the", " ", "boundary", 
     " ", "of", " ", "the", " ", "thrust", " ", "limited", " ", "reachable", 
     " ", "set"}], ",", " ", 
    RowBox[{
     RowBox[{
     "and", " ", "given", " ", "the", " ", "associated", " ", "energy"}], "-", 
     RowBox[{"quadratic", " ", "form"}]}], ",", " ", 
    RowBox[{"computes", " ", "min", " ", "and", " ", "max", " ", "ratios"}], 
    ",", " ", 
    RowBox[{
    "as", " ", "well", " ", "as", " ", "the", " ", "angles", " ", "at", " ", 
     "which", " ", "these", " ", "ratios", " ", "occur"}]}], " ", "**)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"reachablestatsoperator", "[", 
     RowBox[{"b_", ",", " ", "points_", ",", " ", "em_"}], "]"}], " ", ":=", 
    " ", 
    RowBox[{"Module", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"rats", ",", " ", "vecs"}], "}"}], ",", " ", 
      "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{
        RowBox[{"ratio", "[", 
         RowBox[{"pos_", ",", " ", "quad_"}], "]"}], ":=", " ", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"em", " ", "*", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"Eigenvalues", "[", 
                RowBox[{"{", 
                 RowBox[{
                  RowBox[{"Join", "[", 
                   RowBox[{
                    RowBox[{"Join", "[", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"{", 
                    RowBox[{"pos", " ", "//", " ", "Normalize"}], "}"}], " ", 
                    "//", " ", "Transpose"}], ")"}], ".", 
                    RowBox[{"{", 
                    RowBox[{"pos", " ", "//", " ", "Normalize"}], "}"}]}], 
                    ",", " ", 
                    RowBox[{"{", 
                    RowBox[{
                    RowBox[{"{", 
                    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", " ", 
                    RowBox[{"{", 
                    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", " ", 
                    RowBox[{"{", 
                    RowBox[{"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[{"quad", "[", 
                   RowBox[{"[", 
                    RowBox[{
                    RowBox[{"7", ";;", " ", "12"}], ",", " ", 
                    RowBox[{"7", " ", ";;", " ", "12"}]}], "]"}], "]"}]}], 
                 "}"}], "]"}], ")"}], "[", 
              RowBox[{"[", "1", "]"}], "]"}]}], ")"}], " ", "//", " ", 
           "Sqrt"}], ")"}], " ", "/", " ", 
         RowBox[{"Norm", "[", "pos", "]"}]}]}], ";", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"maxminrat", "[", 
         RowBox[{"a_", ",", " ", "v1s_"}], "]"}], " ", ":=", " ", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"MapThread", "[", 
            RowBox[{"ratio", ",", " ", 
             RowBox[{"{", 
              RowBox[{"v1s", ",", " ", 
               RowBox[{"ConstantArray", "[", 
                RowBox[{"a", " ", ",", " ", 
                 RowBox[{"v1s", " ", "//", " ", "Length"}]}], "]"}]}], 
              "}"}]}], "]"}], " ", "//", " ", "Max"}], ",", "   ", 
          RowBox[{
           RowBox[{"MapThread", "[", 
            RowBox[{"ratio", ",", " ", 
             RowBox[{"{", 
              RowBox[{"v1s", ",", " ", 
               RowBox[{"ConstantArray", "[", 
                RowBox[{"a", " ", ",", " ", 
                 RowBox[{"v1s", " ", "//", " ", "Length"}]}], "]"}]}], 
              "}"}]}], "]"}], " ", "//", " ", "Min"}]}], " ", "}"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"maxminvec", "[", 
         RowBox[{"a_", ",", " ", "v1s_"}], "]"}], " ", ":=", " ", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"v1s", "[", 
           RowBox[{"[", "  ", 
            RowBox[{
             RowBox[{"Position", "[", 
              RowBox[{
               RowBox[{"MapThread", "[", 
                RowBox[{"ratio", ",", " ", 
                 RowBox[{"{", 
                  RowBox[{"v1s", ",", " ", 
                   RowBox[{"ConstantArray", "[", 
                    RowBox[{"a", " ", ",", " ", 
                    RowBox[{"v1s", " ", "//", " ", "Length"}]}], "]"}]}], 
                  "}"}]}], "]"}], ",", " ", 
               RowBox[{
                RowBox[{"MapThread", "[", 
                 RowBox[{"ratio", ",", " ", 
                  RowBox[{"{", 
                   RowBox[{"v1s", ",", " ", 
                    RowBox[{"ConstantArray", "[", 
                    RowBox[{"a", " ", ",", " ", 
                    RowBox[{"v1s", " ", "//", " ", "Length"}]}], "]"}]}], 
                   "}"}]}], "]"}], " ", "//", " ", "Max"}]}], "    ", "]"}], 
             "[", 
             RowBox[{"[", 
              RowBox[{"1", ",", " ", "1"}], "]"}], "]"}], " ", "]"}], "]"}], 
          ",", " ", 
          RowBox[{"v1s", "[", 
           RowBox[{"[", "  ", 
            RowBox[{
             RowBox[{"Position", "[", 
              RowBox[{
               RowBox[{"MapThread", "[", 
                RowBox[{"ratio", ",", " ", 
                 RowBox[{"{", 
                  RowBox[{"v1s", ",", " ", 
                   RowBox[{"ConstantArray", "[", 
                    RowBox[{"a", " ", ",", " ", 
                    RowBox[{"v1s", " ", "//", " ", "Length"}]}], "]"}]}], 
                  "}"}]}], "]"}], ",", " ", 
               RowBox[{
                RowBox[{"MapThread", "[", 
                 RowBox[{"ratio", ",", " ", 
                  RowBox[{"{", 
                   RowBox[{"v1s", ",", " ", 
                    RowBox[{"ConstantArray", "[", 
                    RowBox[{"a", " ", ",", " ", 
                    RowBox[{"v1s", " ", "//", " ", "Length"}]}], "]"}]}], 
                   "}"}]}], "]"}], " ", "//", " ", "Min"}]}], "    ", "]"}], 
             "[", 
             RowBox[{"[", 
              RowBox[{"1", ",", " ", "1"}], "]"}], "]"}], " ", "]"}], "]"}]}],
          "}"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"rats", " ", "=", " ", 
        RowBox[{"maxminrat", "[", 
         RowBox[{"b", ",", " ", "points"}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"vecs", " ", "=", " ", 
        RowBox[{"maxminvec", "[", 
         RowBox[{"b", ",", " ", "points"}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{"rats", " ", ",", " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"vecs", "[", 
            RowBox[{"[", 
             RowBox[{"1", ",", " ", 
              RowBox[{"1", ";;", "2"}]}], "]"}], "]"}], ",", " ", 
           RowBox[{"vecs", "[", 
            RowBox[{"[", 
             RowBox[{"2", ",", " ", 
              RowBox[{"1", ";;", "2"}]}], "]"}], "]"}]}], "}"}]}], "}"}]}]}], 
     "\[IndentingNewLine]", "]"}]}], " ", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"*", " ", "\[IndentingNewLine]", 
      RowBox[{"reachablestatsoperatordeb2", "[", 
       RowBox[{"b_", ",", " ", "points_", ",", " ", "em_"}], "]"}]}], " ", ":=",
      " ", 
     RowBox[{"Module", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"rats", ",", " ", "vecs"}], "}"}], ",", " ", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{
         RowBox[{"ratio", "[", 
          RowBox[{"pos_", ",", " ", "quad_"}], "]"}], ":=", " ", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"em", " ", "*", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"Eigenvalues", "[", 
                 RowBox[{"{", 
                  RowBox[{
                   RowBox[{"Join", "[", 
                    RowBox[{
                    RowBox[{"Join", "[", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"{", 
                    RowBox[{"pos", " ", "//", " ", "Normalize"}], "}"}], " ", 
                    "//", " ", "Transpose"}], ")"}], ".", 
                    RowBox[{"{", 
                    RowBox[{"pos", " ", "//", " ", "Normalize"}], "}"}]}], 
                    ",", " ", 
                    RowBox[{"{", 
                    RowBox[{
                    RowBox[{"{", 
                    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", " ", 
                    RowBox[{"{", 
                    RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", " ", 
                    RowBox[{"{", 
                    RowBox[{"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[{"quad", "[", 
                    RowBox[{"[", 
                    RowBox[{
                    RowBox[{"7", ";;", " ", "12"}], ",", " ", 
                    RowBox[{"7", " ", ";;", " ", "12"}]}], "]"}], "]"}]}], 
                  "}"}], "]"}], ")"}], "[", 
               RowBox[{"[", "1", "]"}], "]"}]}], ")"}], " ", "//", " ", 
            "Sqrt"}], ")"}], " ", "/", " ", 
          RowBox[{"Norm", "[", "pos", "]"}]}]}], ";", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"maxminrat", "[", 
          RowBox[{"a_", ",", " ", "v1s_"}], "]"}], " ", ":=", " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{"MapThread", "[", 
             RowBox[{"ratio", ",", " ", 
              RowBox[{"{", 
               RowBox[{"v1s", ",", " ", 
                RowBox[{"ConstantArray", "[", 
                 RowBox[{"a", " ", ",", " ", 
                  RowBox[{"v1s", " ", "//", " ", "Length"}]}], "]"}]}], 
               "}"}]}], "]"}], " ", "//", " ", "Max"}], ",", "   ", 
           RowBox[{
            RowBox[{"MapThread", "[", 
             RowBox[{"ratio", ",", " ", 
              RowBox[{"{", 
               RowBox[{"v1s", ",", " ", 
                RowBox[{"ConstantArray", "[", 
                 RowBox[{"a", " ", ",", " ", 
                  RowBox[{"v1s", " ", "//", " ", "Length"}]}], "]"}]}], 
               "}"}]}], "]"}], " ", "//", " ", "Min"}]}], " ", "}"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"maxminvec", "[", 
          RowBox[{"a_", ",", " ", "v1s_"}], "]"}], " ", ":=", " ", 
         RowBox[{"v1s", "[", 
          RowBox[{"[", "  ", 
           RowBox[{"Position", "[", 
            RowBox[{
             RowBox[{"MapThread", "[", 
              RowBox[{"ratio", ",", " ", 
               RowBox[{"{", 
                RowBox[{"v1s", ",", " ", 
                 RowBox[{"ConstantArray", "[", 
                  RowBox[{"a", " ", ",", " ", 
                   RowBox[{"v1s", " ", "//", " ", "Length"}]}], "]"}]}], 
                "}"}]}], "]"}], ",", " ", 
             RowBox[{
              RowBox[{"MapThread", "[", 
               RowBox[{"ratio", ",", " ", 
                RowBox[{"{", 
                 RowBox[{"v1s", ",", " ", 
                  RowBox[{"ConstantArray", "[", 
                   RowBox[{"a", " ", ",", " ", 
                    RowBox[{"v1s", " ", "//", " ", "Length"}]}], "]"}]}], 
                 "}"}]}], "]"}], " ", "//", " ", "Max"}]}], "    ", "]"}], 
           " ", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
        "\[IndentingNewLine]", 
        RowBox[{"vecs", " ", "=", " ", 
         RowBox[{"maxminvec", "[", 
          RowBox[{"b", ",", " ", "points"}], "]"}]}], ";", 
        "\[IndentingNewLine]", "vecs"}]}], "\[IndentingNewLine]", "]"}]}], 
    " ", "\[IndentingNewLine]", "**)"}], " "}]}]], "Input",
 CellChangeTimes->{{3.922126565314114*^9, 3.922126580771452*^9}, {
   3.922126704468804*^9, 3.922126705947434*^9}, 3.922126901551271*^9, 
   3.9221287324070587`*^9, {3.9221291080143538`*^9, 3.922129119159079*^9}, {
   3.9221291992585163`*^9, 3.922129215343298*^9}, {3.922129258128023*^9, 
   3.9221292896449633`*^9}, {3.922129324289394*^9, 3.9221293491407127`*^9}, {
   3.922129421307498*^9, 3.922129442139203*^9}, {3.922129475190606*^9, 
   3.922129482963978*^9}, {3.922129563707798*^9, 3.9221296621161118`*^9}, {
   3.9221299022067547`*^9, 3.9221299154689493`*^9}, {3.922130026196035*^9, 
   3.922130073548706*^9}, {3.922130127501897*^9, 3.922130128714967*^9}, {
   3.9221302355590467`*^9, 3.922130258260676*^9}, {3.9221303248841333`*^9, 
   3.922130325435637*^9}, {3.9222801711925163`*^9, 3.9222801740848722`*^9}, {
   3.922805176894771*^9, 3.9228051821294937`*^9}, {3.9228052962390738`*^9, 
   3.9228053149612913`*^9}, {3.9228054188591337`*^9, 3.922805419696224*^9}, 
   3.923731495073716*^9, {3.92373575092951*^9, 3.92373575125555*^9}, {
   3.923735802630307*^9, 3.9237358098486433`*^9}, {3.9237358438933887`*^9, 
   3.923735845617166*^9}, {3.923735952599894*^9, 3.923735953595971*^9}, {
   3.923736364681592*^9, 3.923736373905263*^9}, {3.923736665499807*^9, 
   3.923736670141395*^9}, {3.9237368530528383`*^9, 3.92373691075611*^9}, {
   3.923737050541265*^9, 3.923737051272838*^9}, {3.925221832769929*^9, 
   3.925221832999977*^9}, {3.9252224410477858`*^9, 3.9252224417497063`*^9}, {
   3.925224119017887*^9, 3.925224120207304*^9}, {3.925224194821704*^9, 
   3.925224308175753*^9}, {3.925224338581411*^9, 3.925224401138693*^9}, {
   3.925224558765996*^9, 3.92522465806172*^9}, {3.925224701592553*^9, 
   3.925224751420607*^9}, {3.925224790531036*^9, 3.9252248783316727`*^9}, {
   3.925224914043165*^9, 3.925225012937273*^9}, {3.925225243610983*^9, 
   3.925225266052886*^9}, {3.925225303427202*^9, 3.925225303708529*^9}, {
   3.925226094619844*^9, 3.92522610210445*^9}, {3.925228172419969*^9, 
   3.92522820483424*^9}, {3.925229037988138*^9, 3.925229088448476*^9}, {
   3.925230359394208*^9, 3.925230376521163*^9}, {3.925390513633038*^9, 
   3.925390601610162*^9}, {3.925390647928447*^9, 3.925390670514942*^9}, {
   3.925390770669362*^9, 3.925390772512945*^9}, {3.925390844027473*^9, 
   3.9253908462704144`*^9}, {3.925390984141727*^9, 3.925391028836824*^9}, {
   3.925391382136241*^9, 3.9253914070561934`*^9}, {3.925391631443481*^9, 
   3.925391640251116*^9}, {3.925391815375328*^9, 3.9253918572642403`*^9}, {
   3.925392072302026*^9, 3.925392081394278*^9}, {3.925392114805657*^9, 
   3.925392123162113*^9}, {3.925392154790612*^9, 3.92539219201752*^9}, {
   3.925392359577207*^9, 3.9253923794095497`*^9}, {3.925398771223444*^9, 
   3.9253987769295607`*^9}, 3.925398913717326*^9, {3.925399405698003*^9, 
   3.925399428731737*^9}, {3.925400754790189*^9, 3.9254007677067432`*^9}, {
   3.925401389745448*^9, 3.9254013898255157`*^9}, {3.925732381773876*^9, 
   3.9257324073138113`*^9}, {3.9257324778067017`*^9, 3.925732530592381*^9}, {
   3.926683451828351*^9, 3.926683478615221*^9}, 3.926683513636981*^9, {
   3.926683544016572*^9, 3.9266835810251226`*^9}, {3.926683635140427*^9, 
   3.926683654277988*^9}, {3.92668376628743*^9, 3.926683775845339*^9}, 
   3.9266838094619293`*^9, {3.926683974642449*^9, 3.926683978064101*^9}},
 CellLabel->
  "In[445]:=",ExpressionUUID->"c4ceae39-7397-4043-ac31-479218ca2bcb"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"*", " ", "Sample", " ", "parameter", " ", "values"}], " ", 
   "**)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"mu1", " ", "=", " ", 
     RowBox[{"3.986", " ", 
      RowBox[{"10", "^", "14"}]}]}], " ", ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"endtime", " ", "=", " ", "1500"}], ";"}], "\[IndentingNewLine]", 
   
   RowBox[{
    RowBox[{"alpha", " ", "=", " ", "7780000"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"n", " ", "=", " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"mu1", " ", "/", " ", 
        RowBox[{"alpha", "^", "3"}]}], ")"}], "^", 
      RowBox[{"(", 
       RowBox[{"1", "/", "2"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"tmax", "=", 
     RowBox[{"2", "/", "1000"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"comparableenergy", " ", "=", " ", 
     RowBox[{"endtime", " ", 
      RowBox[{"tmax", "^", "2"}]}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.922126565314114*^9, 3.922126580771452*^9}, {
   3.922126704468804*^9, 3.922126705947434*^9}, 3.922126901551271*^9, 
   3.9221287324070587`*^9, {3.9221291080143538`*^9, 3.922129119159079*^9}, {
   3.9221291992585163`*^9, 3.922129215343298*^9}, {3.922129258128023*^9, 
   3.9221292896449633`*^9}, {3.922129324289394*^9, 3.9221293491407127`*^9}, {
   3.922129421307498*^9, 3.922129442139203*^9}, {3.922129475190606*^9, 
   3.922129482963978*^9}, {3.922129563707798*^9, 3.9221296621161118`*^9}, {
   3.9221299022067547`*^9, 3.9221299154689493`*^9}, {3.922130026196035*^9, 
   3.922130073548706*^9}, {3.922130127501897*^9, 3.922130128714967*^9}, {
   3.9221302355590467`*^9, 3.922130258260676*^9}, {3.9221303248841333`*^9, 
   3.922130325435637*^9}, {3.9222801711925163`*^9, 3.9222801740848722`*^9}, {
   3.922805176894771*^9, 3.9228051821294937`*^9}, {3.9228052962390738`*^9, 
   3.9228053149612913`*^9}, {3.9228054188591337`*^9, 3.922805419696224*^9}, 
   3.923731495073716*^9, {3.92373575092951*^9, 3.92373575125555*^9}, {
   3.923735802630307*^9, 3.9237358098486433`*^9}, {3.9237358438933887`*^9, 
   3.923735845617166*^9}, {3.923735952599894*^9, 3.923735953595971*^9}, {
   3.923736364681592*^9, 3.923736373905263*^9}, {3.923736665499807*^9, 
   3.923736670141395*^9}, {3.9237368530528383`*^9, 3.92373691075611*^9}, {
   3.923737050541265*^9, 3.923737051272838*^9}, {3.925221832769929*^9, 
   3.925221832999977*^9}, {3.9252224410477858`*^9, 3.9252224417497063`*^9}, {
   3.925224119017887*^9, 3.925224120207304*^9}, {3.925224194821704*^9, 
   3.925224308175753*^9}, {3.925224338581411*^9, 3.925224401138693*^9}, {
   3.925224558765996*^9, 3.92522465806172*^9}, {3.925224701592553*^9, 
   3.925224751420607*^9}, {3.925224790531036*^9, 3.9252248783316727`*^9}, {
   3.925224914043165*^9, 3.925225012937273*^9}, {3.925225243610983*^9, 
   3.925225266052886*^9}, {3.925225303427202*^9, 3.925225303708529*^9}, {
   3.925226094619844*^9, 3.92522610210445*^9}, {3.925228172419969*^9, 
   3.92522820483424*^9}, {3.925229037988138*^9, 3.925229088448476*^9}, {
   3.925230359394208*^9, 3.925230376521163*^9}, {3.925390513633038*^9, 
   3.925390601610162*^9}, {3.925390647928447*^9, 3.925390670514942*^9}, {
   3.925390770669362*^9, 3.925390772512945*^9}, {3.925390844027473*^9, 
   3.9253908462704144`*^9}, {3.925390984141727*^9, 3.925391028836824*^9}, {
   3.925391382136241*^9, 3.9253914070561934`*^9}, {3.925391631443481*^9, 
   3.925391640251116*^9}, {3.925391815375328*^9, 3.9253918572642403`*^9}, {
   3.925392072302026*^9, 3.925392081394278*^9}, {3.925392114805657*^9, 
   3.925392123162113*^9}, {3.925392154790612*^9, 3.92539219201752*^9}, {
   3.925392359577207*^9, 3.9253923794095497`*^9}, {3.925398771223444*^9, 
   3.9253987769295607`*^9}, 3.925398913717326*^9, {3.925399405698003*^9, 
   3.925399428731737*^9}, {3.925400754790189*^9, 3.9254007677067432`*^9}, {
   3.925401389745448*^9, 3.9254013898255157`*^9}, {3.925732381773876*^9, 
   3.9257324073138113`*^9}, 3.9257324778067017`*^9, {3.925732535345351*^9, 
   3.925732550063822*^9}},
 CellLabel->
  "In[446]:=",ExpressionUUID->"d82cf4ae-f3bf-4a07-b427-7d6864e8cdbf"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "*", " ", "test", " ", "of", " ", "comparison", " ", "statistics"}], " ", 
   "**)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"testbounds", " ", "=", " ", 
     RowBox[{"thrustreachable", "[", 
      RowBox[{"xrf", ",", " ", "endtime", ",", " ", "n", ",", " ", "tmax"}], 
      "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"testbounds", "[", 
    RowBox[{"[", "1", "]"}], "]"}], "\[IndentingNewLine]", 
   RowBox[{"teststats", " ", "=", " ", 
    RowBox[{"reachablestatsoperator", "[", " ", 
     RowBox[{
      RowBox[{"energyform", "[", 
       RowBox[{"1600", ",", " ", "n"}], "]"}], ",", "testbounds", ",", " ", 
      "comparableenergy"}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.922126565314114*^9, 3.922126580771452*^9}, {
   3.922126704468804*^9, 3.922126705947434*^9}, 3.922126901551271*^9, 
   3.9221287324070587`*^9, {3.9221291080143538`*^9, 3.922129119159079*^9}, {
   3.9221291992585163`*^9, 3.922129215343298*^9}, {3.922129258128023*^9, 
   3.9221292896449633`*^9}, {3.922129324289394*^9, 3.9221293491407127`*^9}, {
   3.922129421307498*^9, 3.922129442139203*^9}, {3.922129475190606*^9, 
   3.922129482963978*^9}, {3.922129563707798*^9, 3.9221296621161118`*^9}, {
   3.9221299022067547`*^9, 3.9221299154689493`*^9}, {3.922130026196035*^9, 
   3.922130073548706*^9}, {3.922130127501897*^9, 3.922130128714967*^9}, {
   3.9221302355590467`*^9, 3.922130258260676*^9}, {3.9221303248841333`*^9, 
   3.922130325435637*^9}, {3.9222801711925163`*^9, 3.9222801740848722`*^9}, {
   3.922805176894771*^9, 3.9228051821294937`*^9}, {3.9228052962390738`*^9, 
   3.9228053149612913`*^9}, {3.9228054188591337`*^9, 3.922805419696224*^9}, 
   3.923731495073716*^9, {3.92373575092951*^9, 3.92373575125555*^9}, {
   3.923735802630307*^9, 3.9237358098486433`*^9}, {3.9237358438933887`*^9, 
   3.923735845617166*^9}, {3.923735952599894*^9, 3.923735953595971*^9}, {
   3.923736364681592*^9, 3.923736373905263*^9}, {3.923736665499807*^9, 
   3.923736670141395*^9}, {3.9237368530528383`*^9, 3.92373691075611*^9}, {
   3.923737050541265*^9, 3.923737051272838*^9}, {3.925221832769929*^9, 
   3.925221832999977*^9}, {3.9252224410477858`*^9, 3.9252224417497063`*^9}, {
   3.925224119017887*^9, 3.925224120207304*^9}, {3.925224194821704*^9, 
   3.925224308175753*^9}, {3.925224338581411*^9, 3.925224401138693*^9}, {
   3.925224558765996*^9, 3.92522465806172*^9}, {3.925224701592553*^9, 
   3.925224751420607*^9}, {3.925224790531036*^9, 3.9252248783316727`*^9}, {
   3.925224914043165*^9, 3.925225012937273*^9}, {3.925225243610983*^9, 
   3.925225266052886*^9}, {3.925225303427202*^9, 3.925225303708529*^9}, {
   3.925226094619844*^9, 3.92522610210445*^9}, {3.925228172419969*^9, 
   3.92522820483424*^9}, {3.925229037988138*^9, 3.925229088448476*^9}, {
   3.925230359394208*^9, 3.925230376521163*^9}, {3.925390513633038*^9, 
   3.925390601610162*^9}, {3.925390647928447*^9, 3.925390670514942*^9}, {
   3.925390770669362*^9, 3.925390772512945*^9}, {3.925390844027473*^9, 
   3.9253908462704144`*^9}, {3.925390984141727*^9, 3.925391028836824*^9}, {
   3.925391382136241*^9, 3.9253914070561934`*^9}, {3.925391631443481*^9, 
   3.925391640251116*^9}, {3.925391815375328*^9, 3.9253918572642403`*^9}, {
   3.925392072302026*^9, 3.925392081394278*^9}, {3.925392114805657*^9, 
   3.925392123162113*^9}, {3.925392154790612*^9, 3.92539219201752*^9}, {
   3.925392359577207*^9, 3.9253923794095497`*^9}, {3.925398771223444*^9, 
   3.9253987769295607`*^9}, 3.925398913717326*^9, {3.925399405698003*^9, 
   3.925399428731737*^9}, {3.925400754790189*^9, 3.9254007677067432`*^9}, {
   3.925401389745448*^9, 3.9254013898255157`*^9}, {3.925732381773876*^9, 
   3.9257324073138113`*^9}, 3.9257324778067017`*^9, {3.925732535345351*^9, 
   3.925732579478861*^9}, {3.926164751301188*^9, 3.926164753074662*^9}, {
   3.926589345663061*^9, 3.926589347450185*^9}, {3.926590053174768*^9, 
   3.9265900550693693`*^9}, {3.926683682595344*^9, 3.926683683204872*^9}, {
   3.926683984804964*^9, 3.926683985476653*^9}},
 CellLabel->
  "In[452]:=",ExpressionUUID->"5bac8550-244b-4b24-a50f-7e5196fb9d83"],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate", "izero", 
   "\"Integral and error estimates are 0 on all integration subregions. Try \
increasing the value of the MinRecursion option. If value of integral may be \
0, specify a finite value for the AccuracyGoal option.\"", 2, 452, 77, 
   30967495211780212173, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{
  3.926595031341083*^9, 3.9265951516707287`*^9, 3.9265954474845057`*^9, 
   3.9266834839602957`*^9, 3.926683521690254*^9, 3.92668368723067*^9, 
   3.926683814515615*^9, {3.92668398951153*^9, 3.926684006535088*^9}},
 CellLabel->
  "During evaluation of \
In[452]:=",ExpressionUUID->"eff9712d-9f55-4034-a2be-6f2e3e3fbfad"],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate", "izero", 
   "\"Integral and error estimates are 0 on all integration subregions. Try \
increasing the value of the MinRecursion option. If value of integral may be \
0, specify a finite value for the AccuracyGoal option.\"", 2, 452, 78, 
   30967495211780212173, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{
  3.926595031341083*^9, 3.9265951516707287`*^9, 3.9265954474845057`*^9, 
   3.9266834839602957`*^9, 3.926683521690254*^9, 3.92668368723067*^9, 
   3.926683814515615*^9, {3.92668398951153*^9, 3.92668400654907*^9}},
 CellLabel->
  "During evaluation of \
In[452]:=",ExpressionUUID->"eac3f0c3-802f-4c77-9220-00843369023f"],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate", "izero", 
   "\"Integral and error estimates are 0 on all integration subregions. Try \
increasing the value of the MinRecursion option. If value of integral may be \
0, specify a finite value for the AccuracyGoal option.\"", 2, 452, 79, 
   30967495211780212173, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{
  3.926595031341083*^9, 3.9265951516707287`*^9, 3.9265954474845057`*^9, 
   3.9266834839602957`*^9, 3.926683521690254*^9, 3.92668368723067*^9, 
   3.926683814515615*^9, {3.92668398951153*^9, 3.9266840065677567`*^9}},
 CellLabel->
  "During evaluation of \
In[452]:=",ExpressionUUID->"4f784c11-ffc5-48de-b2a7-267534e8e887"],

Cell[BoxData[
 TemplateBox[{
  "General", "stop", 
   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"NIntegrate\\\", \
\\\"::\\\", \\\"izero\\\"}], \\\"MessageName\\\"]\\) will be suppressed \
during this calculation.\"", 2, 452, 80, 30967495211780212173, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{
  3.926595031341083*^9, 3.9265951516707287`*^9, 3.9265954474845057`*^9, 
   3.9266834839602957`*^9, 3.926683521690254*^9, 3.92668368723067*^9, 
   3.926683814515615*^9, {3.92668398951153*^9, 3.9266840065704*^9}},
 CellLabel->
  "During evaluation of \
In[452]:=",ExpressionUUID->"65cf75d5-3200-4fa4-ac34-291b9a9c1807"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"940.8270021995455`", ",", 
   RowBox[{"-", "2295.3700606773796`"}], ",", "0.`"}], "}"}]], "Output",
 CellChangeTimes->{
  3.9265950408120813`*^9, 3.926595161089231*^9, 3.926595456767788*^9, 
   3.9266834934970617`*^9, 3.926683531026675*^9, 3.926683696531247*^9, 
   3.926683823848234*^9, {3.926683999005316*^9, 3.926684016138186*^9}},
 CellLabel->
  "Out[453]=",ExpressionUUID->"a3aaf8db-05f5-43c5-8a7c-f859e2f64fec"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"1.4950655015545977`", ",", "1.2205414440778182`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "815.5673445233721`"}], ",", "2292.9850514040545`"}], 
      "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"1154.587620745528`", ",", "1612.8072558243039`"}], "}"}]}], 
    "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.9265950408120813`*^9, 3.926595161089231*^9, 3.926595456767788*^9, 
   3.9266834934970617`*^9, 3.926683531026675*^9, 3.926683696531247*^9, 
   3.926683823848234*^9, {3.926683999005316*^9, 3.926684016886485*^9}},
 CellLabel->
  "Out[454]=",ExpressionUUID->"5d0e7b46-eaab-4ee8-8572-a4759784af16"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "*", " ", "Threading", " ", "the", " ", "different", " ", "terminal", " ",
      "times", " ", "through", " ", "our", " ", "energy", " ", "and", " ", 
     "thrust", " ", "limited", " ", "RRS", " ", "computations"}], ",", " ", 
    RowBox[{
    "and", " ", "then", " ", "through", " ", "the", " ", "ratio", " ", 
     "computations"}]}], " ", "**)"}], " ", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"starts", " ", "//", " ", "Clear"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"finaltime", " ", "//", " ", "Clear"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"starts", " ", "=", " ", "1500"}], ";"}], " ", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"finaltime", " ", "=", " ", "1700"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"terminaltimes", " ", "=", " ", 
     RowBox[{"Range", "[", 
      RowBox[{"starts", ",", " ", "finaltime", ",", " ", "50"}], "]"}]}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ns", " ", "=", " ", 
     RowBox[{"ConstantArray", "[", 
      RowBox[{"n", ",", " ", 
       RowBox[{"terminaltimes", " ", "//", " ", "Length"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"tmaxs", " ", "=", " ", 
     RowBox[{"ConstantArray", "[", 
      RowBox[{"tmax", ",", " ", 
       RowBox[{"terminaltimes", " ", "//", " ", "Length"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"xrfs", " ", "=", " ", 
     RowBox[{"ConstantArray", "[", 
      RowBox[{"xrf", ",", " ", 
       RowBox[{"terminaltimes", " ", "//", " ", "Length"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"bounds", " ", "=", " ", 
     RowBox[{"MapThread", "[", 
      RowBox[{"thrustreachable", ",", " ", 
       RowBox[{"{", 
        RowBox[{
        "xrfs", ",", " ", "terminaltimes", ",", " ", "ns", ",", " ", 
         "tmaxs"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"forms", " ", "=", " ", 
     RowBox[{"MapThread", "[", 
      RowBox[{"energyform", ",", " ", 
       RowBox[{"{", 
        RowBox[{"terminaltimes", ",", " ", "ns"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"comparableenergies", " ", "=", " ", 
     RowBox[{"terminaltimes", " ", 
      RowBox[{"tmax", "^", "2"}]}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.922126565314114*^9, 3.922126580771452*^9}, {
   3.922126704468804*^9, 3.922126705947434*^9}, 3.922126901551271*^9, 
   3.9221287324070587`*^9, {3.9221291080143538`*^9, 3.922129119159079*^9}, {
   3.9221291992585163`*^9, 3.922129215343298*^9}, {3.922129258128023*^9, 
   3.9221292896449633`*^9}, {3.922129324289394*^9, 3.9221293491407127`*^9}, {
   3.922129421307498*^9, 3.922129442139203*^9}, {3.922129475190606*^9, 
   3.922129482963978*^9}, {3.922129563707798*^9, 3.9221296621161118`*^9}, {
   3.9221299022067547`*^9, 3.9221299154689493`*^9}, {3.922130026196035*^9, 
   3.922130073548706*^9}, {3.922130127501897*^9, 3.922130128714967*^9}, {
   3.9221302355590467`*^9, 3.922130258260676*^9}, {3.9221303248841333`*^9, 
   3.922130325435637*^9}, {3.9222801711925163`*^9, 3.9222801740848722`*^9}, {
   3.922805176894771*^9, 3.9228051821294937`*^9}, {3.9228052962390738`*^9, 
   3.9228053149612913`*^9}, {3.9228054188591337`*^9, 3.922805419696224*^9}, 
   3.923731495073716*^9, {3.92373575092951*^9, 3.92373575125555*^9}, {
   3.923735802630307*^9, 3.9237358098486433`*^9}, {3.9237358438933887`*^9, 
   3.923735845617166*^9}, {3.923735952599894*^9, 3.923735953595971*^9}, {
   3.923736364681592*^9, 3.923736373905263*^9}, {3.923736665499807*^9, 
   3.923736670141395*^9}, {3.9237368530528383`*^9, 3.92373691075611*^9}, {
   3.923737050541265*^9, 3.923737051272838*^9}, {3.925221832769929*^9, 
   3.925221832999977*^9}, {3.9252224410477858`*^9, 3.9252224417497063`*^9}, {
   3.925224119017887*^9, 3.925224120207304*^9}, {3.925224194821704*^9, 
   3.925224308175753*^9}, {3.925224338581411*^9, 3.925224401138693*^9}, {
   3.925224558765996*^9, 3.92522465806172*^9}, {3.925224701592553*^9, 
   3.925224751420607*^9}, {3.925224790531036*^9, 3.9252248783316727`*^9}, {
   3.925224914043165*^9, 3.925225012937273*^9}, {3.925225243610983*^9, 
   3.925225266052886*^9}, {3.925225303427202*^9, 3.925225303708529*^9}, {
   3.925226094619844*^9, 3.92522610210445*^9}, {3.925228172419969*^9, 
   3.92522820483424*^9}, {3.925229037988138*^9, 3.925229088448476*^9}, {
   3.925230359394208*^9, 3.925230376521163*^9}, {3.925390513633038*^9, 
   3.925390601610162*^9}, {3.925390647928447*^9, 3.925390670514942*^9}, {
   3.925390770669362*^9, 3.925390772512945*^9}, {3.925390844027473*^9, 
   3.9253908462704144`*^9}, {3.925390984141727*^9, 3.925391028836824*^9}, {
   3.925391382136241*^9, 3.9253914070561934`*^9}, {3.925391631443481*^9, 
   3.925391640251116*^9}, 3.925391815375328*^9, {3.925393775414851*^9, 
   3.92539378727888*^9}, {3.925394026278256*^9, 3.925394092562438*^9}, {
   3.925394850436942*^9, 3.925394853524756*^9}, {3.9253959392732363`*^9, 
   3.925395948346939*^9}, 3.925398766759432*^9, 3.9253989236872663`*^9, {
   3.925732624554319*^9, 3.925732636476462*^9}},
 CellLabel->
  "In[455]:=",ExpressionUUID->"8d5d3b56-3d98-4aa2-94a6-05733121195f"],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate", "izero", 
   "\"Integral and error estimates are 0 on all integration subregions. Try \
increasing the value of the MinRecursion option. If value of integral may be \
0, specify a finite value for the AccuracyGoal option.\"", 2, 463, 81, 
   30967495211780212173, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.926595041571376*^9, 3.92659516192593*^9, 
  3.92659545753093*^9, 3.92668401692129*^9},
 CellLabel->
  "During evaluation of \
In[455]:=",ExpressionUUID->"4f611342-2991-4279-bd60-64ef961eabc9"],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate", "izero", 
   "\"Integral and error estimates are 0 on all integration subregions. Try \
increasing the value of the MinRecursion option. If value of integral may be \
0, specify a finite value for the AccuracyGoal option.\"", 2, 463, 82, 
   30967495211780212173, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.926595041571376*^9, 3.92659516192593*^9, 
  3.92659545753093*^9, 3.926684016941003*^9},
 CellLabel->
  "During evaluation of \
In[455]:=",ExpressionUUID->"347a72b1-dd69-49f6-a75e-43ff8fbb3af6"],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate", "izero", 
   "\"Integral and error estimates are 0 on all integration subregions. Try \
increasing the value of the MinRecursion option. If value of integral may be \
0, specify a finite value for the AccuracyGoal option.\"", 2, 463, 83, 
   30967495211780212173, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.926595041571376*^9, 3.92659516192593*^9, 
  3.92659545753093*^9, 3.926684016946014*^9},
 CellLabel->
  "During evaluation of \
In[455]:=",ExpressionUUID->"38e85ed0-d423-4e11-a186-f308288d5025"],

Cell[BoxData[
 TemplateBox[{
  "General", "stop", 
   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"NIntegrate\\\", \
\\\"::\\\", \\\"izero\\\"}], \\\"MessageName\\\"]\\) will be suppressed \
during this calculation.\"", 2, 463, 84, 30967495211780212173, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.926595041571376*^9, 3.92659516192593*^9, 
  3.92659545753093*^9, 3.9266840169495783`*^9},
 CellLabel->
  "During evaluation of \
In[455]:=",ExpressionUUID->"42b5eef7-7bd6-4890-a9a4-8981c303a410"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"*", "Computing", " ", "the", " ", "commparisons"}], " ", "**)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"computedmmratios", " ", "=", " ", 
     RowBox[{"MapThread", "[", 
      RowBox[{"reachablestatsoperator", ",", " ", 
       RowBox[{"{", 
        RowBox[{"forms", ",", " ", "bounds", ",", " ", "comparableenergies"}],
         "}"}]}], " ", "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"maxrats", " ", "=", " ", 
     RowBox[{"computedmmratios", "[", 
      RowBox[{"[", 
       RowBox[{";;", ",", " ", "1", ",", " ", "1"}], "]"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"minrats", " ", "=", " ", 
     RowBox[{"computedmmratios", "[", 
      RowBox[{"[", 
       RowBox[{";;", ",", " ", "1", ",", " ", "2"}], "]"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"maxvecs", " ", "=", " ", 
     RowBox[{"computedmmratios", "[", 
      RowBox[{"[", 
       RowBox[{";;", ",", " ", "2", ",", " ", "1"}], "]"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"minvecs", " ", "=", " ", 
     RowBox[{"computedmmratios", "[", 
      RowBox[{"[", 
       RowBox[{";;", ",", " ", "2", ",", " ", "2"}], "]"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"vectoangle", "[", "k_", "]"}], " ", ":=", " ", 
     RowBox[{"Mod", "[", 
      RowBox[{
       RowBox[{"ArcTan", "[", 
        RowBox[{
         RowBox[{"k", "[", 
          RowBox[{"[", "1", "]"}], "]"}], ",", 
         RowBox[{"k", "[", 
          RowBox[{"[", "2", "]"}], "]"}]}], "]"}], ",", " ", "Pi"}], "]"}]}], 
    ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"maxangles", " ", "=", " ", 
     RowBox[{"Map", "[", 
      RowBox[{"vectoangle", ",", " ", "maxvecs"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"minangles", " ", "=", " ", 
     RowBox[{"Map", "[", 
      RowBox[{"vectoangle", ",", " ", "minvecs"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.925228209161273*^9, 3.9252282623345222`*^9}, {
   3.925228547553074*^9, 3.9252285869691057`*^9}, {3.925229003359222*^9, 
   3.925229019027536*^9}, {3.92522910616368*^9, 3.925229122342781*^9}, {
   3.925229187628474*^9, 3.925229192705195*^9}, {3.9252292535220833`*^9, 
   3.9252293230642138`*^9}, 3.92522939642738*^9, {3.925229456371237*^9, 
   3.925229459443162*^9}, {3.925229793982645*^9, 3.9252298350229692`*^9}, {
   3.9252298956887074`*^9, 3.925229901712088*^9}, {3.925230168532755*^9, 
   3.925230170088011*^9}, {3.9252308603700438`*^9, 3.925230865004015*^9}, {
   3.925231056433466*^9, 3.925231096187351*^9}, {3.925231162748597*^9, 
   3.925231240811164*^9}, {3.9253915974915533`*^9, 3.925391623802257*^9}, {
   3.925392788862295*^9, 3.9253928688039303`*^9}, {3.9253930227047997`*^9, 
   3.925393096173184*^9}, {3.925393708871986*^9, 3.925393715717189*^9}, {
   3.9253937494103193`*^9, 3.92539380770921*^9}, {3.9253940410763283`*^9, 
   3.925394048157661*^9}, {3.925394218562673*^9, 3.925394270437191*^9}, {
   3.9253943364122066`*^9, 3.9253944592534313`*^9}, {3.925394506159732*^9, 
   3.925394539707529*^9}, {3.925394585846653*^9, 3.9253947217617083`*^9}, {
   3.925732625999062*^9, 3.9257326790058823`*^9}, {3.92659480157522*^9, 
   3.926594819657933*^9}, {3.926594849837088*^9, 3.9265948818601294`*^9}, {
   3.92659497999238*^9, 3.9265950048502913`*^9}, {3.926595369479834*^9, 
   3.926595394481021*^9}},
 CellLabel->
  "In[466]:=",ExpressionUUID->"ce85b294-c6a9-4a49-8a9c-f7900fa64e64"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"*", " ", "Plotting", " ", "code"}], " ", "**)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"l1", " ", "=", " ", 
    RowBox[{"ListLinePlot", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"{", 
        RowBox[{"terminaltimes", ",", " ", "maxrats"}], "}"}], " ", "//", " ",
        "Transpose"}], ",", " ", 
      RowBox[{
      "PlotLabel", " ", "->", " ", "\"\<Max Ratios vs Terminal Time\>\""}], 
      ",", " ", 
      RowBox[{"AxesLabel", "\[Rule]", " ", 
       RowBox[{"{", 
        RowBox[{"\"\<Terminal Time (s)\>\"", ",", " ", "\"\<Max Ratio\>\""}], 
        " ", "}"}]}]}], "   ", "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"l2", " ", "=", " ", 
    RowBox[{"ListLinePlot", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"{", 
        RowBox[{"terminaltimes", ",", " ", "minrats"}], "}"}], " ", "//", " ",
        "Transpose"}], ",", " ", 
      RowBox[{
      "PlotLabel", "\[Rule]", " ", "\"\<Min Ratios vs Terminal Time\>\""}], 
      ",", " ", 
      RowBox[{"AxesLabel", "\[Rule]", " ", 
       RowBox[{"{", 
        RowBox[{"\"\<Terminal Time (s)\>\"", ",", " ", "\"\<Min Ratio\>\""}], 
        " ", "}"}]}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"l3", " ", "=", " ", 
    RowBox[{"ListLinePlot", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"{", 
        RowBox[{"terminaltimes", ",", " ", "maxangles"}], "}"}], " ", "//", 
       " ", "Transpose"}], ",", " ", 
      RowBox[{
      "PlotLabel", " ", "->", " ", "\"\<Max Angles vs Terminal Time\>\""}], 
      ",", " ", 
      RowBox[{"AxesLabel", "\[Rule]", " ", 
       RowBox[{"{", 
        RowBox[{
        "\"\<Terminal Time (s)\>\"", ",", " ", "\"\<Max Angle (Rad)\>\""}], 
        " ", "}"}]}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"l4", " ", "=", " ", 
    RowBox[{"ListLinePlot", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"{", 
        RowBox[{"terminaltimes", ",", " ", "minangles"}], "}"}], " ", "//", 
       " ", "Transpose"}], ",", " ", 
      RowBox[{"PlotLabel", "\[Rule]", "\"\<Min Angles vs Terminal Time\>\""}],
       ",", " ", 
      RowBox[{"AxesLabel", "\[Rule]", " ", 
       RowBox[{"{", 
        RowBox[{
        "\"\<Terminal Time (s)\>\"", ",", " ", "\"\<Min Angle (Rad)\>\""}], 
        "}"}]}]}], "]"}]}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"Export", "[", 
    RowBox[{"\"\<maxratt.png\>\"", ",", " ", "l1"}], "]"}], 
   "\[IndentingNewLine]", 
   RowBox[{"Export", "[", 
    RowBox[{"\"\<minratt.png\>\"", ",", " ", "l2"}], "]"}], 
   "\[IndentingNewLine]", 
   RowBox[{"Export", "[", 
    RowBox[{"\"\<maxanglet.png\>\"", ",", " ", "l3"}], "]"}], 
   "\[IndentingNewLine]", 
   RowBox[{"Export", "[", 
    RowBox[{"\"\<minanglet.png\>\"", ",", " ", "l4"}], "]"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.925228209161273*^9, 3.9252282623345222`*^9}, {
   3.925228547553074*^9, 3.9252285869691057`*^9}, {3.925229003359222*^9, 
   3.925229019027536*^9}, {3.92522910616368*^9, 3.925229122342781*^9}, {
   3.925229187628474*^9, 3.925229192705195*^9}, {3.9252292535220833`*^9, 
   3.9252293230642138`*^9}, 3.92522939642738*^9, {3.925229456371237*^9, 
   3.925229459443162*^9}, {3.925229793982645*^9, 3.9252298350229692`*^9}, {
   3.9252298956887074`*^9, 3.925229901712088*^9}, {3.925230168532755*^9, 
   3.925230170088011*^9}, {3.9252308603700438`*^9, 3.925230865004015*^9}, {
   3.925231056433466*^9, 3.925231096187351*^9}, {3.925231162748597*^9, 
   3.925231240811164*^9}, {3.9253915974915533`*^9, 3.925391623802257*^9}, {
   3.925392788862295*^9, 3.9253928688039303`*^9}, {3.9253930227047997`*^9, 
   3.925393096173184*^9}, {3.925393708871986*^9, 3.925393715717189*^9}, {
   3.9253937494103193`*^9, 3.92539380770921*^9}, {3.9253940410763283`*^9, 
   3.925394048157661*^9}, {3.925394218562673*^9, 3.925394270437191*^9}, {
   3.9253943364122066`*^9, 3.9253944592534313`*^9}, {3.925394506159732*^9, 
   3.925394539707529*^9}, {3.925394585846653*^9, 3.9253947217617083`*^9}, {
   3.925732625999062*^9, 3.9257326864852552`*^9}, {3.925734596082385*^9, 
   3.9257346769914703`*^9}, {3.9265510554270563`*^9, 3.9265510644355373`*^9}},
 
 CellLabel->
  "In[474]:=",ExpressionUUID->"79927019-60b3-4e35-99ba-21182e4ed03c"],

Cell[BoxData[
 GraphicsBox[{{}, {{}, {}, 
    {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
      NCache[
       Rational[7, 360], 0.019444444444444445`]], AbsoluteThickness[1.6], 
     LineBox[{{1500., 1.3000790899281465`}, {1550., 1.315302800620411}, {
      1600., 1.3314768606067289`}, {1650., 1.3485986850124496`}, {1700., 
      1.3666702767168}}]}}, {{}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{
    FormBox["\"Terminal Time (s)\"", TraditionalForm], 
    FormBox["\"Max Ratio\"", TraditionalForm]},
  AxesOrigin->{1495.8333333333335`, 1.2963795795509991`},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{
   "OptimizePlotMarkers" -> True, "OptimizePlotMarkers" -> True, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        Identity[
         Part[#, 1]], 
        Identity[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        Identity[
         Part[#, 1]], 
        Identity[
         Part[#, 2]]}& )}},
  PlotLabel->FormBox["\"Max Ratios vs Terminal Time\"", TraditionalForm],
  PlotRange->{{1495.8333333333335`, 1700.}, {1.2963795795509991`, 
   1.3666702767168}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.9265950934175787`*^9, 3.9265952152473717`*^9, 
  3.926595438273725*^9, 3.9265955084257393`*^9, 3.926684069043392*^9},
 CellLabel->
  "Out[474]=",ExpressionUUID->"52363998-c52a-4c02-837e-c2e581d99825"],

Cell[BoxData[
 GraphicsBox[{{}, {{}, {}, 
    {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
      NCache[
       Rational[7, 360], 0.019444444444444445`]], AbsoluteThickness[1.6], 
     LineBox[{{1500., 1.126565280962002}, {1550., 1.1251828701097148`}, {
      1600., 1.1237403185767365`}, {1650., 1.1223380038993862`}, {1700., 
      1.1209945965255026`}}]}}, {{}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{
    FormBox["\"Terminal Time (s)\"", TraditionalForm], 
    FormBox["\"Min Ratio\"", TraditionalForm]},
  AxesOrigin->{1495.8333333333335`, 1.1206851140568082`},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{
   "OptimizePlotMarkers" -> True, "OptimizePlotMarkers" -> True, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        Identity[
         Part[#, 1]], 
        Identity[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        Identity[
         Part[#, 1]], 
        Identity[
         Part[#, 2]]}& )}},
  PlotLabel->FormBox["\"Min Ratios vs Terminal Time\"", TraditionalForm],
  PlotRange->{{1495.8333333333335`, 1700.}, {1.1206851140568082`, 
   1.126565280962002}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.9265950934175787`*^9, 3.9265952152473717`*^9, 
  3.926595438273725*^9, 3.9265955084257393`*^9, 3.9266840690613203`*^9},
 CellLabel->
  "Out[475]=",ExpressionUUID->"8b0cfa6a-f599-4ae9-8ad8-390a0aadaf98"],

Cell[BoxData[
 GraphicsBox[{{}, {{}, {}, 
    {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
      NCache[
       Rational[7, 360], 0.019444444444444445`]], AbsoluteThickness[1.6], 
     LineBox[{{1500., 1.8983725018871254`}, {1550., 1.8714309976496424`}, {
      1600., 1.851134423723824}, {1650., 1.8301588086974332`}, {1700., 
      1.8083517399806868`}}]}}, {{}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{
    FormBox["\"Terminal Time (s)\"", TraditionalForm], 
    FormBox["\"Max Angle (Rad)\"", TraditionalForm]},
  AxesOrigin->{1495.8333333333335`, 1.8033505865414405`},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{
   "OptimizePlotMarkers" -> True, "OptimizePlotMarkers" -> True, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        Identity[
         Part[#, 1]], 
        Identity[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        Identity[
         Part[#, 1]], 
        Identity[
         Part[#, 2]]}& )}},
  PlotLabel->FormBox["\"Max Angles vs Terminal Time\"", TraditionalForm],
  PlotRange->{{1495.8333333333335`, 1700.}, {1.8033505865414405`, 
   1.8983725018871254`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.9265950934175787`*^9, 3.9265952152473717`*^9, 
  3.926595438273725*^9, 3.9265955084257393`*^9, 3.926684069079473*^9},
 CellLabel->
  "Out[476]=",ExpressionUUID->"049120fb-0e9e-4618-a403-5ae8817e81ff"],

Cell[BoxData[
 GraphicsBox[{{}, {{}, {}, 
    {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
      NCache[
       Rational[7, 360], 0.019444444444444445`]], AbsoluteThickness[1.6], 
     LineBox[{{1500., 0.9932904725652137}, {1550., 0.9866285810360867}, {
      1600., 0.9647959962510121}, {1650., 0.9422956556678738}, {1700., 
      0.9367701937414349}}]}}, {{}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{
    FormBox["\"Terminal Time (s)\"", TraditionalForm], 
    FormBox["\"Min Angle (Rad)\"", TraditionalForm]},
  AxesOrigin->{1495.8333333333335`, 0.9336301782512253},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{
   "OptimizePlotMarkers" -> True, "OptimizePlotMarkers" -> True, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        Identity[
         Part[#, 1]], 
        Identity[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        Identity[
         Part[#, 1]], 
        Identity[
         Part[#, 2]]}& )}},
  PlotLabel->FormBox["\"Min Angles vs Terminal Time\"", TraditionalForm],
  PlotRange->{{1495.8333333333335`, 1700.}, {0.9336301782512253, 
   0.9932904725652137}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.9265950934175787`*^9, 3.9265952152473717`*^9, 
  3.926595438273725*^9, 3.9265955084257393`*^9, 3.926684069096758*^9},
 CellLabel->
  "Out[477]=",ExpressionUUID->"b45db974-1449-4b0a-ba6f-228b6b03e56d"],

Cell[BoxData["\<\"maxratt.png\"\>"], "Output",
 CellChangeTimes->{3.9265950934175787`*^9, 3.9265952152473717`*^9, 
  3.926595438273725*^9, 3.9265955084257393`*^9, 3.926684069181739*^9},
 CellLabel->
  "Out[478]=",ExpressionUUID->"0f7de5c2-d415-4b06-921d-fd3c00693231"],

Cell[BoxData["\<\"minratt.png\"\>"], "Output",
 CellChangeTimes->{3.9265950934175787`*^9, 3.9265952152473717`*^9, 
  3.926595438273725*^9, 3.9265955084257393`*^9, 3.9266840692297373`*^9},
 CellLabel->
  "Out[479]=",ExpressionUUID->"75a46f05-e452-438e-aa61-dae41de09c64"],

Cell[BoxData["\<\"maxanglet.png\"\>"], "Output",
 CellChangeTimes->{3.9265950934175787`*^9, 3.9265952152473717`*^9, 
  3.926595438273725*^9, 3.9265955084257393`*^9, 3.9266840693029118`*^9},
 CellLabel->
  "Out[480]=",ExpressionUUID->"d02796a5-7068-4a55-8f49-07f789e44e64"],

Cell[BoxData["\<\"minanglet.png\"\>"], "Output",
 CellChangeTimes->{3.9265950934175787`*^9, 3.9265952152473717`*^9, 
  3.926595438273725*^9, 3.9265955084257393`*^9, 3.926684069382543*^9},
 CellLabel->
  "Out[481]=",ExpressionUUID->"d57b4af6-eb8b-469b-9a41-835d44b89610"]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"endtime", " ", "//", " ", "Clear"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"plotrrs", "[", 
   RowBox[{
   "terminaltime_", ",", " ", "quadform_", ",", " ", "thrustmax_", ",", " ", 
    "n_"}], "]"}], " ", ":=", " ", 
  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[{"ax1length", " ", "=", " ", 
      RowBox[{"Norm", "[", "ax1", "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"ax2length", " ", "=", " ", 
      RowBox[{"Norm", "[", "ax2", "]"}]}], ";", "\[IndentingNewLine]", 
     "\[IndentingNewLine]", 
     RowBox[{"testbounds", " ", "=", " ", 
      RowBox[{"thrustreachable", "[", 
       RowBox[{
       "xrf", ",", " ", "terminaltime", ",", " ", "n", ",", " ", 
        "thrustmax"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"teststats", " ", "=", " ", 
      RowBox[{"reachablestatsoperator", "[", 
       RowBox[{"quadform", ",", "testbounds", ",", " ", 
        RowBox[{"terminaltime", " ", 
         RowBox[{"thrustmax", "^", "2"}]}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", "\[IndentingNewLine]", 
     RowBox[{"plot", " ", "=", " ", 
      RowBox[{"Show", "[", 
       RowBox[{
        RowBox[{"Graphics", "[", 
         RowBox[{"GeometricTransformation", "[", 
          RowBox[{
           RowBox[{"Circle", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"0", ",", "0"}], "}"}], ",", " ", 
             RowBox[{"{", 
              RowBox[{"ax1length", ",", " ", "ax2length"}], " ", "}"}]}], 
            "]"}], ",", " ", 
           RowBox[{"RotationTransform", "[", 
            RowBox[{"ax1", " ", "//", " ", "vectoangle"}], "]"}]}], "]"}], 
         "]"}], ",", " ", 
        RowBox[{"Graphics", "[", 
         RowBox[{"Point", "[", 
          RowBox[{"testbounds", "[", 
           RowBox[{"[", 
            RowBox[{";;", ",", " ", 
             RowBox[{"1", " ", ";;", " ", "2"}]}], "]"}], "]"}], "]"}], "]"}],
         ",", " ", 
        RowBox[{"Graphics", "[", 
         RowBox[{"Arrow", "[", 
          RowBox[{"ax1length", " ", 
           RowBox[{"AnglePath", "[", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"teststats", "[", 
               RowBox[{"[", 
                RowBox[{"2", ",", " ", "1"}], "]"}], "]"}], " ", "//", " ", 
              "vectoangle"}], "}"}], "]"}]}], "]"}], "]"}], ",", " ", 
        RowBox[{"Graphics", "[", 
         RowBox[{"Arrow", "[", 
          RowBox[{"ax2length", " ", 
           RowBox[{"AnglePath", "[", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"teststats", "[", 
               RowBox[{"[", 
                RowBox[{"2", ",", " ", "2"}], "]"}], "]"}], " ", "//", " ", 
              "vectoangle"}], "}"}], "]"}]}], "]"}], "]"}], ",", " ", 
        RowBox[{"Graphics", "[", 
         RowBox[{"Text", "[", 
          RowBox[{
           RowBox[{"Style", "[", 
            RowBox[{"\"\<Maximum Ratio\>\"", ",", "10"}], "]"}], ",", 
           RowBox[{"ax1length", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"teststats", "[", 
               RowBox[{"[", 
                RowBox[{"2", ",", " ", "1"}], "]"}], "]"}], " ", "//", " ", 
              "Normalize"}], ")"}]}]}], "]"}], "]"}], ",", " ", 
        RowBox[{"Graphics", "[", 
         RowBox[{"Text", "[", 
          RowBox[{
           RowBox[{"Style", "[", 
            RowBox[{"\"\<Minimum Ratio\>\"", ",", "10"}], "]"}], ",", 
           RowBox[{"ax2length", " ", 
            RowBox[{"(", 
             RowBox[{
              RowBox[{"teststats", "[", 
               RowBox[{"[", 
                RowBox[{"2", ",", " ", "2"}], "]"}], "]"}], " ", "//", " ", 
              "Normalize"}], ")"}]}]}], "]"}], "]"}], ",", " ", 
        RowBox[{"Axes", " ", "->", " ", "True"}], ",", " ", 
        RowBox[{
        "PlotLabel", " ", "->", " ", 
         "\"\<Energy Reachable and Thrust Limited RRS\>\""}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"StringForm", "[", 
         RowBox[{"\"\<RRSttt``.png\>\"", ",", " ", "terminaltime"}], "]"}], 
        " ", "//", " ", "ToString"}], ",", " ", "plot"}], "]"}]}]}], 
   "\[IndentingNewLine]", "]"}]}]}], "Input",
 CellChangeTimes->{
  3.925229320881812*^9, 3.9252293779140663`*^9, {3.9254013374236717`*^9, 
   3.925401341755189*^9}, {3.9254013740203543`*^9, 3.9254014437877617`*^9}, {
   3.925401479897015*^9, 3.9254014833598957`*^9}, {3.925401530421446*^9, 
   3.925401740584506*^9}, {3.9254018090693903`*^9, 3.9254018294148617`*^9}, {
   3.925401934551457*^9, 3.9254020126347017`*^9}, {3.925402097888081*^9, 
   3.925402174279312*^9}, {3.9254022572805243`*^9, 3.925402266134878*^9}, {
   3.925402304002576*^9, 3.925402365190275*^9}, {3.925402519118141*^9, 
   3.925402520273685*^9}, {3.9254025550305347`*^9, 3.925402563230091*^9}, {
   3.925402595319421*^9, 3.9254026030856647`*^9}, {3.9254026544623613`*^9, 
   3.925402752665216*^9}, {3.925403124442355*^9, 3.925403127791978*^9}, {
   3.925403206423682*^9, 3.9254032548910437`*^9}, {3.925732813848874*^9, 
   3.925732823395372*^9}, {3.925732974302396*^9, 3.925732975237857*^9}, {
   3.9257330428940153`*^9, 3.925733102905794*^9}, {3.925733215568214*^9, 
   3.925733296544441*^9}, {3.9257333623010483`*^9, 3.925733377719996*^9}, {
   3.9257334383481693`*^9, 3.925733445343869*^9}, {3.9257335184801903`*^9, 
   3.9257335340493507`*^9}, {3.9257335885833063`*^9, 3.9257335887665567`*^9}, 
   3.925733763492587*^9, {3.925733811593409*^9, 3.925733814593515*^9}, {
   3.925734020839611*^9, 3.925734021536331*^9}, {3.9257340861485662`*^9, 
   3.925734112655384*^9}, 3.926550663795846*^9, 3.926551037526328*^9, {
   3.926595430311865*^9, 3.926595432133421*^9}},
 CellLabel->
  "In[482]:=",ExpressionUUID->"c099018c-97a1-4afd-94a3-22d93d786922"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.9254015616916523`*^9, 3.9254015907450657`*^9}, 
   3.925401756139267*^9, 3.925734541257496*^9, 3.92573494204535*^9},
 CellLabel->
  "In[484]:=",ExpressionUUID->"007c2eb7-74f7-497e-a06c-ad5f8bbcdc75"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"*", " ", "Saving", " ", "the", " ", "plots"}], " ", "**)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"MapThread", "[", 
   RowBox[{"plotrrs", ",", " ", 
    RowBox[{"{", 
     RowBox[{
     "terminaltimes", ",", " ", "forms", ",", " ", "tmaxs", ",", " ", "ns"}], 
     "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.922130394729796*^9, 3.922130403962502*^9}, {
   3.922180491997245*^9, 3.9221804955055523`*^9}, {3.923731873051819*^9, 
   3.923731880740302*^9}, 3.9237320117676287`*^9, {3.9237320886336603`*^9, 
   3.923732111431079*^9}, {3.923732164899686*^9, 3.923732341876013*^9}, {
   3.9237323729020157`*^9, 3.923732403403569*^9}, {3.9237324351318502`*^9, 
   3.9237325455525303`*^9}, {3.9237326033405037`*^9, 3.923732648442669*^9}, {
   3.923732724273677*^9, 3.923732766189991*^9}, {3.925222551651752*^9, 
   3.925222583598201*^9}, {3.925222628939603*^9, 3.925222631229698*^9}, {
   3.9252228245788717`*^9, 3.92522287172933*^9}, {3.925223089930551*^9, 
   3.9252233081904984`*^9}, {3.925223346030294*^9, 3.925223351813615*^9}, {
   3.925224030673959*^9, 3.925224039679594*^9}, {3.9253949284599943`*^9, 
   3.925394932338594*^9}, {3.9253949747087297`*^9, 3.925394994263081*^9}, {
   3.9253950440999403`*^9, 3.925395044699041*^9}, {3.925395153178814*^9, 
   3.925395173012898*^9}, {3.925395413572029*^9, 3.9253954467682447`*^9}, {
   3.925395673790967*^9, 3.925395737340582*^9}, {3.925395805943186*^9, 
   3.925395847660926*^9}, {3.925395927937634*^9, 3.9253959308566427`*^9}, {
   3.9253963972248096`*^9, 3.9253964599349613`*^9}, {3.925396545734774*^9, 
   3.9253966092903337`*^9}, {3.925396640587646*^9, 3.925396673097431*^9}, {
   3.92539670632568*^9, 3.925396711381776*^9}, {3.925396758933502*^9, 
   3.925396781938702*^9}, {3.925396816469982*^9, 3.925396831575831*^9}, {
   3.925396864949047*^9, 3.925396867591926*^9}, {3.925396915687825*^9, 
   3.925396966788624*^9}, 3.9253970196417828`*^9, {3.92539710177389*^9, 
   3.925397102459566*^9}, {3.9253973800783043`*^9, 3.925397427278467*^9}, 
   3.925397609042029*^9, {3.925397831370624*^9, 3.925397896095368*^9}, 
   3.925397933440997*^9, {3.925397978089177*^9, 3.9253980306253767`*^9}, {
   3.9253980606428747`*^9, 3.925398065827478*^9}, {3.9253980984519053`*^9, 
   3.925398145865192*^9}, {3.9253981816590357`*^9, 3.925398183315136*^9}, {
   3.925398250798353*^9, 3.9253982604022017`*^9}, {3.925398327278713*^9, 
   3.925398344032898*^9}, {3.9253983773474703`*^9, 3.925398397483679*^9}, {
   3.925398526485132*^9, 3.925398533802897*^9}, {3.925398623468104*^9, 
   3.925398718819539*^9}, {3.925398754751903*^9, 3.925398757109392*^9}, {
   3.9253993221663923`*^9, 3.92539936527999*^9}, {3.925399575618786*^9, 
   3.9253996108552713`*^9}, {3.925399680416231*^9, 3.925399741113248*^9}, {
   3.9253997712101717`*^9, 3.9253998198213997`*^9}, {3.9253998631630173`*^9, 
   3.925399891180518*^9}, {3.925399966873437*^9, 3.925399997246299*^9}, {
   3.925400188827908*^9, 3.92540019981277*^9}, {3.925400642266828*^9, 
   3.925400716572877*^9}, {3.925400863484913*^9, 3.925400867239274*^9}, {
   3.925400912926683*^9, 3.925401007919718*^9}, {3.925401051023543*^9, 
   3.9254010796211367`*^9}, {3.925401147955761*^9, 3.925401241325233*^9}, {
   3.925401595644754*^9, 3.925401606745082*^9}, {3.925401720597919*^9, 
   3.925401752822945*^9}, {3.925734863761654*^9, 3.925734877749889*^9}, {
   3.925734908826676*^9, 3.925734950911256*^9}},
 CellLabel->
  "In[485]:=",ExpressionUUID->"6ad10b9c-46d8-4f63-9072-e5fba7ca1875"],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate", "izero", 
   "\"Integral and error estimates are 0 on all integration subregions. Try \
increasing the value of the MinRecursion option. If value of integral may be \
0, specify a finite value for the AccuracyGoal option.\"", 2, 485, 85, 
   30967495211780212173, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.926595093789057*^9, 3.9265952160483923`*^9, 
  3.926595508785144*^9, 3.92668406944382*^9},
 CellLabel->
  "During evaluation of \
In[485]:=",ExpressionUUID->"fe48cc3f-f616-4d7e-837e-a318416c00b4"],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate", "izero", 
   "\"Integral and error estimates are 0 on all integration subregions. Try \
increasing the value of the MinRecursion option. If value of integral may be \
0, specify a finite value for the AccuracyGoal option.\"", 2, 485, 86, 
   30967495211780212173, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.926595093789057*^9, 3.9265952160483923`*^9, 
  3.926595508785144*^9, 3.926684069462262*^9},
 CellLabel->
  "During evaluation of \
In[485]:=",ExpressionUUID->"9fb7fe6b-a913-4517-8490-5502a5a362d6"],

Cell[BoxData[
 TemplateBox[{
  "NIntegrate", "izero", 
   "\"Integral and error estimates are 0 on all integration subregions. Try \
increasing the value of the MinRecursion option. If value of integral may be \
0, specify a finite value for the AccuracyGoal option.\"", 2, 485, 87, 
   30967495211780212173, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.926595093789057*^9, 3.9265952160483923`*^9, 
  3.926595508785144*^9, 3.926684069483116*^9},
 CellLabel->
  "During evaluation of \
In[485]:=",ExpressionUUID->"78388102-b816-4b63-9690-f24a2c79cf1e"],

Cell[BoxData[
 TemplateBox[{
  "General", "stop", 
   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"NIntegrate\\\", \
\\\"::\\\", \\\"izero\\\"}], \\\"MessageName\\\"]\\) will be suppressed \
during this calculation.\"", 2, 485, 88, 30967495211780212173, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.926595093789057*^9, 3.9265952160483923`*^9, 
  3.926595508785144*^9, 3.926684069486535*^9},
 CellLabel->
  "During evaluation of \
In[485]:=",ExpressionUUID->"c57b0dd1-b8e1-4b67-8037-19d0a78c46cf"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"\<\"RRSttt1500.png\"\>", ",", "\<\"RRSttt1550.png\"\>", 
   ",", "\<\"RRSttt1600.png\"\>", ",", "\<\"RRSttt1650.png\"\>", 
   ",", "\<\"RRSttt1700.png\"\>"}], "}"}]], "Output",
 CellChangeTimes->{3.9265951152843447`*^9, 3.9265952663798428`*^9, 
  3.926595557730755*^9, 3.92668411985211*^9},
 CellLabel->
  "Out[485]=",ExpressionUUID->"7d569499-535e-4d28-a510-185955a6a89e"]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.925221703501926*^9, 3.925221703502062*^9}, {
  3.925401767454959*^9, 3.925401802506412*^9}, {3.925402030555917*^9, 
  3.925402036836368*^9}, {3.9254021913999777`*^9, 3.925402203590562*^9}},
 CellLabel->
  "In[486]:=",ExpressionUUID->"c785bbf2-3350-4182-9c0e-5d65684a06da"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.925221705176982*^9, 3.9252217051771517`*^9}},
 CellLabel->
  "In[487]:=",ExpressionUUID->"8659b878-ae33-4e5f-9107-fb888126cfd7"]
},
WindowSize->{Full, Full},
WindowMargins->{{352, Automatic}, {Automatic, 50}},
DockedCells->{},
FrontEndVersion->"12.3 for Mac OS X ARM (64-bit) (July 9, 2021)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"0817f4e8-3072-4fda-8106-aa1fffaead3d"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 775, 18, 73, "Input",ExpressionUUID->"344c9343-693b-47e9-b882-096618076cd6"],
Cell[1358, 42, 257, 4, 34, "Output",ExpressionUUID->"5dc626c2-b9aa-4372-87c6-912baed454f4"]
}, Open  ]],
Cell[CellGroupData[{
Cell[1652, 51, 11332, 286, 289, "Input",ExpressionUUID->"b895b316-708f-4092-a477-0f3a66f14cfd"],
Cell[12987, 339, 1715, 27, 31, "Message",ExpressionUUID->"a6b82077-2ddc-44eb-9426-d2e8355a1a8a"]
}, Open  ]],
Cell[14717, 369, 4952, 111, 535, "Input",ExpressionUUID->"38e6d68a-01ae-469f-8372-72472f6edb77"],
Cell[19672, 482, 22049, 511, 868, "Input",ExpressionUUID->"e8b154ed-e19c-47e3-a98d-7cc40f6a965d"],
Cell[41724, 995, 16404, 360, 559, "Input",ExpressionUUID->"c4ceae39-7397-4043-ac31-479218ca2bcb"],
Cell[58131, 1357, 4156, 73, 157, "Input",ExpressionUUID->"d82cf4ae-f3bf-4a07-b427-7d6864e8cdbf"],
Cell[CellGroupData[{
Cell[62312, 1434, 4139, 66, 94, "Input",ExpressionUUID->"5bac8550-244b-4b24-a50f-7e5196fb9d83"],
Cell[66454, 1502, 700, 14, 28, "Message",ExpressionUUID->"eff9712d-9f55-4034-a2be-6f2e3e3fbfad"],
Cell[67157, 1518, 699, 14, 28, "Message",ExpressionUUID->"eac3f0c3-802f-4c77-9220-00843369023f"],
Cell[67859, 1534, 702, 14, 28, "Message",ExpressionUUID->"4f784c11-ffc5-48de-b2a7-267534e8e887"],
Cell[68564, 1550, 649, 13, 28, "Message",ExpressionUUID->"65cf75d5-3200-4fa4-ac34-291b9a9c1807"],
Cell[69216, 1565, 457, 9, 34, "Output",ExpressionUUID->"a3aaf8db-05f5-43c5-8a7c-f859e2f64fec"],
Cell[69676, 1576, 736, 19, 34, "Output",ExpressionUUID->"5d0e7b46-eaab-4ee8-8572-a4759784af16"]
}, Open  ]],
Cell[CellGroupData[{
Cell[70449, 1600, 5301, 104, 283, "Input",ExpressionUUID->"8d5d3b56-3d98-4aa2-94a6-05733121195f"],
Cell[75753, 1706, 573, 12, 28, "Message",ExpressionUUID->"4f611342-2991-4279-bd60-64ef961eabc9"],
Cell[76329, 1720, 574, 12, 28, "Message",ExpressionUUID->"347a72b1-dd69-49f6-a75e-43ff8fbb3af6"],
Cell[76906, 1734, 574, 12, 28, "Message",ExpressionUUID->"38e85ed0-d423-4e11-a186-f308288d5025"],
Cell[77483, 1748, 527, 11, 28, "Message",ExpressionUUID->"42b5eef7-7bd6-4890-a9a4-8981c303a410"]
}, Open  ]],
Cell[78025, 1762, 3616, 80, 241, "Input",ExpressionUUID->"ce85b294-c6a9-4a49-8a9c-f7900fa64e64"],
Cell[CellGroupData[{
Cell[81666, 1846, 4250, 95, 241, "Input",ExpressionUUID->"79927019-60b3-4e35-99ba-21182e4ed03c"],
Cell[85919, 1943, 1793, 45, 217, "Output",ExpressionUUID->"52363998-c52a-4c02-837e-c2e581d99825"],
Cell[87715, 1990, 1801, 45, 214, "Output",ExpressionUUID->"8b0cfa6a-f599-4ae9-8ad8-390a0aadaf98"],
Cell[89519, 2037, 1807, 45, 210, "Output",ExpressionUUID->"049120fb-0e9e-4618-a403-5ae8817e81ff"],
Cell[91329, 2084, 1801, 45, 211, "Output",ExpressionUUID->"b45db974-1449-4b0a-ba6f-228b6b03e56d"],
Cell[93133, 2131, 268, 4, 34, "Output",ExpressionUUID->"0f7de5c2-d415-4b06-921d-fd3c00693231"],
Cell[93404, 2137, 270, 4, 34, "Output",ExpressionUUID->"75a46f05-e452-438e-aa61-dae41de09c64"],
Cell[93677, 2143, 272, 4, 34, "Output",ExpressionUUID->"d02796a5-7068-4a55-8f49-07f789e44e64"],
Cell[93952, 2149, 270, 4, 34, "Output",ExpressionUUID->"d57b4af6-eb8b-469b-9a41-835d44b89610"]
}, Open  ]],
Cell[94237, 2156, 10371, 245, 493, "Input",ExpressionUUID->"c099018c-97a1-4afd-94a3-22d93d786922"],
Cell[104611, 2403, 248, 4, 30, "Input",ExpressionUUID->"007c2eb7-74f7-497e-a06c-ad5f8bbcdc75"],
Cell[CellGroupData[{
Cell[104884, 2411, 3518, 54, 52, "Input",ExpressionUUID->"6ad10b9c-46d8-4f63-9072-e5fba7ca1875"],
Cell[108405, 2467, 577, 12, 28, "Message",ExpressionUUID->"fe48cc3f-f616-4d7e-837e-a318416c00b4"],
Cell[108985, 2481, 578, 12, 28, "Message",ExpressionUUID->"9fb7fe6b-a913-4517-8490-5502a5a362d6"],
Cell[109566, 2495, 578, 12, 28, "Message",ExpressionUUID->"78388102-b816-4b63-9690-f24a2c79cf1e"],
Cell[110147, 2509, 529, 11, 28, "Message",ExpressionUUID->"c57b0dd1-b8e1-4b67-8037-19d0a78c46cf"],
Cell[110679, 2522, 413, 8, 34, "Output",ExpressionUUID->"7d569499-535e-4d28-a510-185955a6a89e"]
}, Open  ]],
Cell[111107, 2533, 321, 5, 30, "Input",ExpressionUUID->"c785bbf2-3350-4182-9c0e-5d65684a06da"],
Cell[111431, 2540, 177, 3, 30, "Input",ExpressionUUID->"8659b878-ae33-4e5f-9107-fb888126cfd7"]
}
]
*)