PACAGetModes.nb

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

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

(* CreatedBy='Mathematica 11.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    124964,       3107]
NotebookOptionsPosition[    119748,       3011]
NotebookOutlinePosition[    120105,       3027]
CellTagsIndexPosition[    120062,       3024]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"initialSeedPoly", "[", 
    RowBox[{"exp_", ",", "x_"}], "]"}], ":=", 
   RowBox[{"Expand", "[", 
    RowBox[{"exp", "/.", 
     RowBox[{"{", 
      RowBox[{"x", "\[Rule]", 
       RowBox[{"(", 
        RowBox[{"x", "-", 
         RowBox[{
          RowBox[{"Coefficient", "[", 
           RowBox[{"exp", ",", 
            RowBox[{"x", "^", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"Exponent", "[", 
                RowBox[{"exp", ",", "x"}], "]"}], "-", "1"}], ")"}]}]}], 
           "]"}], "/", 
          RowBox[{"Exponent", "[", 
           RowBox[{"exp", ",", "x"}], "]"}]}]}], ")"}]}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"stageTwoSeedPoly", "[", 
    RowBox[{"exp_", ",", "x_"}], "]"}], ":=", 
   RowBox[{"Expand", "[", 
    RowBox[{
     RowBox[{"x", "^", 
      RowBox[{"Exponent", "[", 
       RowBox[{"exp", ",", "x"}], "]"}]}], "+", 
     RowBox[{"Sum", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"-", 
         RowBox[{"Abs", "[", 
          RowBox[{"Coefficient", "[", 
           RowBox[{"exp", ",", 
            RowBox[{"x", "^", "i"}]}], "]"}], "]"}]}], 
        RowBox[{"x", "^", "i"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", 
         RowBox[{
          RowBox[{"Exponent", "[", 
           RowBox[{"exp", ",", "x"}], "]"}], "-", "2"}]}], "}"}]}], "]"}], 
     "-", 
     RowBox[{"Abs", "[", 
      RowBox[{"exp", "-", 
       RowBox[{"Sum", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"Coefficient", "[", 
           RowBox[{"exp", ",", 
            RowBox[{"x", "^", "i"}]}], "]"}], 
          RowBox[{"x", "^", "i"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"i", ",", "1", ",", 
           RowBox[{"Exponent", "[", 
            RowBox[{"exp", ",", "x"}], "]"}]}], "}"}]}], "]"}]}], "]"}]}], 
    "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"initialSeed", "[", 
    RowBox[{"exp_", ",", "x_"}], "]"}], ":=", "\[IndentingNewLine]", 
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"Clear", "[", 
       RowBox[{"n", ",", "z", ",", "exp1", ",", "i", ",", "c1", ",", "flag"}],
        "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"i", "=", "0"}], ";", "\[IndentingNewLine]", 
      RowBox[{"n", " ", "=", 
       RowBox[{"Exponent", "[", 
        RowBox[{"exp", ",", "x"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"c1", "=", 
       RowBox[{"Coefficient", "[", 
        RowBox[{"exp", ",", 
         RowBox[{"x", "^", 
          RowBox[{"(", 
           RowBox[{"n", "-", "1"}], ")"}]}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"exp1", "=", 
       RowBox[{"stageTwoSeedPoly", "[", 
        RowBox[{
         RowBox[{"initialSeedPoly", "[", 
          RowBox[{"exp", ",", "x"}], "]"}], ",", "x"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"While", "[", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"Not", "[", 
         RowBox[{"And", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"exp1", "/.", 
              RowBox[{"x", "\[Rule]", "i"}]}], ")"}], "\[LessEqual]", "0"}], 
           ",", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"exp1", "/.", 
              RowBox[{"x", "\[Rule]", 
               RowBox[{"i", "+", "1"}]}]}], ")"}], "\[GreaterEqual]", "0"}]}],
           "]"}], "]"}], ",", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"i", "=", 
          RowBox[{"i", "+", "1"}]}], ";"}]}], "\[IndentingNewLine]", "]"}], 
      ";", "\[IndentingNewLine]", 
      RowBox[{"Return", "[", 
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"-", "c1"}], "/", "n"}], "+", 
          RowBox[{
           RowBox[{"Round", "[", 
            RowBox[{"i", "+", "1"}], "]"}], 
           RowBox[{"Exp", "[", 
            RowBox[{"I", "*", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"2", 
                RowBox[{"Pi", "/", "n"}], "*", 
                RowBox[{"(", 
                 RowBox[{"k", "-", "1"}], ")"}]}], "+", 
               RowBox[{
                RowBox[{"Pi", "/", "2"}], "/", "n"}]}], ")"}]}], "]"}]}]}], 
         ",", 
         RowBox[{"{", 
          RowBox[{"k", ",", "1", ",", "n"}], "}"}]}], "]"}], "]"}], ";"}], 
     "\[IndentingNewLine]", ",", "1"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"PACARootFinder", "[", 
    RowBox[{"exp_", ",", "x_"}], "]"}], ":=", 
   RowBox[{"If", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"Exponent", "[", 
       RowBox[{"exp", ",", "x"}], "]"}], "\[Equal]", "1"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"-", 
       RowBox[{"(", 
        RowBox[{"exp", "-", 
         RowBox[{"x", " ", 
          RowBox[{"Coefficient", "[", 
           RowBox[{"exp", ",", "x"}], "]"}]}]}], ")"}]}], "/", 
      RowBox[{"Coefficient", "[", 
       RowBox[{"exp", ",", "x"}], "]"}]}], ",", "\[IndentingNewLine]", 
     RowBox[{"Do", "[", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{
        RowBox[{"Clear", "[", 
         RowBox[{"n1", ",", "z", ",", "z1", ",", "i", ",", "flag"}], "]"}], 
        ";", "\[IndentingNewLine]", 
        RowBox[{"n1", "=", 
         RowBox[{"Exponent", "[", 
          RowBox[{"exp", ",", "x"}], "]"}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"z", " ", "=", 
         RowBox[{"N", "[", 
          RowBox[{"initialSeed", "[", 
           RowBox[{"exp", ",", "x"}], "]"}], "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"While", "[", 
         RowBox[{"True", ",", "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"z1", " ", "=", " ", 
            RowBox[{"Table", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"z", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "+", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{"exp", "/.", 
                  RowBox[{"x", "\[Rule]", 
                   RowBox[{"z", "[", 
                    RowBox[{"[", "i", "]"}], "]"}]}]}], ")"}], "/", 
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"(", 
                    RowBox[{"exp", "/.", 
                    RowBox[{"x", "\[Rule]", 
                    RowBox[{"z", "[", 
                    RowBox[{"[", "i", "]"}], "]"}]}]}], ")"}], "*", 
                   RowBox[{"Sum", "[", 
                    RowBox[{
                    RowBox[{"If", "[", 
                    RowBox[{
                    RowBox[{"k", "\[NotEqual]", "i"}], ",", 
                    RowBox[{"1", "/", 
                    RowBox[{"(", 
                    RowBox[{
                    RowBox[{"z", "[", 
                    RowBox[{"[", "i", "]"}], "]"}], "-", 
                    RowBox[{"z", "[", 
                    RowBox[{"[", "k", "]"}], "]"}]}], ")"}]}], ",", "0"}], 
                    "]"}], ",", 
                    RowBox[{"{", 
                    RowBox[{"k", ",", "1", ",", "n1"}], "}"}]}], "]"}]}], "-", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"D", "[", 
                    RowBox[{"exp", ",", "x"}], "]"}], "/.", 
                    RowBox[{"x", "\[Rule]", 
                    RowBox[{"z", "[", 
                    RowBox[{"[", "i", "]"}], "]"}]}]}], ")"}]}], ")"}]}]}], 
              ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", "1", ",", "n1"}], "}"}]}], "]"}]}], ";", 
           "\[IndentingNewLine]", 
           RowBox[{"flag", "=", "True"}], ";", "\[IndentingNewLine]", 
           RowBox[{"For", "[", 
            RowBox[{
             RowBox[{"i", "=", "1"}], ",", 
             RowBox[{"i", "\[LessEqual]", "n1"}], ",", 
             RowBox[{"i", "++"}], ",", 
             RowBox[{"flag", "=", 
              RowBox[{"And", "[", 
               RowBox[{"flag", ",", 
                RowBox[{
                 RowBox[{"Abs", "[", 
                  RowBox[{
                   RowBox[{"z", "[", 
                    RowBox[{"[", "i", "]"}], "]"}], "-", 
                   RowBox[{"z1", "[", 
                    RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "<", 
                 RowBox[{"10", "^", 
                  RowBox[{"(", 
                   RowBox[{"-", "5"}], ")"}]}]}]}], "]"}]}]}], "]"}], ";", 
           "\[IndentingNewLine]", 
           RowBox[{"If", "[", 
            RowBox[{"flag", ",", 
             RowBox[{"Break", "[", "]"}]}], "]"}], ";", "\[IndentingNewLine]", 
           RowBox[{"For", "[", 
            RowBox[{
             RowBox[{"i", "=", "1"}], ",", 
             RowBox[{"i", "\[LessEqual]", "n1"}], ",", 
             RowBox[{"i", "++"}], ",", 
             RowBox[{
              RowBox[{"z", "[", 
               RowBox[{"[", "i", "]"}], "]"}], "=", 
              RowBox[{"z1", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}], ";"}]}], 
         "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
        RowBox[{"Return", "[", "z1", "]"}], ";"}], "\[IndentingNewLine]", ",",
        "1"}], "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"NewPACARootFinder", "[", 
    RowBox[{"exp_", ",", "x_"}], "]"}], ":=", "\[IndentingNewLine]", 
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"Clear", "[", 
       RowBox[{
       "n", ",", "cp", ",", "exp1", ",", "exp2", ",", "exp3", ",", "c0"}], 
       "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"n", "=", 
       RowBox[{"Exponent", "[", 
        RowBox[{"exp", ",", "x"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"cp", "=", 
       RowBox[{"Coefficient", "[", 
        RowBox[{"exp", ",", 
         RowBox[{"x", "^", "n"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"exp1", "=", 
       RowBox[{"N", "[", 
        RowBox[{"Expand", "[", 
         RowBox[{"exp", "/", "cp"}], "]"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"c0", "=", 
       RowBox[{"exp1", "/.", 
        RowBox[{"x", "\[Rule]", "0"}]}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"(*", 
       RowBox[{"Refining", " ", "the", " ", "coefficients"}], "*)"}], 
      "\[IndentingNewLine]", 
      RowBox[{"exp2", "=", 
       RowBox[{
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Abs", "[", "c0", "]"}], "<", "0.001"}], ",", 
          RowBox[{
           RowBox[{
            RowBox[{"Sign", "[", 
             RowBox[{"Re", "[", "c0", "]"}], "]"}], "0.001"}], "+", 
           RowBox[{"I", " ", 
            RowBox[{"Sign", "[", 
             RowBox[{"Im", "[", "c0", "]"}], "]"}], "0.001"}]}], ",", "c0"}], 
         "]"}], "+", 
        RowBox[{"Sum", "[", " ", 
         RowBox[{
          RowBox[{
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{
              RowBox[{"Abs", "[", 
               RowBox[{"Coefficient", "[", 
                RowBox[{"exp1", ",", 
                 RowBox[{"x", "^", "i"}]}], "]"}], "]"}], "<", "0.001"}], ",", 
             RowBox[{
              RowBox[{
               RowBox[{"Sign", "[", 
                RowBox[{"Re", "[", 
                 RowBox[{"Coefficient", "[", 
                  RowBox[{"exp1", ",", 
                   RowBox[{"x", "^", "i"}]}], "]"}], "]"}], "]"}], "0.001"}], 
              "+", 
              RowBox[{"I", " ", 
               RowBox[{"Sign", "[", 
                RowBox[{"Im", "[", 
                 RowBox[{"Coefficient", "[", 
                  RowBox[{"exp1", ",", 
                   RowBox[{"x", "^", "i"}]}], "]"}], "]"}], "]"}], 
               "0.001"}]}], ",", 
             RowBox[{"Coefficient", "[", 
              RowBox[{"exp1", ",", 
               RowBox[{"x", "^", "i"}]}], "]"}]}], " ", "]"}], "*", 
           RowBox[{"x", "^", "i"}]}], ",", 
          RowBox[{"{", 
           RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"cp", "=", 
       RowBox[{"Coefficient", "[", 
        RowBox[{"exp2", ",", 
         RowBox[{"x", "^", "n"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"exp3", "=", 
       RowBox[{"SetPrecision", "[", 
        RowBox[{
         RowBox[{"N", "[", 
          RowBox[{"Expand", "[", 
           RowBox[{"exp2", "/", "cp"}], "]"}], "]"}], ",", "2"}], "]"}]}], 
      ";", "\[IndentingNewLine]", 
      RowBox[{"Return", "[", 
       RowBox[{"PACARootFinder", "[", 
        RowBox[{"exp3", ",", "x"}], "]"}], "]"}]}], "\[IndentingNewLine]", 
     ",", "1"}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.744727208433558*^9, 3.744727216081555*^9}, {
   3.7447273513389416`*^9, 3.7447274616475663`*^9}, {3.7447274948862286`*^9, 
   3.744727631521207*^9}, 3.744727786298504*^9, 3.7447279313542523`*^9, {
   3.7447279913114758`*^9, 3.744728321824041*^9}, {3.7447283556618853`*^9, 
   3.7447284049985685`*^9}, {3.7447284678856754`*^9, 
   3.7447287248889847`*^9}, {3.744728793279539*^9, 3.744729058064641*^9}, {
   3.7447291321308684`*^9, 3.74472913226149*^9}, {3.7447294133684244`*^9, 
   3.7447294728260126`*^9}, {3.744729518248832*^9, 3.744729575640523*^9}, {
   3.744729646801509*^9, 3.7447296671752625`*^9}, {3.744729699227769*^9, 
   3.744729754132455*^9}, {3.7447297845161533`*^9, 3.7447298339148645`*^9}, {
   3.744729872359327*^9, 3.744729876635719*^9}, {3.7447299595476007`*^9, 
   3.74472997676904*^9}, {3.744730028332126*^9, 3.744730033075099*^9}, 
   3.7447301547617216`*^9, 3.7447302332245646`*^9, {3.744730264523342*^9, 
   3.744730265481803*^9}, {3.74478764622153*^9, 3.7447876477086926`*^9}, {
   3.744787679296955*^9, 3.7447876850929155`*^9}, {3.744788576285486*^9, 
   3.744788638854446*^9}, {3.7448068459179335`*^9, 3.7448068687116814`*^9}, {
   3.744806912257843*^9, 3.7448069586225233`*^9}, {3.7448070338655405`*^9, 
   3.744807037340858*^9}, {3.7448087810164776`*^9, 3.744808948945901*^9}, {
   3.7448090637327843`*^9, 3.7448090731764555`*^9}, {3.744809109015911*^9, 
   3.7448091133766575`*^9}, 3.7448091709846973`*^9, {3.744809217715214*^9, 
   3.7448092312385807`*^9}, {3.744809308780244*^9, 3.7448093161710186`*^9}, {
   3.744809546877994*^9, 3.7448095504487453`*^9}, {3.7448097541074905`*^9, 
   3.7448097751546373`*^9}, {3.7448098524034915`*^9, 3.744809903346184*^9}, {
   3.7448099545658402`*^9, 3.744809961015873*^9}, 3.7448100228756695`*^9, {
   3.7448101321324153`*^9, 3.7448102263398247`*^9}, {3.744813838762951*^9, 
   3.7448138879210005`*^9}, {3.744813963367264*^9, 3.7448139642789497`*^9}, 
   3.744814044305109*^9, {3.7448140898493147`*^9, 3.744814144635936*^9}, {
   3.744814279648568*^9, 3.74481428832794*^9}, {3.74481432620319*^9, 
   3.7448143310768094`*^9}, {3.744814439035221*^9, 3.744814633715416*^9}, {
   3.7448146745792785`*^9, 3.7448146792842426`*^9}, {3.7448147434291844`*^9, 
   3.7448148161694813`*^9}, {3.744814886787711*^9, 3.7448148886007423`*^9}, {
   3.7448150569518485`*^9, 3.7448150867311707`*^9}, {3.7448151369313364`*^9, 
   3.744815209641858*^9}, {3.7448152413140917`*^9, 3.744815243806426*^9}, {
   3.744815291088516*^9, 3.7448152922757187`*^9}, {3.744815331105657*^9, 
   3.74481534502633*^9}, {3.744815423518159*^9, 3.7448154345060797`*^9}, {
   3.744815478719944*^9, 3.7448155065343714`*^9}, {3.7448155879204807`*^9, 
   3.744815599321472*^9}, {3.7448161613677063`*^9, 3.7448161908890676`*^9}, {
   3.7448162437366114`*^9, 3.744816250336466*^9}, {3.744816324816804*^9, 
   3.7448164053166513`*^9}, {3.7448165279891653`*^9, 3.744816528226059*^9}, {
   3.7449592246509447`*^9, 3.7449592344422045`*^9}, {3.7449592968030186`*^9, 
   3.744959298395156*^9}, {3.744959356958131*^9, 3.7449593606114817`*^9}, 
   3.7449594025552597`*^9, {3.744987556084035*^9, 3.7449875607841263`*^9}, {
   3.7449877006626797`*^9, 3.744987717147438*^9}, {3.7449881276897063`*^9, 
   3.7449882347084885`*^9}, {3.7449883678122954`*^9, 
   3.7449883774841366`*^9}, {3.747993700260203*^9, 3.7479937011193542`*^9}, {
   3.7479937404218664`*^9, 3.747993741186053*^9}, {3.7479942583388667`*^9, 
   3.747994259693874*^9}, {3.747994303473939*^9, 3.7479943158003445`*^9}, {
   3.7479946777573447`*^9, 3.7479946855407887`*^9}, {3.747994727745816*^9, 
   3.7479947368214784`*^9}, {3.74799539788995*^9, 
   3.747995403791401*^9}},ExpressionUUID->"c618399f-f8df-4a55-9ad5-\
ffac92a4333b"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"NewPACARootFinder", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"5", "*", 
     RowBox[{"x", "^", "2"}]}], "-", "2"}], ",", "x"}], "]"}]], "Input",
 CellChangeTimes->{{3.7479942637547183`*^9, 3.747994282258483*^9}, 
   3.747994322807146*^9, {3.747994698862751*^9, 
   3.7479947029233594`*^9}},ExpressionUUID->"24262446-3bf0-494d-b923-\
b6e34627688b"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"0.6324555320336759`", "\[VeryThinSpace]", "+", 
    RowBox[{"2.0679515313825692`*^-25", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "0.6324555320336759`"}], "-", 
    RowBox[{"2.0679515313825692`*^-25", " ", "\[ImaginaryI]"}]}]}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.7479943014825354`*^9, 3.7479943235315733`*^9}, 
   3.7479944868444586`*^9, {3.7479946931340866`*^9, 3.7479947073702173`*^9}, 
   3.7479947390053167`*^9},ExpressionUUID->"bc7b092e-409d-4e85-88ec-\
d3a80d8900c8"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"NewPACARootFinder", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"5", "*", 
     RowBox[{"x", "^", "6.1"}]}], "-", "2.3"}], ",", "x"}], "]"}]], "Input",
 CellChangeTimes->{{3.7479944881774535`*^9, 
  3.7479944943048816`*^9}},ExpressionUUID->"84759498-767a-4843-818e-\
86b0b9fd7ea5"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"0.8804696668472677`", "\[VeryThinSpace]", "-", 
    RowBox[{"3.579650570602829`*^-18", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.4532594464407327`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.7548395249656574`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "0.41380026635921197`"}], "+", 
    RowBox[{"0.7771719074948474`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "0.879302243209895`"}], "-", 
    RowBox[{"0.04532548206235946`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "0.41380026635922806`"}], "-", 
    RowBox[{"0.7771719074948416`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.45325944644073235`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.7548395249656575`", " ", "\[ImaginaryI]"}]}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.7479946899951105`*^9, 
  3.747994742071562*^9},ExpressionUUID->"57c39d3a-ebb9-49c3-b75f-\
81d574c074b7"]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"xp", "[", 
    RowBox[{"exp_", ",", "x_"}], "]"}], ":=", 
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"Clear", "[", 
       RowBox[{"exp1", ",", "n", ",", "slnroot", ",", "i", ",", "m"}], "]"}], 
      ";", "\[IndentingNewLine]", 
      RowBox[{"exp1", "=", 
       RowBox[{"Numerator", "[", 
        RowBox[{"Together", "[", "exp", "]"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"n", "=", 
       RowBox[{"Exponent", "[", 
        RowBox[{"exp1", ",", "x"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"slnroot", " ", "=", " ", 
       RowBox[{"NewPACARootFinder", "[", 
        RowBox[{"exp1", ",", "x"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"m", " ", "=", " ", "0"}], ";", "\[IndentingNewLine]", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"i", "=", "1"}], ",", 
        RowBox[{"i", "\[LessEqual]", "n"}], ",", 
        RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"r", "=", 
          RowBox[{"Re", "[", 
           RowBox[{"slnroot", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "]"}]}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{"im", "=", 
          RowBox[{"Abs", "[", 
           RowBox[{"Im", "[", 
            RowBox[{"slnroot", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "]"}], "]"}]}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"im", "<", 
             RowBox[{"10", "^", 
              RowBox[{"(", 
               RowBox[{"-", "3"}], ")"}]}]}], "&&", 
            RowBox[{"m", "<", "r"}]}], ",", 
           RowBox[{"m", " ", "=", "r"}]}], "]"}], ";"}]}], 
       "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"Return", "[", "m", "]"}]}], "\[IndentingNewLine]", ",", "1"}], 
    "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"padZero", "[", 
    RowBox[{"S_", ",", "j_"}], "]"}], ":=", "\[IndentingNewLine]", 
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"Clear", "[", 
       RowBox[{"k", ",", "l", ",", "y"}], "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"Length", "[", "S", "]"}], "<", 
         RowBox[{"j", "+", "2"}]}], ",", "\[IndentingNewLine]", 
        RowBox[{"Do", "[", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           RowBox[{"y", "=", "S"}], ";", "\[IndentingNewLine]", 
           RowBox[{"l", "=", 
            RowBox[{"Length", "[", "S", "]"}]}], ";", "\[IndentingNewLine]", 
           RowBox[{"For", "[", 
            RowBox[{
             RowBox[{"k", "=", "1"}], ",", 
             RowBox[{"k", "\[LessEqual]", 
              RowBox[{"j", "+", "2", "-", "l"}]}], ",", 
             RowBox[{"k", "++"}], ",", 
             RowBox[{"AppendTo", "[", 
              RowBox[{"y", ",", "0"}], "]"}]}], "]"}], ";"}], 
          "\[IndentingNewLine]", ",", "1"}], "]"}]}], "\[IndentingNewLine]", 
       "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"Return", "[", "y", "]"}], ";"}], "\[IndentingNewLine]", ",", 
     "1"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"aExpression", "[", 
    RowBox[{
    "n_", ",", "exp_", ",", "x_", ",", "W_", ",", "\[Omega]_", ",", 
     "\[Beta]_"}], "]"}], ":=", "\[IndentingNewLine]", 
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"Clear", "[", 
       RowBox[{
       "S", ",", "T", ",", "U", ",", "k", ",", "root", ",", "dexp", ",", "p", 
        ",", "q", ",", "r", ",", "R"}], "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"root", " ", "=", 
       RowBox[{"xp", "[", 
        RowBox[{"exp", ",", "x"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"dexp", "=", 
       RowBox[{"D", "[", 
        RowBox[{"exp", ",", "x"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"R", "=", 
       RowBox[{"Table", "[", 
        RowBox[{"1", ",", 
         RowBox[{"{", 
          RowBox[{"k", ",", "1", ",", 
           RowBox[{"n", "+", "1"}]}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"S", "=", 
       RowBox[{"N", "[", 
        RowBox[{"CoefficientList", "[", 
         RowBox[{
          RowBox[{"Series", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"x", "^", "4"}], " ", "exp"}], ",", 
            RowBox[{"{", 
             RowBox[{"x", ",", "root", ",", "n"}], "}"}]}], "]"}], ",", "x"}],
          "]"}], "]"}]}], ";", 
      RowBox[{"S1", "=", 
       RowBox[{"padZero", "[", 
        RowBox[{"S", ",", "n"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"T", "=", 
       RowBox[{"N", "[", 
        RowBox[{"CoefficientList", "[", 
         RowBox[{
          RowBox[{"Series", "[", 
           RowBox[{
            RowBox[{
             RowBox[{
              RowBox[{"-", "2"}], " ", "exp", " ", "root", " ", 
              SuperscriptBox["x", "3"]}], "+", 
             RowBox[{"2", " ", "exp", " ", 
              SuperscriptBox["x", "4"]}], "-", 
             RowBox[{"dexp", " ", "root", " ", 
              SuperscriptBox["x", "4"]}], "+", 
             RowBox[{"dexp", " ", 
              SuperscriptBox["x", "5"]}], "-", 
             RowBox[{
             "2", " ", "\[ImaginaryI]", " ", "root", " ", "x", " ", 
              "\[Omega]"}], "+", 
             RowBox[{"2", " ", "\[ImaginaryI]", " ", 
              SuperscriptBox["x", "2"], " ", "\[Omega]"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"x", ",", "root", ",", "n"}], "}"}]}], "]"}], ",", "x"}],
          "]"}], "]"}]}], ";", 
      RowBox[{"T1", "=", 
       RowBox[{"padZero", "[", 
        RowBox[{"T", ",", "n"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"U", "=", 
       RowBox[{"N", "[", 
        RowBox[{"CoefficientList", "[", 
         RowBox[{
          RowBox[{"Series", "[", 
           RowBox[{
            RowBox[{
             RowBox[{
              RowBox[{"-", 
               SuperscriptBox["root", "2"]}], " ", "W"}], "-", 
             RowBox[{
              SuperscriptBox["root", "2"], " ", 
              SuperscriptBox["W", "2"]}], "+", 
             RowBox[{"2", " ", "root", " ", "W", " ", "x"}], "+", 
             RowBox[{"2", " ", "root", " ", 
              SuperscriptBox["W", "2"], " ", "x"}], "-", 
             RowBox[{"W", " ", 
              SuperscriptBox["x", "2"]}], "-", 
             RowBox[{
              SuperscriptBox["W", "2"], " ", 
              SuperscriptBox["x", "2"]}], "+", 
             RowBox[{"dexp", " ", 
              SuperscriptBox["root", "2"], " ", 
              SuperscriptBox["x", "5"], " ", "\[Beta]"}], "-", 
             RowBox[{"2", " ", "dexp", " ", "root", " ", 
              SuperscriptBox["x", "6"], " ", "\[Beta]"}], "+", 
             RowBox[{"dexp", " ", 
              SuperscriptBox["x", "7"], " ", "\[Beta]"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"x", ",", "root", ",", "n"}], "}"}]}], "]"}], ",", "x"}],
          "]"}], "]"}]}], ";", 
      RowBox[{"U1", "=", 
       RowBox[{"padZero", "[", 
        RowBox[{"U", ",", "n"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"q", "=", 
       RowBox[{"S1", "[", 
        RowBox[{"[", "1", "]"}], "]"}]}], ";", 
      RowBox[{"r", "=", 
       RowBox[{"T1", "[", 
        RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"F", "=", 
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{"Table", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"l", 
             RowBox[{"(", 
              RowBox[{"l", "-", "1"}], ")"}], 
             RowBox[{"S1", "[", 
              RowBox[{"[", 
               RowBox[{"k", "-", "l", "+", "1"}], "]"}], "]"}]}], "+", 
            RowBox[{"l", " ", 
             RowBox[{"T1", "[", 
              RowBox[{"[", 
               RowBox[{"k", "-", "l", "+", "1"}], "]"}], "]"}]}], "+", 
            RowBox[{"U1", "[", 
             RowBox[{"[", 
              RowBox[{"k", "-", "l", "+", "1"}], "]"}], "]"}]}], ",", 
           RowBox[{"{", 
            RowBox[{"l", ",", "1", ",", 
             RowBox[{"k", "-", "1"}]}], "}"}]}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"k", ",", "1", ",", 
           RowBox[{"n", "+", "1"}]}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"p", "=", 
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "1"}], "/", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"k", "*", 
             RowBox[{"(", 
              RowBox[{"k", "-", "1"}], ")"}], "q"}], "+", 
            RowBox[{"k", " ", "r"}]}], " ", ")"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"k", ",", "2", ",", 
           RowBox[{"n", "+", "1"}]}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"k", "=", "2"}], ",", 
        RowBox[{"k", "\[LessEqual]", 
         RowBox[{"n", "+", "1"}]}], ",", 
        RowBox[{"k", "++"}], ",", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{
          RowBox[{"R", "[", 
           RowBox[{"[", "k", "]"}], "]"}], "=", " ", 
          RowBox[{"Sum", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"R", "[", 
              RowBox[{"[", "l", "]"}], "]"}], 
             RowBox[{
              RowBox[{"F", "[", 
               RowBox[{"[", "k", "]"}], "]"}], "[", 
              RowBox[{"[", "l", "]"}], "]"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"l", ",", "1", ",", 
              RowBox[{"k", "-", "1"}]}], "}"}]}], "]"}]}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"R", "[", 
           RowBox[{"[", "k", "]"}], "]"}], "=", 
          RowBox[{"Simplify", "[", 
           RowBox[{
            RowBox[{"p", " ", "[", 
             RowBox[{"[", 
              RowBox[{"k", "-", "1"}], "]"}], "]"}], 
            RowBox[{"R", "[", 
             RowBox[{"[", "k", "]"}], "]"}]}], "]"}]}], ";"}]}], 
       "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"Return", "[", "R", "]"}]}], "\[IndentingNewLine]", ",", "1"}], 
    "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"PACAGetModes", "[", 
    RowBox[{"m_", ",", "exp_", ",", "x_", ",", "W_", ",", "\[Beta]_"}], "]"}],
    ":=", "\[IndentingNewLine]", 
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"Clear", "[", 
       RowBox[{
       "exp1", ",", "CoeffList", ",", "RootCalc", ",", "\[Omega]", ",", 
        "exp1", ",", "s", ",", "root", ",", "r"}], "]"}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"root", "=", 
       RowBox[{"xp", "[", 
        RowBox[{"exp", ",", "x"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"CoeffList", "=", 
       RowBox[{"aExpression", "[", 
        RowBox[{
        "m", ",", "exp", ",", "x", ",", "W", ",", "\[Omega]", ",", 
         "\[Beta]"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"RootCalc", "=", 
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{"-", "root"}], ")"}], "^", 
          RowBox[{"(", 
           RowBox[{"n", "-", "1"}], ")"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"n", ",", "1", ",", 
           RowBox[{"m", "+", "1"}]}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"exp1", "=", 
       RowBox[{"Sum", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"CoeffList", "[", 
           RowBox[{"[", "n", "]"}], "]"}], 
          RowBox[{"RootCalc", "[", 
           RowBox[{"[", "n", "]"}], "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"n", ",", "1", ",", 
           RowBox[{"m", "+", "1"}]}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"s", "=", 
       RowBox[{"Numerator", "[", 
        RowBox[{"Together", "[", "exp1", "]"}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"r", "=", 
       RowBox[{"NewPACARootFinder", "[", 
        RowBox[{"s", ",", "\[Omega]"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"Return", "[", "r", "]"}]}], "\[IndentingNewLine]", ",", "1"}], 
    "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.7449521913452673`*^9, 3.7449521926222906`*^9}, {
   3.744952298478713*^9, 3.744952359510972*^9}, {3.744952540645899*^9, 
   3.744952541780607*^9}, {3.74495632732815*^9, 3.7449563873713894`*^9}, {
   3.7449564204422054`*^9, 3.744956465708021*^9}, 3.7449566018008127`*^9, {
   3.744956680334388*^9, 3.744956680592323*^9}, {3.7449567159826*^9, 
   3.744956819817954*^9}, {3.744956895080092*^9, 3.744956935320843*^9}, {
   3.7449569704361954`*^9, 3.7449570310093985`*^9}, {3.7449598090046844`*^9, 
   3.7449598261777163`*^9}, {3.7449598710509853`*^9, 
   3.7449600459635916`*^9}, {3.7449600989602623`*^9, 
   3.7449601450708337`*^9}, {3.7449602805139313`*^9, 
   3.7449603849949074`*^9}, {3.744960422513364*^9, 3.7449604626667967`*^9}, {
   3.7449605564618635`*^9, 3.744960600608215*^9}, {3.744960683715615*^9, 
   3.744960688261647*^9}, {3.7449607504998703`*^9, 3.7449607507063265`*^9}, {
   3.744962079501486*^9, 3.7449620844717884`*^9}, {3.7449621211666656`*^9, 
   3.7449621579450254`*^9}, {3.7449622922449093`*^9, 3.744962296593932*^9}, {
   3.7449624736016874`*^9, 3.7449624835261865`*^9}, {3.744962757954794*^9, 
   3.7449627620580196`*^9}, {3.744963088981325*^9, 3.744963113622847*^9}, {
   3.7449631524776816`*^9, 3.7449631876269264`*^9}, {3.7449632177780695`*^9, 
   3.7449632215225515`*^9}, {3.7449633489689784`*^9, 3.744963356124629*^9}, 
   3.744963397142354*^9, {3.744963665776214*^9, 3.7449636796022425`*^9}, {
   3.7449637140266976`*^9, 3.74496375455307*^9}, {3.7449638191613135`*^9, 
   3.744963821867092*^9}, {3.744963910704813*^9, 3.7449639500782986`*^9}, {
   3.7449640070632195`*^9, 3.744964059219761*^9}, {3.7449642397388444`*^9, 
   3.744964244505695*^9}, {3.7449642756385436`*^9, 3.7449643359373837`*^9}, 
   3.7449847594718065`*^9, {3.744984797904793*^9, 3.7449848333107758`*^9}, {
   3.7449848782590113`*^9, 3.7449848829576163`*^9}, {3.744984920975251*^9, 
   3.7449849210544653`*^9}, {3.7449849730821877`*^9, 
   3.7449849902419386`*^9}, {3.7449850463692293`*^9, 3.744985052597821*^9}, {
   3.744985156529217*^9, 3.74498524934625*^9}, {3.744985304159028*^9, 
   3.74498537536794*^9}, {3.744985409969453*^9, 3.7449854103453226`*^9}, {
   3.744985613635949*^9, 3.7449856323683805`*^9}, {3.744985672710044*^9, 
   3.744985694925222*^9}, {3.744985824001689*^9, 3.7449858649886155`*^9}, {
   3.7449859320495806`*^9, 3.744985934417283*^9}, {3.7449859663466706`*^9, 
   3.744986013936191*^9}, {3.7449863383777556`*^9, 3.7449863654824476`*^9}, {
   3.7449865679406085`*^9, 3.744986581369054*^9}, {3.744986627324888*^9, 
   3.744986640683814*^9}, {3.7449870411032457`*^9, 3.7449870831348896`*^9}, {
   3.7449871153803473`*^9, 3.7449871556052513`*^9}, {3.744987388504115*^9, 
   3.744987395979446*^9}, {3.74498744042052*^9, 3.744987443495529*^9}, {
   3.744987489710492*^9, 3.7449875518667173`*^9}, {3.7449876113492804`*^9, 
   3.7449876116045866`*^9}, {3.744987770922022*^9, 3.744987788678034*^9}, {
   3.744987938323732*^9, 3.744987952161995*^9}, {3.7449880293128595`*^9, 
   3.744988115995302*^9}, {3.7449882461142464`*^9, 3.7449882887715664`*^9}, {
   3.7449883436297865`*^9, 
   3.744988361886529*^9}},ExpressionUUID->"1ee244e2-427e-48a3-b3dd-\
5109a29a98c9"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"start", "=", 
   RowBox[{"AbsoluteTime", "[", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"PACAGetModes", "[", 
  RowBox[{"3", ",", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"1", "-", 
      RowBox[{"2", " ", "M", " ", "x"}], "+", 
      RowBox[{
       RowBox[{"1", "/", 
        RowBox[{"x", "^", "2"}]}], "/", 
       RowBox[{"\[Lambda]", "^", "2"}]}]}], "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"M", "\[Rule]", "6"}], ",", 
       RowBox[{"\[Lambda]", "\[Rule]", "2"}]}], "}"}]}], ")"}], ",", "x", ",",
    "0", ",", 
   RowBox[{"-", "3"}]}], "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"end", "=", 
   RowBox[{"AbsoluteTime", "[", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{"end", "-", "start"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"Clear", "[", 
  RowBox[{"end", ",", "start"}], "]"}]}], "Input",
 CellChangeTimes->{{3.744956490793466*^9, 3.7449565604486055`*^9}, {
   3.744956632787993*^9, 3.744956653299733*^9}, {3.7449568385678916`*^9, 
   3.7449568556584706`*^9}, {3.744957051533181*^9, 3.7449571035202875`*^9}, {
   3.7449571836059713`*^9, 3.7449572341427083`*^9}, {3.7449573352068343`*^9, 
   3.7449574054247503`*^9}, {3.7449591140757027`*^9, 
   3.7449591159235663`*^9}, {3.7449591730565305`*^9, 
   3.7449592069126625`*^9}, {3.7449592605661507`*^9, 
   3.7449592854872437`*^9}, {3.744959377529653*^9, 3.7449593884704*^9}, {
   3.744959437063609*^9, 3.744959494261145*^9}, {3.7449598443796253`*^9, 
   3.744959853525857*^9}, {3.7449598953249383`*^9, 3.7449598972134824`*^9}, {
   3.7449599357401667`*^9, 3.7449599365518723`*^9}, {3.744959989944744*^9, 
   3.744959991392049*^9}, {3.744960051474907*^9, 3.7449600557440677`*^9}, {
   3.7449603961190977`*^9, 3.7449604148846235`*^9}, {3.7449604727056103`*^9, 
   3.744960535168971*^9}, {3.7449606323509073`*^9, 3.7449606526775684`*^9}, {
   3.744960707777443*^9, 3.744960717566272*^9}, {3.7449607599559455`*^9, 
   3.744960761797044*^9}, {3.744962114471135*^9, 3.744962134201274*^9}, {
   3.74496245416081*^9, 3.7449624549521456`*^9}, {3.7449628948918757`*^9, 
   3.744962953477597*^9}, {3.7449631348462453`*^9, 3.7449631397185116`*^9}, 
   3.744963726499408*^9, {3.744963763847467*^9, 3.7449638104766407`*^9}, {
   3.744963852060524*^9, 3.744963891463736*^9}, {3.7449641368774*^9, 
   3.7449641379567986`*^9}, {3.744964359207734*^9, 3.744964380128233*^9}, {
   3.7449644561107235`*^9, 3.744964462227107*^9}, 3.7449645248522487`*^9, {
   3.744984906375869*^9, 3.744984932424016*^9}, {3.7449855367959557`*^9, 
   3.7449855475525303`*^9}, {3.744985655329217*^9, 3.7449856620590897`*^9}, {
   3.7449857078880143`*^9, 3.744985708335929*^9}, {3.7449860237663555`*^9, 
   3.7449860315075693`*^9}, {3.7449860896081944`*^9, 3.744986148621373*^9}, {
   3.7449875834243383`*^9, 3.744987596712726*^9}, {3.744987744875256*^9, 
   3.7449877578126473`*^9}, {3.7449879781563187`*^9, 3.744987979438125*^9}, {
   3.7449883136999283`*^9, 
   3.744988332364383*^9}},ExpressionUUID->"eb71481b-61af-462e-b1c4-\
10bb29f4fad0"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"3.6458673416981062`", "\[VeryThinSpace]", "+", 
    RowBox[{"5.931269349648283`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{
    RowBox[{"-", "3.6458673416981062`"}], "+", 
    RowBox[{"5.931269349648283`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.2024283136726881`", " ", "\[ImaginaryI]"}]}]}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.7449874516723375`*^9, 3.7449874807991023`*^9}, {
   3.7449875733084054`*^9, 3.744987602296988*^9}, {3.744987739575447*^9, 
   3.744987759370823*^9}, {3.7449879674517*^9, 3.7449879847433786`*^9}, {
   3.744988310287976*^9, 3.744988330367915*^9}, 3.7449883966509943`*^9, 
   3.747981898479988*^9, 3.747983682849817*^9, {3.7479937305004635`*^9, 
   3.747993748354534*^9}, 3.7479943490137067`*^9, 
   3.747994746077136*^9},ExpressionUUID->"5aa78c3f-5602-42cb-8b52-\
e65758bc3d21"],

Cell[BoxData["0.014533`5.613900267048052"], "Print",
 CellChangeTimes->{{3.744987451704632*^9, 3.7449874807991023`*^9}, {
   3.7449875733084054`*^9, 3.744987602296988*^9}, {3.744987739590665*^9, 
   3.744987759370823*^9}, {3.7449879674517*^9, 3.7449879847433786`*^9}, {
   3.744988310287976*^9, 3.7449883303837547`*^9}, 3.7449883966509943`*^9, 
   3.7479818984949884`*^9, 3.747983682867354*^9, {3.7479937305004635`*^9, 
   3.747993748354534*^9}, 3.7479943490137067`*^9, 
   3.747994746077136*^9},ExpressionUUID->"91a72e39-aa1d-4f66-a249-\
509d0017dc04"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"PACAPlotModes", "[", 
   RowBox[{"m_", ",", "exp_", ",", "x_", ",", "W_", ",", "\[Beta]_"}], "]"}], 
  ":=", " ", 
  RowBox[{"ListPlot", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"Tooltip", "[", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Re", "[", "#1", "]"}], ",", 
          RowBox[{"Abs", "[", 
           RowBox[{"Im", "[", "#1", "]"}], "]"}]}], "}"}], "]"}], "&"}], 
      ")"}], "/@", 
     RowBox[{"PACAGetModes", "[", 
      RowBox[{"m", ",", "exp", ",", "x", ",", "W", ",", "\[Beta]"}], "]"}]}], 
    ",", 
    RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Re", "[", "\[Omega]", "]"}], ",", 
       RowBox[{"Im", "[", "\[Omega]", "]"}]}], "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.7449595164416637`*^9, 3.7449595637149577`*^9}, {
  3.7449596860351095`*^9, 3.7449596888635116`*^9}, {3.744964109906399*^9, 
  3.7449641274667683`*^9}},
 NumberMarks->False,ExpressionUUID->"c72a2cb8-2837-4f31-a1b2-0030f61b6a56"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PACAPlotModes", "[", 
  RowBox[{"34", ",", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"1", "-", 
      RowBox[{"2", " ", "M", " ", "x"}], "+", 
      RowBox[{
       RowBox[{"1", "/", 
        RowBox[{"x", "^", "2"}]}], "/", 
       RowBox[{"l", "^", "2"}]}]}], "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"M", "\[Rule]", "6"}], ",", 
       RowBox[{"\[Lambda]", "\[Rule]", "2"}]}], "}"}]}], ")"}], ",", "x", ",",
    "1", ",", "1"}], "]"}]], "Input",
 CellChangeTimes->{{3.7449596203501534`*^9, 3.744959620650761*^9}, {
  3.7449596516727095`*^9, 3.7449597288995686`*^9}, {3.7449602066804605`*^9, 
  3.7449602631554923`*^9}, {3.7449607861794558`*^9, 3.744960787613287*^9}, {
  3.7449627749735093`*^9, 3.7449627765384626`*^9}, {3.7449628562430778`*^9, 
  3.744962882769644*^9}, {3.7449629702943077`*^9, 3.7449630300678215`*^9}, {
  3.7449640851245193`*^9, 3.7449640871517353`*^9}, {3.744964386966626*^9, 
  3.744964387364276*^9}, {3.744964446962843*^9, 
  3.744964498850993*^9}},ExpressionUUID->"97a87942-c0a5-4551-bc5a-\
5bc7c895b364"],

