energythrustcomparisondebugging_ratio.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[    162661,       3572]
NotebookOptionsPosition[    157265,       3490]
NotebookOutlinePosition[    157723,       3508]
CellTagsIndexPosition[    157680,       3505]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData["\[AliasDelimiter]"], "Input",
 CellChangeTimes->{3.9371572956486073`*^9},
 CellLabel->
  "In[378]:=",ExpressionUUID->"688a2e9f-7ca3-48dd-ac22-f5d7415eb864"],

Cell[BoxData["\[AliasDelimiter]"], "Output",
 CellChangeTimes->{3.93715826631149*^9},
 CellLabel->
  "Out[378]=",ExpressionUUID->"bf6b0adc-9856-4441-b15f-85fb2df8cafb"]
}, Open  ]],

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}, {3.926597433337858*^9, 
   3.926597444998395*^9}, {3.926689030305884*^9, 3.926689032715295*^9}, 
   3.9266891199795103`*^9},
 CellLabel->
  "In[379]:=",ExpressionUUID->"344c9343-693b-47e9-b882-096618076cd6"],

Cell[BoxData["\<\"/Users/maxzweig/energythrust\"\>"], "Output",
 CellChangeTimes->{3.937158266334857*^9},
 CellLabel->
  "Out[379]=",ExpressionUUID->"c8f475e4-98de-4808-8c45-b977a1e68ae0"]
}, Open  ]],

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[380]:=",ExpressionUUID->"b895b316-708f-4092-a477-0f3a66f14cfd"],

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_", ",", 
     "numpoints_"}], "]"}], " ", ":=", " ", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "points", ",", " ", "tfs", ",", " ", "hs", ",", " ", "As", ",", " ", 
       "ns", ",", " ", "Rrs", ",", " ", "pointsT", ",", " ", "rd", ",", " ", 
       "rdpoints", ",", " ", "rdboundary", ",", " ", "first2"}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"points", " ", "=", " ", 
       RowBox[{"CirclePoints", "[", "numpoints", "]"}]}], ";", 
      "\[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", ",", " ", "numpoints"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"tfs", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"hs", " ", "=", " ", 
       RowBox[{"ConstantArray", "[", 
        RowBox[{"1", ",", " ", "numpoints"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"hs", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"ns", " ", "=", " ", 
       RowBox[{"ConstantArray", "[", 
        RowBox[{"na", ",", " ", "numpoints"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"ns", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"Rrs", " ", "=", " ", 
       RowBox[{"ConstantArray", "[", 
        RowBox[{"1", ",", " ", "numpoints"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Rrs", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"emaxs", " ", "=", " ", 
       RowBox[{"ConstantArray", "[", 
        RowBox[{"emax", ",", " ", "numpoints"}], "]"}]}], ";", 
      "\[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}, {3.926598537937785*^9, 
   3.926598572651369*^9}},
 CellLabel->
  "In[382]:=",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_", ",", " ", "h_"}], "]"}], " ", "=", " ", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{"72", " ", "t", " ", 
              RowBox[{"Cos", "[", 
               RowBox[{"h", " ", "t"}], "]"}]}], "+", 
             FractionBox[
              RowBox[{
               RowBox[{"154", " ", "h", " ", "t"}], "+", 
               RowBox[{"48", " ", 
                SuperscriptBox["h", "3"], " ", 
                SuperscriptBox["t", "3"]}], "-", 
               RowBox[{"384", " ", 
                RowBox[{"Sin", "[", 
                 RowBox[{"h", " ", "t"}], "]"}]}], "-", 
               RowBox[{"27", " ", 
                RowBox[{"Sin", "[", 
                 RowBox[{"2", " ", "h", " ", "t"}], "]"}]}]}], 
              RowBox[{"4", " ", "h"}]]}], ",", 
            RowBox[{
             RowBox[{
              RowBox[{"-", "3"}], " ", "h", " ", 
              SuperscriptBox["t", "2"]}], "+", 
             FractionBox[
              RowBox[{"6", "-", 
               RowBox[{"6", " ", 
                RowBox[{"Cos", "[", 
                 RowBox[{"h", " ", "t"}], "]"}]}]}], "h"]}], ",", "0", ",", 
            FractionBox[
             RowBox[{"8", "+", 
              RowBox[{"12", " ", 
               SuperscriptBox["h", "2"], " ", 
               SuperscriptBox["t", "2"]}], "-", 
              RowBox[{"8", " ", 
               RowBox[{"Cos", "[", 
                RowBox[{"h", " ", "t"}], "]"}]}], "-", 
              RowBox[{"24", " ", "h", " ", "t", " ", 
               RowBox[{"Sin", "[", 
                RowBox[{"h", " ", "t"}], "]"}]}], "+", 
              RowBox[{"9", " ", 
               SuperscriptBox[
                RowBox[{"Sin", "[", 
                 RowBox[{"h", " ", "t"}], "]"}], "2"]}]}], 
             RowBox[{"2", " ", 
              SuperscriptBox["h", "2"]}]], ",", 
            FractionBox[
             RowBox[{
              RowBox[{"23", " ", "h", " ", "t"}], "+", 
              RowBox[{"6", " ", 
               SuperscriptBox["h", "3"], " ", 
               SuperscriptBox["t", "3"]}], "+", 
              RowBox[{"42", " ", "h", " ", "t", " ", 
               RowBox[{"Cos", "[", 
                RowBox[{"h", " ", "t"}], "]"}]}], "-", 
              RowBox[{"56", " ", 
               RowBox[{"Sin", "[", 
                RowBox[{"h", " ", "t"}], "]"}]}], "-", 
              RowBox[{
               FractionBox["9", "2"], " ", 
               RowBox[{"Sin", "[", 
                RowBox[{"2", " ", "h", " ", "t"}], "]"}]}]}], 
             SuperscriptBox["h", "2"]], ",", "0"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{
              RowBox[{"-", "3"}], " ", "h", " ", 
              SuperscriptBox["t", "2"]}], "+", 
             FractionBox[
              RowBox[{"6", "-", 
               RowBox[{"6", " ", 
                RowBox[{"Cos", "[", 
                 RowBox[{"h", " ", "t"}], "]"}]}]}], "h"]}], ",", "t", ",", 
            "0", ",", 
            FractionBox[
             RowBox[{"2", " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                 RowBox[{"-", "h"}], " ", "t"}], "+", 
                RowBox[{"Sin", "[", 
                 RowBox[{"h", " ", "t"}], "]"}]}], ")"}]}], 
             SuperscriptBox["h", "2"]], ",", 
            RowBox[{
             RowBox[{"-", 
              FractionBox[
               RowBox[{"3", " ", 
                SuperscriptBox["t", "2"]}], "2"]}], "+", 
             FractionBox[
              RowBox[{"4", "-", 
               RowBox[{"4", " ", 
                RowBox[{"Cos", "[", 
                 RowBox[{"h", " ", "t"}], "]"}]}]}], 
              SuperscriptBox["h", "2"]]}], ",", "0"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"0", ",", "0", ",", 
            FractionBox[
             RowBox[{
              RowBox[{"2", " ", "h", " ", "t"}], "+", 
              RowBox[{"Sin", "[", 
               RowBox[{"2", " ", "h", " ", "t"}], "]"}]}], 
             RowBox[{"4", " ", "h"}]], ",", "0", ",", "0", ",", 
            FractionBox[
             SuperscriptBox[
              RowBox[{"Sin", "[", 
               RowBox[{"h", " ", "t"}], "]"}], "2"], 
             RowBox[{"2", " ", 
              SuperscriptBox["h", "2"]}]]}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{
            FractionBox[
             RowBox[{"8", "+", 
              RowBox[{"12", " ", 
               SuperscriptBox["h", "2"], " ", 
               SuperscriptBox["t", "2"]}], "-", 
              RowBox[{"8", " ", 
               RowBox[{"Cos", "[", 
                RowBox[{"h", " ", "t"}], "]"}]}], "-", 
              RowBox[{"24", " ", "h", " ", "t", " ", 
               RowBox[{"Sin", "[", 
                RowBox[{"h", " ", "t"}], "]"}]}], "+", 
              RowBox[{"9", " ", 
               SuperscriptBox[
                RowBox[{"Sin", "[", 
                 RowBox[{"h", " ", "t"}], "]"}], "2"]}]}], 
             RowBox[{"2", " ", 
              SuperscriptBox["h", "2"]}]], ",", 
            FractionBox[
             RowBox[{"2", " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                 RowBox[{"-", "h"}], " ", "t"}], "+", 
                RowBox[{"Sin", "[", 
                 RowBox[{"h", " ", "t"}], "]"}]}], ")"}]}], 
             SuperscriptBox["h", "2"]], ",", "0", ",", 
            FractionBox[
             RowBox[{
              RowBox[{"26", " ", "h", " ", "t"}], "-", 
              RowBox[{"32", " ", 
               RowBox[{"Sin", "[", 
                RowBox[{"h", " ", "t"}], "]"}]}], "+", 
              RowBox[{"3", " ", 
               RowBox[{"Sin", "[", 
                RowBox[{"2", " ", "h", " ", "t"}], "]"}]}]}], 
             RowBox[{"4", " ", 
              SuperscriptBox["h", "3"]}]], ",", 
            FractionBox[
             RowBox[{"3", " ", 
              SuperscriptBox[
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "h"}], " ", "t"}], "+", 
                 RowBox[{"Sin", "[", 
                  RowBox[{"h", " ", "t"}], "]"}]}], ")"}], "2"]}], 
             SuperscriptBox["h", "3"]], ",", "0"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{
            FractionBox[
             RowBox[{
              RowBox[{"23", " ", "h", " ", "t"}], "+", 
              RowBox[{"6", " ", 
               SuperscriptBox["h", "3"], " ", 
               SuperscriptBox["t", "3"]}], "+", 
              RowBox[{"42", " ", "h", " ", "t", " ", 
               RowBox[{"Cos", "[", 
                RowBox[{"h", " ", "t"}], "]"}]}], "-", 
              RowBox[{"56", " ", 
               RowBox[{"Sin", "[", 
                RowBox[{"h", " ", "t"}], "]"}]}], "-", 
              RowBox[{
               FractionBox["9", "2"], " ", 
               RowBox[{"Sin", "[", 
                RowBox[{"2", " ", "h", " ", "t"}], "]"}]}]}], 
             SuperscriptBox["h", "2"]], ",", 
            RowBox[{
             RowBox[{"-", 
              FractionBox[
               RowBox[{"3", " ", 
                SuperscriptBox["t", "2"]}], "2"]}], "+", 
             FractionBox[
              RowBox[{"4", "-", 
               RowBox[{"4", " ", 
                RowBox[{"Cos", "[", 
                 RowBox[{"h", " ", "t"}], "]"}]}]}], 
              SuperscriptBox["h", "2"]]}], ",", "0", ",", 
            FractionBox[
             RowBox[{"3", " ", 
              SuperscriptBox[
               RowBox[{"(", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "h"}], " ", "t"}], "+", 
                 RowBox[{"Sin", "[", 
                  RowBox[{"h", " ", "t"}], "]"}]}], ")"}], "2"]}], 
             SuperscriptBox["h", "3"]], ",", 
            FractionBox[
             RowBox[{
              RowBox[{"14", " ", "h", " ", "t"}], "+", 
              RowBox[{"3", " ", 
               SuperscriptBox["h", "3"], " ", 
               SuperscriptBox["t", "3"]}], "+", 
              RowBox[{"24", " ", "h", " ", "t", " ", 
               RowBox[{"Cos", "[", 
                RowBox[{"h", " ", "t"}], "]"}]}], "-", 
              RowBox[{"32", " ", 
               RowBox[{"Sin", "[", 
                RowBox[{"h", " ", "t"}], "]"}]}], "-", 
              RowBox[{"3", " ", 
               RowBox[{"Sin", "[", 
                RowBox[{"2", " ", "h", " ", "t"}], "]"}]}]}], 
             SuperscriptBox["h", "3"]], ",", "0"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"0", ",", "0", ",", 
            FractionBox[
             SuperscriptBox[
              RowBox[{"Sin", "[", 
               RowBox[{"h", " ", "t"}], "]"}], "2"], 
             RowBox[{"2", " ", 
              SuperscriptBox["h", "2"]}]], ",", "0", ",", "0", ",", 
            RowBox[{"-", 
             FractionBox[
              RowBox[{
               RowBox[{
                RowBox[{"-", "2"}], " ", "h", " ", "t"}], "+", 
               RowBox[{"Sin", "[", 
                RowBox[{"2", " ", "h", " ", "t"}], "]"}]}], 
              RowBox[{"4", " ", 
               SuperscriptBox["h", "3"]}]]}]}], "}"}]}], "}"}]}], ";", 
       RowBox[{"(*", "Good", "*)"}], "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"r", "[", 
         RowBox[{"t", ",", " ", "n"}], "]"}], " ", "//", " ", "MatrixForm"}], 
       ";", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"sig", "[", 
         RowBox[{"t_", ",", " ", "n_"}], "]"}], " ", "=", " ", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"4", "-", 
             RowBox[{"3", " ", 
              RowBox[{"Cos", "[", 
               RowBox[{"n", " ", "t"}], "]"}]}]}], ",", "0", ",", "0", ",", 
            FractionBox[
             RowBox[{"Sin", "[", 
              RowBox[{"n", " ", "t"}], "]"}], "n"], ",", 
            FractionBox[
             RowBox[{"2", "-", 
              RowBox[{"2", " ", 
               RowBox[{"Cos", "[", 
                RowBox[{"n", " ", "t"}], "]"}]}]}], "n"], ",", "0"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"6", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{
                RowBox[{"-", "n"}], " ", "t"}], "+", 
               RowBox[{"Sin", "[", 
                RowBox[{"n", " ", "t"}], "]"}]}], ")"}]}], ",", "1", ",", "0",
             ",", 
            FractionBox[
             RowBox[{"2", " ", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"-", "1"}], "+", 
                RowBox[{"Cos", "[", 
                 RowBox[{"n", " ", "t"}], "]"}]}], ")"}]}], "n"], ",", 
            RowBox[{
             RowBox[{
              RowBox[{"-", "3"}], " ", "t"}], "+", 
             FractionBox[
              RowBox[{"4", " ", 
               RowBox[{"Sin", "[", 
                RowBox[{"n", " ", "t"}], "]"}]}], "n"]}], ",", "0"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"0", ",", "0", ",", 
            RowBox[{"Cos", "[", 
             RowBox[{"n", " ", "t"}], "]"}], ",", "0", ",", "0", ",", 
            FractionBox[
             RowBox[{"Sin", "[", 
              RowBox[{"n", " ", "t"}], "]"}], "n"]}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"3", " ", "n", " ", 
             RowBox[{"Sin", "[", 
              RowBox[{"n", " ", "t"}], "]"}]}], ",", "0", ",", "0", ",", 
            RowBox[{"Cos", "[", 
             RowBox[{"n", " ", "t"}], "]"}], ",", 
            RowBox[{"2", " ", 
             RowBox[{"Sin", "[", 
              RowBox[{"n", " ", "t"}], "]"}]}], ",", "0"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"6", " ", "n", " ", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"-", "1"}], "+", 
               RowBox[{"Cos", "[", 
                RowBox[{"n", " ", "t"}], "]"}]}], ")"}]}], ",", "0", ",", "0",
             ",", 
            RowBox[{
             RowBox[{"-", "2"}], " ", 
             RowBox[{"Sin", "[", 
              RowBox[{"n", " ", "t"}], "]"}]}], ",", 
            RowBox[{
             RowBox[{"-", "3"}], "+", 
             RowBox[{"4", " ", 
              RowBox[{"Cos", "[", 
               RowBox[{"n", " ", "t"}], "]"}]}]}], ",", "0"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"0", ",", "0", ",", 
            RowBox[{
             RowBox[{"-", "n"}], " ", 
             RowBox[{"Sin", "[", 
              RowBox[{"n", " ", "t"}], "]"}]}], ",", "0", ",", "0", ",", 
            RowBox[{"Cos", "[", 
             RowBox[{"n", " ", "t"}], "]"}]}], "}"}]}], "}"}]}], ";", " ", 
       RowBox[{"(*", " ", "Good", " ", "*)"}], "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"gamma", "[", 
         RowBox[{"t_", ",", " ", "t0_", ",", " ", "n_"}], "]"}], " ", "=", 
        " ", 
        RowBox[{
         RowBox[{"sig", "[", 
          RowBox[{"t", ",", " ", "n"}], "]"}], ".", " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"sig", "[", 
            RowBox[{"t0", ",", " ", "n"}], "]"}], " ", "//", " ", "Inverse"}],
           ")"}]}]}], ";", "\[IndentingNewLine]", "\[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, {3.9271430218101254`*^9, 3.927143022174961*^9}, {
   3.927145449011572*^9, 3.9271454536186132`*^9}},
 CellLabel->
  "In[383]:=",ExpressionUUID->"e8b154ed-e19c-47e3-a98d-7cc40f6a965d"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"reachablestatsoperator", "[", 
    RowBox[{"points_", ",", "axis1_", ",", "axis2_"}], "]"}], " ", ":=", " ", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"rats", ",", " ", "vecs"}], "}"}], ",", " ", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{
       RowBox[{"ellipseVecFromUnitVec", "[", 
        RowBox[{"uv_", ",", "ax1_", ",", " ", "ax2_"}], "]"}], ":=", 
       RowBox[{
        RowBox[{"1", "/", 
         RowBox[{"Sqrt", "[", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"Normalize", "[", "ax1", "]"}], ".", "uv"}], ")"}], "^",
              "2"}], "/", 
            RowBox[{
             RowBox[{"Norm", "[", "ax1", "]"}], "^", "2"}]}], "+", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"Normalize", "[", "ax2", "]"}], ".", "uv"}], ")"}], "^",
              "2"}], "/", 
            RowBox[{
             RowBox[{"Norm", "[", "ax2", "]"}], "^", "2"}]}]}], "]"}]}], "*", 
        "uv"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"ratio", "[", "pos_", "]"}], ":=", " ", 
       RowBox[{
        RowBox[{"Norm", "[", 
         RowBox[{"ellipseVecFromUnitVec", "[", 
          RowBox[{
           RowBox[{"Normalize", "[", 
            RowBox[{"pos", "[", 
             RowBox[{"[", 
              RowBox[{"1", ";;", "2"}], "]"}], "]"}], "]"}], ",", "axis1", ",",
            " ", "axis2"}], "]"}], " ", "]"}], "/", " ", 
        RowBox[{"Norm", "[", 
         RowBox[{"pos", "[", 
          RowBox[{"[", 
           RowBox[{"1", ";;", "2"}], "]"}], "]"}], "]"}]}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"maxminrat", "[", " ", "v1s_", "]"}], " ", ":=", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{"ratio", "/@", "v1s"}], ")"}], "//", " ", "Max"}], ",", 
         " ", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"ratio", "/@", "v1s"}], ")"}], "//", " ", "Min"}]}], 
        "}"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"toq2q1", "[", "vec_", "]"}], " ", ":=", " ", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"vec", "[", 
              RowBox[{"[", "1", "]"}], "]"}], " ", "<", " ", "0"}], " ", "&&",
             " ", 
            RowBox[{
             RowBox[{"vec", "[", 
              RowBox[{"[", "2", "]"}], "]"}], " ", "<", " ", "0"}]}], ")"}], 
          " ", "||", " ", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{
             RowBox[{"vec", "[", 
              RowBox[{"[", "1", "]"}], "]"}], " ", ">", "0"}], " ", "&&", " ", 
            RowBox[{
             RowBox[{"vec", "[", 
              RowBox[{"[", "2", "]"}], "]"}], " ", "<", " ", "0"}]}], ")"}]}],
          ",", " ", 
         RowBox[{"-", "vec"}], ",", " ", "vec"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"maxminvec", "[", " ", "v1s_", "]"}], " ", ":=", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{"v1s", "[", 
           RowBox[{"[", " ", 
            RowBox[{
             RowBox[{"Position", "[", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"ratio", "/@", "v1s"}], ")"}], ",", " ", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"ratio", "/@", "v1s"}], ")"}], " ", "//", " ", 
                "Max"}]}], "    ", "]"}], "[", 
             RowBox[{"[", 
              RowBox[{"1", ",", "1"}], "]"}], "]"}], "]"}], "]"}], " ", "//", 
          " ", "toq2q1"}], " ", ",", " ", 
         RowBox[{
          RowBox[{"v1s", "[", 
           RowBox[{"[", " ", 
            RowBox[{
             RowBox[{"Position", "[", 
              RowBox[{
               RowBox[{"(", 
                RowBox[{"ratio", "/@", "v1s"}], ")"}], ",", " ", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"ratio", "/@", "v1s"}], ")"}], " ", "//", " ", 
                "Min"}]}], "   ", "]"}], "[", 
             RowBox[{"[", 
              RowBox[{"1", ",", "1"}], "]"}], "]"}], "]"}], "]"}], " ", "//", 
          " ", "toq2q1"}]}], "}"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"rats", " ", "=", " ", 
       RowBox[{"maxminrat", "[", " ", "points", "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"vecs", " ", "=", " ", 
       RowBox[{"maxminvec", "[", " ", "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]"}]], "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.926684011548202*^9, 3.9266840393723965`*^9}, {3.9266841522493362`*^9, 
   3.9266841810590825`*^9}, {3.9266842131791306`*^9, 
   3.9266842547767754`*^9}, {3.9266853734777565`*^9, 3.926685490163378*^9}, {
   3.9266855493388186`*^9, 3.9266855528594637`*^9}, {3.927040824709537*^9, 
   3.9270409233106833`*^9}, {3.927041020179716*^9, 3.927041028909778*^9}, {
   3.9270413163353453`*^9, 3.927041518469928*^9}, {3.927041694006281*^9, 
   3.927041701705345*^9}, {3.927041753282625*^9, 3.927041791423614*^9}, {
   3.927041924969586*^9, 3.927041950906375*^9}, {3.927042111178555*^9, 
   3.9270421542033167`*^9}, 3.927042190719461*^9, {3.927125450144771*^9, 
   3.927125564310648*^9}, {3.9271259672731857`*^9, 3.927125977009632*^9}, 
   3.927126162049254*^9, {3.92712646017297*^9, 3.927126466076231*^9}, 
   3.927197482099365*^9, 3.9271976975801487`*^9, 3.927198272170766*^9, {
   3.927283043326962*^9, 3.927283049838214*^9}, {3.927283080631194*^9, 
   3.9272831488858013`*^9}, {3.927283443333994*^9, 3.927283459169217*^9}, {
   3.927283883848501*^9, 3.927284029456459*^9}, {3.927285690310808*^9, 
   3.927285700849557*^9}, {3.9272857973328323`*^9, 3.9272857990263987`*^9}, {
   3.927286372195821*^9, 3.9272863740914927`*^9}, {3.9272900489768553`*^9, 
   3.9272900611472473`*^9}, {3.927290257761067*^9, 3.927290299442458*^9}, {
   3.927290409286565*^9, 3.927290410906271*^9}, {3.927290858858644*^9, 
   3.927290864614212*^9}, {3.927291221444138*^9, 3.927291243150619*^9}, {
   3.928612698997448*^9, 3.928612741597005*^9}, {3.928686456005248*^9, 
   3.928686457705928*^9}, {3.9370848755315228`*^9, 3.9370848761810703`*^9}, {
   3.941923575763579*^9, 
   3.941923576226108*^9}},ExpressionUUID->"c4ceae39-7397-4043-ac31-\
479218ca2bcb"],

