elliptical_gaussian_correct.nb

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

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

(* CreatedBy='Mathematica 12.1' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     45971,        936]
NotebookOptionsPosition[     42858,        884]
NotebookOutlinePosition[     43345,        902]
CellTagsIndexPosition[     43302,        899]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"amplitude", "[", 
    RowBox[{"t_", ",", " ", "dd_", ",", " ", "\[Omega]\[Omega]_"}], "]"}], ":=", 
   RowBox[{"Module", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"tt", ",", " ", "d", ",", " ", "g"}], "}"}], ",", "  ", 
     RowBox[{
      RowBox[{"tt", "=", "t"}], ";", " ", 
      RowBox[{"d", "=", "dd"}], ";", " ", 
      RowBox[{"g", "=", "\[Omega]\[Omega]"}], ";", "\[IndentingNewLine]", 
      RowBox[{"ret", "=", "0"}], ";", "\[IndentingNewLine]", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"m", "=", "0"}], ",", " ", 
        RowBox[{"m", "\[LessEqual]", " ", 
         RowBox[{
          RowBox[{"Length", "[", "d", "]"}], " ", "-", " ", "2"}]}], ",", " ", 
        RowBox[{"m", "++"}], ",", " ", 
        RowBox[{"ret", "+=", 
         RowBox[{
          RowBox[{"d", "[", 
           RowBox[{"[", 
            RowBox[{"m", "+", "2"}], "]"}], "]"}], "*", 
          RowBox[{"Cos", "[", 
           RowBox[{
            RowBox[{"m", "*", 
             RowBox[{"d", "[", 
              RowBox[{"[", "1", "]"}], "]"}]}], " ", "+", " ", 
            RowBox[{"m", "*", "g", "*", "tt"}]}], "]"}]}]}]}], "]"}], ";", 
      "\[IndentingNewLine]", "ret"}]}], "]"}]}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"sim", "[", 
   RowBox[{"dd_", ",", " ", "scale_"}], "]"}], ":=", " ", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"d", ",", " ", "sc"}], "}"}], ",", "  ", 
    RowBox[{
     RowBox[{"d", "=", "dd"}], ";", " ", 
     RowBox[{"sc", " ", "=", "scale"}], ";", " ", "\[IndentingNewLine]", 
     RowBox[{"\[Phi]", " ", "=", " ", 
      RowBox[{"d", "[", 
       RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"\[Omega]", " ", "=", " ", 
      RowBox[{"2", "*", "\[Pi]", "*", "18"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"U", " ", "=", " ", 
      RowBox[{"2", "\[Pi]", "*", "1"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"a", " ", "=", " ", ".001"}], ";", "\[IndentingNewLine]", 
     RowBox[{"end", " ", "=", " ", 
      RowBox[{
       RowBox[{"\[Pi]", "/", 
        RowBox[{"(", 
         RowBox[{"a", "*", "2", "*", "U"}], ")"}]}], "*", "sc"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"h", "[", 
       RowBox[{"x_", ",", " ", "e_", ",", " ", "\[Omega]\[Omega]_"}], "]"}], ":=",
       " ", 
      RowBox[{"Module", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"xx", ",", " ", "ee", ",", " ", "g"}], "}"}], ",", " ", 
        RowBox[{
         RowBox[{"xx", "=", "x"}], ";", " ", 
         RowBox[{"ee", "=", "e"}], ";", 
         RowBox[{"g", "=", "\[Omega]\[Omega]"}], ";", "\[IndentingNewLine]", 
         RowBox[{"ret", "=", "0"}], ";", "\[IndentingNewLine]", 
         RowBox[{"For", "[", 
          RowBox[{
           RowBox[{"i", "=", "0"}], ",", " ", 
           RowBox[{"i", "\[LessEqual]", "12"}], ",", " ", 
           RowBox[{"i", "++"}], ",", " ", 
           RowBox[{"ret", " ", "+=", " ", 
            RowBox[{"U", "*", 
             RowBox[{"Cos", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"-", "i"}], "*", "\[Omega]", "*", "xx"}], " ", "-", 
               " ", 
               RowBox[{"\[Pi]", "/", "2"}], " ", "+", " ", 
               RowBox[{"a", "*", 
                RowBox[{"amplitude", "[", 
                 RowBox[{"xx", ",", " ", "ee", ",", " ", "g"}], "]"}]}]}], 
              "]"}]}]}]}], "]"}], ";", "\[IndentingNewLine]", "ret"}]}], 
       "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"s", "=", 
      RowBox[{"NDSolve", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"\[ImaginaryI]", "*", 
            RowBox[{
             RowBox[{"y", "'"}], "[", "x", "]"}]}], "\[Equal]", 
           RowBox[{
            RowBox[{"h", "[", 
             RowBox[{"x", ",", " ", "d", ",", " ", "\[Omega]"}], "]"}], "*", 
            RowBox[{"y", "[", "x", "]"}]}]}], " ", ",", 
          RowBox[{
           RowBox[{"y", "[", "0", "]"}], "\[Equal]", 
           RowBox[{"1", "/", 
            RowBox[{"Sqrt", "[", "2", "]"}]}]}]}], "}"}], ",", "y", ",", 
        RowBox[{"{", 
         RowBox[{"x", ",", "0", ",", "end"}], "}"}], ",", " ", 
        RowBox[{"MaxSteps", "\[Rule]", "Infinity"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"up", " ", "=", " ", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Re", "[", 
         RowBox[{
          RowBox[{"Evaluate", "[", 
           RowBox[{
            RowBox[{"y", "[", "end", "]"}], "/.", "s"}], "]"}], "[", 
          RowBox[{"[", "1", "]"}], "]"}], "]"}], ",", " ", 
        RowBox[{"Im", "[", 
         RowBox[{
          RowBox[{"Evaluate", "[", 
           RowBox[{
            RowBox[{"y", "[", "end", "]"}], "/.", "s"}], "]"}], "[", 
          RowBox[{"[", "1", "]"}], "]"}], "]"}]}], "}"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"up", " ", "=", " ", 
      RowBox[{
       RowBox[{"Evaluate", "[", 
        RowBox[{
         RowBox[{"y", "[", "end", "]"}], "/.", "s"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"s", "=", 
      RowBox[{"NDSolve", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"\[ImaginaryI]", "*", 
            RowBox[{
             RowBox[{"y", "'"}], "[", "x", "]"}]}], "\[Equal]", 
           RowBox[{
            RowBox[{"-", 
             RowBox[{"h", "[", 
              RowBox[{"x", ",", " ", "d", ",", " ", "\[Omega]"}], "]"}]}], 
            "*", 
            RowBox[{"y", "[", "x", "]"}]}]}], " ", ",", 
          RowBox[{
           RowBox[{"y", "[", "0", "]"}], "\[Equal]", 
           RowBox[{"1", "/", 
            RowBox[{"Sqrt", "[", "2", "]"}]}]}]}], "}"}], ",", "y", ",", 
        RowBox[{"{", 
         RowBox[{"x", ",", "0", ",", "end"}], "}"}], ",", " ", 
        RowBox[{"MaxSteps", "\[Rule]", "Infinity"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"down", " ", "=", " ", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Re", "[", 
         RowBox[{
          RowBox[{"Evaluate", "[", 
           RowBox[{
            RowBox[{"y", "[", "end", "]"}], "/.", "s"}], "]"}], "[", 
          RowBox[{"[", "1", "]"}], "]"}], "]"}], ",", " ", 
        RowBox[{"Im", "[", 
         RowBox[{
          RowBox[{"Evaluate", "[", 
           RowBox[{
            RowBox[{"y", "[", "end", "]"}], "/.", "s"}], "]"}], "[", 
          RowBox[{"[", "1", "]"}], "]"}], "]"}]}], "}"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"down", " ", "=", " ", 
      RowBox[{
       RowBox[{"Evaluate", "[", 
        RowBox[{
         RowBox[{"y", "[", "end", "]"}], "/.", "s"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"Re", "[", 
      RowBox[{
       RowBox[{"up", "*", 
        RowBox[{"Conjugate", "[", "down", "]"}]}], " ", "+", " ", 
       RowBox[{"down", "*", 
        RowBox[{"Conjugate", "[", "up", "]"}]}]}], "]"}]}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"s", "=", 
   RowBox[{"Import", "[", 
    RowBox[{
    "\"\<~/repos/single_site_addressability/final_notebooks/elliptical2040.\
csv\>\"", ",", "\"\<String\>\""}], "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"data", "=", 
   RowBox[{"ImportString", "[", 
    RowBox[{
     RowBox[{"StringReplace", "[", 
      RowBox[{"s", ",", 
       RowBox[{"\"\< \>\"", "\[Rule]", "\"\<\>\""}]}], "]"}], ",", 
     "\"\<Data\>\""}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.8536743657896433`*^9, 3.8536745959091463`*^9}, {
   3.853674863443431*^9, 3.853674910443141*^9}, {3.853675117457375*^9, 
   3.85367514401513*^9}, {3.8536751991385098`*^9, 3.853675236597302*^9}, {
   3.853675285775385*^9, 3.8536753780140257`*^9}, {3.853675881570705*^9, 
   3.8536758823138943`*^9}, {3.85367731031667*^9, 3.853677313491145*^9}, {
   3.8536776376572227`*^9, 3.853677649067718*^9}, {3.853677791095517*^9, 
   3.85367779136609*^9}, 3.853678322137643*^9, {3.8536783645148783`*^9, 
   3.853678383960507*^9}, {3.853678711141357*^9, 3.853678711461541*^9}, 
   3.85367914713828*^9, {3.853684900570665*^9, 3.853684922849593*^9}, {
   3.853685027168192*^9, 3.8536850317596483`*^9}, 3.853685074920278*^9, {
   3.853685146102857*^9, 3.853685146494647*^9}, {3.853688452240942*^9, 
   3.8536884803584146`*^9}, {3.853688582478001*^9, 3.853688582884939*^9}, {
   3.853688665454718*^9, 3.8536886662537622`*^9}, {3.853689040618313*^9, 
   3.853689047870369*^9}, {3.8536891489687366`*^9, 3.85368917533892*^9}, {
   3.8536892402007217`*^9, 3.853689240567411*^9}, {3.853689350527638*^9, 
   3.853689354866362*^9}, {3.853689387768202*^9, 3.8536894150002193`*^9}, {
   3.853689496087983*^9, 3.8536895044291058`*^9}, {3.8536897063567*^9, 
   3.85368972683926*^9}, {3.8536897869803658`*^9, 3.8536897925199013`*^9}, {
   3.853689873654457*^9, 3.853689901259185*^9}, {3.853689962915304*^9, 
   3.8536899688658247`*^9}, {3.853690002900214*^9, 3.8536900336260757`*^9}, {
   3.8536901020886517`*^9, 3.8536901326870003`*^9}, {3.8536901824630127`*^9, 
   3.853690278973731*^9}, {3.853690359336279*^9, 3.853690393388207*^9}, {
   3.853690432791032*^9, 3.853690433371882*^9}, {3.853690493925581*^9, 
   3.853690503619574*^9}, {3.853690553124702*^9, 3.853690553314909*^9}, {
   3.853690685082426*^9, 3.853690695765847*^9}, {3.85369074385006*^9, 
   3.853690772607856*^9}, {3.853690892077045*^9, 3.8536908932066317`*^9}, {
   3.853690939240532*^9, 3.853690968782983*^9}, {3.8536910539345827`*^9, 
   3.853691094380601*^9}, {3.8536911793746223`*^9, 3.8536912397569847`*^9}, {
   3.8536912851707373`*^9, 3.85369132063616*^9}, {3.8536913613573723`*^9, 
   3.853691364074049*^9}, {3.853691403147917*^9, 3.853691503906432*^9}, {
   3.85369163952085*^9, 3.853691695860112*^9}, {3.8536917479523563`*^9, 
   3.8536917869674673`*^9}, {3.8536923804095373`*^9, 3.853692468089088*^9}, {
   3.853692580119742*^9, 3.853692646205879*^9}, {3.853692740463381*^9, 
   3.853692783699418*^9}, {3.853693019156251*^9, 3.8536930250258007`*^9}, {
   3.853693073468309*^9, 3.8536930890977716`*^9}, {3.8536931580918016`*^9, 
   3.853693168443248*^9}, {3.853694737846693*^9, 3.8536947431690407`*^9}, {
   3.853695897399675*^9, 3.8536958974869432`*^9}, {3.8536959605827293`*^9, 
   3.853695985652339*^9}, {3.853696032158483*^9, 3.8536960756910973`*^9}, {
   3.853696119489854*^9, 3.853696119850593*^9}, {3.85369615272761*^9, 
   3.853696179506233*^9}, {3.853696334500503*^9, 3.853696350704192*^9}, {
   3.853696403162277*^9, 3.8536964702813883`*^9}, {3.853696601777157*^9, 
   3.853696602029481*^9}, {3.853696650329557*^9, 3.8536966515315228`*^9}, {
   3.853696809438369*^9, 3.853696809835968*^9}, {3.853696881462143*^9, 
   3.853696918778885*^9}, {3.853696954269558*^9, 3.853696956667102*^9}, {
   3.853697095148136*^9, 3.85369712355396*^9}, {3.8536972231476583`*^9, 
   3.853697226055709*^9}, 3.853697262587307*^9, {3.85369733708943*^9, 
   3.8536973742628107`*^9}, {3.853697640766515*^9, 3.8536977007575893`*^9}, {
   3.8536978268771067`*^9, 3.8536978271778727`*^9}, {3.853698668659555*^9, 
   3.8536986779291162`*^9}, {3.853699235231475*^9, 3.853699279611374*^9}, {
   3.853699402660597*^9, 3.853699406194165*^9}, {3.853699645123085*^9, 
   3.853699645351893*^9}, {3.853699781200683*^9, 3.853699818486618*^9}, {
   3.853699887496151*^9, 3.853699887629603*^9}, {3.85369996843965*^9, 
   3.853699968586257*^9}, {3.853700022102954*^9, 3.853700034667981*^9}, {
   3.853702355217246*^9, 3.853702376117786*^9}, {3.85370245609006*^9, 
   3.8537024569493933`*^9}, {3.853702490170854*^9, 3.853702490637806*^9}, {
   3.853702547031966*^9, 3.8537025863875723`*^9}, {3.853702700295763*^9, 
   3.85370270757061*^9}, {3.853702986051464*^9, 3.853702989240418*^9}, {
   3.853703032472535*^9, 3.853703033684074*^9}, {3.853703073103657*^9, 
   3.853703127014783*^9}, {3.853703191313239*^9, 3.853703212407836*^9}, {
   3.853703261881538*^9, 3.853703263821433*^9}, {3.853703322957757*^9, 
   3.8537033458717127`*^9}, {3.853703436343142*^9, 3.853703484503354*^9}, {
   3.853708943278688*^9, 3.853709041655178*^9}, {3.853709081494081*^9, 
   3.853709087996724*^9}, {3.8537091333450212`*^9, 3.853709155340178*^9}, {
   3.853709264799343*^9, 3.8537093952662163`*^9}, {3.853709457714802*^9, 
   3.8537094929771137`*^9}, {3.8537096237931423`*^9, 3.853709769567865*^9}, {
   3.853709861406527*^9, 3.853709864741437*^9}, {3.853710014517325*^9, 
   3.853710033061511*^9}, {3.853710091089306*^9, 3.8537101149723053`*^9}, {
   3.853710274612808*^9, 3.853710274889159*^9}, {3.853710488896709*^9, 
   3.85371048949736*^9}, {3.853710665742984*^9, 3.853710694045446*^9}, {
   3.853710725782402*^9, 3.8537107265571632`*^9}, {3.853711015068223*^9, 
   3.8537110152525797`*^9}, {3.853711160650436*^9, 3.853711160888929*^9}, {
   3.853711379255548*^9, 3.853711387414712*^9}, {3.8537114536057034`*^9, 
   3.8537114927291517`*^9}, {3.8537115259667597`*^9, 
   3.8537116637229137`*^9}, {3.853711785524444*^9, 3.8537117858028507`*^9}, {
   3.853711836532168*^9, 3.85371187613792*^9}, {3.8537119091945543`*^9, 
   3.853712028488655*^9}, {3.853713726720662*^9, 3.8537138302878313`*^9}, 
   3.85371456856003*^9, {3.853764659189323*^9, 3.853764665429399*^9}, {
   3.853764725251747*^9, 3.8537647259137583`*^9}, {3.853764784138938*^9, 
   3.85376479036432*^9}, {3.853764820916912*^9, 3.85376483967481*^9}, {
   3.853764879317566*^9, 3.8537649139276133`*^9}, {3.853764966385396*^9, 
   3.8537649665673227`*^9}, {3.853765054050329*^9, 3.853765102986249*^9}, {
   3.853765233214212*^9, 3.85376529044595*^9}, {3.853771465633646*^9, 
   3.853771471447762*^9}, {3.8537731187543716`*^9, 3.853773120271694*^9}, {
   3.853773413662854*^9, 3.853773464502344*^9}, {3.853773528622183*^9, 
   3.853773680426334*^9}, {3.853773711373385*^9, 3.853773731154826*^9}, {
   3.853773772995131*^9, 3.853773990490493*^9}, {3.8537740963284817`*^9, 
   3.85377412555855*^9}, {3.853774166376012*^9, 3.853774166573874*^9}, {
   3.8537743501422777`*^9, 3.853774401971589*^9}, {3.853779469731584*^9, 
   3.8537795802743387`*^9}, {3.853779659992814*^9, 3.8537798078312397`*^9}, {
   3.8537798964006166`*^9, 3.853779947269903*^9}, {3.8537799804958763`*^9, 
   3.853779980893667*^9}, {3.853780139208432*^9, 3.853780184888959*^9}, {
   3.8546343094846077`*^9, 3.8546343378393784`*^9}, {3.854634988285323*^9, 
   3.854634997284748*^9}, {3.854635072241373*^9, 3.854635075192412*^9}, {
   3.854635251101366*^9, 3.854635251482029*^9}, {3.854643634567586*^9, 
   3.854643634963047*^9}, {3.85464429733599*^9, 3.854644297805299*^9}, {
   3.854644336117484*^9, 3.854644371724382*^9}, {3.854644849998622*^9, 
   3.854644870191777*^9}, {3.854645030185088*^9, 3.854645040999422*^9}, {
   3.8546451690610743`*^9, 3.854645172748939*^9}, {3.854645388652811*^9, 
   3.854645436829324*^9}, {3.854645903815405*^9, 3.85464590459067*^9}, 
   3.854646350643763*^9, 3.8546463891233673`*^9, {3.8546464827253923`*^9, 
   3.8546464835609837`*^9}, {3.854646513650613*^9, 3.854646521649877*^9}, {
   3.854646597641128*^9, 3.854646605603918*^9}, {3.854646670256452*^9, 
   3.854646693974844*^9}, {3.854646779319399*^9, 3.854646790941516*^9}, {
   3.854646850534377*^9, 3.854646867947379*^9}, 3.854646991251322*^9, 
   3.854647051538505*^9, 3.854647490014264*^9, 3.8546476997675343`*^9, 
   3.8546477336726723`*^9, 3.854647775979744*^9, {3.854651906288433*^9, 
   3.854651916984133*^9}, {3.8546519599851303`*^9, 3.854651960370017*^9}, {
   3.854652006770667*^9, 3.854652006905072*^9}, {3.854652064377034*^9, 
   3.85465207821658*^9}, 3.854652167456909*^9, {3.8546527054824533`*^9, 
   3.8546527793049994`*^9}, {3.854652909175799*^9, 3.854652911248272*^9}, {
   3.854653042495326*^9, 3.854653042606783*^9}, {3.854653194537889*^9, 
   3.85465322398234*^9}, {3.854653259901227*^9, 3.854653281341522*^9}, {
   3.8546533157234573`*^9, 3.854653319347863*^9}, 3.854653392412509*^9, {
   3.85465343219446*^9, 3.854653593585658*^9}, {3.854653634520596*^9, 
   3.85465370937074*^9}, {3.854653747216405*^9, 3.854653753743925*^9}, {
   3.85465383028835*^9, 3.8546538375746317`*^9}, {3.854653917926985*^9, 
   3.854653935095462*^9}, {3.854654141893053*^9, 3.854654214419475*^9}, {
   3.854654245508334*^9, 3.854654254178976*^9}, {3.8546542849640293`*^9, 
   3.854654292338689*^9}},
 CellLabel->"In[66]:=",ExpressionUUID->"e1f875ed-0c62-4f28-b64a-da4f248d4070"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"DiscretePlot", "[", 
   RowBox[{
    RowBox[{"sim", "[", 
     RowBox[{
      RowBox[{"data", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ",", " ", "k"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"k", ",", " ", "0", ",", " ", "1", ",", " ", ".1"}], "}"}]}], 
   "]"}], " "}]], "Input",
 CellChangeTimes->{{3.8536781476071*^9, 3.853678151684187*^9}, 
   3.853699771852398*^9, {3.8537034462881317`*^9, 3.853703459863743*^9}, {
   3.853703490808908*^9, 3.853703511311789*^9}, {3.853707917758478*^9, 
   3.853707954076024*^9}, {3.854635009804063*^9, 3.854635064635055*^9}, {
   3.854635100980604*^9, 3.8546351511781187`*^9}, {3.854636632549409*^9, 
   3.8546366430116*^9}, 3.854644178881977*^9, {3.854644303161846*^9, 
   3.854644321032157*^9}, {3.85464474859794*^9, 3.8546447487873983`*^9}, 
   3.85464667712156*^9, {3.85464687124504*^9, 3.854646872157959*^9}, {
   3.8546469450791397`*^9, 3.8546469453493834`*^9}, {3.854652529078702*^9, 
   3.854652530639539*^9}, {3.854652800214374*^9, 3.854652801388836*^9}, 
   3.8546532951213427`*^9, {3.854654064265699*^9, 3.854654065352481*^9}, {
   3.8546541202501917`*^9, 3.854654121077506*^9}, {3.854654206945467*^9, 
   3.854654207278981*^9}},
 CellLabel->"In[70]:=",ExpressionUUID->"c1dd616e-4838-483c-95ee-772fe1cdbbf3"],

Cell[BoxData[
 GraphicsBox[{{
    {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.012833333333333334`],
      AbsoluteThickness[1.6], 
     TagBox[{
       {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
        0.012833333333333334`], AbsoluteThickness[1.6], Opacity[0.2], 
        LineBox[{{{0.6000000000000001, -0.056164141920211574`}, {
           0.6000000000000001, 0}}, {{
           0.7000000000000001, -0.32132731379006735`}, {
           0.7000000000000001, 0}}, {{0.8, -0.5629452443761653}, {0.8, 0}}, {{
           0.9, -0.7633818464665422}, {0.9, 0}}, {{1., -0.9079530342524365}, {
           1., 0}}},
         VertexColors->None]}, 
       {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
        0.012833333333333334`], AbsoluteThickness[1.6], Opacity[0.2], 
        LineBox[{{{0., 0.9999999999999998}, {0., 0}}, {{0.1, 
           0.9633469064178873}, {0.1, 0}}, {{0.2, 0.8563024623805333}, {
           0.2, 0}}, {{0.30000000000000004`, 0.6867057569711399}, {
           0.30000000000000004`, 0}}, {{0.4, 0.46696746678932316`}, {
           0.4, 0}}, {{0.5, 0.21316224125045957`}, {0.5, 0}}},
         VertexColors->None]}},
      Annotation[#, "Charting`Private`Tag$17504#1"]& ]}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
      0.012833333333333334`], AbsoluteThickness[1.6], {}, 
      PointBox[{{0., 0.9999999999999998}, {0.1, 0.9633469064178873}, {0.2, 
       0.8563024623805333}, {0.30000000000000004`, 0.6867057569711399}, {0.4, 
       0.46696746678932316`}, {0.5, 0.21316224125045957`}, {
       0.6000000000000001, -0.056164141920211574`}, {
       0.7000000000000001, -0.32132731379006735`}, {
       0.8, -0.5629452443761653}, {0.9, -0.7633818464665422}, {
       1., -0.9079530342524365}}], {}},
     Annotation[#, "Charting`Private`Tag$17504#1"]& ]}, {}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{
   "MessagesHead" -> DiscretePlot, "AxisPadding" -> Scaled[0.02], 
    "DefaultBoundaryStyle" -> Automatic, 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "DefaultPlotStyle" -> {
      Directive[
       RGBColor[0.368417, 0.506779, 0.709798], 
       AbsoluteThickness[1.6]], 
      Directive[
       RGBColor[0.880722, 0.611041, 0.142051], 
       AbsoluteThickness[1.6]], 
      Directive[
       RGBColor[0.560181, 0.691569, 0.194885], 
       AbsoluteThickness[1.6]], 
      Directive[
       RGBColor[0.922526, 0.385626, 0.209179], 
       AbsoluteThickness[1.6]], 
      Directive[
       RGBColor[0.528488, 0.470624, 0.701351], 
       AbsoluteThickness[1.6]], 
      Directive[
       RGBColor[0.772079, 0.431554, 0.102387], 
       AbsoluteThickness[1.6]], 
      Directive[
       RGBColor[0.363898, 0.618501, 0.782349], 
       AbsoluteThickness[1.6]], 
      Directive[
       RGBColor[1, 0.75, 0], 
       AbsoluteThickness[1.6]], 
      Directive[
       RGBColor[0.647624, 0.37816, 0.614037], 
       AbsoluteThickness[1.6]], 
      Directive[
       RGBColor[0.571589, 0.586483, 0.], 
       AbsoluteThickness[1.6]], 
      Directive[
       RGBColor[0.915, 0.3325, 0.2125], 
       AbsoluteThickness[1.6]], 
      Directive[
       RGBColor[0.40082222609352647`, 0.5220066643438841, 0.85], 
       AbsoluteThickness[1.6]], 
      Directive[
       RGBColor[0.9728288904374106, 0.621644452187053, 0.07336199581899142], 
       AbsoluteThickness[1.6]], 
      Directive[
       RGBColor[0.736782672705901, 0.358, 0.5030266573755369], 
       AbsoluteThickness[1.6]], 
      Directive[
       RGBColor[0.28026441037696703`, 0.715, 0.4292089322474965], 
       AbsoluteThickness[1.6]]}, "DomainPadding" -> Scaled[0.02], 
    "PointSizeFunction" -> "SmallPointSize", "RangePadding" -> Scaled[0.05]},
  PlotRange->{{0., 1.}, {-0.9079530342524365, 0.9999999999999998}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.8546350440756693`*^9, 3.854635154875169*^9}, {
   3.854635414080759*^9, 3.854635427643181*^9}, 3.8546355021808987`*^9, 
   3.854636673276359*^9, 3.8546436399520283`*^9, 3.854644210345174*^9, {
   3.854644306777555*^9, 3.8546443458652763`*^9}, 3.854644788704174*^9, 
   3.854644970163806*^9, 3.854646540174344*^9, {3.8546466739355917`*^9, 
   3.854646743175687*^9}, 3.854646784104439*^9, {3.8546468149325333`*^9, 
   3.8546468779653397`*^9}, 3.854646971293462*^9, 3.854647187943108*^9, 
   3.8546475188827467`*^9, 3.854647728367765*^9, 3.854647763094267*^9, 
   3.854651921434339*^9, 3.854651987763356*^9, {3.854652060635347*^9, 
   3.8546521047059183`*^9}, 3.854652139116807*^9, 3.8546524145698757`*^9, 
   3.854652547782436*^9, 3.854652577809824*^9, 3.854652616890245*^9, 
   3.854652781165638*^9, 3.854652862445115*^9, 3.8546530714528313`*^9, 
   3.8546531831101522`*^9, 3.8546532368717823`*^9, {3.854653295440649*^9, 
   3.854653323785205*^9}, 3.854653412194949*^9, {3.854653469552642*^9, 
   3.854653491342918*^9}, {3.8546536023163424`*^9, 3.854653616145735*^9}, {
   3.854653692006339*^9, 3.854653718126913*^9}, {3.8546537556169453`*^9, 
   3.854653765419855*^9}, {3.8546538129180326`*^9, 3.854653842309174*^9}, 
   3.854654066021278*^9, {3.854654121752162*^9, 3.8546541496988497`*^9}, {
   3.854654199576558*^9, 3.854654255314829*^9}, 3.854654286427655*^9, 
   3.8546543211863937`*^9},
 CellLabel->"Out[70]=",ExpressionUUID->"ce7c9e20-61d4-485b-9acf-19df45f99f8d"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"amplitude", "[", 
   RowBox[{"q", ",", " ", 
    RowBox[{"{", 
     RowBox[{"1", ",", " ", "2"}], "}"}], ",", " ", "q"}], "]"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"sim", "[", 
    RowBox[{
     RowBox[{"data", "[", 
      RowBox[{"[", "1", "]"}], "]"}], ",", " ", "1"}], "]"}], "*)"}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.853707931488508*^9, 3.8537079408277197`*^9}, {
   3.854635204674899*^9, 3.854635205057366*^9}, {3.8546353954724627`*^9, 
   3.854635395879298*^9}, {3.854635435048851*^9, 3.8546354351244802`*^9}, {
   3.8546354906854687`*^9, 3.8546354909265013`*^9}, {3.854643989347321*^9, 
   3.85464399729051*^9}, 3.854645025412787*^9, {3.8546528186764402`*^9, 
   3.8546528211402473`*^9}, {3.85465310139575*^9, 3.854653107184914*^9}, {
   3.854653230700836*^9, 3.854653246520028*^9}, {3.85465332939217*^9, 
   3.854653347638596*^9}, {3.854653492742963*^9, 3.854653497469236*^9}, {
   3.8546536206691637`*^9, 3.8546536210681973`*^9}, {3.854653695269279*^9, 
   3.8546536954273977`*^9}, {3.854653886543046*^9, 3.854653891138061*^9}},
 CellLabel->"In[32]:=",ExpressionUUID->"046e7e43-7d67-4e52-a480-5506f5e9ab54"],