Cell[BoxData[
 GraphicsBox[{{}, {{{}, 
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
      0.012833333333333334`], AbsoluteThickness[1.6], PointBox[CompressedData["
1:eJxTTMoPSmViYGBQAmIQXX7g9jbB0hf2m+yt+VUbz9gnvfpmqBH1wj5PODf4
765N9qpOgqlrk57bS53em5K3Y5P9sZd7li6c+dRe4Y9ucOr/U/YmT0rEN3E+
tl/1c3r8t6ar9jMvBcoc7L1vr8lkcnhnxR37U0kfljCdvm6/dOmtdAnhh/Yn
Gicpv3p01v5V3hvzXK3H9lbNiv5TJ2+x/7+97rVy5BMYfz+Mfxyifv9rqHqo
efth5kHt2w+zD+qe/TD3QN27H+ZeeYh/9rNC/ZMN8e9+b6h/6yHhsf8BNDz+
i0s8yNr0fH+qwJ4LV65ct3/vYmH9Qe7Zfsapgd365Q/s776a+lJvzeP9n3nO
L/16+ol99cZY+1UrHuy/qbLiTmrDC3u+J/UrNNJu7hefJ1P06c1r+3XLrBvu
tFzcz6jEt95V9519uFLtQze2ffsjWatmKaZ+sGdguDYlqPq6pYVNuvz6zA/2
kevqQPL24Z8rwfJLn4H120/hhegXWgO2zz5aHWIfMEoZtl0RsnZbdloKyHRI
3gy2375kA8T+NIh77a24IO7d4Az2j/1jqH8+Q/xrLwj1LwCMwBqn
       "]]}, {}}, {
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.7642906733945115, 0.22269118541453065`}}]},
       RowBox[{"{", 
         RowBox[{"0.7642906733945115`", ",", "0.22269118541453065`"}], "}"}]],
      Annotation[#, {0.7642906733945115, 0.22269118541453065`}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.7610054946774876, 0.07316573416490055}}]},
       RowBox[{"{", 
         RowBox[{"0.7610054946774876`", ",", "0.07316573416490055`"}], "}"}]],
      Annotation[#, {0.7610054946774876, 0.07316573416490055}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.7307955717061004, 0.07312669715410572}}]},
       RowBox[{"{", 
         RowBox[{"0.7307955717061004`", ",", "0.07312669715410572`"}], "}"}]],
      Annotation[#, {0.7307955717061004, 0.07312669715410572}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.6750038373059446, 0.21091906129365245`}}]},
       RowBox[{"{", 
         RowBox[{"0.6750038373059446`", ",", "0.21091906129365245`"}], "}"}]],
      Annotation[#, {0.6750038373059446, 0.21091906129365245`}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.5949335535179956, 0.3361183103114863}}]},
       RowBox[{"{", 
         RowBox[{"0.5949335535179956`", ",", "0.3361183103114863`"}], "}"}]],
      Annotation[#, {0.5949335535179956, 0.3361183103114863}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.49302699819171775`, 0.44486850798424143`}}]},
       RowBox[{"{", 
         RowBox[{"0.49302699819171775`", ",", "0.44486850798424143`"}], "}"}]],
      Annotation[#, {0.49302699819171775`, 0.44486850798424143`}, 
       "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.37176576716732923`, 0.5335809734726139}}]},
       RowBox[{"{", 
         RowBox[{"0.37176576716732923`", ",", "0.5335809734726139`"}], "}"}]],
      Annotation[#, {0.37176576716732923`, 0.5335809734726139}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.23348738419464454`, 0.5989290325410057}}]},
       RowBox[{"{", 
         RowBox[{"0.23348738419464454`", ",", "0.5989290325410057`"}], "}"}]],
      Annotation[#, {0.23348738419464454`, 0.5989290325410057}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.08037694150835942, 0.6358813857566473}}]},
       RowBox[{"{", 
         RowBox[{"0.08037694150835942`", ",", "0.6358813857566473`"}], "}"}]],
      Annotation[#, {0.08037694150835942, 0.6358813857566473}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{-0.08037694150835942, 0.6358813857566473}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.08037694150835942`"}], ",", 
           "0.6358813857566473`"}], "}"}]],
      Annotation[#, {-0.08037694150835942, 0.6358813857566473}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{-0.23348738419464451`, 0.5989290325410058}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.23348738419464451`"}], ",", 
           "0.5989290325410058`"}], "}"}]],
      Annotation[#, {-0.23348738419464451`, 0.5989290325410058}, 
       "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{-0.37176576716732923`, 0.5335809734726139}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.37176576716732923`"}], ",", 
           "0.5335809734726139`"}], "}"}]],
      Annotation[#, {-0.37176576716732923`, 0.5335809734726139}, 
       "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], 
        PointBox[{{-0.49302699819171775`, 0.44486850798424143`}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.49302699819171775`"}], ",", 
           "0.44486850798424143`"}], "}"}]],
      Annotation[#, {-0.49302699819171775`, 0.44486850798424143`}, 
       "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{-0.5949335535179956, 0.3361183103114863}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.5949335535179956`"}], ",", "0.3361183103114863`"}],
          "}"}]],
      Annotation[#, {-0.5949335535179956, 0.3361183103114863}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{-0.6750038373059446, 0.21091906129365245`}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.6750038373059446`"}], ",", 
           "0.21091906129365245`"}], "}"}]],
      Annotation[#, {-0.6750038373059446, 0.21091906129365245`}, 
       "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{-0.7307955717060998, 0.07312669715410543}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.7307955717060998`"}], ",", 
           "0.07312669715410543`"}], "}"}]],
      Annotation[#, {-0.7307955717060998, 0.07312669715410543}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{-0.7610054946774886, 0.07316573416490006}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.7610054946774886`"}], ",", 
           "0.07316573416490006`"}], "}"}]],
      Annotation[#, {-0.7610054946774886, 0.07316573416490006}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{-0.7642906733945124, 0.22269118541453192`}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.7642906733945124`"}], ",", 
           "0.22269118541453192`"}], "}"}]],
      Annotation[#, {-0.7642906733945124, 0.22269118541453192`}, 
       "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{-0.740529477771588, 0.37236519226303805`}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.740529477771588`"}], ",", "0.37236519226303805`"}],
          "}"}]],
      Annotation[#, {-0.740529477771588, 0.37236519226303805`}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{-0.6912766605984454, 0.5145490380351988}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.6912766605984454`"}], ",", "0.5145490380351988`"}],
          "}"}]],
      Annotation[#, {-0.6912766605984454, 0.5145490380351988}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{-0.61476846334639, 0.6498974073842007}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.61476846334639`"}], ",", "0.6498974073842007`"}], 
         "}"}]],
      Annotation[#, {-0.61476846334639, 0.6498974073842007}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{-0.5205889928452473, 0.7656735715583195}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.5205889928452473`"}], ",", "0.7656735715583195`"}],
          "}"}]],
      Annotation[#, {-0.5205889928452473, 0.7656735715583195}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{-0.39686027960812076`, 0.8726742008920522}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.39686027960812076`"}], ",", 
           "0.8726742008920522`"}], "}"}]],
      Annotation[#, {-0.39686027960812076`, 0.8726742008920522}, 
       "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{-0.2737342121387013, 0.943026391892147}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.2737342121387013`"}], ",", "0.943026391892147`"}], 
         "}"}]],
      Annotation[#, {-0.2737342121387013, 0.943026391892147}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{-0.11728327756021294`, 1.024690249847785}}]},
       RowBox[{"{", 
         RowBox[{
           RowBox[{"-", "0.11728327756021294`"}], ",", "1.024690249847785`"}],
          "}"}]],
      Annotation[#, {-0.11728327756021294`, 1.024690249847785}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], 
        PointBox[{{4.630954467616832*^-30, 1.0258017756166407`}}]},
       RowBox[{"{", 
         RowBox[{"4.630954467616832`*^-30", ",", "1.0258017756166407`"}], 
         "}"}]],
      Annotation[#, {4.630954467616832*^-30, 1.0258017756166407`}, 
       "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.11728327756161984`, 1.0246902498467614`}}]},
       RowBox[{"{", 
         RowBox[{"0.11728327756161984`", ",", "1.0246902498467614`"}], "}"}]],
      Annotation[#, {0.11728327756161984`, 1.0246902498467614`}, 
       "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.2737342121396103, 0.9430263918915665}}]},
       RowBox[{"{", 
         RowBox[{"0.2737342121396103`", ",", "0.9430263918915665`"}], "}"}]],
      Annotation[#, {0.2737342121396103, 0.9430263918915665}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.5205889928450935, 0.7656735715583908}}]},
       RowBox[{"{", 
         RowBox[{"0.5205889928450935`", ",", "0.7656735715583908`"}], "}"}]],
      Annotation[#, {0.5205889928450935, 0.7656735715583908}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], 
        PointBox[{{3.8941392834202444`*^-24, 2.000000199642076}}]},
       RowBox[{"{", 
         RowBox[{"3.8941392834202444`*^-24", ",", "2.000000199642076`"}], 
         "}"}]],
      Annotation[#, {3.8941392834202444`*^-24, 2.000000199642076}, 
       "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.39686027960742914`, 0.8726742008925741}}]},
       RowBox[{"{", 
         RowBox[{"0.39686027960742914`", ",", "0.8726742008925741`"}], "}"}]],
      Annotation[#, {0.39686027960742914`, 0.8726742008925741}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.6147684633463768, 0.6498974073841233}}]},
       RowBox[{"{", 
         RowBox[{"0.6147684633463768`", ",", "0.6498974073841233`"}], "}"}]],
      Annotation[#, {0.6147684633463768, 0.6498974073841233}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.69127666059841, 0.5145490380352239}}]},
       RowBox[{"{", 
         RowBox[{"0.69127666059841`", ",", "0.5145490380352239`"}], "}"}]],
      Annotation[#, {0.69127666059841, 0.5145490380352239}, "Tooltip"]& ], 
     TagBox[
      TooltipBox[
       {Opacity[0.], PointBox[{{0.7405294777715866, 0.3723651922630334}}]},
       RowBox[{"{", 
         RowBox[{"0.7405294777715866`", ",", "0.3723651922630334`"}], "}"}]],
      Annotation[#, {0.7405294777715866, 0.3723651922630334}, 
       "Tooltip"]& ]}}, {}, {}, {{}, {}}},
  AspectRatio->1,
  Axes->{True, True},
  AxesLabel->{None, None},
  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]],
  ImagePadding->All,
  Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
         Part[#, 1]], 
        (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
         Part[#, 1]], 
        (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{-0.7642906733945124, 0.7642906733945115}, {
    0, 2.000000199642076}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.744959667565563*^9, 3.7449597062316027`*^9}, 
   3.7449598019694495`*^9, {3.744960209027565*^9, 3.7449602642436304`*^9}, 
   3.744960807226019*^9, 3.7449627777514343`*^9, {3.7449628582078238`*^9, 
   3.744962884572164*^9}, {3.7449629717142367`*^9, 3.7449630316580625`*^9}, {
   3.7449640691878014`*^9, 3.7449640894627347`*^9}, 3.7449641304862323`*^9, 
   3.7449643544111176`*^9, 3.7449643907374763`*^9, {3.7449644523224325`*^9, 
   3.744964500839737*^9}, 3.744985514856537*^9, 3.744986042616168*^9, 
   3.747981902765833*^9, 3.7479836892199683`*^9, 3.7479937789104595`*^9, 
   3.747994749940632*^9},ExpressionUUID->"f7c71fb8-c2cc-48be-abc2-\
2953179e0ef4"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"PACATableModes", "[", 
    RowBox[{"m_", ",", "exp_", ",", "x_", ",", "W_", ",", "\[Beta]_"}], "]"}],
    ":=", "\[IndentingNewLine]", 
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"modes", "=", 
       RowBox[{"PACAGetModes", "[", 
        RowBox[{"m", ",", "exp", ",", "x", ",", "W", ",", "\[Beta]"}], 
        "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"rearrange", "=", 
       RowBox[{"Table", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Re", "[", 
            RowBox[{"modes", "[", 
             RowBox[{"[", "k", "]"}], "]"}], "]"}], ",", 
           RowBox[{"Im", "[", 
            RowBox[{"modes", "[", 
             RowBox[{"[", "k", "]"}], "]"}], "]"}]}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"k", ",", "1", ",", "m"}], "}"}]}], "]"}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"Return", "[", 
       RowBox[{"TableForm", "[", 
        RowBox[{"rearrange", ",", 
         RowBox[{"TableHeadings", "\[Rule]", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"Range", "[", "m", "]"}], ",", 
            RowBox[{"{", 
             RowBox[{"\"\<\[Omega]R\>\"", ",", "\"\<\[Omega]I\>\""}], "}"}]}],
            "}"}]}]}], "]"}], "]"}]}], "\[IndentingNewLine]", ",", "1"}], 
    "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.7449648069224763`*^9, 3.744964877351656*^9}, {
  3.7449649861880445`*^9, 3.7449650170569115`*^9}, {3.7449650762186584`*^9, 
  3.744965142009015*^9}, {3.74496517928568*^9, 3.7449652015147104`*^9}, {
  3.744965237043289*^9, 3.74496529481324*^9}, {3.7449654186044755`*^9, 
  3.744965425366977*^9}, {3.7449654715440364`*^9, 
  3.744965525221916*^9}},ExpressionUUID->"1ae650d2-7bf3-474d-befe-\
66e26555971c"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PACATableModes", "[", 
  RowBox[{"25", ",", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"1", "-", 
      RowBox[{"2", " ", "M", " ", "x"}], "+", 
      RowBox[{
       RowBox[{"1", "/", 
        RowBox[{"x", "^", "2"}]}], "/", 
       RowBox[{"\[Lambda]", "^", "2"}]}]}], "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"M", "\[Rule]", "6"}], ",", 
       RowBox[{"\[Lambda]", "\[Rule]", "2"}]}], "}"}]}], ")"}], ",", "x", ",",
    "1", ",", "1"}], "]"}]], "Input",
 CellChangeTimes->{{3.7449652977303715`*^9, 3.7449653322135963`*^9}, {
  3.744965495343505*^9, 
  3.7449655060354857`*^9}},ExpressionUUID->"ef264675-4148-453b-8957-\
a1c8578d0a32"],

Cell[BoxData[
 TagBox[
  TagBox[GridBox[{
     {
      StyleBox["\[Null]",
       ShowStringCharacters->False], 
      TagBox["\<\"\[Omega]R\"\>",
       HoldForm], 
      TagBox["\<\"\[Omega]I\"\>",
       HoldForm]},
     {
      TagBox["1",
       HoldForm], "0.6902298792219161`", 
      RowBox[{"-", "0.1344038733787388`"}]},
     {
      TagBox["2",
       HoldForm], "0.6505809196035197`", "0.04949324192831641`"},
     {
      TagBox["3",
       HoldForm], "0.5665562575457388`", "0.2164382583181933`"},
     {
      TagBox["4",
       HoldForm], 
      RowBox[{"-", "0.09981042579103763`"}], "0.5241363253028833`"},
     {
      TagBox["5",
       HoldForm], "0.44340011600509227`", "0.35770302784140856`"},
     {
      TagBox["6",
       HoldForm], "0.2857860869862677`", "0.4636499178063146`"},
     {
      TagBox["7",
       HoldForm], "1.2106462030595684`*^-38", 
      RowBox[{"-", "2.0000003003034146`"}]},
     {
      TagBox["8",
       HoldForm], 
      RowBox[{"-", "0.2857860869862677`"}], "0.4636499178063146`"},
     {
      TagBox["9",
       HoldForm], 
      RowBox[{"-", "0.44340011600509227`"}], "0.35770302784140856`"},
     {
      TagBox["10",
       HoldForm], "0.09981042579103763`", "0.5241363253028833`"},
     {
      TagBox["11",
       HoldForm], 
      RowBox[{"-", "0.5665562575457389`"}], "0.21643825831819338`"},
     {
      TagBox["12",
       HoldForm], 
      RowBox[{"-", "0.6505809196035199`"}], "0.04949324192831641`"},
     {
      TagBox["13",
       HoldForm], 
      RowBox[{"-", "0.690229879221916`"}], 
      RowBox[{"-", "0.13440387337873902`"}]},
     {
      TagBox["14",
       HoldForm], 
      RowBox[{"-", "0.6853971104820735`"}], 
      RowBox[{"-", "0.32314178756699624`"}]},
     {
      TagBox["15",
       HoldForm], 
      RowBox[{"-", "0.6331338523231823`"}], 
      RowBox[{"-", "0.509630184765819`"}]},
     {
      TagBox["16",
       HoldForm], 
      RowBox[{"-", "0.5389751646494793`"}], 
      RowBox[{"-", "0.6763780456515694`"}]},
     {
      TagBox["17",
       HoldForm], 
      RowBox[{"-", "0.4072312639877149`"}], 
      RowBox[{"-", "0.8254952500867435`"}]},
     {
      TagBox["18",
       HoldForm], 
      RowBox[{"-", "0.2395543240183497`"}], 
      RowBox[{"-", "0.934538011010911`"}]},
     {
      TagBox["19",
       HoldForm], 
      RowBox[{"-", "0.06996399198432866`"}], 
      RowBox[{"-", "1.0237916460017036`"}]},
     {
      TagBox["20",
       HoldForm], "0.06996399198581472`", 
      RowBox[{"-", "1.023791646001167`"}]},
     {
      TagBox["21",
       HoldForm], "0.2395543240189389`", 
      RowBox[{"-", "0.9345380110107225`"}]},
     {
      TagBox["22",
       HoldForm], "0.4072312639879013`", 
      RowBox[{"-", "0.8254952500867611`"}]},
     {
      TagBox["23",
       HoldForm], "0.5389751646495403`", 
      RowBox[{"-", "0.6763780456515457`"}]},
     {
      TagBox["24",
       HoldForm], "0.6331338523231888`", 
      RowBox[{"-", "0.5096301847658269`"}]},
     {
      TagBox["25",
       HoldForm], "0.6853971104820756`", 
      RowBox[{"-", "0.32314178756699347`"}]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxDividers->{
     "Columns" -> {False, True, {False}, False}, "ColumnsIndexed" -> {}, 
      "Rows" -> {False, True, {False}, False}, "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[2.0999999999999996`]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}],
   {OutputFormsDump`HeadedRows, OutputFormsDump`HeadedColumns}],
  Function[BoxForm`e$, 
   TableForm[
   BoxForm`e$, 
    TableHeadings -> {{1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 
      17, 18, 19, 20, 21, 22, 23, 24, 25}, {
      "\[Omega]R", "\[Omega]I"}}]]]], "Output",
 CellChangeTimes->{{3.744965318941342*^9, 3.744965340593396*^9}, 
   3.7449654376314945`*^9, {3.744965485842355*^9, 3.7449655274825697`*^9}, 
   3.744985520339907*^9, 3.744986050206764*^9, 3.747981906515972*^9, 
   3.7479836952990727`*^9, 3.7479937842826023`*^9, 
   3.7479947562063007`*^9},ExpressionUUID->"0d20787a-cbad-483c-ad21-\
53876bc92d8d"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"PACAShowModes", "[", 
   RowBox[{"m_", ",", "exp_", ",", "x_", ",", "W_", ",", "\[Beta]_"}], "]"}], 
  ":=", 
  RowBox[{"Print", "[", 
   RowBox[{
    RowBox[{"Show", "[", 
     RowBox[{
      RowBox[{"PACAPlotModes", "[", 
       RowBox[{"m", ",", "exp", ",", "x", ",", "W", ",", "\[Beta]"}], "]"}], 
      ",", 
      RowBox[{"ImageSize", "\[Rule]", "Large"}]}], "]"}], ",", "\"\<\\t\>\"", 
    ",", 
    RowBox[{"PACATableModes", "[", 
     RowBox[{"m", ",", "exp", ",", "x", ",", "W", ",", "\[Beta]"}], "]"}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.744965837653386*^9, 3.74496584574325*^9}, {
  3.7449659670587206`*^9, 3.744965999663277*^9}, {3.744966103994971*^9, 
  3.7449661107246923`*^9}, {3.7449661580403204`*^9, 3.744966172665733*^9}, {
  3.744966281730274*^9, 
  3.7449663598802547`*^9}},ExpressionUUID->"8d7c1a55-460c-4b70-9151-\
1c2c9bd2d128"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PACAShowModes", "[", 
  RowBox[{"35", ",", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"1", "-", 
      RowBox[{"2", " ", "M", " ", "x"}], "+", 
      RowBox[{
       RowBox[{"1", "/", 
        RowBox[{"x", "^", "2"}]}], "/", 
       RowBox[{"\[Lambda]", "^", "2"}]}]}], "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"M", "\[Rule]", "6"}], ",", 
       RowBox[{"\[Lambda]", "\[Rule]", "2"}]}], "}"}]}], ")"}], ",", "x", ",",
    "1", ",", "1"}], "]"}]], "Input",
 CellChangeTimes->{{3.7449660014521074`*^9, 3.7449660157073283`*^9}, {
  3.7449661139444036`*^9, 3.7449661147320004`*^9}, {3.7449661821146865`*^9, 
  3.7449662025665674`*^9}, {3.744966320876524*^9, 3.7449663216554885`*^9}, {
  3.7449663672924128`*^9, 
  3.7449663680978804`*^9}},ExpressionUUID->"065c127c-0626-4ef4-bd10-\
1723e1c6a173"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
   GraphicsBox[{{}, {{{}, 
       {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
        0.012833333333333334`], AbsoluteThickness[1.6], 
        PointBox[CompressedData["
1:eJxTTMoPSmViYGBQBmIQzfzI9caCLS/s530V3LvK+ri9Cdfh5BWpL+xnB916
yei1yL6oIV1cxOO5vYWKrOcm6d329RwcD/akPLW3Ve2oVO0+Z2/SLLN41vFH
9mq2HpZm367Z/9V81/Gb7b69nI2+4rl5d+199sU33wi7bl9w58xc39aH9vnr
xHQYXc7aH5b7ffJ182P7sNTLZz8pbLFffWDxvBlLnsD4+2F8qPr9MPVQ8/bD
zIPatx9mH9Q9+2Hugbp3vz3UvQUQ/+yH+cce4t/9isEQ/356CA6P/TbQ8MiK
EXGp0Xux//rihTYO/6/Yi01um1b46Nl+/SvTjzEK37OPLztwJGzfk/27ZuUb
Nwc/thfMO/llP/Oj/T8+cFaLWT2331h4XvP+vTv7iwunp6RNemXfr1537XfA
1f0pUbanlxS+tfessf3ot//4/vD51rqqme/tg7awCne/2bDfd/X9lGNzPtgD
o6iB+6zLrlnT3goD2Q6Wm8Hy9k7TIfKFhmD99pJZNmD9Z13A5tv3mEDM504H
22+v0Aax/14p2H32hp8g7vteCHa/vfd8iPtlJoH9Z58C9d9GiP/tj0D9DwCe
0Btx
         "]]}, {}}, {
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.7720493530351294, 0.18150839073124997`}}]},
         RowBox[{"{", 
           RowBox[{"0.7720493530351294`", ",", "0.18150839073124997`"}], 
           "}"}]],
        Annotation[#, {0.7720493530351294, 0.18150839073124997`}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.7624093960260496, 0.03572088222168653}}]},
         RowBox[{"{", 
           RowBox[{"0.7624093960260496`", ",", "0.03572088222168653`"}], 
           "}"}]],
        Annotation[#, {0.7624093960260496, 0.03572088222168653}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.7275486428370768, 0.10589136395459853`}}]},
         RowBox[{"{", 
           RowBox[{"0.7275486428370768`", ",", "0.10589136395459853`"}], 
           "}"}]],
        Annotation[#, {0.7275486428370768, 0.10589136395459853`}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.6685470939267047, 0.2386214106957513}}]},
         RowBox[{"{", 
           RowBox[{"0.6685470939267047`", ",", "0.2386214106957513`"}], "}"}]],
        Annotation[#, {0.6685470939267047, 0.2386214106957513}, "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.5868657289284784, 0.3587775763889617}}]},
         RowBox[{"{", 
           RowBox[{"0.5868657289284784`", ",", "0.3587775763889617`"}], "}"}]],
        Annotation[#, {0.5868657289284784, 0.3587775763889617}, "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.4848011815271603, 0.4628176998202792}}]},
         RowBox[{"{", 
           RowBox[{"0.4848011815271603`", ",", "0.4628176998202792`"}], "}"}]],
        Annotation[#, {0.4848011815271603, 0.4628176998202792}, "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.3646756442021555, 0.5475223619495306}}]},
         RowBox[{"{", 
           RowBox[{"0.3646756442021555`", ",", "0.5475223619495306`"}], "}"}]],
        Annotation[#, {0.3646756442021555, 0.5475223619495306}, "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.22863783505214338`, 0.6098536438884924}}]},
         RowBox[{"{", 
           RowBox[{"0.22863783505214338`", ",", "0.6098536438884924`"}], 
           "}"}]],
        Annotation[#, {0.22863783505214338`, 0.6098536438884924}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.078627753496581, 0.6450923059874375}}]},
         RowBox[{"{", 
           RowBox[{"0.078627753496581`", ",", "0.6450923059874375`"}], "}"}]],
        Annotation[#, {0.078627753496581, 0.6450923059874375}, "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{-0.078627753496581, 0.6450923059874375}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.078627753496581`"}], ",", 
             "0.6450923059874375`"}], "}"}]],
        Annotation[#, {-0.078627753496581, 0.6450923059874375}, "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], 
          PointBox[{{-0.22863783505214338`, 0.6098536438884924}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.22863783505214338`"}], ",", 
             "0.6098536438884924`"}], "}"}]],
        Annotation[#, {-0.22863783505214338`, 0.6098536438884924}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{-0.3646756442021555, 0.5475223619495306}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.3646756442021555`"}], ",", 
             "0.5475223619495306`"}], "}"}]],
        Annotation[#, {-0.3646756442021555, 0.5475223619495306}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{-0.4848011815271603, 0.4628176998202792}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.4848011815271603`"}], ",", 
             "0.4628176998202792`"}], "}"}]],
        Annotation[#, {-0.4848011815271603, 0.4628176998202792}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{-0.5868657289284784, 0.3587775763889617}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.5868657289284784`"}], ",", 
             "0.3587775763889617`"}], "}"}]],
        Annotation[#, {-0.5868657289284784, 0.3587775763889617}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], 
          PointBox[{{-0.6685470939267047, 0.23862141069575135`}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.6685470939267047`"}], ",", 
             "0.23862141069575135`"}], "}"}]],
        Annotation[#, {-0.6685470939267047, 0.23862141069575135`}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], 
          PointBox[{{-0.7275486428370765, 0.10589136395459853`}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.7275486428370765`"}], ",", 
             "0.10589136395459853`"}], "}"}]],
        Annotation[#, {-0.7275486428370765, 0.10589136395459853`}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{-0.7624093960260508, 0.03572088222168746}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.7624093960260508`"}], ",", 
             "0.03572088222168746`"}], "}"}]],
        Annotation[#, {-0.7624093960260508, 0.03572088222168746}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], 
          PointBox[{{-0.7720493530351276, 0.18150839073124725`}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.7720493530351276`"}], ",", 
             "0.18150839073124725`"}], "}"}]],
        Annotation[#, {-0.7720493530351276, 0.18150839073124725`}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{-0.755674488990633, 0.32807928010020987`}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.755674488990633`"}], ",", 
             "0.32807928010020987`"}], "}"}]],
        Annotation[#, {-0.755674488990633, 0.32807928010020987`}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{-0.7151420535998081, 0.469910091341549}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.7151420535998081`"}], ",", 
             "0.469910091341549`"}], "}"}]],
        Annotation[#, {-0.7151420535998081, 0.469910091341549}, "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{-0.6482347338494029, 0.6039443974121561}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.6482347338494029`"}], ",", 
             "0.6039443974121561`"}], "}"}]],
        Annotation[#, {-0.6482347338494029, 0.6039443974121561}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{-0.562957742785729, 0.7258407977201662}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.562957742785729`"}], ",", 
             "0.7258407977201662`"}], "}"}]],
        Annotation[#, {-0.562957742785729, 0.7258407977201662}, "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{-0.4511030109695425, 0.8303710903815457}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.4511030109695425`"}], ",", 
             "0.8303710903815457`"}], "}"}]],
        Annotation[#, {-0.4511030109695425, 0.8303710903815457}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], 
          PointBox[{{-0.33306785532533295`, 0.9201225251698131}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.33306785532533295`"}], ",", 
             "0.9201225251698131`"}], "}"}]],
        Annotation[#, {-0.33306785532533295`, 0.9201225251698131}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], 
          PointBox[{{-0.18552576808267165`, 0.9815851100401954}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.18552576808267165`"}], ",", 
             "0.9815851100401954`"}], "}"}]],
        Annotation[#, {-0.18552576808267165`, 0.9815851100401954}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], 
          PointBox[{{-0.06610936369334411, 1.0382751407984927`}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "0.06610936369334411`"}], ",", 
             "1.0382751407984927`"}], "}"}]],
        Annotation[#, {-0.06610936369334411, 1.0382751407984927`}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], 
          PointBox[{{-5.250899696766081*^-28, 2.0000001484757606`}}]},
         RowBox[{"{", 
           RowBox[{
             RowBox[{"-", "5.250899696766081`*^-28"}], ",", 
             "2.0000001484757606`"}], "}"}]],
        Annotation[#, {-5.250899696766081*^-28, 2.0000001484757606`}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.06610936369334021, 1.0382751407973534`}}]},
         RowBox[{"{", 
           RowBox[{"0.06610936369334021`", ",", "1.0382751407973534`"}], 
           "}"}]],
        Annotation[#, {0.06610936369334021, 1.0382751407973534`}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.18552576808213986`, 0.9815851100459582}}]},
         RowBox[{"{", 
           RowBox[{"0.18552576808213986`", ",", "0.9815851100459582`"}], 
           "}"}]],
        Annotation[#, {0.18552576808213986`, 0.9815851100459582}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.3330678553257485, 0.9201225251687375}}]},
         RowBox[{"{", 
           RowBox[{"0.3330678553257485`", ",", "0.9201225251687375`"}], "}"}]],
        Annotation[#, {0.3330678553257485, 0.9201225251687375}, "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.45110301096939115`, 0.8303710903821333}}]},
         RowBox[{"{", 
           RowBox[{"0.45110301096939115`", ",", "0.8303710903821333`"}], 
           "}"}]],
        Annotation[#, {0.45110301096939115`, 0.8303710903821333}, 
         "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.5629577427859507, 0.7258407977202009}}]},
         RowBox[{"{", 
           RowBox[{"0.5629577427859507`", ",", "0.7258407977202009`"}], "}"}]],
        Annotation[#, {0.5629577427859507, 0.7258407977202009}, "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.6482347338492777, 0.6039443974122859}}]},
         RowBox[{"{", 
           RowBox[{"0.6482347338492777`", ",", "0.6039443974122859`"}], "}"}]],
        Annotation[#, {0.6482347338492777, 0.6039443974122859}, "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.7151420535997803, 0.4699100913415519}}]},
         RowBox[{"{", 
           RowBox[{"0.7151420535997803`", ",", "0.4699100913415519`"}], "}"}]],
        Annotation[#, {0.7151420535997803, 0.4699100913415519}, "Tooltip"]& ], 
       TagBox[
        TooltipBox[
         {Opacity[0.], PointBox[{{0.7556744889906409, 0.3280792801002088}}]},
         RowBox[{"{", 
           RowBox[{"0.7556744889906409`", ",", "0.3280792801002088`"}], "}"}]],
        Annotation[#, {0.7556744889906409, 0.3280792801002088}, 
         "Tooltip"]& ]}}, {}, {}, {{}, {}}},
    AspectRatio->1,
    Axes->{True, True},
    AxesLabel->{None, None},
    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]],
    ImagePadding->All,
    ImageSize->Large,
    Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
          (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
           Part[#, 1]], 
          (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
           Part[#, 2]]}& ), "CopiedValueFunction" -> ({
          (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
           Part[#, 1]], 
          (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
           Part[#, 2]]}& )}},
    PlotRange->{{-0.7720493530351276, 0.7720493530351294}, {
      0, 2.0000001484757606`}},
    PlotRangeClipping->True,
    PlotRangePadding->{{
       Scaled[0.02], 
       Scaled[0.02]}, {
       Scaled[0.02], 
       Scaled[0.05]}},
    Ticks->{Automatic, Automatic}], "\[InvisibleSpace]", "\<\"\\t\"\>", 
   "\[InvisibleSpace]", 
   TagBox[
    TagBox[GridBox[{
       {
        StyleBox["\[Null]",
         ShowStringCharacters->False], 
        TagBox["\<\"\[Omega]R\"\>",
         HoldForm], 
        TagBox["\<\"\[Omega]I\"\>",
         HoldForm]},
       {
        TagBox["1",
         HoldForm], "0.7720493530351294`", 
        RowBox[{"-", "0.18150839073124997`"}]},
       {
        TagBox["2",
         HoldForm], "0.7624093960260496`", 
        RowBox[{"-", "0.03572088222168653`"}]},
       {
        TagBox["3",
         HoldForm], "0.7275486428370768`", "0.10589136395459853`"},
       {
        TagBox["4",
         HoldForm], "0.6685470939267047`", "0.2386214106957513`"},
       {
        TagBox["5",
         HoldForm], "0.5868657289284784`", "0.3587775763889617`"},
       {
        TagBox["6",
         HoldForm], "0.4848011815271603`", "0.4628176998202792`"},
       {
        TagBox["7",
         HoldForm], "0.3646756442021555`", "0.5475223619495306`"},
       {
        TagBox["8",
         HoldForm], "0.22863783505214338`", "0.6098536438884924`"},
       {
        TagBox["9",
         HoldForm], "0.078627753496581`", "0.6450923059874375`"},
       {
        TagBox["10",
         HoldForm], 
        RowBox[{"-", "0.078627753496581`"}], "0.6450923059874375`"},
       {
        TagBox["11",
         HoldForm], 
        RowBox[{"-", "0.22863783505214338`"}], "0.6098536438884924`"},
       {
        TagBox["12",
         HoldForm], 
        RowBox[{"-", "0.3646756442021555`"}], "0.5475223619495306`"},
       {
        TagBox["13",
         HoldForm], 
        RowBox[{"-", "0.4848011815271603`"}], "0.4628176998202792`"},
       {
        TagBox["14",
         HoldForm], 
        RowBox[{"-", "0.5868657289284784`"}], "0.3587775763889617`"},
       {
        TagBox["15",
         HoldForm], 
        RowBox[{"-", "0.6685470939267047`"}], "0.23862141069575135`"},
       {
        TagBox["16",
         HoldForm], 
        RowBox[{"-", "0.7275486428370765`"}], "0.10589136395459853`"},
       {
        TagBox["17",
         HoldForm], 
        RowBox[{"-", "0.7624093960260508`"}], 
        RowBox[{"-", "0.03572088222168746`"}]},
       {
        TagBox["18",
         HoldForm], 
        RowBox[{"-", "0.7720493530351276`"}], 
        RowBox[{"-", "0.18150839073124725`"}]},
       {
        TagBox["19",
         HoldForm], 
        RowBox[{"-", "0.755674488990633`"}], 
        RowBox[{"-", "0.32807928010020987`"}]},
       {
        TagBox["20",
         HoldForm], 
        RowBox[{"-", "0.7151420535998081`"}], 
        RowBox[{"-", "0.469910091341549`"}]},
       {
        TagBox["21",
         HoldForm], 
        RowBox[{"-", "0.6482347338494029`"}], 
        RowBox[{"-", "0.6039443974121561`"}]},
       {
        TagBox["22",
         HoldForm], 
        RowBox[{"-", "0.562957742785729`"}], 
        RowBox[{"-", "0.7258407977201662`"}]},
       {
        TagBox["23",
         HoldForm], 
        RowBox[{"-", "0.4511030109695425`"}], 
        RowBox[{"-", "0.8303710903815457`"}]},
       {
        TagBox["24",
         HoldForm], 
        RowBox[{"-", "0.33306785532533295`"}], 
        RowBox[{"-", "0.9201225251698131`"}]},
       {
        TagBox["25",
         HoldForm], 
        RowBox[{"-", "0.18552576808267165`"}], 
        RowBox[{"-", "0.9815851100401954`"}]},
       {
        TagBox["26",
         HoldForm], 
        RowBox[{"-", "0.06610936369334411`"}], 
        RowBox[{"-", "1.0382751407984927`"}]},
       {
        TagBox["27",
         HoldForm], 
        RowBox[{"-", "5.250899696766081`*^-28"}], 
        RowBox[{"-", "2.0000001484757606`"}]},
       {
        TagBox["28",
         HoldForm], "0.06610936369334021`", 
        RowBox[{"-", "1.0382751407973534`"}]},
       {
        TagBox["29",
         HoldForm], "0.18552576808213986`", 
        RowBox[{"-", "0.9815851100459582`"}]},
       {
        TagBox["30",
         HoldForm], "0.3330678553257485`", 
        RowBox[{"-", "0.9201225251687375`"}]},
       {
        TagBox["31",
         HoldForm], "0.45110301096939115`", 
        RowBox[{"-", "0.8303710903821333`"}]},
       {
        TagBox["32",
         HoldForm], "0.5629577427859507`", 
        RowBox[{"-", "0.7258407977202009`"}]},
       {
        TagBox["33",
         HoldForm], "0.6482347338492777`", 
        RowBox[{"-", "0.6039443974122859`"}]},
       {
        TagBox["34",
         HoldForm], "0.7151420535997803`", 
        RowBox[{"-", "0.4699100913415519`"}]},
       {
        TagBox["35",
         HoldForm], "0.7556744889906409`", 
        RowBox[{"-", "0.3280792801002088`"}]}
      },
      GridBoxAlignment->{
       "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
        "RowsIndexed" -> {}},
      GridBoxDividers->{
       "Columns" -> {False, True, {False}, False}, "ColumnsIndexed" -> {}, 
        "Rows" -> {False, True, {False}, False}, "RowsIndexed" -> {}},
      GridBoxSpacings->{"Columns" -> {
          Offset[0.27999999999999997`], {
           Offset[2.0999999999999996`]}, 
          Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
          Offset[0.2], {
           Offset[0.4]}, 
          Offset[0.2]}, "RowsIndexed" -> {}}],
     {OutputFormsDump`HeadedRows, OutputFormsDump`HeadedColumns}],
    Function[BoxForm`e$, 
     TableForm[
     BoxForm`e$, 
      TableHeadings -> {{1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 
        16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 
        33, 34, 35}, {"\[Omega]R", "\[Omega]I"}}]]]}],
  SequenceForm[
   Graphics[{{}, {{{}, {
        Hue[0.67, 0.6, 0.6], 
        Directive[
         PointSize[0.012833333333333334`], 
         RGBColor[0.368417, 0.506779, 0.709798], 
         AbsoluteThickness[1.6]], 
        Point[CompressedData["
1:eJxTTMoPSmViYGBQBmIQzfzI9caCLS/s530V3LvK+ri9Cdfh5BWpL+xnB916
yei1yL6oIV1cxOO5vYWKrOcm6d329RwcD/akPLW3Ve2oVO0+Z2/SLLN41vFH
9mq2HpZm367Z/9V81/Gb7b69nI2+4rl5d+199sU33wi7bl9w58xc39aH9vnr
xHQYXc7aH5b7ffJ182P7sNTLZz8pbLFffWDxvBlLnsD4+2F8qPr9MPVQ8/bD
zIPatx9mH9Q9+2Hugbp3vz3UvQUQ/+yH+cce4t/9isEQ/356CA6P/TbQ8MiK
EXGp0Xux//rihTYO/6/Yi01um1b46Nl+/SvTjzEK37OPLztwJGzfk/27ZuUb
Nwc/thfMO/llP/Oj/T8+cFaLWT2331h4XvP+vTv7iwunp6RNemXfr1537XfA
1f0pUbanlxS+tfessf3ot//4/vD51rqqme/tg7awCne/2bDfd/X9lGNzPtgD
o6iB+6zLrlnT3goD2Q6Wm8Hy9k7TIfKFhmD99pJZNmD9Z13A5tv3mEDM504H
22+v0Aax/14p2H32hp8g7vteCHa/vfd8iPtlJoH9Z58C9d9GiP/tj0D9DwCe
0Btx
         "]]}, {}}, {
       Tooltip[{
         Opacity[0.], 
         Point[{{0.7720493530351294, 0.18150839073124997`}}]}, {
        0.7720493530351294, 0.18150839073124997`}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.7624093960260496, 0.03572088222168653}}]}, {
        0.7624093960260496, 0.03572088222168653}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.7275486428370768, 0.10589136395459853`}}]}, {
        0.7275486428370768, 0.10589136395459853`}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.6685470939267047, 0.2386214106957513}}]}, {
        0.6685470939267047, 0.2386214106957513}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.5868657289284784, 0.3587775763889617}}]}, {
        0.5868657289284784, 0.3587775763889617}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.4848011815271603, 0.4628176998202792}}]}, {
        0.4848011815271603, 0.4628176998202792}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.3646756442021555, 0.5475223619495306}}]}, {
        0.3646756442021555, 0.5475223619495306}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.22863783505214338`, 0.6098536438884924}}]}, {
        0.22863783505214338`, 0.6098536438884924}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.078627753496581, 0.6450923059874375}}]}, {
        0.078627753496581, 0.6450923059874375}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.078627753496581, 
          0.6450923059874375}}]}, {-0.078627753496581, 0.6450923059874375}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.22863783505214338`, 
          0.6098536438884924}}]}, {-0.22863783505214338`, 
        0.6098536438884924}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.3646756442021555, 
          0.5475223619495306}}]}, {-0.3646756442021555, 0.5475223619495306}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.4848011815271603, 
          0.4628176998202792}}]}, {-0.4848011815271603, 0.4628176998202792}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.5868657289284784, 
          0.3587775763889617}}]}, {-0.5868657289284784, 0.3587775763889617}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.6685470939267047, 
          0.23862141069575135`}}]}, {-0.6685470939267047, 
        0.23862141069575135`}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.7275486428370765, 
          0.10589136395459853`}}]}, {-0.7275486428370765, 
        0.10589136395459853`}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.7624093960260508, 
          0.03572088222168746}}]}, {-0.7624093960260508, 
        0.03572088222168746}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.7720493530351276, 
          0.18150839073124725`}}]}, {-0.7720493530351276, 
        0.18150839073124725`}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.755674488990633, 
          0.32807928010020987`}}]}, {-0.755674488990633, 
        0.32807928010020987`}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.7151420535998081, 
          0.469910091341549}}]}, {-0.7151420535998081, 0.469910091341549}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.6482347338494029, 
          0.6039443974121561}}]}, {-0.6482347338494029, 0.6039443974121561}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.562957742785729, 
          0.7258407977201662}}]}, {-0.562957742785729, 0.7258407977201662}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.4511030109695425, 
          0.8303710903815457}}]}, {-0.4511030109695425, 0.8303710903815457}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.33306785532533295`, 
          0.9201225251698131}}]}, {-0.33306785532533295`, 
        0.9201225251698131}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.18552576808267165`, 
          0.9815851100401954}}]}, {-0.18552576808267165`, 
        0.9815851100401954}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-0.06610936369334411, 
          1.0382751407984927`}}]}, {-0.06610936369334411, 
        1.0382751407984927`}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{-5.250899696766081*^-28, 
          2.0000001484757606`}}]}, {-5.250899696766081*^-28, 
        2.0000001484757606`}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.06610936369334021, 1.0382751407973534`}}]}, {
        0.06610936369334021, 1.0382751407973534`}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.18552576808213986`, 0.9815851100459582}}]}, {
        0.18552576808213986`, 0.9815851100459582}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.3330678553257485, 0.9201225251687375}}]}, {
        0.3330678553257485, 0.9201225251687375}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.45110301096939115`, 0.8303710903821333}}]}, {
        0.45110301096939115`, 0.8303710903821333}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.5629577427859507, 0.7258407977202009}}]}, {
        0.5629577427859507, 0.7258407977202009}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.6482347338492777, 0.6039443974122859}}]}, {
        0.6482347338492777, 0.6039443974122859}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.7151420535997803, 0.4699100913415519}}]}, {
        0.7151420535997803, 0.4699100913415519}], 
       Tooltip[{
         Opacity[0.], 
         Point[{{0.7556744889906409, 0.3280792801002088}}]}, {
        0.7556744889906409, 0.3280792801002088}]}}, {}, {}, {{}, {}}}, {
    ImageSize -> Large, DisplayFunction -> Identity, PlotRangePadding -> {{
        Scaled[0.02], 
        Scaled[0.02]}, {
        Scaled[0.02], 
        Scaled[0.05]}}, AxesOrigin -> {0, 0}, 
     PlotRange -> {{-0.7720493530351276, 0.7720493530351294}, {
       0, 2.0000001484757606`}}, PlotRangeClipping -> True, ImagePadding -> 
     All, DisplayFunction -> Identity, AspectRatio -> 1, Axes -> {True, True},
      AxesLabel -> {None, None}, 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 -> {"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
           (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
            Part[#, 1]], 
           (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
            Part[#, 2]]}& ), "CopiedValueFunction" -> ({
           (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[
            Part[#, 1]], 
           (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[
            Part[#, 2]]}& )}}, 
     PlotRange -> {{-0.7720493530351276, 0.7720493530351294}, {
       0, 2.0000001484757606`}}, PlotRangeClipping -> True, 
     PlotRangePadding -> {{
        Scaled[0.02], 
        Scaled[0.02]}, {
        Scaled[0.02], 
        Scaled[0.05]}}, Ticks -> {Automatic, Automatic}}], "\t", 
   TableForm[CompressedData["
1:eJxTTMoPSmViYGBQBmIQzfzI9caCLS/s530V3LvK+vh+E67DyStSX9jPDrr1
ktFr0f6ihnRxEY/n9hYqsp6bpHfb13NwPNiT8tTeVrWjUrX7nL1Js8ziWccf
2avZeliafbtm/1fzXcdvtvv2cjb6iufm3bX32RfffCPsun3BnTNzfVsf2uev
E9NhdDlrf1ju98nXzY/tw1Ivn/2ksMV+9YHF82YseQLj74fxoer3w9RDzdsP
Mw9q336YfVD37Ie5B+re/fZQ9xZA/LMf5h97iH/3KwZD/PvpITg89ttAwyMr
RsSlRu/F/uuLF9o4/L+yX2xy27TCR8/261+ZfoxR+N7++LIDR8L2Pdm/a1a+
cXPw4/2CeSe/7Gd+tP/HB85qMavn+zcWnte8f+/O/uLC6Slpk17t71evu/Y7
4Or+lCjb00sK3+73rLH96Lf/+P7w+da6qpnv9wdtYRXufrNhv+/q+ynH5nzY
D4yiBu6zLrtmTXsrDGQfsNwMlrd3mg6RLzQE67eXzLIB6z/rAjbfvscEYj53
Oth+e4U2iP33SsHuszf8BHHf90Kw++2950PcLzMJ7D/7FKj/NkL8b38E6n8A
Gcgl8Q==
    "], TableHeadings -> {{1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 
      16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 
      34, 35}, {"\[Omega]R", "\[Omega]I"}}]],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.7449661161497536`*^9, {3.7449661629983263`*^9, 3.7449662032066207`*^9}, {
   3.7449663070080557`*^9, 3.7449663704848843`*^9}, 3.7449855269340243`*^9, 
   3.7479819146235538`*^9, 3.74798370155177*^9, 
   3.747993793220028*^9},ExpressionUUID->"ad046499-cdec-4488-9c59-\
b11f41555054"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"RandomReal", "[", 
  RowBox[{"{", 
   RowBox[{"0.01", ",", "12"}], "}"}], "]"}]], "Input",
 CellChangeTimes->{{3.747983717208945*^9, 
  3.7479837515588474`*^9}},ExpressionUUID->"0840f0cc-db28-4c5c-8d27-\
8a81c6f5c08f"],

Cell[BoxData["8.940970958824291`"], "Output",
 CellChangeTimes->{{3.7479837277359915`*^9, 
  3.7479837541905384`*^9}},ExpressionUUID->"099971d2-eb7b-42e9-8059-\
460fb84e2cb6"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"N", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"(", 
      RowBox[{"\[Mu]", "/", "\[Beta]"}], ")"}], "^", "\[Beta]"}], 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"\[Beta]", "-", "1"}], ")"}], "^", 
     RowBox[{"(", 
      RowBox[{"\[Beta]", "-", "1"}], ")"}]}]}], "]"}], "/.", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"\[Mu]", "\[Rule]", "1.81"}], ",", 
    RowBox[{"\[Beta]", "\[Rule]", "1.51"}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.7479838398909945`*^9, 3.7479838494334145`*^9}, {
  3.747983900811057*^9, 
  3.747983903752104*^9}},ExpressionUUID->"f2c78e73-404c-4c46-abec-\
a8c254e1c93f"],