Cell[CellGroupData[{

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]", "n", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"tmax", "=", 
     RowBox[{"2", "/", "1000"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"comparableenergy", " ", "=", " ", 
     RowBox[{"endtime", " ", 
      RowBox[{"tmax", "^", "2"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"2", 
    RowBox[{"Pi", "/", "n"}]}]}]}]], "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}, {3.9265987875782557`*^9, 3.9265987950193143`*^9}, 
   3.927045210645069*^9, 3.927045974568512*^9, {3.927288023019156*^9, 
   3.9272880292292547`*^9}, {3.928062902574766*^9, 3.928062903277153*^9}, {
   3.928064330718052*^9, 3.928064331218998*^9}, 3.928173579546324*^9, {
   3.9292836992637053`*^9, 3.929283699555962*^9}},
 CellLabel->
  "In[385]:=",ExpressionUUID->"d82cf4ae-f3bf-4a07-b427-7d6864e8cdbf"],

Cell[BoxData["0.0009200242014178264`"], "Output",
 CellChangeTimes->{3.9371582663607597`*^9},
 CellLabel->
  "Out[389]=",ExpressionUUID->"1a121e2d-c5a7-4194-a9e3-db1ccd73b003"],

Cell[BoxData["6829.369594296242`"], "Output",
 CellChangeTimes->{3.93715826636172*^9},
 CellLabel->
  "Out[392]=",ExpressionUUID->"921cd2a4-14e5-4da2-aab2-8f6af6b1b8de"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "*", " ", "test", " ", "of", " ", "comparison", " ", "statistics"}], " ", 
   "**)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"testbounds", " ", "=", " ", 
     RowBox[{"thrustreachable", "[", 
      RowBox[{
      "xrf", ",", " ", "endtime", ",", " ", "n", ",", " ", "tmax", ",", 
       "numsamples"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"testbounds", "[", 
    RowBox[{"[", "1", "]"}], "]"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     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.926598594426783*^9, 3.92659859677666*^9}, 
   3.926598817678035*^9, {3.9270412460335903`*^9, 3.927041250194068*^9}},
 CellLabel->
  "In[393]:=",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, 393, 146, 
   31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.937158266375347*^9},
 CellLabel->
  "During evaluation of \
In[393]:=",ExpressionUUID->"7fa8c881-5945-4648-8706-f0ab69c4167f"],

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, 393, 147, 
   31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.937158266388566*^9},
 CellLabel->
  "During evaluation of \
In[393]:=",ExpressionUUID->"e239c939-241c-451b-a455-84c7eba290ed"],

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, 393, 148, 
   31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.937158266393214*^9},
 CellLabel->
  "During evaluation of \
In[393]:=",ExpressionUUID->"ddc48787-99d5-4410-adf4-c15e469a8fb5"],

Cell[BoxData[
 TemplateBox[{
  "General", "stop", 
   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"NIntegrate\\\", \
\\\"::\\\", \\\"izero\\\"}], \\\"MessageName\\\"]\\) will be suppressed \
during this calculation.\"", 2, 393, 149, 31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.9371582663978043`*^9},
 CellLabel->
  "During evaluation of \
In[393]:=",ExpressionUUID->"6441a518-3f1b-4554-9b0f-ce650c807471"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1660.5911298629444`", ",", 
   RowBox[{"-", "2186.5496343518266`"}], ",", "0.`"}], "}"}]], "Output",
 CellChangeTimes->{3.93715826645767*^9},
 CellLabel->
  "Out[394]=",ExpressionUUID->"cb9fc03e-29f5-4456-9028-92af7216e525"]
}, 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[{"numsamples", "=", "10"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"starts", " ", "//", " ", "Clear"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"finaltime", " ", "//", " ", "Clear"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"starts", " ", "=", " ", "60"}], ";"}], " ", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"finaltime", " ", "=", " ", "6829"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"terminaltimes", " ", "=", " ", 
     RowBox[{"Range", "[", 
      RowBox[{"starts", ",", " ", "finaltime", ",", " ", "60"}], "]"}]}], 
    ";"}], "\[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[{"numpoints", "=", 
    RowBox[{"ConstantArray", "[", 
     RowBox[{"numsamples", ",", 
      RowBox[{"terminaltimes", " ", "//", " ", "Length"}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"bounds", " ", "=", " ", 
     RowBox[{"MapThread", "[", 
      RowBox[{"thrustreachable", ",", " ", 
       RowBox[{"{", 
        RowBox[{
        "xrfs", ",", " ", "terminaltimes", ",", " ", "ns", ",", " ", "tmaxs", 
         ",", "numpoints"}], "}"}]}], "]"}]}], ";"}], "\[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}, {3.9265987251142817`*^9, 
   3.9265987503874435`*^9}, 3.9265988191386824`*^9, {3.9265989024588385`*^9, 
   3.926598941579545*^9}, {3.9265993152642*^9, 3.926599315423295*^9}, {
   3.926602454579197*^9, 3.9266024557932014`*^9}, {3.9266028687714787`*^9, 
   3.9266028694403257`*^9}, {3.927140588986698*^9, 3.927140592128375*^9}, {
   3.928062996667963*^9, 3.928062997992062*^9}, {3.928173582154059*^9, 
   3.9281735904305553`*^9}, {3.9281765453279448`*^9, 3.928176546288322*^9}, {
   3.9281921153947973`*^9, 3.928192116998913*^9}, {3.928606137852817*^9, 
   3.928606137987651*^9}, {3.9286236118921747`*^9, 3.9286236120767927`*^9}, {
   3.929284149903611*^9, 3.9292841500568647`*^9}, 3.929284445251458*^9, 
   3.929284614924912*^9, {3.929483959206221*^9, 3.929483959291903*^9}, {
   3.929832976067952*^9, 3.929832976403141*^9}, {3.9298332492573833`*^9, 
   3.929833250543584*^9}, {3.930227370433585*^9, 3.930227370585516*^9}, 
   3.930230296071521*^9, 3.9305156091435204`*^9, {3.93051793721502*^9, 
   3.930517937351413*^9}, {3.930518107547576*^9, 3.930518107667583*^9}, {
   3.930518597689735*^9, 3.930518597806704*^9}, {3.937158126056772*^9, 
   3.937158126630921*^9}, {3.94138686382432*^9, 
   3.941386865590104*^9}},ExpressionUUID->"8d5d3b56-3d98-4aa2-94a6-\
05733121195f"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",", "10", ",",
    "10"}], "}"}]], "Output",
 CellChangeTimes->{3.9371582664605083`*^9},
 CellLabel->
  "Out[404]=",ExpressionUUID->"4f4fe017-2e4f-427e-aa32-efb1184e1464"],

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, 405, 150, 
   31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.937158266470469*^9},
 CellLabel->
  "During evaluation of \
In[395]:=",ExpressionUUID->"67149ad9-5abe-4451-8c8f-de3b1e27ce8c"],

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, 405, 151, 
   31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.937158266481818*^9},
 CellLabel->
  "During evaluation of \
In[395]:=",ExpressionUUID->"b6a72f8f-877e-4b71-b1eb-cb8bce01e97f"],

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, 405, 152, 
   31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.937158266494276*^9},
 CellLabel->
  "During evaluation of \
In[395]:=",ExpressionUUID->"4eec3367-8fc3-45c8-89b6-79bb53979300"],

Cell[BoxData[
 TemplateBox[{
  "General", "stop", 
   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"NIntegrate\\\", \
\\\"::\\\", \\\"izero\\\"}], \\\"MessageName\\\"]\\) will be suppressed \
during this calculation.\"", 2, 405, 153, 31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.937158266497589*^9},
 CellLabel->
  "During evaluation of \
In[395]:=",ExpressionUUID->"de6d00bf-0dd6-45c6-8ff3-bb644d13a8f7"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"getAxes", "[", 
    RowBox[{"terminaltime_", ",", " ", "quadform_", ",", " ", "thrustmax_"}], 
    "]"}], " ", ":=", " ", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"testbounds", ",", " ", "teststats"}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{
       RowBox[{"energysetaxes", "[", "quad_", "]"}], " ", ":=", 
       "\[IndentingNewLine]", 
       RowBox[{"Eigenvectors", "[", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Join", "[", 
           RowBox[{
            RowBox[{"Join", "[", 
             RowBox[{
              RowBox[{"IdentityMatrix", "[", "2", "]"}], ",", " ", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"{", 
                 RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", 
                RowBox[{"{", 
                 RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}]}], "}"}], 
              ",", " ", "2"}], "]"}], ",", " ", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
               "0", ",", "0", ",", "0", ",", " ", "0", ",", "0", ",", "0"}], 
               "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
               "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}],
               ",", 
              RowBox[{"{", 
               RowBox[{
               "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}],
               ",", 
              RowBox[{"{", 
               RowBox[{
               "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], 
               "}"}]}], "}"}]}], "]"}], ",", " ", 
          RowBox[{"quad", "[", 
           RowBox[{"[", 
            RowBox[{
             RowBox[{"7", ";;", " ", "12"}], ",", " ", 
             RowBox[{"7", " ", ";;", " ", "12"}]}], "]"}], "]"}]}], "}"}], 
        "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"energyseteigenvalues", "[", "quad_", "]"}], " ", ":=", "  ", 
       RowBox[{"Eigenvalues", "[", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Join", "[", 
           RowBox[{
            RowBox[{"Join", "[", 
             RowBox[{
              RowBox[{"IdentityMatrix", "[", "2", "]"}], ",", " ", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"{", 
                 RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", 
                RowBox[{"{", 
                 RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}]}], "}"}], 
              ",", " ", "2"}], "]"}], ",", " ", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
               "0", ",", "0", ",", "0", ",", " ", "0", ",", "0", ",", "0"}], 
               "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
               "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}],
               ",", 
              RowBox[{"{", 
               RowBox[{
               "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}],
               ",", 
              RowBox[{"{", 
               RowBox[{
               "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], 
               "}"}]}], "}"}]}], "]"}], ",", " ", 
          RowBox[{"quad", "[", 
           RowBox[{"[", 
            RowBox[{
             RowBox[{"7", ";;", " ", "12"}], ",", " ", 
             RowBox[{"7", " ", ";;", " ", "12"}]}], "]"}], "]"}]}], "}"}], 
        "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"vecs", " ", "=", " ", 
       RowBox[{"energysetaxes", "[", "quadform", "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"vals", " ", "=", " ", 
       RowBox[{"energyseteigenvalues", "[", "quadform", "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"ax1", " ", "=", " ", 
       RowBox[{
        RowBox[{
         RowBox[{"vecs", "[", 
          RowBox[{"[", "1", "]"}], "]"}], "[", 
         RowBox[{"[", 
          RowBox[{"1", ";;", "2"}], "]"}], "]"}], " ", "*", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"terminaltime", " ", 
            RowBox[{"thrustmax", "^", "2"}], "*", " ", 
            RowBox[{"vals", "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], ")"}], " ", "//", " ", 
          "Sqrt"}], ")"}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"ax2", " ", "=", " ", 
       RowBox[{
        RowBox[{
         RowBox[{"vecs", "[", 
          RowBox[{"[", "2", "]"}], "]"}], "[", 
         RowBox[{"[", 
          RowBox[{"1", ";;", "2"}], "]"}], "]"}], " ", "*", " ", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"terminaltime", " ", 
            RowBox[{"thrustmax", "^", "2"}], "*", " ", 
            RowBox[{"vals", "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], ")"}], " ", "//", " ", 
          "Sqrt"}], ")"}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"ax1", ",", "ax2"}], "}"}]}]}], "]"}]}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"getAxes", "[", 
   RowBox[{"3600", ",", " ", 
    RowBox[{"energyform", "[", 
     RowBox[{"3600", ",", " ", "n"}], "]"}], ",", " ", "tmax"}], "]"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"*", "\[IndentingNewLine]", "axes"}], " ", "=", " ", 
     RowBox[{"getAxes", "[", 
      RowBox[{"3600", ",", " ", 
       RowBox[{"energyform", "[", 
        RowBox[{"3600", ",", " ", "n"}], "]"}], ",", " ", "tmax"}], "]"}]}], ";",
     "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"ellipseVecFromUnitVec", "[", 
      RowBox[{"uv_", ",", "ax1_", ",", " ", "ax2_"}], "]"}], ":=", 
     RowBox[{
      RowBox[{"1", "/", 
       RowBox[{"Sqrt", "[", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"Normalize", "[", "ax1", "]"}], ".", "uv"}], ")"}], "^", 
           "2"}], "/", 
          RowBox[{
           RowBox[{"Norm", "[", "ax1", "]"}], "^", "2"}]}], "+", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"Normalize", "[", "ax2", "]"}], ".", "uv"}], ")"}], "^", 
           "2"}], "/", 
          RowBox[{
           RowBox[{"Norm", "[", "ax2", "]"}], "^", "2"}]}]}], "]"}]}], "*", 
      "uv"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"pos", " ", "=", " ", 
     RowBox[{"ellipseVecFromUnitVec", "[", 
      RowBox[{"delta", ",", " ", 
       RowBox[{"axes", "[", 
        RowBox[{"[", "1", "]"}], "]"}], ",", " ", 
       RowBox[{"axes", "[", 
        RowBox[{"[", "2", "]"}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"pos6", "=", " ", 
     RowBox[{"Maximize", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"Join", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], 
           "}"}], ",", " ", 
          RowBox[{"{", 
           RowBox[{
           "x1", ",", " ", "x2", ",", " ", "x3", ",", "x4", ",", "x5", ",", 
            "x6"}], "}"}]}], "]"}], ".", 
        RowBox[{"energyform", "[", 
         RowBox[{"3600", ",", " ", "n"}], "]"}], ".", 
        RowBox[{"Join", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], 
           "}"}], ",", 
          RowBox[{"{", 
           RowBox[{
           "x1", ",", " ", "x2", ",", " ", "x3", ",", "x4", ",", "x5", ",", 
            "x6"}], "}"}]}], "]"}]}], ",", " ", 
       RowBox[{
        RowBox[{"x1", " ", "==", 
         RowBox[{"pos", "[", 
          RowBox[{"[", "1", "]"}], "]"}]}], " ", "&&", " ", 
        RowBox[{"x2", "==", 
         RowBox[{"pos", "[", 
          RowBox[{"[", "2", "]"}], "]"}]}], " ", "&&", " ", 
        RowBox[{
         RowBox[{
          RowBox[{"Join", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], 
             "}"}], ",", " ", 
            RowBox[{"{", 
             RowBox[{
             "x1", ",", " ", "x2", ",", " ", "x3", ",", "x4", ",", "x5", ",", 
              "x6"}], "}"}]}], "]"}], ".", 
          RowBox[{"energyform", "[", 
           RowBox[{"3600", ",", " ", "n"}], "]"}], ".", 
          RowBox[{"Join", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], 
             "}"}], ",", 
            RowBox[{"{", 
             RowBox[{
             "x1", ",", " ", "x2", ",", " ", "x3", ",", "x4", ",", "x5", ",", 
              "x6"}], "}"}]}], "]"}]}], " ", "==", " ", 
         RowBox[{"3600", " ", 
          RowBox[{"tmax", "^", "2"}]}]}]}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"x1", ",", "x2", ",", "x3", ",", "x4", ",", "x5", ",", "x6"}],
         "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", "pos6"}], 
   "\[IndentingNewLine]", "**)"}]}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.9266862880748672`*^9, 3.92668631696879*^9}, {
  3.9266863733692923`*^9, 3.926686379975095*^9}, {3.9266864127496576`*^9, 
  3.9266864133944426`*^9}, {3.937158038489808*^9, 3.9371580817709837`*^9}, {
  3.9388170605335827`*^9, 3.938817102534504*^9}, {3.9388173196549473`*^9, 
  3.938817323692441*^9}},ExpressionUUID->"a7361479-b25e-4f65-bd51-\
39ad7bd317c6"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"18504.73878058888`", ",", 
     RowBox[{"-", "38556.18796070086`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "9241.209744303138`"}], ",", 
     RowBox[{"-", "4435.245841965427`"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.93715839886477*^9},
 CellLabel->
  "Out[409]=",ExpressionUUID->"58a102f0-0ea4-4016-ad5a-945abb3748d6"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"*", "Computing", " ", "the", " ", "commparisons"}], " ", "**)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"ellipseAxes", "=", 
     RowBox[{"MapThread", "[", 
      RowBox[{"getAxes", ",", " ", 
       RowBox[{"{", 
        RowBox[{"terminaltimes", ",", "forms", ",", " ", "tmaxs"}], "}"}]}], 
      " ", "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"computedmmratios", " ", "=", " ", 
     RowBox[{"MapThread", "[", 
      RowBox[{"reachablestatsoperator", ",", " ", 
       RowBox[{"{", 
        RowBox[{"bounds", ",", " ", 
         RowBox[{"ellipseAxes", "[", 
          RowBox[{"[", 
           RowBox[{";;", ",", "1"}], "]"}], "]"}], ",", 
         RowBox[{"ellipseAxes", "[", 
          RowBox[{"[", 
           RowBox[{";;", ",", "2"}], "]"}], "]"}]}], "}"}]}], " ", "]"}]}], 
    ";"}], "\[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[{"ArcTan", "[", 
      RowBox[{
       RowBox[{"k", "[", 
        RowBox[{"[", "1", "]"}], "]"}], ",", 
       RowBox[{"k", "[", 
        RowBox[{"[", "2", "]"}], "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"reflection", "[", "th_", "]"}], " ", ":=", " ", 
     RowBox[{"ArcTan", "[", 
      RowBox[{"1", " ", "/", " ", 
       RowBox[{"(", 
        RowBox[{"th", " ", "//", " ", "Tan"}], ")"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"toq1q2", "[", "th_", "]"}], ":=", " ", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"th", " ", "<", " ", "0"}], ",", " ", 
       RowBox[{"Pi", " ", "+", " ", "th"}], ",", " ", "th"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"maxangles", " ", "=", " ", 
     RowBox[{"Map", "[", 
      RowBox[{"vectoangle", ",", " ", "maxvecs"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"minangles", " ", "=", " ", 
     RowBox[{"Map", "[", 
      RowBox[{"vectoangle", ",", " ", "minvecs"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"maxangles", " ", "=", " ", 
     RowBox[{"Map", "[", 
      RowBox[{"reflection", ",", " ", "maxangles"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"minangles", " ", "=", " ", 
     RowBox[{"Map", "[", 
      RowBox[{"reflection", ",", " ", "minangles"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"maxangles", " ", "=", " ", 
     RowBox[{"Map", "[", 
      RowBox[{"toq1q2", ",", " ", "maxangles"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"minangles", " ", "=", " ", 
     RowBox[{"Map", "[", 
      RowBox[{"toq1q2", ",", " ", "minangles"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", 
   "\[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.9266003847116704`*^9, 
   3.9266004034045787`*^9}, {3.926600769931049*^9, 3.926600773632177*^9}, 
   3.9266863399664884`*^9, {3.9266863938894014`*^9, 3.926686424014738*^9}, {
   3.9266864822723403`*^9, 3.926686624195409*^9}, {3.9270392942533617`*^9, 
   3.92703929996848*^9}, {3.9270395018491592`*^9, 3.9270395045792713`*^9}, {
   3.927039608371284*^9, 3.927039678486886*^9}, {3.927039832848875*^9, 
   3.927039849208479*^9}, {3.927039890364254*^9, 3.927039891097527*^9}, {
   3.927040205611505*^9, 3.927040208525373*^9}, {3.927040268684758*^9, 
   3.9270403175238113`*^9}, {3.927040365063211*^9, 3.9270403708375607`*^9}, {
   3.927040418058381*^9, 3.927040438069558*^9}, {3.927040496709632*^9, 
   3.927040500805028*^9}, {3.9270405585998898`*^9, 3.927040560001799*^9}, {
   3.9270405908981953`*^9, 3.927040594180419*^9}, {3.9270425635006104`*^9, 
   3.927042600676921*^9}, {3.927042682784376*^9, 3.927042684661235*^9}, {
   3.927042742574403*^9, 3.927042747074641*^9}, {3.9270427895726957`*^9, 
   3.927042875579458*^9}, {3.9270429593933983`*^9, 3.927043007150466*^9}, {
   3.9270431150434504`*^9, 3.927043126322298*^9}, {3.9270431788612127`*^9, 
   3.927043178997695*^9}, {3.927043589473324*^9, 3.927043592606484*^9}, {
   3.927043655476087*^9, 3.927043665783655*^9}, {3.927044122448906*^9, 
   3.92704421249256*^9}, {3.9270443934502583`*^9, 3.9270443981337547`*^9}, {
   3.927044582637257*^9, 3.92704460808392*^9}, {3.927044673232853*^9, 
   3.92704479384916*^9}, {3.927044991363998*^9, 3.927045082539415*^9}, {
   3.927046225863181*^9, 3.927046368036644*^9}, {3.927046398549181*^9, 
   3.927046428510594*^9}, {3.92704651193163*^9, 3.9270465695572863`*^9}, {
   3.927046599812961*^9, 3.927046613087965*^9}, {3.9271257977861433`*^9, 
   3.92712580197996*^9}, {3.927125955362281*^9, 3.927125957956745*^9}, {
   3.927284066104637*^9, 3.9272841420605927`*^9}, {3.927284419433326*^9, 
   3.9272844294092913`*^9}, {3.927284829494843*^9, 3.9272848383676653`*^9}, {
   3.9272848866147127`*^9, 3.927284892527438*^9}, {3.9272849306106787`*^9, 
   3.9272849561509037`*^9}, {3.927285062657856*^9, 3.927285100833256*^9}, {
   3.9272853484219923`*^9, 3.927285361081882*^9}, {3.927285462619608*^9, 
   3.9272854663692493`*^9}, {3.927286047414912*^9, 3.927286053585079*^9}, {
   3.9272865195337973`*^9, 3.9272865206733303`*^9}, {3.9272865556569033`*^9, 
   3.92728656214788*^9}, {3.927289938366933*^9, 3.927289982513041*^9}, {
   3.927290034191708*^9, 3.9272900422640877`*^9}, 3.92729075564371*^9, 
   3.927290885775042*^9, {3.928611476236161*^9, 3.928611490343217*^9}, {
   3.928611530683065*^9, 3.928611573017138*^9}, {3.928612615303673*^9, 
   3.9286126165138807`*^9}, {3.928613314642189*^9, 3.928613402566622*^9}, {
   3.930228792296812*^9, 3.930228814992293*^9}, {3.9302289235763397`*^9, 
   3.9302289314566936`*^9}, {3.93023001914375*^9, 3.930230019887063*^9}, {
   3.930230118532915*^9, 3.930230145698024*^9}},
 CellLabel->
  "In[410]:=",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[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{"terminaltimes", "[", 
             RowBox[{"[", 
              RowBox[{"3", ";;"}], "]"}], "]"}], " ", "/", " ", 
            RowBox[{"(", 
             RowBox[{"2", " ", 
              RowBox[{"Pi", " ", "/", " ", "n"}]}], ")"}]}], ",", " ", 
           RowBox[{"maxrats", "[", 
            RowBox[{"[", 
             RowBox[{"3", ";;"}], "]"}], "]"}]}], "}"}], " ", "//", " ", 
         "Transpose"}], ",", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{"terminaltimes", "[", 
             RowBox[{"[", 
              RowBox[{"3", ";;"}], "]"}], "]"}], " ", "/", " ", 
            RowBox[{"(", 
             RowBox[{"2", " ", 
              RowBox[{"Pi", " ", "/", " ", "n"}]}], ")"}]}], ",", " ", 
           RowBox[{"minrats", "[", 
            RowBox[{"[", 
             RowBox[{"3", ";;"}], "]"}], "]"}]}], "}"}], " ", "//", " ", 
         "Transpose"}]}], "}"}], ",", " ", 
      RowBox[{"AxesLabel", "\[Rule]", " ", 
       RowBox[{"{", 
        RowBox[{"\"\<Terminal Time (Revs)\>\"", ",", " ", "\"\<Ratio\>\""}], 
        " ", "}"}]}], " ", ",", " ", 
      RowBox[{"PlotLegends", "\[Rule]", 
       RowBox[{"LineLegend", "[", 
        RowBox[{"{", 
         RowBox[{"\"\<Maximum Ratio\>\"", ",", "\"\<Minimum Ratio\>\""}], 
         "}"}], "]"}]}]}], "]"}]}], "  ", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"*", "l2"}], " ", "=", " ", 
     RowBox[{"ListLinePlot", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"terminaltimes", "[", 
           RowBox[{"[", 
            RowBox[{"3", ";;"}], "]"}], "]"}], ",", " ", 
          RowBox[{"minrats", "[", 
           RowBox[{"[", 
            RowBox[{"3", ";;"}], "]"}], "]"}]}], "}"}], " ", "//", " ", 
        "Transpose"}], ",", " ", 
       RowBox[{"AxesLabel", "\[Rule]", " ", 
        RowBox[{"{", 
         RowBox[{"\"\<Terminal Time (s)\>\"", ",", " ", "\"\<Ratio\>\""}], 
         " ", "}"}]}]}], "]"}]}], " ", "**)"}], "\[IndentingNewLine]", 
   RowBox[{"l3", " ", "=", " ", 
    RowBox[{"ListLinePlot", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{"terminaltimes", "[", 
             RowBox[{"[", 
              RowBox[{"2", ";;"}], "]"}], "]"}], " ", "/", " ", 
            RowBox[{"(", 
             RowBox[{"2", " ", 
              RowBox[{"Pi", " ", "/", " ", "n"}]}], ")"}]}], ",", " ", 
           RowBox[{"maxangles", "[", 
            RowBox[{"[", 
             RowBox[{"2", ";;"}], "]"}], "]"}]}], "}"}], " ", "//", " ", 
         "Transpose"}], ",", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{"terminaltimes", "[", 
             RowBox[{"[", 
              RowBox[{"2", ";;"}], "]"}], "]"}], " ", "/", " ", 
            RowBox[{"(", 
             RowBox[{"2", " ", 
              RowBox[{"Pi", " ", "/", " ", "n"}]}], ")"}]}], ",", " ", 
           RowBox[{"minangles", "[", 
            RowBox[{"[", 
             RowBox[{"2", ";;"}], "]"}], "]"}]}], "}"}], " ", "//", " ", 
         "Transpose"}]}], "}"}], ",", " ", 
      RowBox[{"AxesLabel", "\[Rule]", " ", 
       RowBox[{"{", 
        RowBox[{"\"\<Terminal Time (Revs)\>\"", ",", 
         " ", "\"\<Angle (Rad)\>\""}], " ", "}"}]}], ",", " ", 
      RowBox[{"PlotLegends", "\[Rule]", 
       RowBox[{"LineLegend", "[", 
        RowBox[{"{", 
         RowBox[{"\"\<Maximum Angle\>\"", ",", "\"\<Minimum Angle\>\""}], 
         "}"}], "]"}]}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"*", "l4"}], " ", "=", " ", 
     RowBox[{"ListLinePlot", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"terminaltimes", "[", 
           RowBox[{"[", 
            RowBox[{"2", ";;"}], "]"}], "]"}], ",", " ", 
          RowBox[{"minangles", "[", 
           RowBox[{"[", 
            RowBox[{"2", ";;"}], "]"}], "]"}]}], "}"}], " ", "//", " ", 
        "Transpose"}], ",", " ", 
       RowBox[{"AxesLabel", "\[Rule]", " ", 
        RowBox[{"{", 
         RowBox[{"\"\<Terminal Time (s)\>\"", ",", 
          " ", "\"\<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}, {3.9266867763785243`*^9, 3.926686805070363*^9}, {
   3.9266868770836587`*^9, 3.9266868829006505`*^9}, {3.927197778322975*^9, 
   3.927197797921386*^9}, {3.9271978532101517`*^9, 3.927197872468264*^9}, {
   3.927197971557757*^9, 3.9271980013953543`*^9}, {3.927198039834714*^9, 
   3.927198041024774*^9}, {3.927198110103691*^9, 3.9271981811740303`*^9}, 
   3.927198218490815*^9, {3.9272841638683*^9, 3.927284207050194*^9}, {
   3.9272843911804533`*^9, 3.927284391524678*^9}, {3.927284626031417*^9, 
   3.9272846263771048`*^9}, {3.927286094551633*^9, 3.92728609907767*^9}, {
   3.927289889080112*^9, 3.927289958532634*^9}, {3.928063064553461*^9, 
   3.9280630933901663`*^9}, {3.928576976981681*^9, 3.9285770212574787`*^9}, {
   3.929285006220785*^9, 3.9292850065975857`*^9}, {3.929285036960163*^9, 
   3.9292850429134197`*^9}, {3.929285101050829*^9, 3.929285137349538*^9}, {
   3.929537653675254*^9, 3.929537683703032*^9}, {3.929538507210573*^9, 
   3.929538532172489*^9}, {3.9295411204416237`*^9, 3.929541177621839*^9}, {
   3.929541277820133*^9, 3.9295413920920763`*^9}, {3.929541495851563*^9, 
   3.929541511872821*^9}, {3.929542527303109*^9, 3.929542566478361*^9}},
 CellLabel->
  "In[425]:=",ExpressionUUID->"79927019-60b3-4e35-99ba-21182e4ed03c"],

Cell[BoxData[
 TemplateBox[{
   GraphicsBox[{{}, {{{}, {}, {
        Hue[0.67, 0.6, 0.6], 
        Directive[
         PointSize[
          NCache[
           Rational[1, 72], 0.013888888888888888`]], 
         RGBColor[0.368417, 0.506779, 0.709798], 
         AbsoluteThickness[1.6]], 
        LineBox[CompressedData["
1:eJw1k3s41FkcxkU3dpFEN6KtrK5K2CT7kpQumi1UE7ZIF7WZ9IRCtq1kK9vu
PpKu203aLiSVbUmkXLqNCTO5DaYZjcsYmt/vuOWy9un8/jjPec7znPM93+/7
ft7JgYK1W7W1tLQEg+v/vcInLdWt7zx8MjsF+Qc1WOi4b9Ls/qvYNFQcM/4X
DY4vM9r2+GAKFjzQDnCJ1eDL/dtYFSEcbXNCA6uglUQYfQ9eFVV6kt+59/dR
oufTP/mMBoF5i/M1zg/BlwV97r3I1cvE9rMTh/GSNUgX1unG5DzGll7f0/13
uPpZkI1bdzU9U4OBmug1us458Fu0/O+efO6/XJwZk90ZKNLgJolftaI2D5eU
ojJJLff/M6xQKVZ9r9ZAGO3bbjjpOY6ZyZOj+7l+XmBEbpJkjyEDojPjtNi/
AFp66Y4dlgztrxCnw+46N9gyMD/Z9d2FS0XweCi3G+PO0H6LwbID1hvWM1gy
uqh6s/Qljjq5Z8QGM7T/V1g+4PE6NIrBrnOJP1uZv8H0AxsvaOIZOs9btL1e
Z1t/iUGCZdAUlZ8QBUW75pmkMXS+EkyXGvgLchhIvc7JN6wWobvf20T0mqHz
iiDYKcwwq2QQmnTokVPIO0R5ZN+JUjF0/nco0HRdnKPLYmj19rhJp0oxfdfw
RfxpLNWjDAafCw1zXFicncTjD0krg6W33QELX5bqU44LvzomrdjHYlagw0zF
23IM716vYxDPUr3E+PzOLwjXWOTeMO8rbBXjmX6+OCuTpfpJcNk5c8qWlyy8
moaW3NJ/DzMbY+nkKpbq+R4x/i3OykYWylmqK/GzKxA4J6TjFmGpvhWYGlhZ
uk2LIHJP2V6BZyXi5JLllrqE6l2JnY3qDqkhgcHDrCVrd1chsmTkiLtjCNW/
Cu6bBRfPjyW41nnV1P63akQmZ9jkjCPUjxokRA50u00lsHc63jg2tQbtEwJG
5bkS6o8UEfw5qwz8CV7G7MnqeSPF3EiHWUbhhPpViyN5yoyMeAL//PXxUlUt
VjdeaZZdIdS/OkgSFpr+lUHQddjCZvPweqiuj19a94xQP+ux6J6j+K6QIH+q
8dk1Z+qxL2mZj1Elof7WQzH6zkIzGUF84fAhblYyzPcX2sk+Euq3DCMdeEcj
mgnW7egJtsuUwZvn4dXdQqj/Mojbs4yPDJ4t9NSl05Z+QLmTkdK2iVAePsA0
RcfIuoGg6Y7MaazkA/aa3oo4XEcoH3Kk+yXNC6wgeOApTh65TY5pZvHfKEsI
5UWOJOZcytxCgoNtxfo9RA6bwy4r/bMJ5UeBJ/uViSfuESz780l4S6wCKxMF
kvJrhPKkwLGbhsPWJxIYzU+vqzFpAM9E+/aMOEL5asANh6632/cTVJVf9xDe
aIB0+/H2qcGE8taAilfF+hF8guTwpPu59h+RWmOVELqcUP4+gpWIwqwcCULG
nZxwv+Ajer2yHFKsCeVRiQ2CtE8Dg3wsyIo5cs1HiaXed0PiBnn4wqcS398I
W6NrR6Dtt1eV0KAEL6E4qnKQjy+8NoIX2mnc50nwpm+rT2xYI8LF/ISoDYTy
2wj+rry4HwMIEi/zn4YPa0K7pZH3/R0cz01Q/NvafiiEYJOr57c7Epuws/h8
WvVeju8mpLS96hKFEUyXu/zBn9aMyNxZ5qERHO/N0I5zNHk7eGaO2nWveNSM
OlGWV0c4x38zMk4GnOnfR5BjZR24yL0F+OlCdmcol4cW8PHgIrOb4FjxxNez
xS2webpGpRPM5UOFnhS38Qu3EPyw09DOYqsKag8t+5t+XF5UsHdha319CCZ8
rXNpFFHhebCb7Q5PLj+t6P1Q8k/1EgJFascw7dhWHA/v31jgxOWpFYsnewXM
tyVI4zWHMGPU0K8/FTTfmsuXGqZ+erxSc4L9n6TvFclqIK6Vb2bM5U0NbS17
ffORBK4J71wkdm1gTmm71vWyNH9tUIfbWu7+xOIr+4JbRS/a8DRVNVPUwOI/
pEmbHg==
         "]]}, {
        Hue[0.9060679774997897, 0.6, 0.6], 
        Directive[
         PointSize[
          NCache[
           Rational[1, 72], 0.013888888888888888`]], 
         RGBColor[0.880722, 0.611041, 0.142051], 
         AbsoluteThickness[1.6]], 
        LineBox[CompressedData["
1:eJw1lHk01WkcxiXJUiiZOi0XxxJtZhpL0zKPrX1KtmISUjeVuEhqCmmRlJwR
IulIqZksRWE6Ji1SKGvWuJfLvZbL5fb+jMjpam6n+/7xO7/znvN93/f7fZ7P
8+r7cJzYigoKChzZ9+3f6vog106aimQNvnf2cQarfwlhLZ/MQMOhiJNWoQxi
Ns468CT8Huw4Lk8tQhh8r8+CWJTskR3EwHj/1tGasIfoS3aQpAXQ/fnId7KV
aPkx8HlhW8qsK0CDi/SnWb70vCLsb/h9SrYPg7yaTtWIkido2x2ySbKHnl8M
J5u4a4wbg6/cMEfVdSWQOLkaVzvT+56jM61UmLqdwV+jsb9t6XgB9knO8VOb
6f0voS99Zn3VnkFN2O6PmqxX6BzO9la2pv2Ugc2W5iiuZTA6dUli057XUN5Q
m1K4ivb3BiF3zRQcLRksujxudeNmOV6ejapQtaD9VkAocbabK1vbzy5v9+ZV
wiyuzPqRJe3/LV5+1SvTXs3A73rSaeNFVWAHtU7+upLOU42i1Vk2ZSYMEvT2
G4g9apCTcDDqPxadrxYai71vdsxhwHO+LnDbXoeR1QbKV9XovHXwzL3fzFJg
EJQcWbgmoB5zTFr+SRol8vnrkRD5d/f0QQKldt9oVtx7SE6IjaL5RK5HA8z4
wlGjZoIUloP7lAcNqGflXhp7R+T6NMLUrlpHq5RgmY/lUmF1I37s11sQ9YTI
9WrC3rd9E/4PCZ7fXSR9M9SEqUsKioT3iFy/Zti27xQrpxM4i5Rq789sgUNq
c++nZCLXswWxkVvPceMJ+paJb8Uub0Vw2kRsTyyR69sK7q1zFutjCE4GNgRz
tn3AwiiD4yujiVzvDzCaMRjac4FAo6DY3sm/DZ552JB4kcj1b4Nub/wVjuy8
22MZP1hcaUfX85Z76VeJ3A8uPKT6UbY3CCzWxPTPzeXCIG1p3GFZ/9/94WHZ
o2durAKCyojA4okqHlJvmhm4lxG5Xx2QjM+rVG8h2FO6K5Yn7kDgY3VFlpjI
/etE/yaXEs5UBuNndc28lfmQuLw61baQ+skH9+PWxJUy3koNtVMcr/HhG+zo
fcee+svHwDXBwNp1DGLfKE+xM+6CSt1MiYIl9bsL/hMBF8gKBjsPThwyL+pC
sWlVkpYJ5b8Lav/qNnrpM9BVG35vtKEbRrae9fz5jJyHbtxXEXxJlPElyu5a
M7e5G5t3PMoO16B5EaAj+AgrTYXB421NmSoHBOidoRZNFBk5LwLsrf9D47SU
IFxSMXNiVAAfztiH9eOUHyFGDPy1No4QbIx/GjoYJYT/0ZjOc8OUJyGSwicj
pSKCWT/ndXJ1euA9PwWFPZSvHjh5Hg7L6SJoa7yzqeZuD/rztQuEPMpbD/hN
que92ggyQ5Pzn1v0wq3dLI/VQvnrxZmz+0wMGwkC5l2en/+6F1l+mRrB9ZTH
PqhrzSxRqyVYVRxx7rZrHxbfCan6WEX57IO9qtWCxbI8KHoEixN6+mBZeLQx
r5Ly2o9XWanH4ysIqqRs16hj/RjSsuNXlVN++2HIjt/uK1snpbs/C50mQqL1
Gw67nPIsQonl5IFv9V422xYfTBLBvdO+OaOC8i2Clce+vkHZfaYC6z/djQag
2ZLikv+O8j4Av92vAj9XE4ycN/+8pXAAC0wjV1TVUf4HMBKta7xSNn+JsYnP
2vWDSB89c928leZhEEdCXI7xuAQXKha8W940iNpTqg5W3TQfYlSGhnpu6yfY
cVjTXJcthnPYruClEpoXMbRzBB28TwTzZ0y9qTUqRurlVHbQV5qfITgdE+mM
yPgR5n6aphg1hCmfMowDtRl5noYw9DlAdbMsDw8cBgJG5gwj/WG5xvR59D0c
hpJpumGzrP4E4bUIM4dR/OVpUIEmI8/bMPTvKendUmdgk1Bv3WwuwcWwLK/k
6fT9lICrM9svRZY/dYvX98vLJMgeE1zKkPX3P03rnV4=
         "]]}}}, {{}, {}}}, {
    DisplayFunction -> Identity, DisplayFunction -> Identity, AspectRatio -> 
     NCache[GoldenRatio^(-1), 0.6180339887498948], Axes -> {True, True}, 
     AxesLabel -> {
       FormBox["\"Terminal Time (Revs)\"", TraditionalForm], 
       FormBox["\"Ratio\"", TraditionalForm]}, 
     AxesOrigin -> {0, 1.0801738368476572`}, 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]]}& )}}, 
     PlotRange -> {{0, 0.9927709880663839}, {1.0801738368476572`, 
       1.370110822031127}}, PlotRangeClipping -> True, PlotRangePadding -> {{
        Scaled[0.02], 
        Scaled[0.02]}, {
        Scaled[0.05], 
        Scaled[0.05]}}, Ticks -> {Automatic, Automatic}}], 
   FormBox[
    FormBox[
     TemplateBox[{"\"Maximum Ratio\"", "\"Minimum Ratio\""}, "LineLegend", 
      DisplayFunction -> (FormBox[
        StyleBox[
         StyleBox[
          PaneBox[
           TagBox[
            GridBox[{{
               TagBox[
                GridBox[{{
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #}, {
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    RGBColor[0.880722, 0.611041, 0.142051], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    RGBColor[0.880722, 0.611041, 0.142051], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #2}}, 
                 GridBoxAlignment -> {
                  "Columns" -> {Center, Left}, "Rows" -> {{Baseline}}}, 
                 AutoDelete -> False, 
                 GridBoxDividers -> {
                  "Columns" -> {{False}}, "Rows" -> {{False}}}, 
                 GridBoxItemSize -> {"Columns" -> {{All}}, "Rows" -> {{All}}},
                  GridBoxSpacings -> {
                  "Columns" -> {{0.5}}, "Rows" -> {{0.8}}}], "Grid"]}}, 
             GridBoxAlignment -> {"Columns" -> {{Left}}, "Rows" -> {{Top}}}, 
             AutoDelete -> False, 
             GridBoxItemSize -> {
              "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
             GridBoxSpacings -> {"Columns" -> {{1}}, "Rows" -> {{0}}}], 
            "Grid"], Alignment -> Left, AppearanceElements -> None, 
           ImageMargins -> {{5, 5}, {5, 5}}, ImageSizeAction -> 
           "ResizeToFit"], LineIndent -> 0, StripOnInput -> False], {
         FontFamily -> "Arial"}, Background -> Automatic, StripOnInput -> 
         False], TraditionalForm]& ), 
      InterpretationFunction :> (RowBox[{"LineLegend", "[", 
         RowBox[{
           RowBox[{"{", 
             RowBox[{
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"PointSize", "[", 
                    FractionBox["1", "72"], "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.24561133333333335`, 0.3378526666666667, 
                    0.4731986666666667], FrameTicks -> None, PlotRangePadding -> 
                    None, ImageSize -> 
                    Dynamic[{
                    Automatic, 
                    1.35 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                    Magnification])}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.368417`", ",", "0.506779`", ",", "0.709798`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.368417, 0.506779, 0.709798]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.368417, 0.506779, 0.709798], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}], 
               ",", 
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"PointSize", "[", 
                    FractionBox["1", "72"], "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.880722, 0.611041, 0.142051], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.587148, 0.40736066666666665`, 0.09470066666666668], 
                    FrameTicks -> None, PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 
                    1.35 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                    Magnification])}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.880722`", ",", "0.611041`", ",", "0.142051`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.880722, 0.611041, 0.142051]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.880722, 0.611041, 0.142051], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}]}],
              "}"}], ",", 
           RowBox[{"{", 
             RowBox[{#, ",", #2}], "}"}], ",", 
           RowBox[{"LegendMarkers", "\[Rule]", 
             RowBox[{"{", 
               RowBox[{
                 RowBox[{"{", 
                   RowBox[{"False", ",", "Automatic"}], "}"}], ",", 
                 RowBox[{"{", 
                   RowBox[{"False", ",", "Automatic"}], "}"}]}], "}"}]}], ",", 
           RowBox[{"Joined", "\[Rule]", 
             RowBox[{"{", 
               RowBox[{"True", ",", "True"}], "}"}]}], ",", 
           RowBox[{"LabelStyle", "\[Rule]", 
             RowBox[{"{", "}"}]}], ",", 
           RowBox[{"LegendLayout", "\[Rule]", "\"Column\""}]}], "]"}]& ), 
      Editable -> True], TraditionalForm], TraditionalForm]},
  "Legended",
  DisplayFunction->(GridBox[{{
      TagBox[
       ItemBox[
        PaneBox[
         TagBox[#, "SkipImageSizeLevel"], Alignment -> {Center, Baseline}, 
         BaselinePosition -> Baseline], DefaultBaseStyle -> "Labeled"], 
       "SkipImageSizeLevel"], 
      ItemBox[#2, DefaultBaseStyle -> "LabeledLabel"]}}, 
    GridBoxAlignment -> {"Columns" -> {{Center}}, "Rows" -> {{Center}}}, 
    AutoDelete -> False, GridBoxItemSize -> Automatic, 
    BaselinePosition -> {1, 1}]& ),
  Editable->True,
  InterpretationFunction->(RowBox[{"Legended", "[", 
     RowBox[{#, ",", 
       RowBox[{"Placed", "[", 
         RowBox[{#2, ",", "After"}], "]"}]}], "]"}]& )]], "Output",
 CellChangeTimes->{3.937158399088071*^9},
 CellLabel->
  "Out[425]=",ExpressionUUID->"d7200a32-ec83-4fa2-b57f-ea6b6ff8d6c1"],

Cell[BoxData[
 TemplateBox[{
   GraphicsBox[{{}, {{{}, {}, {
        Hue[0.67, 0.6, 0.6], 
        Directive[
         PointSize[
          NCache[
           Rational[1, 72], 0.013888888888888888`]], 
         RGBColor[0.368417, 0.506779, 0.709798], 
         AbsoluteThickness[1.6]], 
        LineBox[CompressedData["
1:eJw1k3s41FkYx2eamd9zkiKXahVJNbroakZXvcctRVOGnVZbSUKMlCSpLaLU
07pspRK77UZIFyJ0mZrIXSGEboYZg8HM0GWnotS2j+OP85znPOf2vt/v5zvN
a4+rzygajRb0Y/w/L18WYjLv2xloH3uyUX/6ALwUZGfZDSUDc0mZMmPmAAzv
p0BrlZYwcc4AnHIc73vvSAbM4NlmJSweOX8d0qbdlHJWDQDb21lTe/gWxMd9
sP7KG7mfC2YRnnF8rwHwKrItfm+dDxyr32uifht57w5Y8OubS5MGIKe2bXS4
+B6c3qIu4YhH3hcBZ5eB8zUDOv7ecpg/2loMu1YqcuaPpuPh/wohJNDJ+th3
Gr6qiV3n1FoEl3kvJig/0fDw/4/B136KQ/J7Gq49vPmtjkkJvL18vSu/n4aH
6ykFg9ulhi4/1hrGnHNNW8tAXTNz3b53NDxcXzmYz987wegjDRvHfF7y56UK
qLw4Rc0douHheivBqKRRN4eiY3u9ijeekirQWuzDD9Cj4+H6n0Djl7xsB1M6
Dkg6H8E2roYlt849+WnhSD818FlYaV5kQ8cJpt7TVVtq4YhWmJWuYKS/Z/BM
OHqRzJ+OJW5Jcvf1dfDdI7tG6yid9FsHZgH3q3wS6Xhv4tGCFbvrofubJlWW
ziD910OAbnhhyxUGZr7ZedIkvgEeRS29YJfGIHo8hwGUtl2WwcAXTTZsomc/
hzHs0ltRNxhEn0bgngrS/5jLwBZeVnM7ahpBUatM1hYxiF5NwB6VOSm5hIEL
042HytVNcD1Pb+3+WgbRrxmsW1e3HnjNwG49zGfXxr4AyVlRWLiCQfR8AWsF
B876aRhYYaG6HDvvJbBF5tsnMZlE35dQbbfqdrA+Ex8Keh68h/cKfq7ILPGY
ziR6v4IUc6lRuSUTj8sX2bsGvgb39PTN5+yZRP/XwNP2mHpXwMSpn1ImcOPe
QNby3BrjnUziRwtEHcyLMROzMHfFqe6JWS0gX6nYsbGQRfyRgN0fjMyCYhau
Cg8SDVZLgIsOB5lUsIhfrZA+eaFzcDULby3+JVaiaoX5tuZmFxpYxL82EPgb
ZwS9YuHPUVMXeFJSYHuExKikLOKnFEq9PYs03SxcPEP/Iv+CFBINI0MOvWMR
f6UwxehYqnCQhWPLKbodWwahMjPpXQZF/JbBkz6epWAshTf6Dfpz7shglk2L
0moiRfyXQdkX70jXaRSeqtXXMHN1O6yrexjy91yK8NAOwvi7HD0rCvfckK2Y
2NwOEQ/qWBmYInzIwVrhp3J1pnAerykN+cqhOGx9kt5GivAiB67Mxb3Nk8JH
+ivHDmrkYHc0jZYfQBF+OkDHpKMrJpTCjmcehiqjOwA/iKrziqQITx1gdiO6
mxNL4fGWOW0thp3gGb2j4NsFivDVCTpx3TriFAq/bryypja9E0zjuZMDb1KE
t05Y0CdcqX2XwmmhibmF3C7I/dhrePExRfjrgjuN+yLGVVN496QYo9yyLviL
czMysJkiPCqg0EmPnyel8FJR+LFUgQIaOLct7ecjwqcCNs952rZgAcKjtgSr
EjoVUOzyVXRkISK8doNt+dcKs8UIVw/5CKL3d8PJOT76hhxE+O2GMe1955yt
ED7/z6ZHoaweYOqIDR4tRYTnHvjsYmrgswLhbTY8c7/zPdD+7/r7y1YhwncP
xD3MFHNtEJ4tx6c3zeyFzLlzPdzsEeG9Fxyaqnhxjgh/OM4ZcCrohXyx/7VW
J0T47wXH22vUtusRFrNnea10UIJJSIqygI9IHpTATvUutRQgfKJy8tN5TUpI
6pJL7rkjkg8VWPjzW222IOwi1OFM9VFBpuVQ+tNtiORFBTnLFjm57UDYSJtx
SVejgp6rmbrNvojkRw3vg6MC+EKEO7I+skZFqyHf4bqmOBCRPKmhnJ+5Z9Ze
hLM39O7+YNAHDXk1DyNCEMlXHyho+dLyAwiHvZO86Ejrg8tr0NCXQ4jkrQ+8
1yVtMA5H2CahHjdz+uH7wbSFsyMRyV8/SPp/HW9yHOEx3LJrFaX9UMRPPjR4
AuH/AO9lPBA=
         "]]}, {
        Hue[0.9060679774997897, 0.6, 0.6], 
        Directive[
         PointSize[
          NCache[
           Rational[1, 72], 0.013888888888888888`]], 
         RGBColor[0.880722, 0.611041, 0.142051], 
         AbsoluteThickness[1.6]], 
        LineBox[CompressedData["
1:eJw1k3tcjFkcxrvM+04SYeReySVybdckaXlOalEot1DMalMql4rI5dONFVq1
aWXLuhbltjG60Ka7SlaT1Ix0VWZqukyTZGawYu3nff1xPueP8zvn/H7P833M
PAPWeutoaWkFfl3/7wtt9pnM/hyHtAz+LjlfjJeud9LsB/5E34jWk69Wi8Gc
J4HuGVXPCxIjatnw7dmhqUhpmO9xJPFb/S3sIZ+MdPPFMPdaoaoMuQtLeXSc
U8u3+/cQWuj4rPZfMTwLlxS/XZSJm0qDmjcGEva9+xCMsf/3nKEEwspXg8Ly
srF48oGZp798ez8HA3YZgqUvxfjSGLJm0KI8XJmckNkR/+2/ArQlCit8LMW4
rope6dRciPr2BVmxKTXs/0XwS//z4JzealSGbH5jaPIIdXp/mznrVbP9lMAo
ZY11pKIKKt0Z8RJBKXYqXUiu5TO2vzL4bauKSOZVwPjUe+vzFx9DRx5/LeZI
OdtvOZJ3Xu0bF14ChxGPGzyanqD6eEPMxlGFbP//gF8WcL0k7T52njsbbm5c
AbNDb6MsaBU7jwhz5jrHWburcGai12TFlkos0kxRyW6p2PmeoaQ4yWaJRoWm
deekm5yrcPtSeKPDYjU7bxXSbQtqi8LU2JMQkWXr/xztqz37tj9Qs/M/h2OK
6GF6uxqcBp8TJr9V40ipwKFbX8PqUYMGY6fLvEkaJJq4uGnfqYHFfosXLRYa
Vh8x0i1Gr7hrqsEsz/kzZSIxMj88DzmspWH1ksDDSUdQX65GQYrxQFmPBO4n
nGwG71ez+r2A1Zi9N7jaaqzr5Dy7OaQWSr2sxnxfFatnLURBn1323XgH+SzF
lejZL9Gcr+fnmtfP6vsSSy0/mDumvsXhwJq9AavqkLQ7sLLNq4/Vuw5DU63a
L3zqxdDMHIe1u+sxavngoAkaBat/PVys8lP1qjuRrEkaZRXTgPHeXyJP3mxn
/WgEXbRB+0GiFFa2UR2j0xoh6F9heimrhfWnCU3e63oGvtY9CQvM+VjRhKjS
aUKfTbWsX80wrjUQjTCugaB4Y3STohl9bwIMYysqWf9eYZyXrkZ/XTneHzWd
60G3oDhD33drZjHrZwsiwyIMLdMeongKL3HNHy3Ys9Vi/bR+XcL424JDWsFG
ZkM5JLqM1rY3b0VOyKE1/tM5hPG7FW8mH9TiLeGQDb4f/fj3W5FmPI/H28wh
jP+tmCjO3703iENM9ZXVU5e+RsUsIyn/FIcwPLzGTMGq735K5pDO2622o1+8
RoDc63NfNocwfEgxZOC4pL+SQzJWSa7pbZfCe1LRFB8ZhzC8SLF8n5GtywcO
Ce0tH/JRJcVHt8ald4dQhOFHhrCHvut/NaPIsrjc4O5IGWx0pSulfIowPMng
4bBalbuMIsPnCV81GrUhXvKwZ6w7RRi+2hDx62Hhp50UqRdfXV6Z0ga3w+Fx
PqEUYXhrwyRlTt6W3yhyLTjhXoFVO4b+FbKw+RJFGP7aofaokfXcoYj/mFPj
7pW245OfNDAmnyIMj3IIc+f65YsosiAn7JdkVzn4M6pmH22kCMOnHI+LTvfX
dVFEZ8texZk2OUwy1MVl7ynC8NoBb9G0A440TSoGvF0j93cg8Wp6riePJgy/
HfgxZv4w3kSanL3slh9MdcLATpHrPosmDM+deFcwvMFmAU222q2a5nu2E6Rm
sbPQniYM353wunDrUZ4zTSyk5LTb1C5UPPUX/uxGE4b3Lgh/r3c5v40m/cf4
H5yyurArtCvJbzdNGP67MHKh/o6nwTTJM5/u+cOP3Qi0WzKnIJwmTB664WSs
bW1/kibHy8c/nS3pxvTPRbGup2nC5EMBXp1I710CTVbvMOSbeitw4ZH3wgmX
acLkRYEx2alRT1NoMs5A9+IwlQLbgoLnfrlNEyY/PSgoLLTLvEcTWZqa0ons
QfTmBQXy+zRh8tSDeY59zYEuXHLHpcu/f6QStvMuhORs5BImX0qI1ldfKt/K
JQf7mmpl15RYO7bU57wPlzB5UyLW3fXd9wFcYnfmOXnB78VdC2tBbDCXMPnr
xR5nq2xhKJcMtiq9+bikF8MOxo9NOMYl/wFDFWZO
         "]]}}}, {{}, {}}}, {
    DisplayFunction -> Identity, DisplayFunction -> Identity, AspectRatio -> 
     NCache[GoldenRatio^(-1), 0.6180339887498948], Axes -> {True, True}, 
     AxesLabel -> {
       FormBox["\"Terminal Time (Revs)\"", TraditionalForm], 
       FormBox["\"Angle (Rad)\"", TraditionalForm]}, AxesOrigin -> {0, 0}, 
     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]]}& )}}, 
     PlotRange -> {{0, 0.9927709880663839}, {0, 3.0659074381203197`}}, 
     PlotRangeClipping -> True, PlotRangePadding -> {{
        Scaled[0.02], 
        Scaled[0.02]}, {
        Scaled[0.02], 
        Scaled[0.05]}}, Ticks -> {Automatic, Automatic}}], 
   FormBox[
    FormBox[
     TemplateBox[{"\"Maximum Angle\"", "\"Minimum Angle\""}, "LineLegend", 
      DisplayFunction -> (FormBox[
        StyleBox[
         StyleBox[
          PaneBox[
           TagBox[
            GridBox[{{
               TagBox[
                GridBox[{{
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #}, {
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    RGBColor[0.880722, 0.611041, 0.142051], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    RGBColor[0.880722, 0.611041, 0.142051], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #2}}, 
                 GridBoxAlignment -> {
                  "Columns" -> {Center, Left}, "Rows" -> {{Baseline}}}, 
                 AutoDelete -> False, 
                 GridBoxDividers -> {
                  "Columns" -> {{False}}, "Rows" -> {{False}}}, 
                 GridBoxItemSize -> {"Columns" -> {{All}}, "Rows" -> {{All}}},
                  GridBoxSpacings -> {
                  "Columns" -> {{0.5}}, "Rows" -> {{0.8}}}], "Grid"]}}, 
             GridBoxAlignment -> {"Columns" -> {{Left}}, "Rows" -> {{Top}}}, 
             AutoDelete -> False, 
             GridBoxItemSize -> {
              "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
             GridBoxSpacings -> {"Columns" -> {{1}}, "Rows" -> {{0}}}], 
            "Grid"], Alignment -> Left, AppearanceElements -> None, 
           ImageMargins -> {{5, 5}, {5, 5}}, ImageSizeAction -> 
           "ResizeToFit"], LineIndent -> 0, StripOnInput -> False], {
         FontFamily -> "Arial"}, Background -> Automatic, StripOnInput -> 
         False], TraditionalForm]& ), 
      InterpretationFunction :> (RowBox[{"LineLegend", "[", 
         RowBox[{
           RowBox[{"{", 
             RowBox[{
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"PointSize", "[", 
                    FractionBox["1", "72"], "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.24561133333333335`, 0.3378526666666667, 
                    0.4731986666666667], FrameTicks -> None, PlotRangePadding -> 
                    None, ImageSize -> 
                    Dynamic[{
                    Automatic, 
                    1.35 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                    Magnification])}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.368417`", ",", "0.506779`", ",", "0.709798`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.368417, 0.506779, 0.709798]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.368417, 0.506779, 0.709798], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}], 
               ",", 
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"PointSize", "[", 
                    FractionBox["1", "72"], "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.880722, 0.611041, 0.142051], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.587148, 0.40736066666666665`, 0.09470066666666668], 
                    FrameTicks -> None, PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 
                    1.35 (CurrentValue["FontCapHeight"]/AbsoluteCurrentValue[
                    Magnification])}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.880722`", ",", "0.611041`", ",", "0.142051`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.880722, 0.611041, 0.142051]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.880722, 0.611041, 0.142051], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}]}],
              "}"}], ",", 
           RowBox[{"{", 
             RowBox[{#, ",", #2}], "}"}], ",", 
           RowBox[{"LegendMarkers", "\[Rule]", 
             RowBox[{"{", 
               RowBox[{
                 RowBox[{"{", 
                   RowBox[{"False", ",", "Automatic"}], "}"}], ",", 
                 RowBox[{"{", 
                   RowBox[{"False", ",", "Automatic"}], "}"}]}], "}"}]}], ",", 
           RowBox[{"Joined", "\[Rule]", 
             RowBox[{"{", 
               RowBox[{"True", ",", "True"}], "}"}]}], ",", 
           RowBox[{"LabelStyle", "\[Rule]", 
             RowBox[{"{", "}"}]}], ",", 
           RowBox[{"LegendLayout", "\[Rule]", "\"Column\""}]}], "]"}]& ), 
      Editable -> True], TraditionalForm], TraditionalForm]},
  "Legended",
  DisplayFunction->(GridBox[{{
      TagBox[
       ItemBox[
        PaneBox[
         TagBox[#, "SkipImageSizeLevel"], Alignment -> {Center, Baseline}, 
         BaselinePosition -> Baseline], DefaultBaseStyle -> "Labeled"], 
       "SkipImageSizeLevel"], 
      ItemBox[#2, DefaultBaseStyle -> "LabeledLabel"]}}, 
    GridBoxAlignment -> {"Columns" -> {{Center}}, "Rows" -> {{Center}}}, 
    AutoDelete -> False, GridBoxItemSize -> Automatic, 
    BaselinePosition -> {1, 1}]& ),
  Editable->True,
  InterpretationFunction->(RowBox[{"Legended", "[", 
     RowBox[{#, ",", 
       RowBox[{"Placed", "[", 
         RowBox[{#2, ",", "After"}], "]"}]}], "]"}]& )]], "Output",
 CellChangeTimes->{3.937158399121942*^9},
 CellLabel->
  "Out[426]=",ExpressionUUID->"b276e45d-279a-4b09-8d1b-dd8859af208d"],

Cell[BoxData["\<\"maxratt.png\"\>"], "Output",
 CellChangeTimes->{3.937158399423337*^9},
 CellLabel->
  "Out[427]=",ExpressionUUID->"ba9ad951-a5a6-4746-8554-176cd2dea916"],

Cell[BoxData["\<\"minratt.png\"\>"], "Output",
 CellChangeTimes->{3.937158399456213*^9},
 CellLabel->
  "Out[428]=",ExpressionUUID->"2146ba53-3a5d-4377-ab4d-3da8f8fc97e4"],

Cell[BoxData["\<\"maxanglet.png\"\>"], "Output",
 CellChangeTimes->{3.937158399554262*^9},
 CellLabel->
  "Out[429]=",ExpressionUUID->"a246ddc7-5e3d-4eb5-b03a-8fb63440638e"],

Cell[BoxData["\<\"minanglet.png\"\>"], "Output",
 CellChangeTimes->{3.937158399588705*^9},
 CellLabel->
  "Out[430]=",ExpressionUUID->"8b0f12ef-dcd7-4c45-86a9-d08792f44b46"]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"endtime", " ", "//", " ", "Clear"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  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]", 
      RowBox[{
       RowBox[{"ellipseVecFromUnitVec", "[", 
        RowBox[{"uv_", ",", "ax1_", ",", " ", "ax2_"}], "]"}], ":=", 
       RowBox[{
        RowBox[{"1", "/", 
         RowBox[{"Sqrt", "[", 
          RowBox[{
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"Normalize", "[", "ax1", "]"}], ".", "uv"}], ")"}], "^",
              "2"}], "/", 
            RowBox[{
             RowBox[{"Norm", "[", "ax1", "]"}], "^", "2"}]}], "+", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"Normalize", "[", "ax2", "]"}], ".", "uv"}], ")"}], "^",
              "2"}], "/", 
            RowBox[{
             RowBox[{"Norm", "[", "ax2", "]"}], "^", "2"}]}]}], "]"}]}], "*", 
        "uv"}]}], ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
      RowBox[{"testbounds", " ", "=", " ", 
       RowBox[{"thrustreachable", "[", 
        RowBox[{
        "xrf", ",", " ", "terminaltime", ",", " ", "n", ",", " ", "thrustmax",
          ",", "numsamples"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"teststats", " ", "=", " ", 
       RowBox[{"reachablestatsoperator", "[", 
        RowBox[{"testbounds", ",", "ax1", ",", "ax2"}], " ", "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Print", "[", "testbounds", "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"Print", "[", "teststats", "]"}], ";", "\[IndentingNewLine]", 
      "\[IndentingNewLine]", 
      RowBox[{"flipTransformation", "=", 
       RowBox[{"ReflectionTransform", "[", 
        RowBox[{"{", 
         RowBox[{"1", ",", 
          RowBox[{"-", "1"}]}], "}"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"plot", " ", "=", "  ", 
       RowBox[{"Legended", "[", 
        RowBox[{
         RowBox[{"Show", "[", 
          RowBox[{
           RowBox[{"Graphics", "[", 
            RowBox[{"GeometricTransformation", "[", 
             RowBox[{
              RowBox[{"GeometricTransformation", "[", 
               RowBox[{
                RowBox[{"GeometricTransformation", "[", 
                 RowBox[{
                  RowBox[{"Circle", "[", 
                   RowBox[{
                    RowBox[{"{", 
                    RowBox[{"0", ",", "0"}], "}"}], ",", " ", 
                    RowBox[{"{", 
                    RowBox[{"ax1length", ",", " ", "ax2length"}], " ", 
                    "}"}]}], "]"}], ",", " ", 
                  RowBox[{"RotationTransform", "[", 
                   RowBox[{"ax1", " ", "//", " ", "vectoangle"}], "]"}]}], 
                 "]"}], ",", " ", "flipTransformation"}], "]"}], ",", " ", 
              RowBox[{"ScalingTransform", "[", 
               RowBox[{"{", 
                RowBox[{"0.001", ",", " ", "0.001"}], "}"}], "]"}]}], "]"}], 
            "]"}], ",", " ", 
           RowBox[{"Graphics", "[", 
            RowBox[{"GeometricTransformation", "[", 
             RowBox[{
              RowBox[{"GeometricTransformation", "[", 
               RowBox[{
                RowBox[{"Point", "[", 
                 RowBox[{"testbounds", "[", 
                  RowBox[{"[", 
                   RowBox[{";;", ",", " ", 
                    RowBox[{"1", " ", ";;", " ", "2"}]}], "]"}], "]"}], "]"}],
                 ",", "flipTransformation"}], "]"}], ",", 
              RowBox[{"ScalingTransform", "[", 
               RowBox[{"{", 
                RowBox[{"0.001", ",", " ", "0.001"}], "}"}], "]"}]}], "]"}], 
            "]"}], ",", " ", 
           RowBox[{"Graphics", "[", 
            RowBox[{"{", 
             RowBox[{"Red", ",", " ", 
              RowBox[{"GeometricTransformation", "[", 
               RowBox[{
                RowBox[{"Arrow", "[", 
                 RowBox[{"{", 
                  RowBox[{
                   RowBox[{"{", 
                    RowBox[{"0", ",", "0"}], "}"}], ",", 
                   RowBox[{"0.001", 
                    RowBox[{"ellipseVecFromUnitVec", "[", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"teststats", "[", 
                    RowBox[{"[", 
                    RowBox[{"2", ",", " ", "1"}], "]"}], "]"}], "/", 
                    RowBox[{"Norm", "[", 
                    RowBox[{"teststats", "[", 
                    RowBox[{"[", 
                    RowBox[{"2", ",", " ", "1"}], "]"}], "]"}], "]"}]}], ",", 
                    "ax1", ",", "ax2"}], "]"}]}]}], "}"}], " ", "]"}], ",", 
                " ", "flipTransformation"}], "]"}]}], "}"}], "]"}], ",", " ", 
           RowBox[{"Graphics", "[", 
            RowBox[{"{", 
             RowBox[{"Red", ",", " ", 
              RowBox[{"GeometricTransformation", "[", 
               RowBox[{
                RowBox[{"Arrow", "[", 
                 RowBox[{"{", 
                  RowBox[{
                   RowBox[{"{", 
                    RowBox[{"0", ",", "0"}], "}"}], ",", 
                   RowBox[{
                    RowBox[{"-", "0.001"}], " ", 
                    RowBox[{"ellipseVecFromUnitVec", "[", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"teststats", "[", 
                    RowBox[{"[", 
                    RowBox[{"2", ",", " ", "1"}], "]"}], "]"}], "/", 
                    RowBox[{"Norm", "[", 
                    RowBox[{"teststats", "[", 
                    RowBox[{"[", 
                    RowBox[{"2", ",", " ", "1"}], "]"}], "]"}], "]"}]}], ",", 
                    "ax1", ",", "ax2"}], "]"}]}]}], "}"}], " ", "]"}], ",", 
                " ", "flipTransformation"}], "]"}]}], "}"}], "]"}], ",", " ", 
           RowBox[{"Graphics", "[", 
            RowBox[{"{", 
             RowBox[{"Blue", ",", " ", 
              RowBox[{"GeometricTransformation", "[", 
               RowBox[{
                RowBox[{"Arrow", "[", 
                 RowBox[{"{", 
                  RowBox[{
                   RowBox[{"{", 
                    RowBox[{"0", ",", "0"}], "}"}], ",", 
                   RowBox[{"0.001", 
                    RowBox[{"ellipseVecFromUnitVec", "[", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"teststats", "[", 
                    RowBox[{"[", 
                    RowBox[{"2", ",", " ", "2"}], "]"}], "]"}], "/", 
                    RowBox[{"Norm", "[", 
                    RowBox[{"teststats", "[", 
                    RowBox[{"[", 
                    RowBox[{"2", ",", " ", "2"}], "]"}], "]"}], "]"}]}], ",", 
                    "ax1", ",", "ax2"}], "]"}]}]}], "}"}], " ", "]"}], ",", 
                " ", "flipTransformation"}], "]"}]}], "}"}], "]"}], ",", 
           "\[IndentingNewLine]", 
           RowBox[{"Graphics", "[", 
            RowBox[{"{", 
             RowBox[{"Blue", ",", " ", 
              RowBox[{"GeometricTransformation", "[", 
               RowBox[{
                RowBox[{"Arrow", "[", 
                 RowBox[{"{", 
                  RowBox[{
                   RowBox[{"{", 
                    RowBox[{"0", ",", "0"}], "}"}], ",", 
                   RowBox[{
                    RowBox[{"-", "0.001"}], " ", 
                    RowBox[{"ellipseVecFromUnitVec", "[", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"teststats", "[", 
                    RowBox[{"[", 
                    RowBox[{"2", ",", " ", "2"}], "]"}], "]"}], "/", 
                    RowBox[{"Norm", "[", 
                    RowBox[{"teststats", "[", 
                    RowBox[{"[", 
                    RowBox[{"2", ",", " ", "2"}], "]"}], "]"}], "]"}]}], ",", 
                    "ax1", ",", "ax2"}], "]"}]}]}], "}"}], " ", "]"}], ",", 
                " ", "flipTransformation"}], "]"}]}], "}"}], "]"}], ",", " ", 
           RowBox[{"Axes", " ", "->", " ", "True"}], ",", " ", 
           RowBox[{"ImagePadding", " ", "->", "30"}], ",", "   ", 
           RowBox[{
           "PlotLabel", " ", "->", 
            " ", "\"\<Energy Reachable and Thrust Limited RRS\>\""}], ",", 
           " ", 
           RowBox[{"PlotLegends", " ", "->", " ", 
            RowBox[{"{", 
             
             RowBox[{"\"\<Energy Limited\>\"", ",", 
              " ", "\"\<Thrust Limited\>\"", ",", " ", "\"\<Maximum Ratio\>\"",
               ",", " ", "\"\<Minimum Ratio\>\""}], " ", "}"}]}], ",", 
           RowBox[{"AxesLabel", " ", "->", " ", 
            RowBox[{"{", 
             RowBox[{"\"\<y\>\"", ",", " ", "\"\<x\>\""}], "}"}]}]}], "]"}], ",",
          " ", 
         RowBox[{"Placed", "[", 
          RowBox[{
           RowBox[{"LineLegend", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"Blue", ",", "Red"}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{"\"\<Minimum Ratio\>\"", ",", "\"\<Maximum Ratio\>\""}],
               "}"}]}], "]"}], ",", 
           RowBox[{"{", 
            RowBox[{"0.1", ",", "0.001"}], "}"}]}], "]"}]}], "]"}]}], " ", ";",
       "\[IndentingNewLine]", 
      RowBox[{"Export", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"StringForm", "[", 
          RowBox[{"\"\<RRSttt``.png\>\"", ",", " ", "terminaltime"}], "]"}], 
         " ", "//", " ", "ToString"}], ",", " ", "plot"}], "]"}]}]}], 
    "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]"}], "\[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.9265987712503443`*^9, 3.926598771674278*^9}, 3.9265988205016613`*^9, {
   3.92660091422416*^9, 3.926600918860721*^9}, {3.9266011504446893`*^9, 
   3.9266012301876326`*^9}, {3.9266016137751765`*^9, 
   3.9266016595420837`*^9}, {3.9266018504395485`*^9, 
   3.9266020268418493`*^9}, {3.9266024396534133`*^9, 3.926602443353138*^9}, {
   3.9266025737007895`*^9, 3.926602605291918*^9}, {3.9266027351799917`*^9, 
   3.926602807590149*^9}, {3.9266063470284295`*^9, 3.9266063571262865`*^9}, {
   3.9266066976790695`*^9, 3.9266067142850857`*^9}, {3.926606833694897*^9, 
   3.9266068396557217`*^9}, {3.9266069267239604`*^9, 3.926606927495826*^9}, {
   3.926607058935914*^9, 3.926607062600194*^9}, {3.926607103594391*^9, 
   3.9266071153492184`*^9}, {3.926607883781412*^9, 3.9266078872597804`*^9}, {
   3.926607952324091*^9, 3.9266079730465164`*^9}, {3.9266086640699444`*^9, 
   3.926608680903387*^9}, {3.926608818591425*^9, 3.926608821287423*^9}, {
   3.9266093619853697`*^9, 3.9266094062442007`*^9}, {3.9266841012580876`*^9, 
   3.926684131079776*^9}, {3.9266859391965756`*^9, 3.926685963836713*^9}, {
   3.927128493182778*^9, 3.927128511607149*^9}, {3.927136583714054*^9, 
   3.927136591418873*^9}, {3.927136623911099*^9, 3.927136635420822*^9}, {
   3.9271366740198507`*^9, 3.927136695320702*^9}, {3.927194721917962*^9, 
   3.927194748142041*^9}, {3.9271948359318743`*^9, 3.9271948473790283`*^9}, {
   3.92719495702982*^9, 3.927194965414035*^9}, {3.927195075077073*^9, 
   3.927195100265931*^9}, {3.92719518082654*^9, 3.9271951823286133`*^9}, {
   3.927195236951477*^9, 3.927195264410843*^9}, {3.927195307672089*^9, 
   3.927195389170733*^9}, {3.927195490252577*^9, 3.927195508039556*^9}, {
   3.9271955901228247`*^9, 3.927195628873547*^9}, {3.9271957013711567`*^9, 
   3.927195708022299*^9}, {3.9281736445675364`*^9, 3.928173645083197*^9}, {
   3.928173681549839*^9, 3.92817370269184*^9}, {3.928605106573565*^9, 
   3.928605107274837*^9}, {3.928606780432006*^9, 3.9286068075705347`*^9}, {
   3.928606852582657*^9, 3.928606881121591*^9}, {3.928606934175016*^9, 
   3.9286070032889013`*^9}, {3.928607360476371*^9, 3.928607502244911*^9}, {
   3.928608173061832*^9, 3.92860820983323*^9}, {3.928608409877784*^9, 
   3.92860844724023*^9}, {3.9286088558499393`*^9, 3.9286088599072104`*^9}, {
   3.928609086356181*^9, 3.928609091032599*^9}, {3.9286091338672657`*^9, 
   3.9286091339804983`*^9}, 3.92861012881286*^9, {3.928610570126472*^9, 
   3.9286105713918324`*^9}, {3.928611790008417*^9, 3.928611882406075*^9}, {
   3.928612485909569*^9, 3.928612541994854*^9}, {3.928612780878188*^9, 
   3.9286128005032454`*^9}, {3.928612907282607*^9, 3.928612967871997*^9}, {
   3.928614906349309*^9, 3.928614974157658*^9}, {3.928616736422194*^9, 
   3.92861676286455*^9}, {3.928616820435934*^9, 3.928616831656404*^9}, {
   3.928616926680834*^9, 3.928616984946121*^9}, {3.928617024240378*^9, 
   3.928617087675295*^9}, {3.9286179916642447`*^9, 3.928617993208478*^9}, 
   3.92861851528334*^9, 3.9286185580985737`*^9, {3.928618734480589*^9, 
   3.928618741122864*^9}, {3.928618774659534*^9, 3.928618809236704*^9}, 
   3.928618846953122*^9, 3.928618889276992*^9, {3.9286189345287657`*^9, 
   3.9286189375115347`*^9}, {3.928618997059957*^9, 3.928619000433593*^9}, {
   3.928619669200429*^9, 3.928619669966879*^9}, {3.928619710562914*^9, 
   3.928619717510044*^9}, {3.928622219075193*^9, 3.928622228914513*^9}, {
   3.928622950171659*^9, 3.928622954449678*^9}, {3.928623226164837*^9, 
   3.9286232289002657`*^9}, {3.9413868750741453`*^9, 
   3.941386892075762*^9}},ExpressionUUID->"c099018c-97a1-4afd-94a3-\
22d93d786922"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"plotrrs", "[", 
  RowBox[{
   RowBox[{"terminaltimes", "[", 
    RowBox[{"[", 
     RowBox[{"-", "1"}], "]"}], "]"}], ",", 
   RowBox[{"forms", "[", 
    RowBox[{"[", 
     RowBox[{"-", "1"}], "]"}], "]"}], ",", 
   RowBox[{"tmaxs", "[", 
    RowBox[{"[", 
     RowBox[{"-", "1"}], "]"}], "]"}], ",", 
   RowBox[{"ns", "[", 
    RowBox[{"[", 
     RowBox[{"-", "1"}], "]"}], "]"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.9254015616916523`*^9, 3.9254015907450657`*^9}, 
   3.925401756139267*^9, 3.925734541257496*^9, 3.92573494204535*^9, 
   3.926685712212511*^9},
 CellLabel->
  "In[433]:=",ExpressionUUID->"007c2eb7-74f7-497e-a06c-ad5f8bbcdc75"],

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, 433, 154, 
   31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.937158399630115*^9},
 CellLabel->
  "During evaluation of \
In[433]:=",ExpressionUUID->"468f85bc-ed4f-429e-9355-eba8ca92155c"],

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, 433, 155, 
   31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.937158399643194*^9},
 CellLabel->
  "During evaluation of \
In[433]:=",ExpressionUUID->"86cdf52e-8c35-4a46-8797-e8c0e721e0d4"],

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, 433, 156, 
   31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.937158399655025*^9},
 CellLabel->
  "During evaluation of \
In[433]:=",ExpressionUUID->"9a62e896-a71a-4637-acb1-a95bd08b753c"],

Cell[BoxData[
 TemplateBox[{
  "General", "stop", 
   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"NIntegrate\\\", \
\\\"::\\\", \\\"izero\\\"}], \\\"MessageName\\\"]\\) will be suppressed \
during this calculation.\"", 2, 433, 157, 31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.937158399659396*^9},
 CellLabel->
  "During evaluation of \
In[433]:=",ExpressionUUID->"21fd6517-6823-45fd-86f5-bd480fd2a5c8"],

Cell[BoxData["\<\"RRSttt6780.png\"\>"], "Output",
 CellChangeTimes->{3.937158399779728*^9},
 CellLabel->
  "Out[433]=",ExpressionUUID->"20e66daf-b3ad-4c0e-a911-d521396237c3"]
}, Open  ]],

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[434]:=",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, 434, 158, 
   31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.937158399806857*^9},
 CellLabel->
  "During evaluation of \
In[434]:=",ExpressionUUID->"5bafd69b-cf5f-4a8c-b9c0-3c2ce69523c3"],

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, 434, 159, 
   31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.937158399810587*^9},
 CellLabel->
  "During evaluation of \
In[434]:=",ExpressionUUID->"1925c5e2-697b-4b29-96ad-576763e087ab"],

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, 434, 160, 
   31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.937158399812841*^9},
 CellLabel->
  "During evaluation of \
In[434]:=",ExpressionUUID->"2b2a687a-d60f-4f69-b115-6c7185aeca23"],

Cell[BoxData[
 TemplateBox[{
  "General", "stop", 
   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"NIntegrate\\\", \
\\\"::\\\", \\\"izero\\\"}], \\\"MessageName\\\"]\\) will be suppressed \
during this calculation.\"", 2, 434, 161, 31036657779721760665, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.937158399815235*^9},
 CellLabel->
  "During evaluation of \
In[434]:=",ExpressionUUID->"43715a38-d119-4602-ab69-6239d981658e"],

Cell[BoxData["$Aborted"], "Output",
 CellChangeTimes->{3.937158407974533*^9},
 CellLabel->
  "Out[434]=",ExpressionUUID->"82439d97-3f56-4bd1-a006-8f616d5a5f20"]
}, 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[229]:=",ExpressionUUID->"c785bbf2-3350-4182-9c0e-5d65684a06da"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", 
  "\[IndentingNewLine]", "\[IndentingNewLine]", 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->CompressedData["
1:eJxTTMoPSmViYGAQB2IQfW2lquYRu7eOVqshtFrASqV0p7eOlQsfK4Po12Ef
LbKB9OywC44g+vvaUj8Q/eVsij+IfuH4NhhEa0x8HQqiLX7mZoLVsSwA04ql
PPkg+pHxt0IQ7fFnfxmIPmrNPRlEN1fUzQTR52LkFoJo5o5tJ0G04Zs9p0B0
Zq/VJRBd68N9FUQ7ZSY9BtG5n5++A9EGZ8+dywHSTfM23QTRz36wvAXRM62l
3oFoo91Wv0D0rhNO/0F0VBILby6Q3qCwjB9EH7GJUgDRS/ICVEB0Q5C4Nog+
eUMSTO/6vNoORN9iZbMH0VHhKp4gWn6yGphetqQ4GER3OLwG01pd25rzgbTu
Roc2EN1Qqz0dRIfJdYJpAKg8sS8=
  "],
 CellLabel->
  "In[230]:=",ExpressionUUID->"8659b878-ae33-4e5f-9107-fb888126cfd7"],

Cell[BoxData["  "], "Input",
 CellChangeTimes->{{3.930518070611895*^9, 3.9305180720724163`*^9}},
 CellLabel->
  "In[231]:=",ExpressionUUID->"5a22d7f4-b0be-4536-885b-55cc46f8052d"]
},
WindowSize->{Full, Full},
WindowMargins->{{0, Automatic}, {-22, Automatic}},
DockedCells->{},
TaggingRules-><|"TryRealOnly" -> False|>,
FrontEndVersion->"14.2 for Mac OS X ARM (64-bit) (December 26, 2024)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"9647f2d6-1130-4a31-9fb7-76e3ae5585a8"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 170, 3, 29, "Input",ExpressionUUID->"688a2e9f-7ca3-48dd-ac22-f5d7415eb864"],
Cell[753, 27, 168, 3, 33, "Output",ExpressionUUID->"bf6b0adc-9856-4441-b15f-85fb2df8cafb"]
}, Open  ]],
Cell[CellGroupData[{
Cell[958, 35, 900, 20, 70, "Input",ExpressionUUID->"344c9343-693b-47e9-b882-096618076cd6"],
Cell[1861, 57, 188, 3, 33, "Output",ExpressionUUID->"c8f475e4-98de-4808-8c45-b977a1e68ae0"]
}, Open  ]],
Cell[2064, 63, 11332, 286, 279, "Input",ExpressionUUID->"b895b316-708f-4092-a477-0f3a66f14cfd"],
Cell[13399, 351, 5072, 115, 519, "Input",ExpressionUUID->"38e6d68a-01ae-469f-8372-72472f6edb77"],
Cell[18474, 468, 18658, 457, 700, "Input",ExpressionUUID->"e8b154ed-e19c-47e3-a98d-7cc40f6a965d"],
Cell[37135, 927, 9960, 207, 233, "Input",ExpressionUUID->"c4ceae39-7397-4043-ac31-479218ca2bcb"],
Cell[CellGroupData[{
Cell[47120, 1138, 4600, 79, 213, "Input",ExpressionUUID->"d82cf4ae-f3bf-4a07-b427-7d6864e8cdbf"],
Cell[51723, 1219, 176, 3, 33, "Output",ExpressionUUID->"1a121e2d-c5a7-4194-a9e3-db1ccd73b003"],
Cell[51902, 1224, 169, 3, 33, "Output",ExpressionUUID->"921cd2a4-14e5-4da2-aab2-8f6af6b1b8de"]
}, Open  ]],
Cell[CellGroupData[{
Cell[52108, 1232, 4266, 69, 111, "Input",ExpressionUUID->"5bac8550-244b-4b24-a50f-7e5196fb9d83"],
Cell[56377, 1303, 508, 11, 26, "Message",ExpressionUUID->"7fa8c881-5945-4648-8706-f0ab69c4167f"],
Cell[56888, 1316, 508, 11, 26, "Message",ExpressionUUID->"e239c939-241c-451b-a455-84c7eba290ed"],
Cell[57399, 1329, 508, 11, 26, "Message",ExpressionUUID->"ddc48787-99d5-4410-adf4-c15e469a8fb5"],
Cell[57910, 1342, 461, 10, 26, "Message",ExpressionUUID->"6441a518-3f1b-4554-9b0f-ce650c807471"],
Cell[58374, 1354, 264, 6, 33, "Output",ExpressionUUID->"cb9fc03e-29f5-4456-9028-92af7216e525"]
}, Open  ]],
Cell[CellGroupData[{
Cell[58675, 1365, 6765, 122, 315, "Input",ExpressionUUID->"8d5d3b56-3d98-4aa2-94a6-05733121195f"],
Cell[65443, 1489, 1490, 22, 74, "Output",ExpressionUUID->"4f4fe017-2e4f-427e-aa32-efb1184e1464"],
Cell[66936, 1513, 508, 11, 26, "Message",ExpressionUUID->"67149ad9-5abe-4451-8c8f-de3b1e27ce8c"],
Cell[67447, 1526, 508, 11, 26, "Message",ExpressionUUID->"b6a72f8f-877e-4b71-b1eb-cb8bce01e97f"],
Cell[67958, 1539, 508, 11, 26, "Message",ExpressionUUID->"4eec3367-8fc3-45c8-89b6-79bb53979300"],
Cell[68469, 1552, 459, 10, 26, "Message",ExpressionUUID->"de6d00bf-0dd6-45c6-8ff3-bb644d13a8f7"]
}, Open  ]],
Cell[CellGroupData[{
Cell[68965, 1567, 9400, 244, 437, "Input",ExpressionUUID->"a7361479-b25e-4f65-bd51-39ad7bd317c6"],
Cell[78368, 1813, 426, 12, 33, "Output",ExpressionUUID->"58a102f0-0ea4-4016-ad5a-945abb3748d6"]
}, Open  ]],
Cell[78809, 1828, 8171, 157, 478, "Input",ExpressionUUID->"ce85b294-c6a9-4a49-8a9c-f7900fa64e64"],
Cell[CellGroupData[{
Cell[87005, 1989, 7741, 173, 274, "Input",ExpressionUUID->"79927019-60b3-4e35-99ba-21182e4ed03c"],
Cell[94749, 2164, 16074, 327, 188, "Output",ExpressionUUID->"d7200a32-ec83-4fa2-b57f-ea6b6ff8d6c1"],
Cell[110826, 2493, 16221, 329, 190, "Output",ExpressionUUID->"b276e45d-279a-4b09-8d1b-dd8859af208d"],
Cell[127050, 2824, 171, 3, 33, "Output",ExpressionUUID->"ba9ad951-a5a6-4746-8554-176cd2dea916"],
Cell[127224, 2829, 171, 3, 33, "Output",ExpressionUUID->"2146ba53-3a5d-4377-ab4d-3da8f8fc97e4"],
Cell[127398, 2834, 173, 3, 33, "Output",ExpressionUUID->"a246ddc7-5e3d-4eb5-b03a-8fb63440638e"],
Cell[127574, 2839, 173, 3, 33, "Output",ExpressionUUID->"8b0f12ef-dcd7-4c45-86a9-d08792f44b46"]
}, Open  ]],
Cell[127762, 2845, 19697, 420, 744, "Input",ExpressionUUID->"c099018c-97a1-4afd-94a3-22d93d786922"],
Cell[CellGroupData[{
Cell[147484, 3269, 672, 19, 29, "Input",ExpressionUUID->"007c2eb7-74f7-497e-a06c-ad5f8bbcdc75"],
Cell[148159, 3290, 508, 11, 26, "Message",ExpressionUUID->"468f85bc-ed4f-429e-9355-eba8ca92155c"],
Cell[148670, 3303, 508, 11, 26, "Message",ExpressionUUID->"86cdf52e-8c35-4a46-8797-e8c0e721e0d4"],
Cell[149181, 3316, 508, 11, 26, "Message",ExpressionUUID->"9a62e896-a71a-4637-acb1-a95bd08b753c"],
Cell[149692, 3329, 459, 10, 26, "Message",ExpressionUUID->"21fd6517-6823-45fd-86f5-bd480fd2a5c8"],
Cell[150154, 3341, 174, 3, 33, "Output",ExpressionUUID->"20e66daf-b3ad-4c0e-a911-d521396237c3"]
}, Open  ]],
Cell[CellGroupData[{
Cell[150365, 3349, 3518, 54, 49, "Input",ExpressionUUID->"6ad10b9c-46d8-4f63-9072-e5fba7ca1875"],
Cell[153886, 3405, 508, 11, 26, "Message",ExpressionUUID->"5bafd69b-cf5f-4a8c-b9c0-3c2ce69523c3"],
Cell[154397, 3418, 508, 11, 26, "Message",ExpressionUUID->"1925c5e2-697b-4b29-96ad-576763e087ab"],
Cell[154908, 3431, 508, 11, 26, "Message",ExpressionUUID->"2b2a687a-d60f-4f69-b115-6c7185aeca23"],
Cell[155419, 3444, 459, 10, 26, "Message",ExpressionUUID->"43715a38-d119-4602-ab69-6239d981658e"],
Cell[155881, 3456, 160, 3, 33, "Output",ExpressionUUID->"82439d97-3f56-4bd1-a006-8f616d5a5f20"]
}, Open  ]],
Cell[156056, 3462, 321, 5, 29, "Input",ExpressionUUID->"c785bbf2-3350-4182-9c0e-5d65684a06da"],
Cell[156380, 3469, 699, 14, 151, "Input",ExpressionUUID->"8659b878-ae33-4e5f-9107-fb888126cfd7"],
Cell[157082, 3485, 179, 3, 29, "Input",ExpressionUUID->"5a22d7f4-b0be-4536-885b-55cc46f8052d"]
}
]
*)