Spectrum.nb

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

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

(* CreatedBy='Mathematica 12.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     41466,        990]
NotebookOptionsPosition[     39996,        959]
NotebookOutlinePosition[     40359,        975]
CellTagsIndexPosition[     40316,        972]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "Determine", " ", "the", " ", "universal", " ", "function", " ", "g", " ", 
    "using", " ", "a", " ", "collocation", " ", "method", " ", "thus", " ", 
    "reducing", " ", "the", " ", "infinite", " ", "dimensional", " ", 
    "function", " ", "space", " ", "to", " ", "a", " ", "finite", " ", 
    "vector", " ", "whose", " ", "elements", " ", "are", " ", "g", " ", 
    "evaluated", " ", "at", " ", "its", " ", "Chebyshev", " ", 
    RowBox[{"nodes", ".", " ", "We"}], " ", "then", " ", "using", " ", 
    RowBox[{"Arnoldi", "'"}], "s", " ", "method", " ", "to", " ", "find", " ",
     "the", " ", "spectrum", " ", "of", " ", "the", " ", "linearized", " ", 
    "period", " ", "doubling", " ", 
    RowBox[{"transformation", "."}]}], " ", "*)"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "define", " ", "the", " ", "Chebyshev", " ", "nodes", " ", "and", " ", 
    "period", " ", "doubling", " ", "transformation"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"n", "=", "20"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"x", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Cos", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"2", "*", "i"}], "-", "1"}], ")"}], "*", "Pi"}], ")"}], 
         "/", 
         RowBox[{"(", 
          RowBox[{"2", "*", "n"}], ")"}]}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"\[Phi]", "[", 
      RowBox[{"f_", ",", "x_"}], "]"}], ":=", 
     RowBox[{
      RowBox[{
       RowBox[{"f", "[", "1", "]"}], "*", 
       RowBox[{"f", "[", "x", "]"}]}], "-", 
      RowBox[{"f", "[", 
       RowBox[{"f", "[", 
        RowBox[{
         RowBox[{"f", "[", "1", "]"}], "*", "x"}], "]"}], "]"}]}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "construct", " ", "f", " ", "from", " ", "its", " ", "Chebyshev", " ", 
     "coefficients"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"gc", "[", 
      RowBox[{"c_", ",", "y_"}], "]"}], ":=", 
     RowBox[{"Sum", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"c", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "*", 
        RowBox[{"ChebyshevT", "[", 
         RowBox[{
          RowBox[{"i", "-", "1"}], ",", "y"}], "]"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "get", " ", "the", " ", "coefficients", " ", "from", " ", "the", " ", 
     "values", " ", "of", " ", "g", " ", "at", " ", "the", " ", "Chebyshev", 
     " ", "nodes", " ", "using", " ", "the", " ", "discrete", " ", 
     "Chebyshev", " ", "transform"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"coeff", "[", "g_", "]"}], ":=", 
     RowBox[{"Module", "[", 
      RowBox[{
       RowBox[{"{", "c", "}"}], ",", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{
         RowBox[{"getCoeff", "[", "j_", "]"}], ":=", 
         RowBox[{
          RowBox[{"2", "/", "n"}], "*", 
          RowBox[{"Sum", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"g", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "*", 
             RowBox[{"ChebyshevT", "[", 
              RowBox[{
               RowBox[{"j", "-", "1"}], ",", 
               RowBox[{"x", "[", 
                RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], ",", 
            RowBox[{"{", 
             RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}]}], ";", 
        "\[IndentingNewLine]", "\[IndentingNewLine]", 
        RowBox[{"c", "=", 
         RowBox[{"Table", "[", 
          RowBox[{
           RowBox[{"getCoeff", "[", "i", "]"}], ",", 
           RowBox[{"{", 
            RowBox[{"i", ",", "2", ",", "n"}], "}"}]}], "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"PrependTo", "[", 
         RowBox[{"c", ",", 
          RowBox[{
           RowBox[{"getCoeff", "[", "1", "]"}], "/", "2"}]}], "]"}], ";", 
        "\[IndentingNewLine]", "c"}]}], "\[IndentingNewLine]", "]"}]}], ";"}],
    "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "construct", " ", "g", " ", "from", " ", "its", " ", "Chebyshev", " ", 
     "coefficients"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"expand", "[", "g_", "]"}], ":=", 
     RowBox[{"Module", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"c", ",", "f"}], "}"}], ",", "\[IndentingNewLine]", 
       "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"c", "=", 
         RowBox[{"coeff", "[", "g", "]"}]}], ";", "\[IndentingNewLine]", 
        "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"f", "[", "y_", "]"}], ":=", 
         RowBox[{"gc", "[", 
          RowBox[{"c", ",", "y"}], "]"}]}], ";", "\[IndentingNewLine]", 
        "f"}]}], "\[IndentingNewLine]", "]"}]}], ";"}], "\[IndentingNewLine]",
    "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "construct", " ", "g", " ", "from", " ", "its", " ", "Chebyshev", " ", 
     "coefficients", " ", "and", " ", "evaluate", " ", "it", " ", "at", " ", 
     "the", " ", "Chebyshev", " ", "nodes"}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"H", "[", "g_", "]"}], ":=", 
     RowBox[{"\[Phi]", "[", 
      RowBox[{
       RowBox[{"expand", "[", "g", "]"}], ",", "x"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "perturb", " ", "the", " ", "nth", " ", "component", " ", "of", " ", 
     "vector", " ", "x", " ", "by", " ", "\[Epsilon]"}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"perturb", "[", 
      RowBox[{"x_", ",", "n_", ",", "\[Epsilon]_"}], "]"}], ":=", 
     RowBox[{"Module", "[", 
      RowBox[{
       RowBox[{"{", "y", "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"y", "=", "x"}], ";", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"y", "[", 
          RowBox[{"[", "n", "]"}], "]"}], "=", 
         RowBox[{
          RowBox[{"y", "[", 
           RowBox[{"[", "n", "]"}], "]"}], "+", "\[Epsilon]"}]}], ";", 
        "\[IndentingNewLine]", "y"}]}], "\[IndentingNewLine]", "]"}]}], ";"}],
    "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Epsilon]", "=", 
     RowBox[{"10", "^", 
      RowBox[{"(", 
       RowBox[{"-", "7"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "approximate", " ", "the", " ", "Jacobian", " ", "using", " ", "divided", 
     " ", "differences"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"diff", "[", "g_", "]"}], ":=", 
     RowBox[{"Transpose", "[", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"N", "[", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"H", "[", 
             RowBox[{"perturb", "[", 
              RowBox[{"g", ",", "i", ",", "\[Epsilon]"}], "]"}], "]"}], "-", 
            RowBox[{"H", "[", "g", "]"}]}], ")"}], "/", "\[Epsilon]"}], "]"}],
         ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "initialize", " ", "g", " ", "with", " ", "a", " ", "reasonable", " ", 
     "estimate"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"g", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{"1", "-", 
         RowBox[{"#1", "*", 
          RowBox[{"#2", "^", "2"}]}]}], "&"}], "[", 
       RowBox[{"1.4", ",", "x"}], "]"}], "}"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"it", "=", "5"}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"perform", " ", 
     RowBox[{"Newton", "'"}], "s", " ", "method"}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"k", "=", "1"}], ",", 
     RowBox[{"k", "\[LessEqual]", "it"}], ",", 
     RowBox[{"k", "++"}], ",", "\[IndentingNewLine]", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"A", "=", 
       RowBox[{"diff", "[", 
        RowBox[{"g", "[", 
         RowBox[{"[", "k", "]"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
      "\[IndentingNewLine]", 
      RowBox[{"AppendTo", "[", 
       RowBox[{"g", ",", 
        RowBox[{
         RowBox[{"g", "[", 
          RowBox[{"[", "k", "]"}], "]"}], "-", 
         RowBox[{
          RowBox[{"Inverse", "[", "A", "]"}], ".", 
          RowBox[{"H", "[", 
           RowBox[{"g", "[", 
            RowBox[{"[", "k", "]"}], "]"}], "]"}]}]}]}], "]"}], ";", 
      "\[IndentingNewLine]", "\[IndentingNewLine]", 
      RowBox[{"sup", "=", 
       RowBox[{"Max", "[", 
        RowBox[{"Abs", "[", 
         RowBox[{
          RowBox[{"g", "[", 
           RowBox[{"[", 
            RowBox[{"k", "+", "1"}], "]"}], "]"}], "-", 
          RowBox[{"g", "[", 
           RowBox[{"[", "k", "]"}], "]"}]}], "]"}], "]"}]}], ";", 
      "\[IndentingNewLine]", "\[IndentingNewLine]", 
      RowBox[{"(*", " ", 
       RowBox[{
       "we", " ", "have", " ", "quadratic", " ", "convergence", " ", "to", 
        " ", "the", " ", "universal", " ", "function"}], " ", "*)"}], 
      "\[IndentingNewLine]", 
      RowBox[{"Print", "[", "sup", "]"}], ";"}]}], "\[IndentingNewLine]", 
    "]"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "\[Alpha]", " ", "can", " ", "now", " ", "be", " ", "determined", " ", 
     "from", " ", "g"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{"\[Alpha]", "=", 
    RowBox[{
     RowBox[{"-", "1"}], "/", 
     RowBox[{
      RowBox[{"expand", "[", 
       RowBox[{"g", "[", 
        RowBox[{"[", 
         RowBox[{"it", "+", "1"}], "]"}], "]"}], "]"}], "[", "1", "]"}]}]}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"ListPlot", "[", 
    RowBox[{
     RowBox[{"Transpose", "[", 
      RowBox[{"{", 
       RowBox[{"x", ",", 
        RowBox[{"g", "[", 
         RowBox[{"[", 
          RowBox[{"it", "+", "1"}], "]"}], "]"}]}], "}"}], "]"}], ",", 
     RowBox[{"AxesLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"\"\<x\>\"", ",", "\"\<y\>\""}], "}"}]}], ",", 
     RowBox[{"ImageSize", "\[Rule]", "Large"}]}], "]"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"Max", "[", 
    RowBox[{"Abs", "[", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"gc", "[", 
         RowBox[{
          RowBox[{"coeff", "[", 
           RowBox[{"g", "[", 
            RowBox[{"[", 
             RowBox[{"it", "+", "1"}], "]"}], "]"}], "]"}], ",", "#"}], "]"}],
         "-", 
        RowBox[{"gc", "[", 
         RowBox[{
          RowBox[{"coeff", "[", 
           RowBox[{"g", "[", 
            RowBox[{"[", "it", "]"}], "]"}], "]"}], ",", "#"}], "]"}]}], 
       "&"}], "[", 
      RowBox[{"Range", "[", 
       RowBox[{
        RowBox[{"-", "1"}], ",", "1", ",", "0.001"}], "]"}], "]"}], "]"}], 
    "]"}]}]}]], "Input",
 CellChangeTimes->{{3.79507484313233*^9, 3.7950748553789206`*^9}, {
   3.7950749833517003`*^9, 3.7950750871374397`*^9}, {3.7950751543048778`*^9, 
   3.795075175235882*^9}, {3.7950752349004583`*^9, 3.795075239073605*^9}, {
   3.795075280741591*^9, 3.7950753925975075`*^9}, {3.7950755727252216`*^9, 
   3.795075580430415*^9}, {3.7950756145600576`*^9, 3.7950756146218853`*^9}, {
   3.7950756726694927`*^9, 3.7950756728012304`*^9}, {3.7950760225875883`*^9, 
   3.795076049927848*^9}, {3.7950761359707704`*^9, 3.7950761381698008`*^9}, {
   3.7950763697024994`*^9, 3.7950763815878253`*^9}, {3.7950764869451475`*^9, 
   3.795076540775589*^9}, {3.7950765809535465`*^9, 3.795076583432449*^9}, {
   3.7950766135078797`*^9, 3.7950766460631905`*^9}, {3.795077075273151*^9, 
   3.795077087358836*^9}, {3.7950771532129507`*^9, 3.7950772642225556`*^9}, {
   3.7950773849055185`*^9, 3.795077466834695*^9}, {3.795077968043283*^9, 
   3.795077981764943*^9}, {3.7950780298136735`*^9, 3.7950781295155516`*^9}, {
   3.795078195770855*^9, 3.7950783289358044`*^9}, {3.795078624418812*^9, 
   3.7950786432298927`*^9}, {3.7950819660096893`*^9, 3.795082086317921*^9}, {
   3.7950821206413975`*^9, 3.795082131448152*^9}, {3.7950821692285414`*^9, 
   3.7950821695448885`*^9}, {3.7950822146454487`*^9, 3.795082218662984*^9}, {
   3.7950822866917872`*^9, 3.7950822974769964`*^9}, 3.7950823403254566`*^9, {
   3.7950824113091855`*^9, 3.7950824174755692`*^9}, {3.7950824712901983`*^9, 
   3.795082473672692*^9}, {3.7950825526943283`*^9, 3.795082577622246*^9}, {
   3.7950826183243475`*^9, 3.7950826183623247`*^9}, {3.7950827856703124`*^9, 
   3.7950828658362265`*^9}, {3.7950829100359926`*^9, 3.795082926558465*^9}, {
   3.7950829774398203`*^9, 3.7950829844265685`*^9}, {3.795083150927009*^9, 
   3.7950831522170496`*^9}, {3.7950832241896343`*^9, 3.795083229904158*^9}, {
   3.7950834276136537`*^9, 3.795083431210474*^9}, {3.79508346142162*^9, 
   3.7950834783725386`*^9}, {3.7950835194472675`*^9, 3.795083547901086*^9}, {
   3.7950835932379894`*^9, 3.79508365498869*^9}, {3.7950837126931114`*^9, 
   3.7950837139018555`*^9}, {3.795083776328318*^9, 3.7950837854879336`*^9}, {
   3.7950839456042943`*^9, 3.7950840050683985`*^9}, {3.795084063100556*^9, 
   3.7950840639865894`*^9}, {3.795084099551612*^9, 3.795084113804994*^9}, {
   3.7950841521376586`*^9, 3.7950841988962708`*^9}, {3.7950842313157406`*^9, 
   3.7950842510554314`*^9}, {3.7950843644495544`*^9, 
   3.7950844336027966`*^9}, {3.795151561931676*^9, 3.7951515868152013`*^9}, {
   3.7951517794024434`*^9, 3.795151788608837*^9}, 3.7951518594575033`*^9, {
   3.7951519280082645`*^9, 3.795152007417027*^9}, {3.7951520438336964`*^9, 
   3.7951520441149473`*^9}, {3.7951520787932634`*^9, 
   3.7951521034094706`*^9}, {3.7951521637342405`*^9, 3.795152171741838*^9}, {
   3.7951522141056128`*^9, 3.795152224946638*^9}, {3.7951523202439384`*^9, 
   3.7951523826631093`*^9}, {3.7951524822649364`*^9, 3.795152566839775*^9}, {
   3.7951526027278543`*^9, 3.7951526291801605`*^9}, {3.7951526825176272`*^9, 
   3.7951528644930058`*^9}, {3.7951529527132215`*^9, 
   3.7951529996438236`*^9}, {3.795153053290446*^9, 3.795153057011496*^9}, {
   3.7951531384227505`*^9, 3.7951531552109175`*^9}, {3.795153186442442*^9, 
   3.7951532449540253`*^9}, {3.7951532812330613`*^9, 
   3.7951533141520777`*^9}, {3.7951533570623913`*^9, 
   3.7951533949282265`*^9}, {3.7951534291288056`*^9, 
   3.7951534361280785`*^9}, {3.795153483653082*^9, 3.7951534886865997`*^9}, {
   3.7951535416340933`*^9, 3.795153579233449*^9}, {3.7951536522492733`*^9, 
   3.7951536952044697`*^9}, {3.7951537270114594`*^9, 
   3.7951537275230913`*^9}, {3.7951537716003118`*^9, 
   3.7951538392762823`*^9}, {3.7951538844914365`*^9, 3.795153884764741*^9}, {
   3.7951539187718153`*^9, 3.7951539272681074`*^9}, {3.7951539623703113`*^9, 
   3.7951540028042197`*^9}, {3.7951540910464096`*^9, 
   3.7951540911062517`*^9}, {3.7951541276246443`*^9, 
   3.7951541277253466`*^9}, {3.7951541764820642`*^9, 
   3.7951543697924013`*^9}, {3.7951544091691737`*^9, 3.795154417519845*^9}, {
   3.795154520446731*^9, 3.7951545632732954`*^9}, {3.7951546471650443`*^9, 
   3.795154681184123*^9}, {3.7951547137520785`*^9, 3.7951548810728827`*^9}, 
   3.795154920212308*^9, {3.7951549867693896`*^9, 3.795154991591526*^9}, {
   3.79515503193968*^9, 3.795155055090786*^9}, {3.7951550968013077`*^9, 
   3.7951551404993114`*^9}, {3.7951551874460745`*^9, 
   3.7951551948622384`*^9}, {3.7951553943331065`*^9, 3.795155415903221*^9}, {
   3.7951554867173862`*^9, 3.7951555634677935`*^9}, {3.7951556030021305`*^9, 
   3.7951556076829014`*^9}, {3.795155687302083*^9, 3.7951556901770735`*^9}, 
   3.7951557270092125`*^9, {3.795155921063341*^9, 3.7951560002793517`*^9}, {
   3.795156074879015*^9, 3.795156263352704*^9}, {3.7951564748669243`*^9, 
   3.795156474929647*^9}, {3.7951565203308496`*^9, 3.795156582211229*^9}, {
   3.7951566163473077`*^9, 3.7951566654209223`*^9}, {3.7951567235693145`*^9, 
   3.795156723754304*^9}, 3.795156762630292*^9, {3.7951568179630904`*^9, 
   3.7951568180255833`*^9}, {3.79515686148129*^9, 3.7951568615750227`*^9}, {
   3.7951569171319323`*^9, 3.7951569172256765`*^9}, {3.79515696458521*^9, 
   3.79515696843573*^9}, {3.795157563464004*^9, 3.7951575748934336`*^9}, {
   3.795157643022173*^9, 3.7951576431601276`*^9}, {3.7951576774466343`*^9, 
   3.795157762068041*^9}, {3.7951579618311567`*^9, 3.7951579970167227`*^9}, {
   3.795158053878817*^9, 3.7951581305431385`*^9}, {3.795159031162404*^9, 
   3.795159065678258*^9}, {3.795159474403389*^9, 3.7951595120326147`*^9}, {
   3.7951596610977917`*^9, 3.7951597037006807`*^9}, {3.7951598038663783`*^9, 
   3.7951598039826765`*^9}, {3.795159857953082*^9, 3.7951598581312323`*^9}, {
   3.7951602332841234`*^9, 3.7951602532765875`*^9}, {3.7951602891296487`*^9, 
   3.7951603195836835`*^9}, {3.7951603938603835`*^9, 3.795160436524789*^9}, {
   3.7951604847850733`*^9, 3.7951605485510163`*^9}, {3.7951606238188396`*^9, 
   3.7951606499658723`*^9}, {3.7951607048714905`*^9, 3.795160811603306*^9}, {
   3.7951608884827614`*^9, 3.7951610643187327`*^9}, {3.795161097964632*^9, 
   3.7951611343646*^9}, {3.795161170560156*^9, 3.7951612038035135`*^9}, {
   3.795161251184374*^9, 3.795161313497105*^9}, {3.795162316215767*^9, 
   3.7951624847592525`*^9}, {3.795162564072455*^9, 3.795162593575457*^9}, {
   3.795162639478147*^9, 3.795162679870129*^9}, {3.7951628066886296`*^9, 
   3.795162829794921*^9}, {3.7951629183617134`*^9, 3.7951629869034786`*^9}, {
   3.7951630242664604`*^9, 3.795163044097721*^9}, {3.7951631020351753`*^9, 
   3.7951631393192763`*^9}, {3.7951632019178452`*^9, 
   3.7951632051186657`*^9}, {3.7951634458506007`*^9, 
   3.7951635168905897`*^9}, {3.7951635599213743`*^9, 
   3.7951635784202733`*^9}, {3.795163665819278*^9, 3.7951637498991656`*^9}, {
   3.7951637882100344`*^9, 3.795163789767317*^9}, 3.795163857126602*^9, {
   3.7951871776118526`*^9, 3.795187214682426*^9}, {3.7951872467711787`*^9, 
   3.7951872469554405`*^9}, {3.7951873109384947`*^9, 
   3.7951873112541866`*^9}, {3.7951873565261593`*^9, 
   3.7951874203650837`*^9}, {3.795238060250222*^9, 3.79523806147488*^9}, 
   3.7952380986424694`*^9, {3.795238237621951*^9, 3.7952382694983196`*^9}, {
   3.7952383877217827`*^9, 3.7952384446327944`*^9}, {3.795238480783843*^9, 
   3.7952384809115014`*^9}, {3.7952385227201552`*^9, 
   3.7952385249465327`*^9}, {3.7952386503516774`*^9, 3.795238751449489*^9}, {
   3.795238951666271*^9, 3.7952390364334397`*^9}, {3.795239067006586*^9, 
   3.7952390670755773`*^9}, {3.7952391455146604`*^9, 
   3.7952391789329634`*^9}, {3.7952394929532595`*^9, 3.795239560692937*^9}, {
   3.79523960815443*^9, 3.7952398281319785`*^9}, {3.7952398721883707`*^9, 
   3.795239889958288*^9}, {3.795239962233492*^9, 3.7952399883630958`*^9}, {
   3.7952400313269434`*^9, 3.795240038466095*^9}, {3.7952419726238284`*^9, 
   3.7952419818918524`*^9}, {3.7952420393135853`*^9, 3.795242059880535*^9}, {
   3.7952430089972944`*^9, 3.7952430616854534`*^9}, {3.795243208644087*^9, 
   3.7952432089601707`*^9}, {3.79524324051865*^9, 3.7952433064674163`*^9}, {
   3.7952434383548703`*^9, 3.7952434692393928`*^9}, {3.7952435238023205`*^9, 
   3.79524352745129*^9}, 3.7952437263198395`*^9, {3.795243802916382*^9, 
   3.7952438139286203`*^9}, {3.7952439653660507`*^9, 3.7952440550148*^9}, {
   3.7952441250966215`*^9, 3.795244153698166*^9}, {3.7952443475278473`*^9, 
   3.795244351106698*^9}, {3.7952448160896397`*^9, 3.7952448808415613`*^9}, {
   3.7952509392093143`*^9, 3.7952509581080513`*^9}, {3.7952511789554853`*^9, 
   3.7952512065657196`*^9}, {3.795253963425579*^9, 3.7952540114345093`*^9}, 
   3.795411451355878*^9, {3.7954117446902065`*^9, 3.7954117448477426`*^9}, {
   3.7954148177742553`*^9, 3.7954148178221197`*^9}, {3.7954148700593157`*^9, 
   3.795414870121151*^9}, {3.7954962784909105`*^9, 3.7954963987312517`*^9}, {
   3.795496455519002*^9, 3.795496477000684*^9}, {3.7954965156797915`*^9, 
   3.795496674200368*^9}, {3.7954967088148108`*^9, 3.7954967600721774`*^9}, {
   3.7954968760435934`*^9, 3.795496917844624*^9}, {3.795496953532051*^9, 
   3.7954970714309826`*^9}, {3.795497177193751*^9, 3.795497186610553*^9}, {
   3.7954972296370893`*^9, 3.7954973695992947`*^9}, {3.7954974157284946`*^9, 
   3.795497477145604*^9}, {3.795497515543601*^9, 3.7954975899054832`*^9}, {
   3.7954976266352806`*^9, 3.7954977183469863`*^9}, {3.795498062836174*^9, 
   3.7954980629424725`*^9}, {3.795498098512146*^9, 3.7954982694993443`*^9}, {
   3.795677065795061*^9, 3.7956771311405783`*^9}, {3.795677913985016*^9, 
   3.795678008820324*^9}, {3.7956780398853035`*^9, 3.7956780519739504`*^9}, {
   3.798701957718305*^9, 3.798701957891305*^9}, {3.798702297761305*^9, 
   3.7987023689333053`*^9}},
 CellLabel->"In[29]:=",ExpressionUUID->"67f4eb6e-cb07-41ea-a26c-6b540f6cca2f"],

Cell[CellGroupData[{

Cell[BoxData["0.035037481074667365`"], "Print",
 CellChangeTimes->{3.801035192154355*^9},
 CellLabel->
  "During evaluation of \
In[29]:=",ExpressionUUID->"9a0af35c-a984-4d37-b7bd-bab199dc1e16"],

Cell[BoxData["0.0017315591321610535`"], "Print",
 CellChangeTimes->{3.8010351941065345`*^9},
 CellLabel->
  "During evaluation of \
In[29]:=",ExpressionUUID->"e630289a-91e6-4473-9361-a43a0f6c8bec"],

Cell[BoxData["2.544808761995032`*^-6"], "Print",
 CellChangeTimes->{3.8010351957488556`*^9},
 CellLabel->
  "During evaluation of \
In[29]:=",ExpressionUUID->"5811509e-6f0a-4069-bf54-438db5594d3d"],

Cell[BoxData["6.870315427676132`*^-11"], "Print",
 CellChangeTimes->{3.8010351973155317`*^9},
 CellLabel->
  "During evaluation of \
In[29]:=",ExpressionUUID->"b9d3f30e-81ca-41ea-83cc-27a69ab72820"],

Cell[BoxData["2.5979218776228663`*^-14"], "Print",
 CellChangeTimes->{3.8010351989447775`*^9},
 CellLabel->
  "During evaluation of \
In[29]:=",ExpressionUUID->"b2a20b85-e415-47df-8dae-630607fdb942"]
}, Open  ]],

Cell[BoxData["2.5029078750957545`"], "Output",
 CellChangeTimes->{3.8010351990106015`*^9},
 CellLabel->"Out[42]=",ExpressionUUID->"f1adf1d1-a6c9-4eca-a469-3c4a6a9e7445"],

Cell[BoxData[
 GraphicsBox[{{}, {
    {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.012833333333333334`],
      AbsoluteThickness[1.6], PointBox[CompressedData["
1:eJxTTMoPSmViYGAQAWIQvT5BTW//s/f2LJsaMlcJ39x/Tvp763LZ9/ZLXE3b
Z+tc3e/Wa3Q+e8Jb+5PX376WtTm934F1KdsVj9f2Dssub9tft3y/8HvbJ3ah
L+zfOVWvWsFxxP7dOwY39xNP7J8KKV24JnzDfq7e1xd+Ox7Yr1LuDj33+pF9
8sp1yx7V3rCX9Y9+mvbqhf2Z31OWdj48a+/xi+PU9Ii39p59lkZPRbfYe6b9
5Pq/6T2Mv98Hyoeq3+8EVQ81b78U1DyoffuPQO2Dume/hzDEPVD37t/vDHEv
1D/7ty+G+Afq3/1xUP9Cw2N/LTQ8oOG1/9tGSHgBAGQzssU=
      "]]}, {
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
      0.012833333333333334`], AbsoluteThickness[1.6]}, {}}, {
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
      0.012833333333333334`], AbsoluteThickness[1.6]}, {}}}, {{}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{
    FormBox["\"x\"", TraditionalForm], 
    FormBox["\"y\"", 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]],
  ImageSize->Large,
  Method->{
   "OptimizePlotMarkers" -> True, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{-0.996917333733128, 0.996917333733128}, {-0.39182529738425864`,
    0.9906001295414142}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.801035199121419*^9},
 CellLabel->"Out[43]=",ExpressionUUID->"8a5215eb-b618-4019-a25d-96ab0f2c803e"],

Cell[BoxData["2.6145752229922437`*^-14"], "Output",
 CellChangeTimes->{3.801035199460019*^9},
 CellLabel->"Out[44]=",ExpressionUUID->"31aeec02-fed0-4709-88d9-4c632efbaa59"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"define", " ", "g"}], ",", " ", 
    RowBox[{
     RowBox[{
     "the", " ", "linearized", " ", "period", " ", "doubling", " ", 
      "transformation", " ", "and", " ", "\[IndentingNewLine]", "operator", 
      " ", "the", " ", "vector"}], "-", 
     RowBox[{"valued", " ", "operator", " ", "M"}]}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"\[Gamma]", "[", "y_", "]"}], ":=", 
     RowBox[{"gc", "[", 
      RowBox[{
       RowBox[{"coeff", "[", 
        RowBox[{"g", "[", 
         RowBox[{"[", 
          RowBox[{"it", "+", "1"}], "]"}], "]"}], "]"}], ",", "y"}], "]"}]}], 
    ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"assuming", " ", "\[Alpha]", " ", "is", " ", "constant"}], " ", 
    "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"dF", "[", 
      RowBox[{"h_", ",", "y_"}], "]"}], ":=", 
     RowBox[{
      RowBox[{"-", "\[Alpha]"}], "*", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"\[Gamma]", "'"}], "[", 
          RowBox[{"\[Gamma]", "[", 
           RowBox[{
            RowBox[{"-", "y"}], "/", "\[Alpha]"}], "]"}], "]"}], "*", 
         RowBox[{"h", "[", 
          RowBox[{
           RowBox[{"-", "y"}], "/", "\[Alpha]"}], "]"}]}], "+", 
        RowBox[{"h", "[", 
         RowBox[{"\[Gamma]", "[", 
          RowBox[{
           RowBox[{"-", "y"}], "/", "\[Alpha]"}], "]"}], "]"}]}], ")"}]}]}], 
    ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"dF2", "[", 
      RowBox[{"h_", ",", "y_"}], "]"}], ":=", 
     RowBox[{
      RowBox[{"dF", "[", 
       RowBox[{"h", ",", "y"}], "]"}], "-", 
      RowBox[{"\[Alpha]", "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"\[Gamma]", "'"}], "[", 
           RowBox[{"\[Gamma]", "[", 
            RowBox[{
             RowBox[{"-", "y"}], "/", "\[Alpha]"}], "]"}], "]"}], "*", 
          RowBox[{
           RowBox[{"\[Gamma]", "'"}], "[", 
           RowBox[{
            RowBox[{"-", "y"}], "/", "\[Alpha]"}], "]"}], "*", "y"}], "+", 
         RowBox[{"\[Alpha]", "*", 
          RowBox[{"\[Gamma]", "[", 
           RowBox[{"\[Gamma]", "[", 
            RowBox[{
             RowBox[{"-", "y"}], "/", "\[Alpha]"}], "]"}], "]"}]}]}], ")"}], 
       "*", 
       RowBox[{"h", "[", "1", "]"}]}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"M", "[", "g_", "]"}], ":=", 
     RowBox[{"dF", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"gc", "[", 
         RowBox[{
          RowBox[{"coeff", "[", "g", "]"}], ",", "#"}], "]"}], "&"}], ",", 
       "x"}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "determine", " ", "the", " ", "largest", " ", "eigenvalue", " ", "of", 
     " ", "dF", " ", "by", " ", "the", " ", "power", " ", "method"}], " ", 
    "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"b", "=", 
     RowBox[{"g", "[", 
      RowBox[{"[", 
       RowBox[{"it", "+", "1"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"u", "=", 
     RowBox[{"Nest", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"M", "[", "#", "]"}], "/", 
         RowBox[{"Norm", "[", "#", "]"}]}], "&"}], ",", "b", ",", "30"}], 
      "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Delta]", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"M", "[", "u", "]"}], ".", "u"}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{"u", ".", "u"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"NumberForm", "[", 
    RowBox[{"\[Delta]", ",", "8"}], "]"}], "\[IndentingNewLine]", "\n", 
   RowBox[{"(*", " ", 
    RowBox[{"determine", " ", "the", " ", 
     RowBox[{"(", 
      RowBox[{"m", "+", "1"}], ")"}], "th", " ", "largest", " ", 
     "eigenvalues", " ", "using", " ", 
     RowBox[{"Arnoldi", "'"}], "s", " ", "method"}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"m", "=", 
     RowBox[{"n", "+", "1"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"q", "=", 
     RowBox[{"ConstantArray", "[", 
      RowBox[{"0", ",", 
       RowBox[{"{", 
        RowBox[{"n", ",", 
         RowBox[{"m", "+", "1"}]}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"q", "[", 
      RowBox[{"[", 
       RowBox[{"All", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"Normalize", "[", "b", "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"h", "=", 
     RowBox[{"ConstantArray", "[", 
      RowBox[{"0", ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"m", "+", "1"}], ",", 
         RowBox[{"m", "+", "1"}]}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"j", "=", "1"}], ",", 
      RowBox[{"j", "\[LessEqual]", "m"}], ",", 
      RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"v", "=", 
        RowBox[{"M", "[", 
         RowBox[{"q", "[", 
          RowBox[{"[", 
           RowBox[{"All", ",", "j"}], "]"}], "]"}], "]"}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"i", "=", "1"}], ",", 
         RowBox[{"i", "\[LessEqual]", "j"}], ",", 
         RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           RowBox[{"h", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", 
           RowBox[{"v", ".", 
            RowBox[{"q", "[", 
             RowBox[{"[", 
              RowBox[{"All", ",", "i"}], "]"}], "]"}]}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"v", "=", 
           RowBox[{"v", "-", 
            RowBox[{
             RowBox[{"h", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}], "*", 
             RowBox[{"q", "[", 
              RowBox[{"[", 
               RowBox[{"All", ",", "i"}], "]"}], "]"}]}]}]}], ";"}]}], 
        "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"h", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"j", "+", "1"}], ",", "j"}], "]"}], "]"}], "=", 
        RowBox[{"Norm", "[", "v", "]"}]}], ";", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"q", "[", 
         RowBox[{"[", 
          RowBox[{"All", ",", 
           RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "=", 
        RowBox[{"v", "/", 
         RowBox[{"h", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"j", "+", "1"}], ",", "j"}], "]"}], "]"}]}]}], ";"}]}], 
     "\[IndentingNewLine]", "]"}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"NumberForm", "[", 
    RowBox[{
     RowBox[{"MatrixForm", "[", 
      RowBox[{"N", "[", 
       RowBox[{"Round", "[", 
        RowBox[{
         RowBox[{"Eigenvalues", "[", "h", "]"}], ",", 
         RowBox[{"12", "^", 
          RowBox[{"(", 
           RowBox[{"-", "8"}], ")"}]}]}], "]"}], "]"}], "]"}], ",", "12"}], 
    "]"}], "\n", 
   RowBox[{"(*", 
    RowBox[{"NumberForm", "[", 
     RowBox[{
      RowBox[{"MatrixForm", "[", "q", "]"}], ",", "4"}], "]"}], 
    "*)"}]}]}]], "Input",
 CellChangeTimes->{{3.7952513880475388`*^9, 3.7952514322309895`*^9}, {
   3.7952515020605*^9, 3.7952516460833726`*^9}, {3.7952518229453335`*^9, 
   3.7952519428604364`*^9}, {3.7952520042471123`*^9, 
   3.7952520042627335`*^9}, {3.7952523072510586`*^9, 
   3.7952523577281466`*^9}, {3.7952524875274553`*^9, 
   3.7952525160807095`*^9}, {3.795252609948347*^9, 3.795252700292244*^9}, {
   3.795253175645121*^9, 3.7952533097918406`*^9}, {3.795253358280711*^9, 
   3.7952534381938443`*^9}, {3.795253503749679*^9, 3.7952535109029503`*^9}, {
   3.7952536608850875`*^9, 3.7952536812321873`*^9}, {3.795253724064703*^9, 
   3.7952538029168086`*^9}, {3.7952538578998146`*^9, 3.7952539461328692`*^9}, 
   3.7952541424372845`*^9, {3.795254172503744*^9, 3.7952542396294656`*^9}, {
   3.7952542737292027`*^9, 3.7952542904600906`*^9}, {3.7952543572416477`*^9, 
   3.7952543744237213`*^9}, {3.7952544793067*^9, 3.7952544935203457`*^9}, {
   3.7952545523216157`*^9, 3.7952545769145265`*^9}, {3.7952546222722025`*^9, 
   3.7952546378379145`*^9}, {3.7952546974280214`*^9, 3.79525476695671*^9}, {
   3.7952548242147207`*^9, 3.795254870735849*^9}, {3.795254912709693*^9, 
   3.7952550156940317`*^9}, {3.79525509662376*^9, 3.795255244611929*^9}, {
   3.7952552796397066`*^9, 3.7952553273498254`*^9}, {3.79525548614886*^9, 
   3.7952556180704627`*^9}, {3.7952556712653246`*^9, 
   3.7952558267845287`*^9}, {3.7952561136688976`*^9, 3.795256147306993*^9}, {
   3.7952562005054784`*^9, 3.7952562744289513`*^9}, {3.795256327690587*^9, 
   3.7952563608643303`*^9}, {3.7952565456368437`*^9, 
   3.7952566011644654`*^9}, {3.795256699671376*^9, 3.7952567680622663`*^9}, {
   3.7952567987772007`*^9, 3.7952568086936655`*^9}, {3.795256908422246*^9, 
   3.7952570720436087`*^9}, {3.7952572164124346`*^9, 3.795257221694317*^9}, {
   3.7952572779569416`*^9, 3.7952573450257244`*^9}, {3.795257510965373*^9, 
   3.7952575449079957`*^9}, 3.7954114568073025`*^9, {3.795411504648447*^9, 
   3.7954115366463385`*^9}, {3.795411852039313*^9, 3.79541185230956*^9}, {
   3.7954119321122723`*^9, 3.795411940751202*^9}, {3.795412036225009*^9, 
   3.795412037512561*^9}, {3.7954122375300784`*^9, 3.795412241420683*^9}, {
   3.7954123152862425`*^9, 3.7954123531719785`*^9}, {3.795412385013912*^9, 
   3.7954125005234947`*^9}, {3.7954137680730295`*^9, 
   3.7954141072728004`*^9}, {3.7954141384544573`*^9, 3.795414265723445*^9}, {
   3.7954143060855455`*^9, 3.795414306139427*^9}, {3.7954143416475024`*^9, 
   3.795414346375886*^9}, {3.7954144138754563`*^9, 3.795414413943267*^9}, {
   3.7954144550015326`*^9, 3.7954144930102234`*^9}, {3.795414569275995*^9, 
   3.7954145734179263`*^9}, {3.79541462833334*^9, 3.7954146678108263`*^9}, {
   3.79541470331079*^9, 3.7954147244380627`*^9}, {3.7954965779368043`*^9, 
   3.795496585553521*^9}, {3.7954977239467487`*^9, 3.7954978800370817`*^9}, {
   3.7954979505324893`*^9, 3.795498000191965*^9}, 3.7954980520854297`*^9, {
   3.795498379411517*^9, 3.795498379742875*^9}, {3.7954984306456504`*^9, 
   3.795498540231958*^9}, {3.7954985714334655`*^9, 3.795498608651112*^9}, {
   3.7954986855581717`*^9, 3.795498731905563*^9}, {3.7954987864842606`*^9, 
   3.7954987877327995`*^9}, {3.7954988400651793`*^9, 
   3.7954988720211086`*^9}, {3.795678836333595*^9, 3.7956789505349555`*^9}, {
   3.798702138865305*^9, 3.7987021496163054`*^9}},
 CellLabel->"In[45]:=",ExpressionUUID->"b817af9b-0534-42fd-ac41-07af471fc2f2"],