Cell[BoxData["0.9326103855513533`"], "Output",
 CellChangeTimes->{{3.747983850060918*^9, 3.747983854838538*^9}, 
   3.7479839042704363`*^9},ExpressionUUID->"051e7292-9655-48a5-9767-\
3e521503b376"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"RandomReal", "[", 
   RowBox[{"{", 
    RowBox[{"0.01", ",", 
     RowBox[{"Re", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{"\[Mu]", "/", "\[Beta]"}], ")"}], "^", "\[Beta]"}], 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"\[Beta]", "-", "1"}], ")"}], "^", 
        RowBox[{"(", 
         RowBox[{"\[Beta]", "-", "1"}], ")"}]}]}], "]"}]}], "}"}], "]"}], "/.", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"\[Mu]", "\[Rule]", "1.81"}], ",", 
    RowBox[{"\[Beta]", "\[Rule]", "1.51"}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.7479838224647427`*^9, 3.7479838279163547`*^9}, {
  3.7479838636350718`*^9, 
  3.7479838690868015`*^9}},ExpressionUUID->"e7096b41-9ce0-44bd-843f-\
1d92a76971d5"],

Cell[BoxData[
 TemplateBox[{
  "RandomReal","unifr",
   "\"The endpoints specified by \\!\\(\\*RowBox[{\\\"{\\\", \
RowBox[{\\\"0.01`\\\", \\\",\\\", RowBox[{\\\"Re\\\", \\\"[\\\", \
RowBox[{SuperscriptBox[RowBox[{\\\"(\\\", RowBox[{RowBox[{\\\"-\\\", \
\\\"1\\\"}], \\\"+\\\", \\\"\[Beta]\\\"}], \\\")\\\"}], RowBox[{RowBox[{\\\"-\
\\\", \\\"1\\\"}], \\\"+\\\", \\\"\[Beta]\\\"}]], \\\" \\\", \
SuperscriptBox[RowBox[{\\\"(\\\", FractionBox[\\\"\[Mu]\\\", \
\\\"\[Beta]\\\"], \\\")\\\"}], \\\"\[Beta]\\\"]}], \\\"]\\\"}]}], \\\"}\\\"}]\
\\) for the endpoints of the uniform distribution range are not real \
valued.\"",2,69,17,24562841768407440302,"Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{
  3.7479937984564376`*^9},ExpressionUUID->"54d56fcb-ee82-4243-9816-\
890adaa836c5"],

Cell[BoxData["0.09276982682101342`"], "Output",
 CellChangeTimes->{3.7479838291179447`*^9, 3.7479838697720957`*^9, 
  3.7479839619344497`*^9, 
  3.7479937984564376`*^9},ExpressionUUID->"95c899de-bc64-4269-936c-\
3e966294b4a9"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Sqrt", "[", 
  RowBox[{"PACAGetModes", "[", 
   RowBox[{"55", ",", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"1", "-", 
       RowBox[{"\[Mu]", " ", 
        RowBox[{"x", "^", "2"}]}], "+", 
       RowBox[{"Q", " ", 
        RowBox[{"x", "^", 
         RowBox[{"(", 
          RowBox[{"2", "\[Beta]"}], ")"}]}]}]}], "/.", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"\[Mu]", "\[Rule]", "2"}], ",", 
        RowBox[{"\[Beta]", "\[Rule]", 
         RowBox[{"3", "/", "2"}]}], ",", 
        RowBox[{"Q", "\[Rule]", "0.2"}]}], "}"}]}], ")"}], ",", "x", ",", "1",
     ",", "1"}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.747982024220202*^9, 3.747982032422662*^9}, {
  3.747982154123617*^9, 3.747982243956674*^9}, {3.7479823088810353`*^9, 
  3.7479823564203053`*^9}, {3.747982415676242*^9, 3.747982435086567*^9}, {
  3.7479837593576775`*^9, 3.7479838017939105`*^9}, {3.7479839410600233`*^9, 
  3.747983942407189*^9}, {3.74798398061915*^9, 3.747983980885952*^9}, {
  3.7479840430733995`*^9, 3.7479840588615837`*^9}, {3.747984106407897*^9, 
  3.7479841814651413`*^9}, {3.7479842221489573`*^9, 3.7479842238881087`*^9}, {
  3.7479938923343253`*^9, 3.7479938932226086`*^9}, {3.747995484006228*^9, 
  3.7479954875459585`*^9}},ExpressionUUID->"d8244529-2680-459b-880d-\
8cad8d9b37f6"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"0.9241242584892384`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.07109607374406075`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.9193044339584702`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.016908698472529887`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.9115607573009974`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.03721306192329284`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.9008693228718249`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.0911737667607265`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.887293732116698`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.14484069502716143`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.8708714483024131`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.19805391920253784`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.851636832073173`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.250685249486089`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.8296464797375369`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.30262147796400246`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.8049619797166121`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.3537572025652447`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.7776354895172917`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.4040133628776468`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.7476958163729499`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.4533579137364896`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.7150975630477112`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.5018269535854365`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.6795603915381734`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.5495185452876544`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.6402169707823141`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.5962453579896956`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.5962453579896956`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.6402169707823141`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.5495185452876544`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.6795603915381734`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.5018269535854365`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.7150975630477112`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.4533579137364896`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.7476958163729499`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.4040133628776468`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.7776354895172917`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.3537572025652447`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.8049619797166121`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.30262147796400246`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.8296464797375369`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.250685249486089`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.851636832073173`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.19805391920253787`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.8708714483024131`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.14484069502716146`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.887293732116698`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.09117376676072637`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.9008693228718252`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.03721306192329309`", "\[VeryThinSpace]", "+", 
    RowBox[{"0.9115607573009973`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.016908698472530233`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.9193044339584695`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.07109607374405949`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.9241242584892387`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.1250902502668529`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.9260652644103214`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.17875103689584984`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.9249081956207398`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.2322572882197549`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.9208954020598769`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.2847664487456848`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.9142974985552061`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.3365565956194252`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.9038949296265509`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.3886498307436505`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.8917551319521437`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.4368758327775464`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.8769046137415837`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.4872623725727371`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.8560521419301602`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.5356423687910865`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.8405308984166744`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.5761264862667601`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.811506857671763`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.6347023509532584`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.7878275937344071`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.6562971300421209`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.7668489262423376`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"1.0000000694428484`", "\[VeryThinSpace]", "-", 
    RowBox[{"1.0000000694428484`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.7327429577432387`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.7327429577432387`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.7115367136507559`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.7115367136507557`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.7668489262565812`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.6562971300543624`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.787827593735698`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.6347023509543694`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.8115068576712351`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.5761264862656844`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.84053089841662`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.5356423687930425`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.8560521419303804`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.48726237257254534`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.8769046137422254`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.436875832777978`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.8917551319516017`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.3886498307438504`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.9038949296265637`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.33655659561957674`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.9142974985552577`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.28476644874567864`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.9208954020598439`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.2322572882197859`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.9249081956207482`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.17875103689582947`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.9260652644103193`", "\[VeryThinSpace]", "-", 
    RowBox[{"0.1250902502668469`", " ", "\[ImaginaryI]"}]}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.747982244782739*^9, {3.747982334034235*^9, 3.747982363949025*^9}, {
   3.7479824241309457`*^9, 3.7479824396704755`*^9}, 3.7479837924873505`*^9, {
   3.7479839463981886`*^9, 3.7479839823801146`*^9}, {3.7479840451355*^9, 
   3.747984061793767*^9}, 3.747984146156666*^9, 3.7479841829134383`*^9, 
   3.7479842260327015`*^9, 3.7479937544338183`*^9, 3.74799380321996*^9, 
   3.7479938969573603`*^9, 3.7479947722274876`*^9, {3.747995464539677*^9, 
   3.747995493583935*^9}},ExpressionUUID->"f534f25e-6850-4081-9114-\
1151693a6acf"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"N", "[", 
  RowBox[{"3", "/", "2"}], "]"}]], "Input",
 CellChangeTimes->{{3.747995117324941*^9, 
  3.747995128411484*^9}},ExpressionUUID->"036f8525-af8c-48e1-a4f3-\
f1af2ac215c3"],