Cell[BoxData["2"], "Output",
 CellChangeTimes->{{3.854653867389082*^9, 3.854653892108285*^9}, 
   3.854653927645269*^9, {3.854654030217084*^9, 3.8546540456313467`*^9}},
 CellLabel->"Out[32]=",ExpressionUUID->"75e85079-2191-4263-84b8-56957772077e"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"sim", "[", 
  RowBox[{
   RowBox[{"data", "[", 
    RowBox[{"[", "1", "]"}], "]"}], ",", " ", "1"}], "]"}]], "Input",
 CellChangeTimes->{{3.854653820882175*^9, 3.854653820978272*^9}},
 CellLabel->"In[12]:=",ExpressionUUID->"030bf7d1-d5cc-4700-8c6f-56623baef820"],

Cell[BoxData[
 TemplateBox[{
  "NDSolve", "underdet", 
   "\"There are more dependent variables, \\!\\(\\*RowBox[{\\\"{\\\", \
RowBox[{RowBox[{\\\"h\\\", \\\"[\\\", RowBox[{\\\"x\\\", \\\",\\\", \
RowBox[{\\\"{\\\", RowBox[{\\\"0.`\\\", \\\",\\\", \
\\\"0.9998838246301789`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \
\\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\
\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\"}], \\\"}\\\"}]}], \\\"]\\\"}], \\\",\\\", \
RowBox[{\\\"y\\\", \\\"[\\\", \\\"x\\\", \\\"]\\\"}]}], \\\"}\\\"}]\\), than \
equations, so the system is underdetermined.\"", 2, 12, 1, 
   32214179777300378726, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.854653821405376*^9},
 CellLabel->
  "During evaluation of \
In[12]:=",ExpressionUUID->"5dde1a46-0109-4fe9-b8a6-6dda15494276"],

Cell[BoxData[
 TemplateBox[{
  "ReplaceAll", "reps", 
   "\"\\!\\(\\*RowBox[{\\\"{\\\", RowBox[{\\\"NDSolve\\\", \\\"[\\\", \
RowBox[{RowBox[{\\\"{\\\", RowBox[{RowBox[{RowBox[{\\\"\[ImaginaryI]\\\", \
\\\" \\\", RowBox[{SuperscriptBox[\\\"y\\\", \\\"\[Prime]\\\", \
Rule[MultilineFunction, None]], \\\"[\\\", \\\"x\\\", \\\"]\\\"}]}], \\\"\
\[Equal]\\\", RowBox[{RowBox[{\\\"h\\\", \\\"[\\\", RowBox[{\\\"x\\\", \
\\\",\\\", RowBox[{\\\"{\\\", RowBox[{\\\"0.`\\\", \\\",\\\", \
\\\"0.9998838246301789`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \
\\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\
\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\"}], \\\"}\\\"}]}], \\\"]\\\"}], \\\" \\\", \
RowBox[{\\\"y\\\", \\\"[\\\", \\\"x\\\", \\\"]\\\"}]}]}], \\\",\\\", \
RowBox[{RowBox[{\\\"y\\\", \\\"[\\\", \\\"0\\\", \\\"]\\\"}], \
\\\"\[Equal]\\\", FractionBox[\\\"1\\\", SqrtBox[\\\"2\\\"]]}]}], \
\\\"}\\\"}], \\\",\\\", \\\"y\\\", \\\",\\\", RowBox[{\\\"{\\\", \
RowBox[{\\\"x\\\", \\\",\\\", \\\"0\\\", \\\",\\\", \\\"25.`\\\"}], \
\\\"}\\\"}], \\\",\\\", RowBox[{\\\"MaxSteps\\\", \\\"\[Rule]\\\", \\\"\
\[Infinity]\\\"}]}], \\\"]\\\"}], \\\"}\\\"}]\\) is neither a list of \
replacement rules nor a valid dispatch table, and so cannot be used for \
replacing.\"", 2, 12, 2, 32214179777300378726, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.854653821424201*^9},
 CellLabel->
  "During evaluation of \
In[12]:=",ExpressionUUID->"c5d7ea00-89d8-43ab-acf2-9f3a183b62ab"],

Cell[BoxData[
 TemplateBox[{
  "ReplaceAll", "reps", 
   "\"\\!\\(\\*RowBox[{\\\"{\\\", RowBox[{\\\"NDSolve\\\", \\\"[\\\", \
RowBox[{RowBox[{\\\"{\\\", RowBox[{RowBox[{RowBox[{\\\"\[ImaginaryI]\\\", \
\\\" \\\", RowBox[{SuperscriptBox[\\\"y\\\", \\\"\[Prime]\\\", \
Rule[MultilineFunction, None]], \\\"[\\\", \\\"x\\\", \\\"]\\\"}]}], \\\"\
\[Equal]\\\", RowBox[{RowBox[{\\\"h\\\", \\\"[\\\", RowBox[{\\\"x\\\", \
\\\",\\\", RowBox[{\\\"{\\\", RowBox[{\\\"0.`\\\", \\\",\\\", \
\\\"0.9998838246301789`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \
\\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\
\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\"}], \\\"}\\\"}]}], \\\"]\\\"}], \\\" \\\", \
RowBox[{\\\"y\\\", \\\"[\\\", \\\"x\\\", \\\"]\\\"}]}]}], \\\",\\\", \
RowBox[{RowBox[{\\\"y\\\", \\\"[\\\", \\\"0\\\", \\\"]\\\"}], \
\\\"\[Equal]\\\", FractionBox[\\\"1\\\", SqrtBox[\\\"2\\\"]]}]}], \
\\\"}\\\"}], \\\",\\\", \\\"y\\\", \\\",\\\", RowBox[{\\\"{\\\", \
RowBox[{\\\"x\\\", \\\",\\\", \\\"0\\\", \\\",\\\", \\\"25.`\\\"}], \
\\\"}\\\"}], \\\",\\\", RowBox[{\\\"MaxSteps\\\", \\\"\[Rule]\\\", \\\"\
\[Infinity]\\\"}]}], \\\"]\\\"}], \\\"}\\\"}]\\) is neither a list of \
replacement rules nor a valid dispatch table, and so cannot be used for \
replacing.\"", 2, 12, 3, 32214179777300378726, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.854653821431278*^9},
 CellLabel->
  "During evaluation of \
In[12]:=",ExpressionUUID->"6810c86a-ccc6-4a9e-93d6-3f66c9f41536"],

Cell[BoxData[
 TemplateBox[{
  "NDSolve", "underdet", 
   "\"There are more dependent variables, \\!\\(\\*RowBox[{\\\"{\\\", \
RowBox[{RowBox[{\\\"h\\\", \\\"[\\\", RowBox[{\\\"x\\\", \\\",\\\", \
RowBox[{\\\"{\\\", RowBox[{\\\"0.`\\\", \\\",\\\", \
\\\"0.9998838246301789`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \
\\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\
\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \
\\\"0.`\\\", \\\",\\\", \\\"0.`\\\"}], \\\"}\\\"}]}], \\\"]\\\"}], \\\",\\\", \
RowBox[{\\\"y\\\", \\\"[\\\", \\\"x\\\", \\\"]\\\"}]}], \\\"}\\\"}]\\), than \
equations, so the system is underdetermined.\"", 2, 12, 4, 
   32214179777300378726, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.854653821436235*^9},
 CellLabel->
  "During evaluation of \
In[12]:=",ExpressionUUID->"4ec2a5f9-a0b7-4c40-8498-d490a8bb3d63"],

Cell[BoxData[
 TemplateBox[{
  "ReplaceAll", "reps", 
   "\"\\!\\(\\*RowBox[{\\\"{\\\", RowBox[{\\\"NDSolve\\\", \\\"[\\\", \
RowBox[{RowBox[{\\\"{\\\", RowBox[{RowBox[{RowBox[{\\\"\[ImaginaryI]\\\", \
\\\" \\\", RowBox[{SuperscriptBox[\\\"y\\\", \\\"\[Prime]\\\", \
Rule[MultilineFunction, None]], \\\"[\\\", \\\"x\\\", \\\"]\\\"}]}], \\\"\
\[Equal]\\\", RowBox[{RowBox[{\\\"-\\\", RowBox[{\\\"h\\\", \\\"[\\\", \
RowBox[{\\\"x\\\", \\\",\\\", RowBox[{\\\"{\\\", RowBox[{\\\"0.`\\\", \\\",\\\
\", \\\"0.9998838246301789`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\
\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \
\\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\", \\\",\
\\\", \\\"0.`\\\", \\\",\\\", \\\"0.`\\\"}], \\\"}\\\"}]}], \\\"]\\\"}]}], \\\
\" \\\", RowBox[{\\\"y\\\", \\\"[\\\", \\\"x\\\", \\\"]\\\"}]}]}], \\\",\\\", \
RowBox[{RowBox[{\\\"y\\\", \\\"[\\\", \\\"0\\\", \\\"]\\\"}], \
\\\"\[Equal]\\\", FractionBox[\\\"1\\\", SqrtBox[\\\"2\\\"]]}]}], \
\\\"}\\\"}], \\\",\\\", \\\"y\\\", \\\",\\\", RowBox[{\\\"{\\\", \
RowBox[{\\\"x\\\", \\\",\\\", \\\"0\\\", \\\",\\\", \\\"25.`\\\"}], \
\\\"}\\\"}], \\\",\\\", RowBox[{\\\"MaxSteps\\\", \\\"\[Rule]\\\", \\\"\
\[Infinity]\\\"}]}], \\\"]\\\"}], \\\"}\\\"}]\\) is neither a list of \
replacement rules nor a valid dispatch table, and so cannot be used for \
replacing.\"", 2, 12, 5, 32214179777300378726, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.854653821440637*^9},
 CellLabel->
  "During evaluation of \
In[12]:=",ExpressionUUID->"9072c424-c21c-497a-bdd6-1c426d896f3d"],

Cell[BoxData[
 TemplateBox[{
  "General", "stop", 
   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"ReplaceAll\\\", \
\\\"::\\\", \\\"reps\\\"}], \\\"MessageName\\\"]\\) will be suppressed during \
this calculation.\"", 2, 12, 6, 32214179777300378726, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.854653821445375*^9},
 CellLabel->
  "During evaluation of \
In[12]:=",ExpressionUUID->"9dc58a0f-8345-46e7-8ae5-e5e293563f68"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"2", " ", 
    SuperscriptBox[
     RowBox[{"Re", "[", 
      RowBox[{"y", "[", "25.`", "]"}], "]"}], "2"]}], ",", 
   RowBox[{"2", " ", 
    SuperscriptBox[
     RowBox[{"Im", "[", 
      RowBox[{"y", "[", "25.`", "]"}], "]"}], "2"]}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.853707941064587*^9, 3.854635205303454*^9, 3.854635253704734*^9, {
   3.854635392319038*^9, 3.854635435442675*^9}, {3.8546354913209953`*^9, 
   3.854635495128039*^9}, 3.8546439047329273`*^9, {3.854643991726913*^9, 
   3.854643998033065*^9}, 3.854645127403673*^9, 3.85464525697789*^9, 
   3.85464551797893*^9, 3.854646052299553*^9, 3.85464638600403*^9, 
   3.854652821573058*^9, 3.854653107531478*^9, 3.854653186885393*^9, 
   3.854653236907749*^9, {3.8546533273839903`*^9, 3.8546533478678226`*^9}, {
   3.854653395905274*^9, 3.854653414512145*^9}, {3.854653493563828*^9, 
   3.8546534978174257`*^9}, {3.854653618231825*^9, 3.8546536214391413`*^9}, {
   3.85465369565034*^9, 3.854653725096239*^9}, 3.854653767355098*^9, 
   3.8546538214501333`*^9},
 CellLabel->"Out[12]=",ExpressionUUID->"9d24c5f3-86f1-4a4d-9a28-ba83a3a16b7f"]
}, Open  ]],

Cell[BoxData["y"], "Input",
 CellChangeTimes->{{3.8546532482474337`*^9, 3.854653252775475*^9}, {
  3.8546533403755608`*^9, 3.854653342694601*^9}},
 CellLabel->"In[17]:=",ExpressionUUID->"0f8480c7-0ad3-4df5-8536-2d1912be031e"],

Cell[BoxData["y"], "Output",
 CellChangeTimes->{{3.854653249799222*^9, 3.854653253004085*^9}, {
  3.854653340529853*^9, 3.854653343002096*^9}},
 CellLabel->"Out[17]=",ExpressionUUID->"8ad523e0-fc6d-4367-9cb3-4bbd6111bdc9"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.8546531137622223`*^9, 
  3.8546531143364983`*^9}},ExpressionUUID->"4163ba75-88b1-461b-9e8a-\
35864418bfe5"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"sim", "[", 
  RowBox[{"w", ",", " ", "q"}], "]"}]], "Input",
 CellChangeTimes->{{3.854653359068454*^9, 3.8546533633582783`*^9}},
 CellLabel->"In[19]:=",ExpressionUUID->"569104bd-3e2c-4653-834b-18b13f1d6d12"],

Cell[BoxData[
 TemplateBox[{
  "Part", "partd", 
   "\"Part specification \\!\\(\\*RowBox[{\\\"w\\\", \\\"\[LeftDoubleBracket]\
\\\", \\\"1\\\", \\\"\[RightDoubleBracket]\\\"}]\\) is longer than depth of \
object.\"", 2, 19, 1, 32214176670846220627, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.854653363718771*^9},
 CellLabel->
  "During evaluation of \
In[19]:=",ExpressionUUID->"0574253a-9b71-4337-ad0c-494e482f2f68"],

Cell[BoxData[
 TemplateBox[{
  "NDSolve", "ndnl", 
   "\"Endpoint \\!\\(\\*RowBox[{\\\"25.`\\\", \\\" \\\", \\\"q\\\"}]\\) in \
\\!\\(\\*RowBox[{\\\"{\\\", RowBox[{\\\"x\\\", \\\",\\\", \\\"0.`\\\", \
\\\",\\\", RowBox[{\\\"25.`\\\", \\\" \\\", \\\"q\\\"}]}], \\\"}\\\"}]\\) is \
not a real number.\"", 2, 19, 2, 32214176670846220627, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.854653363725963*^9},
 CellLabel->
  "During evaluation of \
In[19]:=",ExpressionUUID->"b1232753-40c4-4ab6-97d8-25eb7ae68415"],

Cell[BoxData[
 TemplateBox[{
  "ReplaceAll", "reps", 
   "\"\\!\\(\\*RowBox[{\\\"{\\\", RowBox[{\\\"NDSolve\\\", \\\"[\\\", \
RowBox[{RowBox[{\\\"{\\\", RowBox[{RowBox[{RowBox[{\\\"\[ImaginaryI]\\\", \
\\\" \\\", RowBox[{SuperscriptBox[\\\"y\\\", \\\"\[Prime]\\\", \
Rule[MultilineFunction, None]], \\\"[\\\", \\\"x\\\", \\\"]\\\"}]}], \\\"\
\[Equal]\\\", RowBox[{\\\"2\\\", \\\" \\\", \\\"\[Pi]\\\", \\\" \\\", \
RowBox[{\\\"Cos\\\", \\\"[\\\", RowBox[{RowBox[{\\\"1.5707963267948966`\\\", \
\\\"\[VeryThinSpace]\\\"}], \\\"+\\\", RowBox[{\\\"Times\\\", \\\"[\\\", \
RowBox[{\\\"\[LeftSkeleton]\\\", \\\"3\\\", \\\"\[RightSkeleton]\\\"}], \\\"]\
\\\"}]}], \\\"]\\\"}], \\\" \\\", RowBox[{\\\"y\\\", \\\"[\\\", \\\"x\\\", \\\
\"]\\\"}]}]}], \\\",\\\", RowBox[{RowBox[{\\\"y\\\", \\\"[\\\", \\\"0\\\", \\\
\"]\\\"}], \\\"\[Equal]\\\", FractionBox[\\\"1\\\", SqrtBox[\\\"2\\\"]]}]}], \
\\\"}\\\"}], \\\",\\\", \\\"y\\\", \\\",\\\", RowBox[{\\\"{\\\", \
RowBox[{\\\"x\\\", \\\",\\\", \\\"0\\\", \\\",\\\", RowBox[{\\\"25.`\\\", \
\\\" \\\", \\\"q\\\"}]}], \\\"}\\\"}], \\\",\\\", RowBox[{\\\"MaxSteps\\\", \
\\\"\[Rule]\\\", \\\"\[Infinity]\\\"}]}], \\\"]\\\"}], \\\"}\\\"}]\\) is \
neither a list of replacement rules nor a valid dispatch table, and so cannot \
be used for replacing.\"", 2, 19, 3, 32214176670846220627, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.854653363729953*^9},
 CellLabel->
  "During evaluation of \
In[19]:=",ExpressionUUID->"eadf93bb-856d-46c9-94df-e3a34d59230f"],

Cell[BoxData[
 TemplateBox[{
  "ReplaceAll", "reps", 
   "\"\\!\\(\\*RowBox[{\\\"{\\\", RowBox[{\\\"NDSolve\\\", \\\"[\\\", \
RowBox[{RowBox[{\\\"{\\\", RowBox[{RowBox[{RowBox[{\\\"\[ImaginaryI]\\\", \
\\\" \\\", RowBox[{SuperscriptBox[\\\"y\\\", \\\"\[Prime]\\\", \
Rule[MultilineFunction, None]], \\\"[\\\", \\\"x\\\", \\\"]\\\"}]}], \\\"\
\[Equal]\\\", RowBox[{\\\"2\\\", \\\" \\\", \\\"\[Pi]\\\", \\\" \\\", \
RowBox[{\\\"Cos\\\", \\\"[\\\", RowBox[{RowBox[{\\\"1.5707963267948966`\\\", \
\\\"\[VeryThinSpace]\\\"}], \\\"+\\\", RowBox[{\\\"Times\\\", \\\"[\\\", \
RowBox[{\\\"\[LeftSkeleton]\\\", \\\"3\\\", \\\"\[RightSkeleton]\\\"}], \\\"]\
\\\"}]}], \\\"]\\\"}], \\\" \\\", RowBox[{\\\"y\\\", \\\"[\\\", \\\"x\\\", \\\
\"]\\\"}]}]}], \\\",\\\", RowBox[{RowBox[{\\\"y\\\", \\\"[\\\", \\\"0\\\", \\\
\"]\\\"}], \\\"\[Equal]\\\", FractionBox[\\\"1\\\", SqrtBox[\\\"2\\\"]]}]}], \
\\\"}\\\"}], \\\",\\\", \\\"y\\\", \\\",\\\", RowBox[{\\\"{\\\", \
RowBox[{\\\"x\\\", \\\",\\\", \\\"0\\\", \\\",\\\", RowBox[{\\\"25.`\\\", \
\\\" \\\", \\\"q\\\"}]}], \\\"}\\\"}], \\\",\\\", RowBox[{\\\"MaxSteps\\\", \
\\\"\[Rule]\\\", \\\"\[Infinity]\\\"}]}], \\\"]\\\"}], \\\"}\\\"}]\\) is \
neither a list of replacement rules nor a valid dispatch table, and so cannot \
be used for replacing.\"", 2, 19, 4, 32214176670846220627, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.8546533637356043`*^9},
 CellLabel->
  "During evaluation of \
In[19]:=",ExpressionUUID->"277aeac2-465b-4814-b700-068368b6dd81"],

Cell[BoxData[
 TemplateBox[{
  "NDSolve", "ndnl", 
   "\"Endpoint \\!\\(\\*RowBox[{\\\"25.`\\\", \\\" \\\", \\\"q\\\"}]\\) in \
\\!\\(\\*RowBox[{\\\"{\\\", RowBox[{\\\"x\\\", \\\",\\\", \\\"0.`\\\", \
\\\",\\\", RowBox[{\\\"25.`\\\", \\\" \\\", \\\"q\\\"}]}], \\\"}\\\"}]\\) is \
not a real number.\"", 2, 19, 5, 32214176670846220627, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.854653363739897*^9},
 CellLabel->
  "During evaluation of \
In[19]:=",ExpressionUUID->"5c3930c3-3d1f-4670-b7e2-81de78256651"],

Cell[BoxData[
 TemplateBox[{
  "ReplaceAll", "reps", 
   "\"\\!\\(\\*RowBox[{\\\"{\\\", RowBox[{\\\"NDSolve\\\", \\\"[\\\", \
RowBox[{RowBox[{\\\"{\\\", RowBox[{RowBox[{RowBox[{\\\"\[ImaginaryI]\\\", \
\\\" \\\", RowBox[{SuperscriptBox[\\\"y\\\", \\\"\[Prime]\\\", \
Rule[MultilineFunction, None]], \\\"[\\\", \\\"x\\\", \\\"]\\\"}]}], \\\"\
\[Equal]\\\", RowBox[{RowBox[{\\\"-\\\", \\\"2\\\"}], \\\" \\\", \
\\\"\[Pi]\\\", \\\" \\\", RowBox[{\\\"Cos\\\", \\\"[\\\", \
RowBox[{RowBox[{\\\"1.5707963267948966`\\\", \\\"\[VeryThinSpace]\\\"}], \
\\\"+\\\", RowBox[{\\\"Times\\\", \\\"[\\\", RowBox[{\\\"\[LeftSkeleton]\\\", \
\\\"3\\\", \\\"\[RightSkeleton]\\\"}], \\\"]\\\"}]}], \\\"]\\\"}], \\\" \\\", \
RowBox[{\\\"y\\\", \\\"[\\\", \\\"x\\\", \\\"]\\\"}]}]}], \\\",\\\", \
RowBox[{RowBox[{\\\"y\\\", \\\"[\\\", \\\"0\\\", \\\"]\\\"}], \
\\\"\[Equal]\\\", FractionBox[\\\"1\\\", SqrtBox[\\\"2\\\"]]}]}], \
\\\"}\\\"}], \\\",\\\", \\\"y\\\", \\\",\\\", RowBox[{\\\"{\\\", \
RowBox[{\\\"x\\\", \\\",\\\", \\\"0\\\", \\\",\\\", RowBox[{\\\"25.`\\\", \
\\\" \\\", \\\"q\\\"}]}], \\\"}\\\"}], \\\",\\\", RowBox[{\\\"MaxSteps\\\", \
\\\"\[Rule]\\\", \\\"\[Infinity]\\\"}]}], \\\"]\\\"}], \\\"}\\\"}]\\) is \
neither a list of replacement rules nor a valid dispatch table, and so cannot \
be used for replacing.\"", 2, 19, 6, 32214176670846220627, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.854653363743662*^9},
 CellLabel->
  "During evaluation of \
In[19]:=",ExpressionUUID->"b8db4c9f-4045-4d6b-a8a8-45a2a8a73d7d"],

Cell[BoxData[
 TemplateBox[{
  "General", "stop", 
   "\"Further output of \\!\\(\\*StyleBox[RowBox[{\\\"ReplaceAll\\\", \
\\\"::\\\", \\\"reps\\\"}], \\\"MessageName\\\"]\\) will be suppressed during \
this calculation.\"", 2, 19, 7, 32214176670846220627, "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{3.854653363747204*^9},
 CellLabel->
  "During evaluation of \
In[19]:=",ExpressionUUID->"09fd4127-62aa-4063-a66c-af2a68a96903"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"2", " ", 
    SuperscriptBox[
     RowBox[{"Re", "[", 
      RowBox[{"y", "[", 
       RowBox[{"25.`", " ", "q"}], "]"}], "]"}], "2"]}], ",", 
   RowBox[{"2", " ", 
    SuperscriptBox[
     RowBox[{"Im", "[", 
      RowBox[{"y", "[", 
       RowBox[{"25.`", " ", "q"}], "]"}], "]"}], "2"]}]}], "}"}]], "Output",
 CellChangeTimes->{3.85465336375082*^9},
 CellLabel->"Out[19]=",ExpressionUUID->"ad25f795-0fd8-47cd-b91b-34aa40e3163a"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{"1", ",", "2"}], "}"}], "[", 
  RowBox[{"[", "1", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.8546539553868513`*^9, 3.8546539743208933`*^9}},
 CellLabel->"In[26]:=",ExpressionUUID->"8e7099af-2f0e-42b9-b08f-a259929f9fcf"],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{{3.8546539648636913`*^9, 3.854653974538732*^9}},
 CellLabel->"Out[26]=",ExpressionUUID->"cb44b216-dca9-41a3-9e78-4701ecb72469"]
}, Open  ]]
},
WindowSize->{1920, 1043},
WindowMargins->{{-1920, Automatic}, {0, Automatic}},
TaggingRules->{
 "WelcomeScreenSettings" -> {"FEStarting" -> False}, "TryRealOnly" -> False},
FrontEndVersion->"12.1 for Linux x86 (64-bit) (June 9, 2020)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"cb6bd9fe-87b4-49f3-80ae-29f9d55e3abd"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 16409, 317, 559, "Input",ExpressionUUID->"e1f875ed-0c62-4f28-b64a-da4f248d4070"],
Cell[CellGroupData[{
Cell[16992, 341, 1302, 24, 30, "Input",ExpressionUUID->"c1dd616e-4838-483c-95ee-772fe1cdbbf3"],
Cell[18297, 367, 5991, 127, 235, "Output",ExpressionUUID->"ce7c9e20-61d4-485b-9acf-19df45f99f8d"]
}, Open  ]],
Cell[CellGroupData[{
Cell[24325, 499, 1220, 23, 97, "Input",ExpressionUUID->"046e7e43-7d67-4e52-a480-5506f5e9ab54"],
Cell[25548, 524, 247, 3, 34, "Output",ExpressionUUID->"75e85079-2191-4263-84b8-56957772077e"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25832, 532, 286, 6, 30, "Input",ExpressionUUID->"030bf7d1-d5cc-4700-8c6f-56623baef820"],
Cell[26121, 540, 929, 17, 22, "Message",ExpressionUUID->"5dde1a46-0109-4fe9-b8a6-6dda15494276"],
Cell[27053, 559, 1600, 26, 44, "Message",ExpressionUUID->"c5d7ea00-89d8-43ab-acf2-9f3a183b62ab"],
Cell[28656, 587, 1600, 26, 44, "Message",ExpressionUUID->"6810c86a-ccc6-4a9e-93d6-3f66c9f41536"],
Cell[30259, 615, 929, 17, 22, "Message",ExpressionUUID->"4ec2a5f9-a0b7-4c40-8498-d490a8bb3d63"],
Cell[31191, 634, 1621, 26, 44, "Message",ExpressionUUID->"9072c424-c21c-497a-bdd6-1c426d896f3d"],
Cell[32815, 662, 454, 10, 22, "Message",ExpressionUUID->"9dc58a0f-8345-46e7-8ae5-e5e293563f68"],
Cell[33272, 674, 1157, 23, 37, "Output",ExpressionUUID->"9d24c5f3-86f1-4a4d-9a28-ba83a3a16b7f"]
}, Open  ]],
Cell[34444, 700, 225, 3, 30, "Input",ExpressionUUID->"0f8480c7-0ad3-4df5-8536-2d1912be031e"],
Cell[34672, 705, 222, 3, 34, "Output",ExpressionUUID->"8ad523e0-fc6d-4367-9cb3-4bbd6111bdc9"],
Cell[34897, 710, 156, 3, 30, "Input",ExpressionUUID->"4163ba75-88b1-461b-9e8a-35864418bfe5"],
Cell[CellGroupData[{
Cell[35078, 717, 231, 4, 30, "Input",ExpressionUUID->"569104bd-3e2c-4653-834b-18b13f1d6d12"],
Cell[35312, 723, 447, 10, 22, "Message",ExpressionUUID->"0574253a-9b71-4337-ad0c-494e482f2f68"],
Cell[35762, 735, 532, 11, 22, "Message",ExpressionUUID->"b1232753-40c4-4ab6-97d8-25eb7ae68415"],
Cell[36297, 748, 1516, 24, 44, "Message",ExpressionUUID->"eadf93bb-856d-46c9-94df-e3a34d59230f"],
Cell[37816, 774, 1518, 24, 44, "Message",ExpressionUUID->"277aeac2-465b-4814-b700-068368b6dd81"],
Cell[39337, 800, 532, 11, 22, "Message",ExpressionUUID->"5c3930c3-3d1f-4670-b7e2-81de78256651"],
Cell[39872, 813, 1539, 25, 44, "Message",ExpressionUUID->"b8db4c9f-4045-4d6b-a8a8-45a2a8a73d7d"],
Cell[41414, 840, 454, 10, 22, "Message",ExpressionUUID->"09fd4127-62aa-4063-a66c-af2a68a96903"],
Cell[41871, 852, 483, 14, 37, "Output",ExpressionUUID->"ad25f795-0fd8-47cd-b91b-34aa40e3163a"]
}, Open  ]],
Cell[CellGroupData[{
Cell[42391, 871, 273, 6, 30, "Input",ExpressionUUID->"8e7099af-2f0e-42b9-b08f-a259929f9fcf"],
Cell[42667, 879, 175, 2, 34, "Output",ExpressionUUID->"cb44b216-dca9-41a3-9e78-4701ecb72469"]
}, Open  ]]
}
]
*)