Cell[BoxData[
 TagBox[
  InterpretationBox["\<\"4.6692016\"\>",
   4.6692016094613695`,
   AutoDelete->True],
  NumberForm[#, 8]& ]], "Output",
 CellChangeTimes->{3.8010352027560863`*^9},
 CellLabel->
  "Out[52]//NumberForm=",ExpressionUUID->"aa0ea510-9a16-43a0-b0ff-\
4b53ff682161"],

Cell[BoxData[
 TagBox[
  TagBox[
   RowBox[{"(", "\[NoBreak]", 
    TagBox[GridBox[{
       {
        InterpretationBox["\<\"4.66920160946\"\>",
         4.669201609456417,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"-2.50290787727\"\>",
         -2.5029078772692688`,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"0.999999993023\"\>",
         0.9999999930229588,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"-0.399535265334\"\>",
         -0.3995352653337132,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"0.159628467069\"\>",
         0.15962846706851447`,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"-0.123652714743\"\>",
         -0.12365271474253639`,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"-0.0637772334383\"\>",
         -0.06377723343832757,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"-0.0573071719779\"\>",
         -0.05730717197785089,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"0.0254810614078\"\>",
         0.025481061407786065`,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"-0.010180600804\"\>",
         -0.010180600803993294`,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"-0.0101275311031\"\>",
         -0.010127531103091421`,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"0.007403847721\"\>",
         0.0074038477209969426`,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"0.00401937807139\"\>",
         0.004019378071386555,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"-0.00187943581673\"\>",
         -0.0018794358167283473`,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"-0.00162513894545\"\>",
         -0.0016251389454494361`,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"-0.000259359877496\"\>",
         -0.00025935987749580857`,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"-0.0000416157249633\"\>",
         -0.00004161572496332495,
         AutoDelete->True]},
       {
        InterpretationBox[
         RowBox[{"\<\"-7.16076993194\"\>", "\[Times]", 
          SuperscriptBox["10", "\<\"-6\"\>"]}],
         -7.160769931936824*^-6,
         AutoDelete->True]},
       {
        InterpretationBox[
         RowBox[{"\<\"-1.21167948507\"\>", "\[Times]", 
          SuperscriptBox["10", "\<\"-6\"\>"]}],
         -1.2116794850727785`*^-6,
         AutoDelete->True]},
       {
        InterpretationBox[
         RowBox[{"\<\"-6.0467690234\"\>", "\[Times]", 
          SuperscriptBox["10", "\<\"-8\"\>"]}],
         -6.046769023395824*^-8,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"0.\"\>",
         0.,
         AutoDelete->True]},
       {
        InterpretationBox["\<\"0.\"\>",
         0.,
         AutoDelete->True]}
      },
      GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}},
      GridBoxSpacings->{"Columns" -> {
          Offset[0.27999999999999997`], {
           Offset[0.5599999999999999]}, 
          Offset[0.27999999999999997`]}, "Rows" -> {
          Offset[0.2], {
           Offset[0.4]}, 
          Offset[0.2]}}],
     Column], "\[NoBreak]", ")"}],
   Function[BoxForm`e$, 
    MatrixForm[BoxForm`e$]]],
  NumberForm[#, 12]& ]], "Output",
 CellChangeTimes->{3.801035205131137*^9},
 CellLabel->
  "Out[58]//NumberForm=",ExpressionUUID->"8da90538-4d39-4ab6-8b42-\
e9a05f8a2d24"]
}, Open  ]]
},
WindowSize->{1680, 968},
WindowMargins->{{-8, Automatic}, {Automatic, -8}},
Magnification->1.5,
FrontEndVersion->"12.0 for Microsoft Windows (64-bit) (April 8, 2019)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 21227, 442, 1987, "Input",ExpressionUUID->"67f4eb6e-cb07-41ea-a26c-6b540f6cca2f"],
Cell[CellGroupData[{
Cell[21832, 468, 194, 4, 34, "Print",ExpressionUUID->"9a0af35c-a984-4d37-b7bd-bab199dc1e16"],
Cell[22029, 474, 197, 4, 34, "Print",ExpressionUUID->"e630289a-91e6-4473-9361-a43a0f6c8bec"],
Cell[22229, 480, 197, 4, 36, "Print",ExpressionUUID->"5811509e-6f0a-4069-bf54-438db5594d3d"],
Cell[22429, 486, 198, 4, 36, "Print",ExpressionUUID->"b9d3f30e-81ca-41ea-83cc-27a69ab72820"],
Cell[22630, 492, 199, 4, 36, "Print",ExpressionUUID->"b2a20b85-e415-47df-8dae-630607fdb942"]
}, Open  ]],
Cell[22844, 499, 169, 2, 49, "Output",ExpressionUUID->"f1adf1d1-a6c9-4eca-a469-3c4a6a9e7445"],
Cell[23016, 503, 2004, 50, 575, "Output",ExpressionUUID->"8a5215eb-b618-4019-a25d-96ab0f2c803e"],
Cell[25023, 555, 172, 2, 49, "Output",ExpressionUUID->"31aeec02-fed0-4709-88d9-4c632efbaa59"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25232, 562, 10890, 266, 1030, "Input",ExpressionUUID->"b817af9b-0534-42fd-ac41-07af471fc2f2"],
Cell[36125, 830, 283, 9, 68, "Output",ExpressionUUID->"aa0ea510-9a16-43a0-b0ff-4b53ff682161"],
Cell[36411, 841, 3569, 115, 653, "Output",ExpressionUUID->"8da90538-4d39-4ab6-8b42-e9a05f8a2d24"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)