Cell[BoxData["1.5`"], "Output",
 CellChangeTimes->{{3.7479951212428923`*^9, 
  3.747995129018222*^9}},ExpressionUUID->"63b737d1-4d98-4deb-91b5-\
f69ac2a77399"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"(", 
  RowBox[{
   RowBox[{"1", "-", 
    RowBox[{"rp", " ", "x"}], "-", 
    RowBox[{"2", " ", 
     RowBox[{
      RowBox[{"rp", "^", "3"}], "/", 
      RowBox[{"R1", "^", "2"}]}], "x"}], "-", 
    RowBox[{
     RowBox[{
      RowBox[{"Q", "^", "2"}], "/", "rp"}], " ", "x"}], "+", 
    RowBox[{
     RowBox[{"Q", "^", "2"}], " ", 
     RowBox[{"x", "^", "2"}]}], "+", 
    RowBox[{
     RowBox[{"1", "/", 
      RowBox[{"x", "^", "2"}]}], "/", 
     RowBox[{"R1", "^", "2"}]}]}], "/.", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"rp", "\[Rule]", "50"}], ",", 
     RowBox[{"R1", "\[Rule]", "1"}], ",", 
     RowBox[{"Q", "\[Rule]", "12"}]}], "}"}]}], ")"}]], "Input",ExpressionUUID\
->"775f6539-6c69-4d87-a122-2a96c8318c01"],

Cell[BoxData[
 RowBox[{"1", "+", 
  FractionBox["1", 
   SuperscriptBox["x", "2"]], "-", 
  FractionBox[
   RowBox[{"6251322", " ", "x"}], "25"], "+", 
  RowBox[{"144", " ", 
   SuperscriptBox["x", "2"]}]}]], "Output",
 CellChangeTimes->{
  3.7491234832980843`*^9},ExpressionUUID->"35bcbd14-256f-4d13-946b-\
f4d0504c87ab"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"PACAGetModes", "[", 
  RowBox[{"3", ",", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"1", "-", 
      RowBox[{"rp", " ", "x"}], "-", 
      RowBox[{"2", " ", 
       RowBox[{
        RowBox[{"rp", "^", "3"}], "/", 
        RowBox[{"R1", "^", "2"}]}], "x"}], "-", 
      RowBox[{
       RowBox[{
        RowBox[{"Q", "^", "2"}], "/", "rp"}], " ", "x"}], "+", 
      RowBox[{
       RowBox[{"Q", "^", "2"}], " ", 
       RowBox[{"x", "^", "2"}]}], "+", 
      RowBox[{
       RowBox[{"1", "/", 
        RowBox[{"x", "^", "2"}]}], "/", 
       RowBox[{"R1", "^", "2"}]}]}], "/.", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"rp", "\[Rule]", "1"}], ",", 
       RowBox[{"R1", "\[Rule]", "0.5"}], ",", 
       RowBox[{"Q", "\[Rule]", "20"}]}], "}"}]}], ")"}], ",", "x", ",", "0", 
   ",", "0"}], "]"}]], "Input",
 CellChangeTimes->{{3.749119260681198*^9, 3.749119262009544*^9}, {
  3.7491212605588923`*^9, 3.749121301693878*^9}, {3.7491213322502203`*^9, 
  3.74912142164489*^9}, {3.7491229303013153`*^9, 3.749122955388872*^9}, {
  3.749123505121518*^9, 3.749123537212034*^9}, {3.749123568624596*^9, 
  3.749123591908461*^9}, {3.749123623690123*^9, 3.749123696088193*^9}, {
  3.749123792619606*^9, 3.74912379338035*^9}, {3.7491238411707573`*^9, 
  3.749123841783115*^9}},ExpressionUUID->"e2d58c15-91cb-4102-910f-\
1d1ed600f73d"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "7289.331566999999`"}], "-", 
    RowBox[{"11500.486477`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"7289.331566999999`", "\[VeryThinSpace]", "-", 
    RowBox[{"11500.486477`", " ", "\[ImaginaryI]"}]}], ",", 
   RowBox[{"0.`", "\[VeryThinSpace]", "-", 
    RowBox[{"1423.1970469999999`", " ", "\[ImaginaryI]"}]}]}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.7491213412406054`*^9, 3.749121382302932*^9}, {
   3.7491214159978113`*^9, 3.749121422303488*^9}, {3.7491229312833853`*^9, 
   3.749122964904944*^9}, {3.749123506679833*^9, 3.74912353760471*^9}, {
   3.749123569628264*^9, 3.749123593381288*^9}, {3.749123633366984*^9, 
   3.749123697026827*^9}, 3.74912379565432*^9, 
   3.749123842136842*^9},ExpressionUUID->"fd52eb57-e57c-438b-98f8-\
364e2d3c9a20"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Table", "[", 
  RowBox[{"r", ",", 
   RowBox[{"{", 
    RowBox[{"r", ",", "1", ",", "3", ",", "1"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.749124409273217*^9, 
  3.749124418753748*^9}},ExpressionUUID->"cddda6b2-3a1d-4c2e-a70f-\
e942836dd46c"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1", ",", "2", ",", "3"}], "}"}]], "Output",
 CellChangeTimes->{
  3.749124419123248*^9},ExpressionUUID->"40399bd2-5e28-4901-8b7b-\
bae32bf453cb"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Table", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"50", "x0"}], ",", 
     RowBox[{
      RowBox[{"Im", "[", 
       RowBox[{"PACAGetModes", "[", 
        RowBox[{"3", ",", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"1", "-", 
            RowBox[{"rp", " ", "x"}], "-", 
            RowBox[{"2", " ", 
             RowBox[{
              RowBox[{"rp", "^", "3"}], "/", 
              RowBox[{"R1", "^", "2"}]}], "x"}], "-", 
            RowBox[{
             RowBox[{
              RowBox[{"Q", "^", "2"}], "/", "rp"}], " ", "x"}], "+", 
            RowBox[{
             RowBox[{"Q", "^", "2"}], " ", 
             RowBox[{"x", "^", "2"}]}], "+", 
            RowBox[{
             RowBox[{"1", "/", 
              RowBox[{"x", "^", "2"}]}], "/", 
             RowBox[{"R1", "^", "2"}]}]}], "/.", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"rp", "\[Rule]", "x0"}], ",", 
             RowBox[{"R1", "\[Rule]", "0.5"}], ",", 
             RowBox[{"Q", "\[Rule]", "10"}]}], "}"}]}], ")"}], ",", "x", ",", 
         "0", ",", "0"}], "]"}], "]"}], "[", 
      RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"x0", ",", "2", ",", "4"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.7491244039335413`*^9, 3.7491244616046343`*^9}, {
  3.7491246295122004`*^9, 3.749124638844488*^9}, {3.749124696376362*^9, 
  3.7491247444188347`*^9}},ExpressionUUID->"79c81e40-ae6f-42e6-afc0-\
14977f409d67"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"100", ",", 
     RowBox[{"-", "4540.38578`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"150", ",", 
     RowBox[{"-", "116242.230172`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"200", ",", 
     RowBox[{"-", "2.496466896205`*^6"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.749124433860552*^9, 3.7491244630570593`*^9}, {
   3.7491246355177603`*^9, 3.749124640107574*^9}, {3.749124704385647*^9, 
   3.749124727236231*^9}, 
   3.749124766317575*^9},ExpressionUUID->"977c75ee-6749-4eed-8eb0-\
ba1e8289e95e"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Table", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"50", "x0"}], ",", 
     RowBox[{
      RowBox[{"Im", "[", 
       RowBox[{"PACAGetModes", "[", 
        RowBox[{"3", ",", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"1", "-", 
            RowBox[{"rp", " ", "x"}], "-", 
            RowBox[{"2", " ", 
             RowBox[{
              RowBox[{"rp", "^", "3"}], "/", 
              RowBox[{"R1", "^", "2"}]}], "x"}], "-", 
            RowBox[{
             RowBox[{
              RowBox[{"Q", "^", "2"}], "/", "rp"}], " ", "x"}], "+", 
            RowBox[{
             RowBox[{"Q", "^", "2"}], " ", 
             RowBox[{"x", "^", "2"}]}], "+", 
            RowBox[{
             RowBox[{"1", "/", 
              RowBox[{"x", "^", "2"}]}], "/", 
             RowBox[{"R1", "^", "2"}]}]}], "/.", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"rp", "\[Rule]", "x0"}], ",", 
             RowBox[{"R1", "\[Rule]", "0.5"}], ",", 
             RowBox[{"Q", "\[Rule]", "20"}]}], "}"}]}], ")"}], ",", "x", ",", 
         "0", ",", "0"}], "]"}], "]"}], "[", 
      RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"x0", ",", "2", ",", "4"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.7491244750741653`*^9, 3.749124482901102*^9}, {
  3.749124774541111*^9, 
  3.749124776599297*^9}},ExpressionUUID->"75e118a9-8e4e-46bf-b7c5-\
67abd15ddbaf"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"100", ",", 
     RowBox[{"-", "1356.4742019999999`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"150", ",", 
     RowBox[{"-", "6015.878108`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"200", ",", 
     RowBox[{"-", "63847.211651`"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.749124475925248*^9, 3.749124485248795*^9}, 
   3.749124777703553*^9},ExpressionUUID->"312ffab3-47fd-474c-975c-\
5026a59a07cc"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Table", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"50", "x0"}], ",", 
     RowBox[{
      RowBox[{"Im", "[", 
       RowBox[{"PACAGetModes", "[", 
        RowBox[{"3", ",", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"1", "-", 
            RowBox[{"rp", " ", "x"}], "-", 
            RowBox[{"2", " ", 
             RowBox[{
              RowBox[{"rp", "^", "3"}], "/", 
              RowBox[{"R1", "^", "2"}]}], "x"}], "-", 
            RowBox[{
             RowBox[{
              RowBox[{"Q", "^", "2"}], "/", "rp"}], " ", "x"}], "+", 
            RowBox[{
             RowBox[{"Q", "^", "2"}], " ", 
             RowBox[{"x", "^", "2"}]}], "+", 
            RowBox[{
             RowBox[{"1", "/", 
              RowBox[{"x", "^", "2"}]}], "/", 
             RowBox[{"R1", "^", "2"}]}]}], "/.", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"rp", "\[Rule]", "x0"}], ",", 
             RowBox[{"R1", "\[Rule]", "0.5"}], ",", 
             RowBox[{"Q", "\[Rule]", "25"}]}], "}"}]}], ")"}], ",", "x", ",", 
         "0", ",", "0"}], "]"}], "]"}], "[", 
      RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"x0", ",", "2", ",", "4"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.749124516062315*^9, 3.749124546491434*^9}, {
  3.7491247851044083`*^9, 
  3.749124787209592*^9}},ExpressionUUID->"58389a1e-8818-478c-b173-\
c29318ab8956"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"100", ",", 
     RowBox[{"-", "1175.738355`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"150", ",", 
     RowBox[{"-", "2554.120918`"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"200", ",", 
     RowBox[{"-", "23034.468945999997`"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.749124523079372*^9, 3.749124547134993*^9}, 
   3.749124787821666*^9},ExpressionUUID->"02f42751-523d-45fd-8812-\
950d05c2a15d"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"N", "[", 
  RowBox[{"100", 
   RowBox[{"Sqrt", "[", "10", "]"}]}], "]"}]], "Input",
 CellChangeTimes->{
  3.749124522699463*^9, {3.749125559287343*^9, 
   3.74912556818382*^9}},ExpressionUUID->"b0ab274b-0b39-4a3b-9750-\
accd7c647701"],

Cell[BoxData["316.22776601683796`"], "Output",
 CellChangeTimes->{
  3.7491255688755417`*^9},ExpressionUUID->"650b202f-6804-4ef5-888b-\
ee0d3a8a309a"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Block", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"{", "}"}], ",", "\[IndentingNewLine]", 
   RowBox[{"Show", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"ListLinePlot", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"100", ",", 
          RowBox[{"-", "4540.385"}]}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"10", "^", "2.5"}], ",", 
          RowBox[{"-", "116242.23"}]}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"1000", ",", 
          RowBox[{
           RowBox[{"-", "2.496"}], "*", 
           RowBox[{"10", "^", "6"}]}]}], "}"}]}], "}"}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"ListLinePlot", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"100", ",", 
           RowBox[{"-", "1356.4742019999999`"}]}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"10", "^", "2.5"}], ",", 
           RowBox[{"-", "6015.878108`"}]}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"1000", ",", 
           RowBox[{"-", "63847.211651`"}]}], "}"}]}], "}"}], ",", 
       RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"ListLinePlot", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"100", ",", 
           RowBox[{"-", "1175.738355`"}]}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"10", "^", "2.5"}], ",", 
           RowBox[{"-", "2554.120918`"}]}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"1000", ",", 
           RowBox[{"-", "23034.468945999997`"}]}], "}"}]}], "}"}], ",", 
       RowBox[{"PlotStyle", "\[Rule]", "Green"}]}], "]"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"Automatic", ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", 
           RowBox[{"10", "^", "6"}]}], ",", "0"}], "}"}]}], "}"}]}], ",", 
     RowBox[{"AxesLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"\"\<r+\>\"", ",", "\"\<Im[\[Omega]]\>\""}], "}"}]}]}], 
    "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", "]"}]], "Input",
 CellChangeTimes->{{3.749123700054109*^9, 3.749123728744789*^9}, {
  3.749123772162978*^9, 3.749123789614767*^9}, {3.7491238631347923`*^9, 
  3.7491241076507807`*^9}, {3.7491244876006393`*^9, 3.749124499757166*^9}, {
  3.749124557351367*^9, 3.749124558442968*^9}, {3.749124626205068*^9, 
  3.7491246526738253`*^9}, {3.749124822312414*^9, 3.749125100791327*^9}, {
  3.749125625045897*^9, 
  3.7491256298718243`*^9}},ExpressionUUID->"4743a23a-d09a-4741-a00e-\
530acbeea330"],

Cell[BoxData[
 GraphicsBox[{{{}, {{}, {}, 
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
      0.019444444444444445`], AbsoluteThickness[1.6], 
      LineBox[{{100., -4540.385}, {316.22776601683796`, -116242.23}, {
       1000., -2.496*^6}}]}}, {}, {}, {}, {}}, {{}, {{}, {}, 
     {RGBColor[1, 0, 0], PointSize[0.019444444444444445`], AbsoluteThickness[
      1.6], LineBox[{{100., -1356.4742019999999`}, {
       316.22776601683796`, -6015.878108}, {
       1000., -63847.211651}}]}}, {}, {}, {}, {}}, {{}, {{}, {}, 
     {RGBColor[0, 1, 0], PointSize[0.019444444444444445`], AbsoluteThickness[
      1.6], LineBox[{{100., -1175.738355}, {
       316.22776601683796`, -2554.120918}, {
       1000., -23034.468945999997`}}]}}, {}, {}, {}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{
    FormBox["\"r+\"", TraditionalForm], 
    FormBox["\"Im[\[Omega]]\"", 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]],
  ImagePadding->All,
  Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{Automatic, {-1000000, 0}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.02]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.749124396659466*^9, 3.74912450225943*^9, 3.749124560151183*^9, 
   3.749124654535553*^9, 3.749124902406542*^9, 3.749124940565027*^9, {
   3.749125024618675*^9, 3.749125052067357*^9}, 3.7491251046115437`*^9, 
   3.749125630965082*^9},ExpressionUUID->"19a24819-f318-4064-9223-\
7cad86a9a80a"]
}, Open  ]]
},
WindowSize->{808, 751},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"11.2 for Mac OS X x86 (32-bit, 64-bit Kernel) (September \
10, 2017)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 16471, 389, 1144, "Input",ExpressionUUID->"c618399f-f8df-4a55-9ad5-ffac92a4333b"],
Cell[CellGroupData[{
Cell[17054, 413, 371, 9, 30, "Input",ExpressionUUID->"24262446-3bf0-494d-b923-b6e34627688b"],
Cell[17428, 424, 553, 12, 37, "Output",ExpressionUUID->"bc7b092e-409d-4e85-88ec-d3a80d8900c8"]
}, Open  ]],
Cell[CellGroupData[{
Cell[18018, 441, 302, 8, 30, "Input",ExpressionUUID->"84759498-767a-4843-818e-86b0b9fd7ea5"],
Cell[18323, 451, 951, 21, 60, "Output",ExpressionUUID->"57c39d3a-ebb9-49c3-b75f-81d574c074b7"]
}, Open  ]],
Cell[19289, 475, 15535, 365, 1324, "Input",ExpressionUUID->"1ee244e2-427e-48a3-b3dd-5109a29a98c9"],
Cell[CellGroupData[{
Cell[34849, 844, 3065, 57, 115, "Input",ExpressionUUID->"eb71481b-61af-462e-b1c4-10bb29f4fad0"],
Cell[37917, 903, 918, 18, 34, "Output",ExpressionUUID->"5aa78c3f-5602-42cb-8b52-e65758bc3d21"],
Cell[38838, 923, 553, 8, 24, "Print",ExpressionUUID->"91a72e39-aa1d-4f66-a249-509d0017dc04"]
}, Open  ]],
Cell[39406, 934, 1107, 29, 30, "Input",ExpressionUUID->"c72a2cb8-2837-4f31-a1b2-0030f61b6a56"],
Cell[CellGroupData[{
Cell[40538, 967, 1072, 24, 30, "Input",ExpressionUUID->"97a87942-c0a5-4551-bc5a-5bc7c895b364"],
Cell[41613, 993, 13164, 306, 385, "Output",ExpressionUUID->"f7c71fb8-c2cc-48be-abc2-2953179e0ef4"]
}, Open  ]],
Cell[54792, 1302, 1817, 44, 136, "Input",ExpressionUUID->"1ae650d2-7bf3-474d-befe-66e26555971c"],
Cell[CellGroupData[{
Cell[56634, 1350, 678, 19, 30, "Input",ExpressionUUID->"ef264675-4148-453b-8957-a1c8578d0a32"],
Cell[57315, 1371, 4318, 138, 461, "Output",ExpressionUUID->"0d20787a-cbad-483c-ad21-53876bc92d8d"]
}, Open  ]],
Cell[61648, 1512, 899, 22, 73, "Input",ExpressionUUID->"8d7c1a55-460c-4b70-9151-1c2c9bd2d128"],
Cell[CellGroupData[{
Cell[62572, 1538, 833, 21, 30, "Input",ExpressionUUID->"065c127c-0626-4ef4-bd10-1723e1c6a173"],
Cell[63408, 1561, 28570, 730, 1231, "Print",ExpressionUUID->"ad046499-cdec-4488-9c59-b11f41555054"]
}, Open  ]],
Cell[CellGroupData[{
Cell[92015, 2296, 242, 6, 30, "Input",ExpressionUUID->"0840f0cc-db28-4c5c-8d27-8a81c6f5c08f"],
Cell[92260, 2304, 175, 3, 34, "Output",ExpressionUUID->"099971d2-eb7b-42e9-8059-460fb84e2cb6"]
}, Open  ]],
Cell[CellGroupData[{
Cell[92472, 2312, 638, 19, 30, "Input",ExpressionUUID->"f2c78e73-404c-4c46-abec-a8c254e1c93f"],
Cell[93113, 2333, 197, 3, 34, "Output",ExpressionUUID->"051e7292-9655-48a5-9767-3e521503b376"]
}, Open  ]],
Cell[CellGroupData[{
Cell[93347, 2341, 762, 22, 30, "Input",ExpressionUUID->"e7096b41-9ce0-44bd-843f-1d92a76971d5"],
Cell[94112, 2365, 803, 15, 44, "Message",ExpressionUUID->"54d56fcb-ee82-4243-9816-890adaa836c5"],
Cell[94918, 2382, 226, 4, 34, "Output",ExpressionUUID->"95c899de-bc64-4269-936c-3e966294b4a9"]
}, Open  ]],
Cell[CellGroupData[{
Cell[95181, 2391, 1315, 29, 30, "Input",ExpressionUUID->"d8244529-2680-459b-880d-8cad8d9b37f6"],
Cell[96499, 2422, 7599, 122, 182, "Output",ExpressionUUID->"f534f25e-6850-4081-9114-1151693a6acf"]
}, Open  ]],
Cell[CellGroupData[{
Cell[104135, 2549, 203, 5, 30, "Input",ExpressionUUID->"036f8525-af8c-48e1-a4f3-f1af2ac215c3"],
Cell[104341, 2556, 159, 3, 34, "Output",ExpressionUUID->"63b737d1-4d98-4deb-91b5-f69ac2a77399"]
}, Open  ]],
Cell[CellGroupData[{
Cell[104537, 2564, 752, 24, 30, "Input",ExpressionUUID->"775f6539-6c69-4d87-a122-2a96c8318c01"],
Cell[105292, 2590, 322, 10, 51, "Output",ExpressionUUID->"35bcbd14-256f-4d13-946b-f4d0504c87ab"]
}, Open  ]],
Cell[CellGroupData[{
Cell[105651, 2605, 1352, 34, 73, "Input",ExpressionUUID->"e2d58c15-91cb-4102-910f-1d1ed600f73d"],
Cell[107006, 2641, 828, 17, 34, "Output",ExpressionUUID->"fd52eb57-e57c-438b-98f8-364e2d3c9a20"]
}, Open  ]],
Cell[CellGroupData[{
Cell[107871, 2663, 276, 7, 30, "Input",ExpressionUUID->"cddda6b2-3a1d-4c2e-a70f-e942836dd46c"],
Cell[108150, 2672, 185, 5, 34, "Output",ExpressionUUID->"40399bd2-5e28-4901-8b7b-bae32bf453cb"]
}, Open  ]],
Cell[CellGroupData[{
Cell[108372, 2682, 1500, 40, 115, "Input",ExpressionUUID->"79c81e40-ae6f-42e6-afc0-14977f409d67"],
Cell[109875, 2724, 592, 16, 37, "Output",ExpressionUUID->"977c75ee-6749-4eed-8eb0-ba1e8289e95e"]
}, Open  ]],
Cell[CellGroupData[{
Cell[110504, 2745, 1448, 40, 115, "Input",ExpressionUUID->"75e118a9-8e4e-46bf-b7c5-67abd15ddbaf"],
Cell[111955, 2787, 489, 14, 34, "Output",ExpressionUUID->"312ffab3-47fd-474c-975c-5026a59a07cc"]
}, Open  ]],
Cell[CellGroupData[{
Cell[112481, 2806, 1448, 40, 115, "Input",ExpressionUUID->"58389a1e-8818-478c-b173-c29318ab8956"],
Cell[113932, 2848, 488, 14, 34, "Output",ExpressionUUID->"02f42751-523d-45fd-8812-950d05c2a15d"]
}, Open  ]],
Cell[CellGroupData[{
Cell[114457, 2867, 258, 7, 30, "Input",ExpressionUUID->"b0ab274b-0b39-4a3b-9750-accd7c647701"],
Cell[114718, 2876, 150, 3, 34, "Output",ExpressionUUID->"650b202f-6804-4ef5-888b-ee0d3a8a309a"]
}, Open  ]],
Cell[CellGroupData[{
Cell[114905, 2884, 2747, 72, 241, "Input",ExpressionUUID->"4743a23a-d09a-4741-a00e-530acbeea330"],
Cell[117655, 2958, 2077, 50, 225, "Output",ExpressionUUID->"19a24819-f318-4064-9223-7cad86a9a80a"]
}, Open  ]]
}
]
*)