Schwinger_minspace.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[     95118,       2706]
NotebookOptionsPosition[     88130,       2490]
NotebookOutlinePosition[     88469,       2505]
CellTagsIndexPosition[     88426,       2502]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"zz", "[", 
   RowBox[{"\[Mu]_", ",", "L_"}], "]"}], ":=", 
  RowBox[{"NIntegrate", "[", 
   RowBox[{
    FractionBox["1", 
     RowBox[{
      RowBox[{"Exp", "[", 
       RowBox[{"\[Mu]", " ", "L", " ", 
        RowBox[{"Cosh", "[", "\[Tau]", "]"}]}], "]"}], "-", "1"}]], ",", 
    RowBox[{"{", 
     RowBox[{"\[Tau]", ",", 
      RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}], ",", 
    RowBox[{"WorkingPrecision", "\[Rule]", "20"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.734625227795106*^9, 3.734625268395537*^9}, {
  3.734625298733664*^9, 3.734625305289575*^9}, {3.734625672281455*^9, 
  3.734625676692836*^9}, {3.734625736256278*^9, 3.734625739865184*^9}, {
  3.734625794000999*^9, 3.734625796036824*^9}, {3.73462584244699*^9, 
  3.734625846370076*^9}, {3.7346258988160057`*^9, 3.734625903919779*^9}},
 Background->RGBColor[1, 0.85, 0.85]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"zmu2", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"i", ",", 
      RowBox[{"Exp", "[", 
       RowBox[{"-", 
        RowBox[{"zz", "[", 
         RowBox[{"1", ",", "i"}], "]"}]}], "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "40"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.7346252008079357`*^9, 3.7346252078615*^9}, {
   3.734625609219603*^9, 3.73462563408287*^9}, {3.734625724850992*^9, 
   3.7346257295000343`*^9}, {3.734628018915842*^9, 3.734628031393264*^9}, 
   3.738847450919568*^9},
 Background->RGBColor[1, 0.85, 0.85]],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "1", ",", "0.30949779042229077648937195979247824721`19.93077444073849"}], 
    "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "2", ",", "0.77652330650203599092319571050841159446`20.59700202203226"}], 
    "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "3", ",", "0.93046262762166402687789783417237859527`21.142225186033297"}],
     "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "4", ",", "0.97763707537731313601821028217692805777`21.645569487819557"}],
     "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "5", ",", "0.99260949078047660980083135596574682883`22.129715844305665"}],
     "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "6", ",", "0.99751070418457723300282059399706490536`22.603382385735262"}],
     "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "7", ",", "0.9991502170995876652370194232087135817`23.070507418780508"}], 
    "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "8", ",", "0.99970703150575089580329351178236612017`23.533115455976002"}],
     "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "9", ",", "0.99989823361422771183161657784456479061`23.99237354992688"}], 
    "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "10", ",", 
     "0.99996443935937202124765312316035250972`24.449022701769746"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
    "11", ",", "0.9999875138886077727418223870019549419`24.903570083581457"}],
     "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "12", ",", "0.99999559833967497811870139512346771222`25.35638251907123"}],
     "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "13", ",", 
     "0.99999844308993994220968118823837702879`25.807736137052878"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
    "14", ",", "0.99999944772566160006549349834867340089`26.25784501601173"}],
     "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "15", ",", "0.99999980360924698717717204835266283484`26.70687892201969"}],
     "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "16", ",", "0.99999993001176358033120216112899325253`27.15497493460725"}],
     "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "17", ",", 
     "0.99999997501067152548259214655274386699`27.602245408864555"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
    "18", ",", 
     "0.99999999106249326874907882510338565035`28.048783616162666"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
    "19", ",", "0.99999999679865741854233989541644235466`28.49466784803615"}],
     "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "20", ",", 
     "0.9999999988517524359133591546641659999060805405320317319404`28.\
939964466979973"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "21", ",", 
     "0.999999999587646405866226715292890218668769758862735044122`29.\
384730215621076"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "22", ",", 
     "0.9999999998517529676999697759719262107137370757231265813133`29.\
82901399181766"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "23", ",", 
     "0.999999999946649097790503480030605764127713163825387588836`30.\
272858231873798"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "24", ",", 
     "0.9999999999807823624380041295057254428994301487326223609303`30.\
71630000151499"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "25", ",", 
     "0.9999999999930716768755295687607825215850805501314573194501`31.\
159371865807042"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "26", ",", 
     "0.9999999999975002452040239277040859104193578198698687692854`31.\
602102589718328"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "27", ",", 
     "0.9999999999990974270937609844325776408036314112832433442562`32.\
04451770742935"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "28", ",", 
     "0.9999999999996738930826222571775207918757917852522182790959`32.\
48663998884464"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "29", ",", 
     "0.9999999999998821009854241344976016777198651120836058204768`32.\
92848982481203"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "30", ",", 
     "0.9999999999999573504500707369542799036548650948304299353604`33.\
370085547480485"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "31", ",", 
     "0.9999999999999845632346889445118784498442705175060855791172`33.\
81144369847786"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "32", ",", 
     "0.9999999999999944098849624760079289819629848717367822158483`34.\
25257925478839"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "33", ",", 
     "0.9999999999999979746775352281610994239556285376337642765416`34.\
693505820092916"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "34", ",", 
     "0.999999999999999265884806828697825163595110691429559322097`35.\
13423578772234"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "35", ",", 
     "0.999999999999999733792970171410545950557504219818472712106`35.\
57478048013312"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "36", ",", 
     "0.9999999999999999034283323377805109749579001189897362485282`36.\
01515026885106"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "37", ",", 
     "0.9999999999999999649534459906487538280857603354315590001541`36.\
4553546780767"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "38", ",", 
     "0.9999999999999999872767656822619048791966870430074178575145`36.\
89540247455265"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "39", ",", 
     "0.9999999999999999953794012724326160781736218965813227149636`37.\
33530174582239"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "40", ",", 
     "0.9999999999999999983214277799797195694895134852935399431824`37.\
77505996863608"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.7346256365924253`*^9, 3.734625683023443*^9, 3.734625730812706*^9, 
   3.734625800118856*^9, {3.734625849978868*^9, 3.73462585399599*^9}, 
   3.734625908011836*^9, 3.734628072447546*^9, 3.7346744425769053`*^9, 
   3.7348458747150183`*^9, 3.734872154808625*^9, 3.7348735267010107`*^9, 
   3.73546492451991*^9, 3.735493256958192*^9, 3.735572420185206*^9, {
   3.738847447165751*^9, 3.738847451910727*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Export", "[", 
  RowBox[{"\"\<ztest\>\"", ",", "zmu2", ",", "\"\<Table\>\"", ",", 
   RowBox[{"\"\<FieldSeparators\>\"", "\[Rule]", "\"\< \>\""}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.733466055078895*^9, 3.73346609592938*^9}, {
   3.733466126079791*^9, 3.7334661821602592`*^9}, {3.733466314852124*^9, 
   3.7334663211859827`*^9}, 3.733466363501338*^9, {3.7334664173982286`*^9, 
   3.733466421357209*^9}, {3.733466570172194*^9, 3.733466588097348*^9}, {
   3.7334667917191973`*^9, 3.733466792200686*^9}, {3.7334668966417933`*^9, 
   3.7334669053113937`*^9}, {3.735464964750519*^9, 3.735464965987082*^9}}],

Cell[BoxData["\<\"ztest\"\>"], "Output",
 CellChangeTimes->{3.733466906679335*^9, 3.733496903071602*^9, 
  3.733807981576309*^9, 3.735464932603867*^9, 3.7354649677740707`*^9, 
  3.735572429227949*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Directory", "[", "]"}]], "Input",
 CellChangeTimes->{{3.7355724544467154`*^9, 3.7355724569238234`*^9}}],

Cell[BoxData["\<\"/home/raynor/ZeroMode/th0/m1\"\>"], "Output",
 CellChangeTimes->{3.735572457626676*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"CreateDirectory", "[", "\"\<ZeroMode\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.7334680669865437`*^9, 3.733468075007784*^9}}],

Cell[BoxData["\<\"/home/raynor/ZeroMode/th0/m4/ZeroMode\"\>"], "Output",
 CellChangeTimes->{3.733468075498551*^9, 3.733496906586548*^9, 
  3.733807985324498*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"SetDirectory", "[", "\"\</home/raynor/ZeroMode\>\"", "]"}]], "Input",\

 CellChangeTimes->{{3.733468077146529*^9, 3.7334680900860367`*^9}, {
  3.733856816688077*^9, 3.733856818869628*^9}, {3.734625928932616*^9, 
  3.7346259330046883`*^9}, {3.734851488401503*^9, 3.734851502260173*^9}}],

Cell[BoxData["\<\"/home/raynor/ZeroMode\"\>"], "Output",
 CellChangeTimes->{{3.7334680872116117`*^9, 3.733468090417677*^9}, 
   3.733477997767555*^9, 3.73349690936483*^9, {3.733502160253498*^9, 
   3.733502165909774*^9}, 3.733807986303892*^9, 3.733856819830125*^9, {
   3.7346259210535707`*^9, 3.734625933457984*^9}, 3.734674448582515*^9, 
   3.73484591148062*^9, {3.734851490152409*^9, 3.734851504591165*^9}, 
   3.7354650334755774`*^9, 3.73884746441096*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"SetDirectory", "[", 
  RowBox[{"ParentDirectory", "[", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.733468100470571*^9, 3.733468123291256*^9}}],

Cell[BoxData["\<\"/home\"\>"], "Output",
 CellChangeTimes->{
  3.733468125714919*^9, 3.733502163480158*^9, {3.733827539853858*^9, 
   3.7338275424488163`*^9}, 3.73385680999263*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"ThetaList", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"\[Pi]", "/", "12"}], ",", 
     RowBox[{"\[Pi]", "/", "6"}], ",", 
     RowBox[{"\[Pi]", "/", "4"}], ",", 
     RowBox[{"\[Pi]", "/", "3"}], ",", 
     RowBox[{"5", 
      RowBox[{"\[Pi]", "/", "12"}]}], ",", 
     RowBox[{"\[Pi]", "/", "2"}], ",", 
     RowBox[{"7", 
      RowBox[{"\[Pi]", "/", "12"}]}], ",", 
     RowBox[{"2", 
      RowBox[{"\[Pi]", "/", "3"}]}], ",", 
     RowBox[{"3", 
      RowBox[{"\[Pi]", "/", "4"}]}], ",", 
     RowBox[{"5", 
      RowBox[{"\[Pi]", "/", "6"}]}], ",", 
     RowBox[{"11", 
      RowBox[{"\[Pi]", "/", "12"}]}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"massList", "=", 
   RowBox[{"{", 
    RowBox[{
    "0.01", ",", "0.02", ",", "0.04", ",", "0.06", ",", "0.07", ",", "0.08", 
     ",", "0.09", ",", "0.1", ",", "0.2", ",", "0.4", ",", "0.7", ",", 
     "1.0"}], "}"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.733468439400272*^9, 3.733468471171887*^9}, {
   3.733499548803596*^9, 3.733499555780487*^9}, {3.734845917700878*^9, 
   3.7348459393996887`*^9}, 3.7348514341706247`*^9, {3.734872176892109*^9, 
   3.734872181056135*^9}, {3.7348722113246813`*^9, 3.7348722224949837`*^9}, {
   3.734872258749413*^9, 3.734872330929001*^9}, {3.734872535426152*^9, 
   3.7348725736816683`*^9}},
 Background->RGBColor[1, 0.9, 0.8]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"ThetaList", "[", 
       RowBox[{"[", "i", "]"}], "]"}], ",", 
      RowBox[{"Cos", "[", 
       RowBox[{"ThetaList", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "11"}], "}"}]}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.733468474971056*^9, 3.733468482391274*^9}, {
   3.733468515255392*^9, 3.733468568812063*^9}, 3.73487218828376*^9, {
   3.734872252076455*^9, 3.734872252661229*^9}, {3.734872314426136*^9, 
   3.734872314560972*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {{}, 
    {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.012833333333333334`],
      AbsoluteThickness[1.6], 
     PointBox[{{0.26179938779914946`, 0.9659258262890683}, {
      0.5235987755982989, 0.8660254037844386}, {0.7853981633974483, 
      0.7071067811865476}, {1.0471975511965979`, 0.5}, {1.3089969389957472`, 
      0.25881904510252074`}, {1.5707963267948966`, 0.}, {
      1.8325957145940461`, -0.25881904510252074`}, {
      2.0943951023931957`, -0.5}, {2.356194490192345, -0.7071067811865476}, {
      2.6179938779914944`, -0.8660254037844386}, {
      2.879793265790644, -0.9659258262890683}}]}, {}}, {}, {}, {{}, {}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  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, 2.879793265790644}, {-0.9659258262890683, 
    0.9659258262890683}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{
  3.733468569475623*^9, 3.733496918221566*^9, 3.734625945255578*^9, 
   3.734674456796854*^9, {3.7348721838118*^9, 3.734872189175617*^9}, {
   3.734872224869413*^9, 3.734872332802452*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"MinCosMat", "[", 
   RowBox[{"l1_", ",", "l2_", ",", "\[Mu]L_"}], "]"}], ":=", 
  RowBox[{"Sum", "[", 
   RowBox[{
    RowBox[{
     FractionBox[
      RowBox[{
       RowBox[{"Factorial", "[", "l1", "]"}], 
       RowBox[{"Factorial", "[", "l2", "]"}]}], 
      RowBox[{
       RowBox[{"Factorial", "[", "n", "]"}], 
       RowBox[{"Factorial", "[", 
        RowBox[{"l1", "-", "n"}], "]"}], 
       RowBox[{"Factorial", "[", 
        RowBox[{"l2", "-", "n"}], "]"}]}]], 
     RowBox[{"(", 
      SuperscriptBox[
       RowBox[{"(", 
        RowBox[{"-", 
         FractionBox[
          RowBox[{"2", "\[Pi]"}], "\[Mu]L"]}], ")"}], 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{"l1", "+", "l2"}], ")"}], "/", "2"}], "-", "n"}]], ")"}]}], 
    ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "0", ",", 
      RowBox[{"Min", "[", 
       RowBox[{"l1", ",", "l2"}], "]"}]}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.733408779723528*^9, 3.733408785401739*^9}, {
   3.733408900170649*^9, 3.73340902501005*^9}, {3.733409365614296*^9, 
   3.733409383570045*^9}, {3.73340953119342*^9, 3.733409544655534*^9}, 
   3.733412443451433*^9, {3.7334126718248987`*^9, 3.733412675558169*^9}, 
   3.733412708241153*^9, {3.733412740553808*^9, 3.733412752230485*^9}, 
   3.733471231697412*^9},
 Background->RGBColor[1, 0.85, 0.85]],

Cell[BoxData[
 RowBox[{
  RowBox[{"MinSinMat", "[", 
   RowBox[{"l1_", ",", "l2_", ",", "\[Mu]L_"}], "]"}], ":=", 
  RowBox[{
   SqrtBox[
    FractionBox[
     RowBox[{"2", "\[Pi]"}], "\[Mu]L"]], 
   RowBox[{"Sum", "[", 
    RowBox[{
     RowBox[{
      FractionBox[
       RowBox[{
        RowBox[{"Factorial", "[", "l1", "]"}], 
        RowBox[{"Factorial", "[", "l2", "]"}]}], 
       RowBox[{
        RowBox[{"Factorial", "[", "n", "]"}], 
        RowBox[{"Factorial", "[", 
         RowBox[{"l1", "-", "n"}], "]"}], 
        RowBox[{"Factorial", "[", 
         RowBox[{"l2", "-", "n"}], "]"}]}]], 
      RowBox[{"(", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"-", 
          FractionBox[
           RowBox[{"2", "\[Pi]"}], "\[Mu]L"]}], ")"}], 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{"l1", "+", "l2", "-", "1"}], ")"}], "/", "2"}], "-", 
         "n"}]], ")"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", "0", ",", 
       RowBox[{"Min", "[", 
        RowBox[{"l1", ",", "l2"}], "]"}]}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.733471002189913*^9, 3.733471002495908*^9}, {
  3.733471360533267*^9, 3.733471403337625*^9}},
 Background->RGBColor[1, 0.85, 0.85]],

Cell[BoxData[
 RowBox[{
  RowBox[{"Nmax", "=", "64"}], ";"}]], "Input",
 CellChangeTimes->{{3.733469356899249*^9, 3.733469365088695*^9}, {
   3.733474955346488*^9, 3.733474955486104*^9}, {3.733475221300501*^9, 
   3.733475268030541*^9}, {3.733475570840164*^9, 3.733475571518869*^9}, 
   3.733481005747081*^9, {3.73349840829177*^9, 3.733498408436243*^9}, {
   3.7338274860309343`*^9, 3.7338274868564157`*^9}, 3.7338292043685493`*^9, {
   3.733856831374539*^9, 3.733856831841073*^9}, 3.733892221053185*^9, {
   3.734851449748386*^9, 3.734851449986897*^9}, {3.734872616598342*^9, 
   3.734872617118101*^9}, {3.734873052965027*^9, 3.734873053020329*^9}, {
   3.734873196708798*^9, 3.7348731972523212`*^9}, {3.7348735749277554`*^9, 
   3.734873575585565*^9}, 3.735643080941182*^9, {3.738847496638112*^9, 
   3.738847503183823*^9}, {3.738860676796172*^9, 3.738860676890194*^9}},
 Background->RGBColor[1, 0.9, 0.8]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"e0p", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"2", "i"}], "+", 
      RowBox[{"1", "/", "2"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"H0p", "=", 
   RowBox[{"DiagonalMatrix", "[", "e0p", "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.733413615293117*^9, 3.733413630929587*^9}, 
   3.733413669330497*^9, {3.733469377975369*^9, 3.733469378732826*^9}, {
   3.733469994777296*^9, 3.733470014999022*^9}},
 Background->RGBColor[1, 0.85, 0.85]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"e0o", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"2", "i"}], "+", "1", "+", 
      RowBox[{"1", "/", "2"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"H0o", "=", 
   RowBox[{"DiagonalMatrix", "[", "e0o", "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.733470005568169*^9, 3.7334700173655863`*^9}},
 Background->RGBColor[1, 0.85, 0.85]],

Cell[BoxData[
 RowBox[{
  RowBox[{"SetUpMatricesToFile", "[", 
   RowBox[{"L_", ",", "fname_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"Mp", ",", "Mo", ",", "MSh"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"Mp", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"MinCosMat", "[", 
         RowBox[{
          RowBox[{"2", "l1"}], ",", 
          RowBox[{"2", "l2"}], ",", "L"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"l1", ",", "0", ",", "Nmax"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"l2", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"M2p", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Chop", "[", 
         RowBox[{"N", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Mp", "[", 
             RowBox[{"[", 
              RowBox[{
               RowBox[{"i", "+", "1"}], ",", 
               RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "/", 
            SqrtBox[
             RowBox[{
              RowBox[{"Factorial", "[", 
               RowBox[{"2", "i"}], "]"}], 
              RowBox[{"Factorial", "[", 
               RowBox[{"2", "j"}], "]"}]}]]}], ",", "20"}], "]"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{"fname", "<>", "\"\<mee\>\""}], ",", "M2p", ",", 
       "\"\<Table\>\"", ",", 
       RowBox[{"\"\<FieldSeparators\>\"", "\[Rule]", "\"\< \>\""}]}], "]"}], 
     ";", "\[IndentingNewLine]", 
     RowBox[{"Mo", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"MinCosMat", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"2", "l1"}], "+", "1"}], ",", 
          RowBox[{
           RowBox[{"2", "l2"}], "+", "1"}], ",", "L"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"l1", ",", "0", ",", "Nmax"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"l2", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"M2o", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Chop", "[", 
         RowBox[{"N", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Mo", "[", 
             RowBox[{"[", 
              RowBox[{
               RowBox[{"i", "+", "1"}], ",", 
               RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "/", 
            SqrtBox[
             RowBox[{
              RowBox[{"Factorial", "[", 
               RowBox[{
                RowBox[{"2", "i"}], "+", "1"}], "]"}], 
              RowBox[{"Factorial", "[", 
               RowBox[{
                RowBox[{"2", "j"}], "+", "1"}], "]"}]}]]}], ",", "20"}], 
          "]"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{"fname", "<>", "\"\<moo\>\""}], ",", "M2o", ",", 
       "\"\<Table\>\"", ",", 
       RowBox[{"\"\<FieldSeparators\>\"", "\[Rule]", "\"\< \>\""}]}], "]"}], 
     ";", "\[IndentingNewLine]", 
     RowBox[{"MSh", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"MinSinMat", "[", 
         RowBox[{
          RowBox[{"2", "l1"}], ",", 
          RowBox[{
           RowBox[{"2", "l2"}], "+", "1"}], ",", "L"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"l1", ",", "0", ",", "Nmax"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"l2", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"M2Sh", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Chop", "[", 
         RowBox[{"N", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"MSh", "[", 
             RowBox[{"[", 
              RowBox[{
               RowBox[{"i", "+", "1"}], ",", 
               RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "/", 
            SqrtBox[
             RowBox[{
              RowBox[{"Factorial", "[", 
               RowBox[{"2", "i"}], "]"}], 
              RowBox[{"Factorial", "[", 
               RowBox[{
                RowBox[{"2", "j"}], "+", "1"}], "]"}]}]]}], ",", "20"}], 
          "]"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{"fname", "<>", "\"\<meo\>\""}], ",", "M2Sh", ",", 
       "\"\<Table\>\"", ",", 
       RowBox[{"\"\<FieldSeparators\>\"", "\[Rule]", "\"\< \>\""}]}], "]"}], 
     ";"}]}], "\[IndentingNewLine]", "]"}]}]], "Input",
 CellChangeTimes->{{3.733498532540843*^9, 3.733498612324437*^9}, {
  3.733498679380898*^9, 3.733498682829219*^9}, {3.73349875098172*^9, 
  3.733498752005642*^9}, {3.733498866903062*^9, 3.733498867178609*^9}, {
  3.734628111915007*^9, 3.734628117254986*^9}, {3.734674250322586*^9, 
  3.734674258030781*^9}, {3.734674321164205*^9, 3.7346743257462807`*^9}, {
  3.73546507681435*^9, 3.735465104595914*^9}, {3.735465261038196*^9, 
  3.735465303964891*^9}, {3.735465537098815*^9, 3.735465537415514*^9}, {
  3.735465741502631*^9, 3.735465778321465*^9}, {3.73549318014614*^9, 
  3.735493185250997*^9}, {3.735493228668056*^9, 3.735493240963909*^9}},
 Background->RGBColor[1, 1, 0.85]],

Cell[BoxData[
 RowBox[{
  RowBox[{"SetUpMatrices", "[", "L_", "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"Mp", ",", "Mo", ",", "MSh"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"Mp", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"MinCosMat", "[", 
         RowBox[{
          RowBox[{"2", "l1"}], ",", 
          RowBox[{"2", "l2"}], ",", "L"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"l1", ",", "0", ",", "Nmax"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"l2", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"M2p", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"zmu2", "[", 
          RowBox[{"[", 
           RowBox[{"L", ",", "2"}], "]"}], "]"}], "*", 
         RowBox[{"N", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Mp", "[", 
             RowBox[{"[", 
              RowBox[{
               RowBox[{"i", "+", "1"}], ",", 
               RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "/", 
            SqrtBox[
             RowBox[{
              RowBox[{"Factorial", "[", 
               RowBox[{"2", "i"}], "]"}], 
              RowBox[{"Factorial", "[", 
               RowBox[{"2", "j"}], "]"}]}]]}], ",", "20"}], "]"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Mo", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"MinCosMat", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"2", "l1"}], "+", "1"}], ",", 
          RowBox[{
           RowBox[{"2", "l2"}], "+", "1"}], ",", "L"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"l1", ",", "0", ",", "Nmax"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"l2", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"M2o", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"zmu2", "[", 
          RowBox[{"[", 
           RowBox[{"L", ",", "2"}], "]"}], "]"}], "*", 
         RowBox[{"N", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"Mo", "[", 
             RowBox[{"[", 
              RowBox[{
               RowBox[{"i", "+", "1"}], ",", 
               RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "/", 
            SqrtBox[
             RowBox[{
              RowBox[{"Factorial", "[", 
               RowBox[{
                RowBox[{"2", "i"}], "+", "1"}], "]"}], 
              RowBox[{"Factorial", "[", 
               RowBox[{
                RowBox[{"2", "j"}], "+", "1"}], "]"}]}]]}], ",", "20"}], 
          "]"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"MSh", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"MinSinMat", "[", 
         RowBox[{
          RowBox[{"2", "l1"}], ",", 
          RowBox[{
           RowBox[{"2", "l2"}], "+", "1"}], ",", "L"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"l1", ",", "0", ",", "Nmax"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"l2", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"M2Sh", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"zmu2", "[", 
          RowBox[{"[", 
           RowBox[{"L", ",", "2"}], "]"}], "]"}], "*", 
         RowBox[{"N", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"MSh", "[", 
             RowBox[{"[", 
              RowBox[{
               RowBox[{"i", "+", "1"}], ",", 
               RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "/", 
            SqrtBox[
             RowBox[{
              RowBox[{"Factorial", "[", 
               RowBox[{"2", "i"}], "]"}], 
              RowBox[{"Factorial", "[", 
               RowBox[{
                RowBox[{"2", "j"}], "+", "1"}], "]"}]}]]}], ",", "20"}], 
          "]"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], ";"}]}], 
   "\[IndentingNewLine]", "]"}]}]], "Input",
 CellChangeTimes->{{3.733498532540843*^9, 3.733498612324437*^9}, {
  3.733498679380898*^9, 3.733498682829219*^9}, {3.73349875098172*^9, 
  3.733498752005642*^9}, {3.733498866903062*^9, 3.733498867178609*^9}, {
  3.734628111915007*^9, 3.734628117254986*^9}, {3.734674250322586*^9, 
  3.734674258030781*^9}, {3.734674321164205*^9, 3.7346743257462807`*^9}, {
  3.738861519176379*^9, 3.738861520593285*^9}},
 Background->RGBColor[1, 1, 0.85]],

Cell[BoxData["First"], "Input",
 CellChangeTimes->{{3.7334986382777767`*^9, 3.7334986498551683`*^9}, {
   3.7334988116597548`*^9, 3.7334988159156237`*^9}, 3.733498869604169*^9, {
   3.734850804112082*^9, 3.734850810486835*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SetUpMatrices", "[", "1", "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.733498624828764*^9, 3.73349863702061*^9}, {
   3.733498755744734*^9, 3.733498756085314*^9}, 3.7334988773506937`*^9}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "m"}], " ", 
    RowBox[{"mass", ":", " ", "\[Theta]"}]}], "\[Rule]", "\[Pi]"}], " ", 
  "*)"}]], "Input",
 CellChangeTimes->{{3.733499935520453*^9, 3.733499943833845*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"CalculateMinSpace", "[", 
   RowBox[{"m_", ",", "Nmax_", ",", "fname_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "veigp", ",", "vveigp", ",", "vecp", ",", "evalvecp", ",", "veigo", ",", 
      "vveigo", ",", "veco", ",", "evalveco", ",", "chmp", ",", "chmo", ",", 
      "shm", ",", "test"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"veigp", "=", 
      RowBox[{"Sort", "[", 
       RowBox[{"Eigenvalues", "[", 
        RowBox[{"H0p", "+", 
         RowBox[{"m", "*", "M2p"}]}], "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{"fname", "<>", "\"\<evalp\>\""}], ",", 
       RowBox[{"Take", "[", 
        RowBox[{"veigp", ",", "11"}], "]"}], ",", "\"\<Table\>\"", ",", 
       RowBox[{"\"\<FieldSeparators\>\"", "\[Rule]", "\"\< \>\""}]}], "]"}], 
     ";", "\[IndentingNewLine]", 
     RowBox[{"vveigp", "=", 
      RowBox[{"Eigenvectors", "[", 
       RowBox[{"H0p", "+", 
        RowBox[{"m", "*", "M2p"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"vecp", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"Sign", "[", 
          RowBox[{"vveigp", "[", 
           RowBox[{"[", 
            RowBox[{
             RowBox[{"Nmax", "+", "1", "-", "i"}], ",", "1"}], "]"}], "]"}], 
          "]"}], 
         RowBox[{"vveigp", "[", 
          RowBox[{"[", 
           RowBox[{"Nmax", "+", "1", "-", "i"}], "]"}], "]"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "0", ",", "10"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"evalvecp", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"vecp", "[", 
            RowBox[{"[", "i", "]"}], "]"}], ".", 
           RowBox[{"(", 
            RowBox[{"H0p", "+", 
             RowBox[{"m", "*", "M2p"}]}], ")"}], ".", 
           RowBox[{"vecp", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], ",", 
          RowBox[{"vecp", "[", 
           RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "11"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"evalvecp", "=", 
      RowBox[{"Sort", "[", 
       RowBox[{"evalvecp", ",", 
        RowBox[{
         RowBox[{
          RowBox[{"#1", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "<", 
          RowBox[{"#2", "[", 
           RowBox[{"[", "1", "]"}], "]"}]}], "&"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", "\[IndentingNewLine]", 
     RowBox[{"veigo", "=", 
      RowBox[{"Sort", "[", 
       RowBox[{"Eigenvalues", "[", 
        RowBox[{"H0o", "+", 
         RowBox[{"m", "*", "M2o"}]}], "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{"fname", "<>", "\"\<evalo\>\""}], ",", 
       RowBox[{"Take", "[", 
        RowBox[{"veigo", ",", "11"}], "]"}], ",", "\"\<Table\>\"", ",", 
       RowBox[{"\"\<FieldSeparators\>\"", "\[Rule]", "\"\< \>\""}]}], "]"}], 
     ";", "\[IndentingNewLine]", 
     RowBox[{"vveigo", "=", 
      RowBox[{"Eigenvectors", "[", 
       RowBox[{"H0o", "+", 
        RowBox[{"m", "*", "M2o"}]}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"veco", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"Sign", "[", 
          RowBox[{"vveigo", "[", 
           RowBox[{"[", 
            RowBox[{
             RowBox[{"Nmax", "+", "1", "-", "i"}], ",", "1"}], "]"}], "]"}], 
          "]"}], 
         RowBox[{"vveigo", "[", 
          RowBox[{"[", 
           RowBox[{"Nmax", "+", "1", "-", "i"}], "]"}], "]"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "0", ",", "10"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"evalveco", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"veco", "[", 
            RowBox[{"[", "i", "]"}], "]"}], ".", 
           RowBox[{"(", 
            RowBox[{"H0o", "+", 
             RowBox[{"m", "*", "M2o"}]}], ")"}], ".", 
           RowBox[{"veco", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], ",", 
          RowBox[{"veco", "[", 
           RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "11"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"evalveco", "=", 
      RowBox[{"Sort", "[", 
       RowBox[{"evalveco", ",", 
        RowBox[{
         RowBox[{
          RowBox[{"#1", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "<", 
          RowBox[{"#2", "[", 
           RowBox[{"[", "1", "]"}], "]"}]}], "&"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", "\[IndentingNewLine]", 
     RowBox[{"chmp", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"evalvecp", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "2"}], "]"}], "]"}], ".", 
         RowBox[{"(", "M2p", ")"}], ".", 
         RowBox[{"evalvecp", "[", 
          RowBox[{"[", 
           RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "11"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "1", ",", "11"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{"fname", "<>", "\"\<chmp\>\""}], ",", "chmp", ",", 
       "\"\<Table\>\"", ",", 
       RowBox[{"\"\<FieldSeparators\>\"", "\[Rule]", "\"\< \>\""}]}], "]"}], 
     ";", "\[IndentingNewLine]", 
     RowBox[{"chmo", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"evalveco", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "2"}], "]"}], "]"}], ".", 
         RowBox[{"(", "M2o", ")"}], ".", 
         RowBox[{"evalveco", "[", 
          RowBox[{"[", 
           RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "11"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "1", ",", "11"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{"fname", "<>", "\"\<chmo\>\""}], ",", "chmo", ",", 
       "\"\<Table\>\"", ",", 
       RowBox[{"\"\<FieldSeparators\>\"", "\[Rule]", "\"\< \>\""}]}], "]"}], 
     ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
     RowBox[{"shm", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"evalvecp", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "2"}], "]"}], "]"}], ".", 
         RowBox[{"(", "M2Sh", ")"}], ".", 
         RowBox[{"evalveco", "[", 
          RowBox[{"[", 
           RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "11"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "1", ",", "11"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{"fname", "<>", "\"\<shm\>\""}], ",", "shm", ",", 
       "\"\<Table\>\"", ",", 
       RowBox[{"\"\<FieldSeparators\>\"", "\[Rule]", "\"\< \>\""}]}], "]"}], 
     ";"}]}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"test", "=", 
      RowBox[{"Sort", "[", 
       RowBox[{"Eigenvalues", "[", 
        RowBox[{
         RowBox[{"DiagonalMatrix", "[", 
          RowBox[{"Take", "[", 
           RowBox[{"veigp", ",", "11"}], "]"}], "]"}], "-", 
         RowBox[{"m", "*", "chmp"}]}], "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{"fname", "<>", " ", "\"\<test\>\""}], ",", "test", ",", 
       "\"\<Table\>\"", ",", 
       RowBox[{"\"\<FieldSeparators\>\"", "\[Rule]", "\"\< \>\""}]}], "]"}], 
     ";"}], "*)"}], "\[IndentingNewLine]", "]"}]}]], "Input",
 CellChangeTimes->{{3.7334689571986647`*^9, 3.7334689806026773`*^9}, {
   3.7334691336340513`*^9, 3.733469173908971*^9}, {3.733469240638956*^9, 
   3.73346924811758*^9}, {3.733469296493482*^9, 3.733469320399613*^9}, 
   3.733469375000362*^9, {3.7334695933702374`*^9, 3.7334696131724033`*^9}, {
   3.7334697388092337`*^9, 3.733469791580422*^9}, 3.733469913801075*^9, {
   3.733469951487596*^9, 3.733469979932344*^9}, {3.733470027805374*^9, 
   3.7334700512205877`*^9}, {3.7334719302327833`*^9, 3.73347195210026*^9}, {
   3.733472220782546*^9, 3.7334722889050417`*^9}, {3.733472327175179*^9, 
   3.7334723273535967`*^9}, {3.733472581323669*^9, 3.733472641951748*^9}, {
   3.73347310127671*^9, 3.7334731989739656`*^9}, {3.733473738574994*^9, 
   3.733473764626265*^9}, {3.733473865843103*^9, 3.7334738726111307`*^9}, {
   3.733473909954258*^9, 3.733473931597665*^9}, {3.733474321584899*^9, 
   3.7334743888975554`*^9}, {3.733474421671094*^9, 3.733474495230625*^9}, {
   3.733474726702187*^9, 3.733474753091535*^9}, {3.733474792180997*^9, 
   3.733474813503304*^9}, {3.7334749106999693`*^9, 3.733474935630713*^9}, {
   3.733475020533498*^9, 3.733475034669263*^9}, {3.7334750829147463`*^9, 
   3.733475181181925*^9}, 3.733475211798917*^9, {3.733475582045281*^9, 
   3.733475597406219*^9}, {3.7334758935423813`*^9, 3.733475927495494*^9}, {
   3.7334765460714817`*^9, 3.733476557179707*^9}, {3.733476671860952*^9, 
   3.7334766771928864`*^9}, {3.73347681085137*^9, 3.733476817172117*^9}, {
   3.733476929470673*^9, 3.7334769307610197`*^9}, {3.7334770660837917`*^9, 
   3.73347707133383*^9}, {3.7334771035150013`*^9, 3.733477117861649*^9}, {
   3.7334771976796722`*^9, 3.733477202780642*^9}, {3.733477304707505*^9, 
   3.733477307776815*^9}, {3.733477360664933*^9, 3.733477402942931*^9}, {
   3.733477528366222*^9, 3.733477532706668*^9}, {3.733477720721971*^9, 
   3.73347780054713*^9}, {3.733477839993087*^9, 3.7334779626980658`*^9}, {
   3.733478115140514*^9, 3.733478126623897*^9}, {3.733478192968296*^9, 
   3.7334781933943367`*^9}, {3.7334821506002903`*^9, 3.733482153438175*^9}, {
   3.733497252712491*^9, 3.733497255855525*^9}, {3.733498239862603*^9, 
   3.73349825487364*^9}, {3.733498295579524*^9, 3.733498296005011*^9}, 
   3.733498611231659*^9, 3.733498942887816*^9, {3.7334989900728703`*^9, 
   3.733498991053515*^9}, {3.733499224047634*^9, 3.7334992415768557`*^9}, 
   3.733499636956987*^9, {3.733892205824765*^9, 3.7338922086728163`*^9}, {
   3.734674214391148*^9, 3.734674224677894*^9}, {3.7348506953874683`*^9, 
   3.734850768056313*^9}, {3.734850835847754*^9, 3.7348508439144707`*^9}, {
   3.734851536040038*^9, 3.7348515397657537`*^9}},
 Background->RGBColor[0.87, 0.94, 1]],

Cell[BoxData[
 RowBox[{
  RowBox[{"ThetaBl", "=", 
   RowBox[{
    RowBox[{"DiagonalMatrix", "[", 
     RowBox[{"Table", "[", 
      RowBox[{
       SqrtBox[
        RowBox[{
         RowBox[{"2", "j"}], "+", "1"}]], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "0", ",", "Nmax"}], "}"}]}], "]"}], "]"}], "+", 
    RowBox[{"DiagonalMatrix", "[", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        SqrtBox[
         RowBox[{"2", "j"}]], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "1", ",", "Nmax"}], "}"}]}], "]"}], ",", 
      RowBox[{"-", "1"}]}], "]"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.733479911316565*^9, 3.733479974823197*^9}, {
  3.7334800443341722`*^9, 3.733480044714157*^9}, {3.73348010471761*^9, 
  3.7334801213013697`*^9}},
 Background->RGBColor[1, 0.85, 0.85]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ThetaBl", "//", "MatrixForm"}]], "Input",
 CellChangeTimes->{{3.7356430863264914`*^9, 3.7356430948671618`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1", "0", "0", "0", "0", "0", "0", "0", "0"},
     {
      SqrtBox["2"], 
      SqrtBox["3"], "0", "0", "0", "0", "0", "0", "0"},
     {"0", "2", 
      SqrtBox["5"], "0", "0", "0", "0", "0", "0"},
     {"0", "0", 
      SqrtBox["6"], 
      SqrtBox["7"], "0", "0", "0", "0", "0"},
     {"0", "0", "0", 
      RowBox[{"2", " ", 
       SqrtBox["2"]}], "3", "0", "0", "0", "0"},
     {"0", "0", "0", "0", 
      SqrtBox["10"], 
      SqrtBox["11"], "0", "0", "0"},
     {"0", "0", "0", "0", "0", 
      RowBox[{"2", " ", 
       SqrtBox["3"]}], 
      SqrtBox["13"], "0", "0"},
     {"0", "0", "0", "0", "0", "0", 
      SqrtBox["14"], 
      SqrtBox["15"], "0"},
     {"0", "0", "0", "0", "0", "0", "0", "4", 
      SqrtBox["17"]}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{{3.735643088760585*^9, 3.735643095489522*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"CalculateMinSpace", "[", 
   RowBox[{"m_", ",", "L_", ",", "\[Theta]_", ",", "Nmax_", ",", "fname_"}], 
   "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "H", ",", "veig", ",", "vveig", ",", "vec", ",", "evalvec", ",", "chm", 
      ",", "shm"}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"H", "=", 
      RowBox[{"ArrayFlatten", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"H0p", "+", 
            RowBox[{"m", "*", "M2p"}]}], ",", 
           RowBox[{"\[Theta]", 
            SqrtBox[
             FractionBox["L", 
              RowBox[{"8", "\[Pi]"}]]], "ThetaBl"}]}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"\[Theta]", 
            SqrtBox[
             FractionBox["L", 
              RowBox[{"8", "\[Pi]"}]]], 
            RowBox[{"Transpose", "[", "ThetaBl", "]"}]}], ",", 
           RowBox[{"H0o", "+", 
            RowBox[{"m", "*", "M2o"}]}]}], "}"}]}], "}"}], "]"}]}], " ", ";", 
     "\[IndentingNewLine]", 
     RowBox[{"veig", "=", 
      RowBox[{"Sort", "[", 
       RowBox[{"Eigenvalues", "[", "H", "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{"fname", "<>", "\"\<eval\>\""}], ",", 
       RowBox[{"Take", "[", 
        RowBox[{"veig", ",", "21"}], "]"}], ",", "\"\<Table\>\"", ",", 
       RowBox[{"\"\<FieldSeparators\>\"", "\[Rule]", "\"\< \>\""}]}], "]"}], 
     ";", "\[IndentingNewLine]", 
     RowBox[{"vveig", "=", 
      RowBox[{"Eigenvectors", "[", "H", "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"vec", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"Sign", "[", 
          RowBox[{"vveig", "[", 
           RowBox[{"[", 
            RowBox[{
             RowBox[{"Nmax", "+", "1", "-", "i"}], ",", "1"}], "]"}], "]"}], 
          "]"}], 
         RowBox[{"vveig", "[", 
          RowBox[{"[", 
           RowBox[{"Nmax", "+", "1", "-", "i"}], "]"}], "]"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "0", ",", "20"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"evalvec", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"vec", "[", 
            RowBox[{"[", "i", "]"}], "]"}], ".", 
           RowBox[{"(", "H", ")"}], ".", 
           RowBox[{"vec", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], ",", 
          RowBox[{"vec", "[", 
           RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"evalvec", "=", 
      RowBox[{"Sort", "[", 
       RowBox[{"evalvec", ",", 
        RowBox[{
         RowBox[{
          RowBox[{"#1", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "<", 
          RowBox[{"#2", "[", 
           RowBox[{"[", "1", "]"}], "]"}]}], "&"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", "\[IndentingNewLine]", 
     RowBox[{"chm", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"evalvec", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "2"}], "]"}], "]"}], ".", 
         RowBox[{"ArrayFlatten", "[", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"M2p", ",", 
              RowBox[{"Table", "[", 
               RowBox[{"0", ",", 
                RowBox[{"{", 
                 RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}], ",", 
                RowBox[{"{", 
                 RowBox[{"j", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], 
             "}"}], ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"Table", "[", 
               RowBox[{"0", ",", 
                RowBox[{"{", 
                 RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}], ",", 
                RowBox[{"{", 
                 RowBox[{"j", ",", "0", ",", "Nmax"}], "}"}]}], "]"}], ",", 
              "M2o"}], "}"}]}], "}"}], "]"}], ".", 
         RowBox[{"evalvec", "[", 
          RowBox[{"[", 
           RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "21"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{"fname", "<>", "\"\<chmp\>\""}], ",", "chm", ",", 
       "\"\<Table\>\"", ",", 
       RowBox[{"\"\<FieldSeparators\>\"", "\[Rule]", "\"\< \>\""}]}], "]"}], 
     ";", "\[IndentingNewLine]", "\[IndentingNewLine]", 
     RowBox[{"shm", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"evalvec", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "2"}], "]"}], "]"}], ".", 
         RowBox[{"ArrayFlatten", "[", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"Table", "[", 
               RowBox[{"0", ",", 
                RowBox[{"{", 
                 RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}], ",", 
                RowBox[{"{", 
                 RowBox[{"j", ",", "0", ",", "Nmax"}], "}"}]}], "]"}], ",", 
              "M2Sh"}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"Transpose", "[", "M2Sh", "]"}], ",", 
              RowBox[{"Table", "[", 
               RowBox[{"0", ",", 
                RowBox[{"{", 
                 RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}], ",", 
                RowBox[{"{", 
                 RowBox[{"j", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], 
             "}"}]}], "}"}], "]"}], ".", 
         RowBox[{"evalvec", "[", 
          RowBox[{"[", 
           RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "21"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Export", "[", 
      RowBox[{
       RowBox[{"fname", "<>", "\"\<shm\>\""}], ",", "shm", ",", 
       "\"\<Table\>\"", ",", 
       RowBox[{"\"\<FieldSeparators\>\"", "\[Rule]", "\"\< \>\""}]}], "]"}], 
     ";"}]}], "\[IndentingNewLine]", "]"}]}]], "Input",
 CellChangeTimes->{{3.7334689571986647`*^9, 3.7334689806026773`*^9}, {
   3.7334691336340513`*^9, 3.733469173908971*^9}, {3.733469240638956*^9, 
   3.73346924811758*^9}, {3.733469296493482*^9, 3.733469320399613*^9}, 
   3.733469375000362*^9, {3.7334695933702374`*^9, 3.7334696131724033`*^9}, {
   3.7334697388092337`*^9, 3.733469791580422*^9}, 3.733469913801075*^9, {
   3.733469951487596*^9, 3.733469979932344*^9}, {3.733470027805374*^9, 
   3.7334700512205877`*^9}, {3.7334719302327833`*^9, 3.73347195210026*^9}, {
   3.733472220782546*^9, 3.7334722889050417`*^9}, {3.733472327175179*^9, 
   3.7334723273535967`*^9}, {3.733472581323669*^9, 3.733472641951748*^9}, {
   3.73347310127671*^9, 3.7334731989739656`*^9}, {3.733473738574994*^9, 
   3.733473764626265*^9}, {3.733473865843103*^9, 3.7334738726111307`*^9}, {
   3.733473909954258*^9, 3.733473931597665*^9}, {3.733474321584899*^9, 
   3.7334743888975554`*^9}, {3.733474421671094*^9, 3.733474495230625*^9}, {
   3.733474726702187*^9, 3.733474753091535*^9}, {3.733474792180997*^9, 
   3.733474813503304*^9}, {3.7334749106999693`*^9, 3.733474935630713*^9}, {
   3.733475020533498*^9, 3.733475034669263*^9}, {3.7334750829147463`*^9, 
   3.733475181181925*^9}, 3.733475211798917*^9, {3.733475582045281*^9, 
   3.733475597406219*^9}, {3.7334758935423813`*^9, 3.733475927495494*^9}, {
   3.7334765460714817`*^9, 3.733476557179707*^9}, {3.733476671860952*^9, 
   3.7334766771928864`*^9}, {3.73347681085137*^9, 3.733476817172117*^9}, {
   3.733476929470673*^9, 3.7334769307610197`*^9}, {3.7334770660837917`*^9, 
   3.73347707133383*^9}, {3.7334771035150013`*^9, 3.733477117861649*^9}, {
   3.7334771976796722`*^9, 3.733477202780642*^9}, {3.733477304707505*^9, 
   3.733477307776815*^9}, {3.733477360664933*^9, 3.733477402942931*^9}, {
   3.733477528366222*^9, 3.733477532706668*^9}, {3.733477720721971*^9, 
   3.73347780054713*^9}, {3.733477839993087*^9, 3.7334779626980658`*^9}, {
   3.733478115140514*^9, 3.733478126623897*^9}, {3.733478192968296*^9, 
   3.7334781933943367`*^9}, {3.733478244974875*^9, 3.733478247666279*^9}, {
   3.733480188546425*^9, 3.733480232576995*^9}, {3.733480283498962*^9, 
   3.733480364915161*^9}, {3.733480400556437*^9, 3.733480543033394*^9}, {
   3.733480594740961*^9, 3.733480633172027*^9}, {3.733480805343444*^9, 
   3.733480858403584*^9}, {3.7334809150265417`*^9, 3.733480949795664*^9}, {
   3.7334811037593203`*^9, 3.7334811064929647`*^9}, {3.7334812894755783`*^9, 
   3.7334813022270203`*^9}, {3.733481360133251*^9, 3.73348136348*^9}, {
   3.7334813940474358`*^9, 3.733481439128265*^9}, {3.733481483616824*^9, 
   3.733481502295176*^9}, {3.733481813481244*^9, 3.733481819981791*^9}, {
   3.7334820322630873`*^9, 3.733482034777286*^9}, {3.7334983804607077`*^9, 
   3.733498385597549*^9}, {3.733499036779974*^9, 3.733499055282302*^9}, {
   3.7334992847952213`*^9, 3.7334992885360813`*^9}, 3.733499655411745*^9, {
   3.7334997131157227`*^9, 3.733499713910386*^9}},
 Background->RGBColor[0.88, 1, 0.88]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"SetUpMatrices", "[", "1", "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"veigo", "=", 
   RowBox[{"Sort", "[", 
    RowBox[{"Eigenvalues", "[", 
     RowBox[{"H0o", "+", 
      RowBox[{"m", "*", "M2o"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vveigo", "=", 
   RowBox[{"Sort", "[", 
    RowBox[{"Eigenvectors", "[", 
     RowBox[{"H0o", "+", 
      RowBox[{"m", "*", "M2o"}]}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"veco", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"Sign", "[", 
       RowBox[{"vveigo", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"Nmax", "+", "1", "-", "i"}], ",", "1"}], "]"}], "]"}], 
       "]"}], 
      RowBox[{"vveigo", "[", 
       RowBox[{"[", 
        RowBox[{"Nmax", "+", "1", "-", "i"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "0", ",", "10"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"evalveco", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{"veco", "[", 
         RowBox[{"[", "i", "]"}], "]"}], ".", 
        RowBox[{"(", 
         RowBox[{"H0o", "+", 
          RowBox[{"m", "*", "M2o"}]}], ")"}], ".", 
        RowBox[{"veco", "[", 
         RowBox[{"[", "i", "]"}], "]"}]}], ",", 
       RowBox[{"veco", "[", 
        RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "11"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"evalveco", "=", 
   RowBox[{"Sort", "[", 
    RowBox[{"evalveco", ",", 
     RowBox[{
      RowBox[{
       RowBox[{"#1", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "<", 
       RowBox[{"#2", "[", 
        RowBox[{"[", "1", "]"}], "]"}]}], "&"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.733826743465952*^9, 3.733826750734221*^9}, {
  3.7338269070481033`*^9, 3.733826917404687*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"CreateDirectory", "[", "\"\<test\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.734626208721593*^9, 3.734626214204255*^9}, {
  3.734846019613875*^9, 3.73484602055336*^9}}],

Cell[BoxData["\<\"/home/raynor/ZeroMode/test\"\>"], "Output",
 CellChangeTimes->{3.7346262402390823`*^9, 3.734846021547538*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"SetDirectory", "[", "\"\<test\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.734846034798349*^9, 3.7348460418942747`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"DirectoryQ", "[", "\"\<th0\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.734626242867621*^9, 3.734626275696384*^9}}],

Cell[BoxData["True"], "Output",
 CellChangeTimes->{{3.7346262503809767`*^9, 3.734626276363102*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"n", "=", "16"}], ",", " ", 
   RowBox[{"n", "\[LessEqual]", "128"}], ",", 
   RowBox[{"n", "=", 
    RowBox[{"n", "*", "2"}]}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"Nmax", "=", "n"}], ";", "\[IndentingNewLine]", 
     RowBox[{"e0p", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"2", "i"}], "+", 
         RowBox[{"1", "/", "2"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"H0p", "=", 
      RowBox[{"DiagonalMatrix", "[", "e0p", "]"}]}], ";", "\n", 
     RowBox[{"e0o", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"2", "i"}], "+", "1", "+", 
         RowBox[{"1", "/", "2"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"H0o", "=", 
      RowBox[{"DiagonalMatrix", "[", "e0o", "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Print", "[", 
      RowBox[{
       RowBox[{"SetUpMatrices", "[", "1", "]"}], "//", "Timing"}], "]"}], ";",
      "\[IndentingNewLine]", 
     RowBox[{"dname", "=", 
      RowBox[{"\"\<th0m1c\>\"", "<>", 
       RowBox[{"ToString", "[", "n", "]"}]}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"Print", "[", 
      RowBox[{
       RowBox[{"CalculateMinSpace", "[", 
        RowBox[{"2.0", ",", "Nmax", ",", "dname"}], "]"}], "//", "Timing"}], 
      "]"}]}], "\[IndentingNewLine]", "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.7338080362839413`*^9, 3.7338080757041903`*^9}, {
   3.7338081461158257`*^9, 3.7338081494737463`*^9}, {3.733808186116248*^9, 
   3.733808189004231*^9}, {3.7338275021002502`*^9, 3.7338275181384478`*^9}, {
   3.734846081882058*^9, 3.734846155022966*^9}, {3.734846272783415*^9, 
   3.73484628352283*^9}, {3.734846405388006*^9, 3.734846423379661*^9}, {
   3.734846560457932*^9, 3.7348465607973537`*^9}, {3.734846832463129*^9, 
   3.734846874325816*^9}, {3.734847333179236*^9, 3.734847338282959*^9}, {
   3.734848465002112*^9, 3.7348484766704893`*^9}, 3.7348510578627987`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.271956`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.734846876042982*^9, 3.7348473389476433`*^9, 
  3.73484848448643*^9, 3.734850773946789*^9, 3.734850850429861*^9, 
  3.7348510589060173`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.175359`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.734846876042982*^9, 3.7348473389476433`*^9, 
  3.73484848448643*^9, 3.734850773946789*^9, 3.734850850429861*^9, 
  3.7348510589766073`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"2.605742`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.734846876042982*^9, 3.7348473389476433`*^9, 
  3.73484848448643*^9, 3.734850773946789*^9, 3.734850850429861*^9, 
  3.734851060767062*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.310773`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.734846876042982*^9, 3.7348473389476433`*^9, 
  3.73484848448643*^9, 3.734850773946789*^9, 3.734850850429861*^9, 
  3.7348510608931513`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"14.663299`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.734846876042982*^9, 3.7348473389476433`*^9, 
  3.73484848448643*^9, 3.734850773946789*^9, 3.734850850429861*^9, 
  3.734851074501718*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1.072614`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.734846876042982*^9, 3.7348473389476433`*^9, 
  3.73484848448643*^9, 3.734850773946789*^9, 3.734850850429861*^9, 
  3.7348510748134947`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"110.994081`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.734846876042982*^9, 3.7348473389476433`*^9, 
  3.73484848448643*^9, 3.734850773946789*^9, 3.734850850429861*^9, 
  3.734851185060668*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"3.835531`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.734846876042982*^9, 3.7348473389476433`*^9, 
  3.73484848448643*^9, 3.734850773946789*^9, 3.734850850429861*^9, 
  3.7348511860936213`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Directory", "[", "]"}]], "Input",
 CellChangeTimes->{{3.7346278711250553`*^9, 3.7346279085062017`*^9}}],

Cell[BoxData["\<\"/home/raynor/th2/m10\"\>"], "Output",
 CellChangeTimes->{{3.734627873816434*^9, 3.7346279095961657`*^9}, 
   3.734873165530587*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"SetDirectory", "[", "\"\</home/raynor/ZeroMode\>\"", "]"}]], "Input",\

 CellChangeTimes->{{3.734873166845841*^9, 3.7348731881641893`*^9}}],

Cell[BoxData["\<\"/home/raynor/ZeroMode\"\>"], "Output",
 CellChangeTimes->{{3.7348731773939734`*^9, 3.734873212006113*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"\[Theta]", " ", "f\[UDoubleDot]gg\[EAcute]s"}], " ", 
  "*)"}]], "Input",
 CellChangeTimes->{{3.734872683473981*^9, 3.734872687031226*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"l", "=", "1"}], ",", " ", 
   RowBox[{"l", "\[LessEqual]", "20"}], ",", 
   RowBox[{"l", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{
      RowBox[{"SetUpMatricesToFile", "[", 
       RowBox[{"l", ",", 
        RowBox[{"\"\<l\>\"", "<>", 
         RowBox[{"ToString", "[", "l", "]"}]}]}], "]"}], "//", "Timing"}], 
     "]"}], "\[IndentingNewLine]", "}"}]}], "\[IndentingNewLine]", 
  "]"}]], "Input",
 CellChangeTimes->{{3.7354658394167967`*^9, 3.735465881110634*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"107.001327`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.735493442207893*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"102.150879`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.735493544295514*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"127.633905`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.735493671855125*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"117.797904`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.735493789584116*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"134.686941`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.735493924193112*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"126.478853`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.7354940506012087`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"132.700737`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.73549418323279*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"119.236392`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.735494302404683*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"134.293602`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.735494436630611*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"130.192046`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.73549456675592*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"136.924141`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.7354947036129313`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"130.172409`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.735494833712471*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"138.042852`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.735494971679244*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"133.572094`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.735495105174902*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"140.678854`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.735495245777837*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"131.351955`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.7354953770666523`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"144.241119`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.735495521234425*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"136.263441`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.7354956574256277`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"141.812079`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.735495799167768*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"136.346024`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7354659913592873`*^9, 3.735495935447262*^9}]
}, Open  ]]
}, Open  ]],

Cell[BoxData[
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"l", "=", "1"}], ",", " ", 
   RowBox[{"l", "\[LessEqual]", "20"}], ",", 
   RowBox[{"l", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"Print", "[", 
      RowBox[{
       RowBox[{"SetUpMatrices", "[", "l", "]"}], "//", "Timing"}], "]"}], 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"j", "=", "1"}], ",", " ", 
       RowBox[{"j", "\[LessEqual]", "11"}], ",", 
       RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"If", " ", "[", 
          RowBox[{
           RowBox[{"!", 
            RowBox[{"DirectoryQ", "[", 
             RowBox[{"\"\<th\>\"", "<>", 
              RowBox[{"ToString", "[", "j", "]"}]}], "]"}]}], ",", 
           RowBox[{"CreateDirectory", "[", 
            RowBox[{"\"\<th\>\"", "<>", 
             RowBox[{"ToString", "[", "j", "]"}]}], "]"}]}], "]"}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{"SetDirectory", "[", 
          RowBox[{"\"\<th\>\"", "<>", 
           RowBox[{"ToString", "[", "j", "]"}]}], "]"}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"For", "[", 
           RowBox[{
            RowBox[{"i", "=", "1"}], ",", " ", 
            RowBox[{"i", "\[LessEqual]", "12"}], ",", " ", 
            RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
            RowBox[{"{", "\[IndentingNewLine]", 
             RowBox[{
              RowBox[{"If", "[", 
               RowBox[{
                RowBox[{"!", 
                 RowBox[{"DirectoryQ", "[", 
                  RowBox[{"\"\<m\>\"", "<>", 
                   RowBox[{"ToString", "[", "i", "]"}]}], "]"}]}], ",", 
                RowBox[{"CreateDirectory", "[", 
                 RowBox[{"\"\<m\>\"", "<>", 
                  RowBox[{"ToString", "[", "i", "]"}]}], "]"}]}], "]"}], ";", 
              "\[IndentingNewLine]", 
              RowBox[{"SetDirectory", "[", 
               RowBox[{"\"\<m\>\"", "<>", 
                RowBox[{"ToString", "[", "i", "]"}]}], "]"}], ";", 
              "\[IndentingNewLine]", 
              RowBox[{"dname", "=", 
               RowBox[{"\"\<th\>\"", "<>", 
                RowBox[{"ToString", "[", "j", "]"}], "<>", "\"\<m\>\"", "<>", 
                
                RowBox[{"ToString", "[", "i", "]"}], "<>", "\"\<L\>\"", "<>", 
                
                RowBox[{"ToString", "[", "l", "]"}]}]}], ";", 
              "\[IndentingNewLine]", 
              RowBox[{"Print", "[", 
               RowBox[{
                RowBox[{"CalculateMinSpace", "[", 
                 RowBox[{
                  RowBox[{"massList", "[", 
                   RowBox[{"[", "i", "]"}], "]"}], ",", "l", ",", 
                  RowBox[{"ThetaList", "[", 
                   RowBox[{"[", "j", "]"}], "]"}], ",", "Nmax", ",", 
                  "dname"}], "]"}], "//", "Timing"}], "]"}], ";", 
              "\[IndentingNewLine]", 
              RowBox[{"SetDirectory", "[", 
               RowBox[{"ParentDirectory", "[", "]"}], "]"}]}], 
             "\[IndentingNewLine]", "}"}]}], "]"}], "\[IndentingNewLine]", 
          RowBox[{"SetDirectory", "[", 
           RowBox[{"ParentDirectory", "[", "]"}], "]"}]}]}], 
        "\[IndentingNewLine]", "}"}]}], "]"}]}], "\[IndentingNewLine]", 
    "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.734872700915489*^9, 3.734872755616193*^9}, {
  3.734872797747842*^9, 3.734872943377454*^9}, {3.734873011661449*^9, 
  3.734873028971923*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Directory", "[", "]"}]], "Input",
 CellChangeTimes->{{3.738861570760494*^9, 3.738861573352564*^9}}],

Cell[BoxData["\<\"/home/raynor/ZeroMode/th00L40/m8\"\>"], "Output",
 CellChangeTimes->{3.738861573822193*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"SetDirectory", "[", 
  RowBox[{"ParentDirectory", "[", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.738861580277215*^9, 3.738861590504023*^9}}],

Cell[BoxData["\<\"/home/raynor/ZeroMode\"\>"], "Output",
 CellChangeTimes->{{3.738861591142027*^9, 3.738861595077642*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"SetDirectory", "[", "\"\<th00L40_64\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.738847721485182*^9, 3.73884772341916*^9}, {
  3.7388477690374117`*^9, 3.738847776019092*^9}, {3.7388615343429537`*^9, 
  3.7388615362925673`*^9}}],

Cell[BoxData["\<\"/home/raynor/ZeroMode/th00L40_64\"\>"], "Output",
 CellChangeTimes->{3.738847723968088*^9, 3.7388477767071657`*^9, 
  3.7388615650223083`*^9, 3.738861596767001*^9}]
}, Open  ]],

Cell[BoxData["S"], "Input",
 CellChangeTimes->{{3.738847739614253*^9, 3.738847741250499*^9}}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{"parity", " ", "figyelembe", " ", "v\[EAcute]tel"}], ",", " ", 
   RowBox[{
    RowBox[{"legut\[OAcute]bb", " ", "\[Theta]"}], "=", "0"}]}], " ", 
  "*)"}]], "Input",
 CellChangeTimes->{{3.734872766390217*^9, 3.734872772901494*^9}, {
   3.734872920177289*^9, 3.734872923270636*^9}, 3.738847604253708*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"l", "=", "21"}], ",", " ", 
   RowBox[{"l", "\[LessEqual]", "40"}], ",", 
   RowBox[{"l", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"Print", "[", 
      RowBox[{
       RowBox[{"SetUpMatrices", "[", "l", "]"}], "//", "Timing"}], "]"}], 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "6"}], ",", " ", 
       RowBox[{"i", "\[LessEqual]", "6"}], ",", " ", 
       RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{"!", 
            RowBox[{"DirectoryQ", "[", 
             RowBox[{"\"\<m\>\"", "<>", 
              RowBox[{"ToString", "[", "i", "]"}]}], "]"}]}], ",", 
           RowBox[{"CreateDirectory", "[", 
            RowBox[{"\"\<m\>\"", "<>", 
             RowBox[{"ToString", "[", "i", "]"}]}], "]"}]}], "]"}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{"SetDirectory", "[", 
          RowBox[{"\"\<m\>\"", "<>", 
           RowBox[{"ToString", "[", "i", "]"}]}], "]"}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{"dname", "=", 
          RowBox[{"\"\<th0m\>\"", "<>", 
           RowBox[{"ToString", "[", "i", "]"}], "<>", "\"\<L\>\"", "<>", 
           RowBox[{"ToString", "[", "l", "]"}]}]}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{"Print", "[", 
          RowBox[{
           RowBox[{"CalculateMinSpace", "[", 
            RowBox[{
             RowBox[{
              RowBox[{"massList", "[", 
               RowBox[{"[", "i", "]"}], "]"}], "*", "l"}], ",", "Nmax", ",", 
             "dname"}], "]"}], "//", "Timing"}], "]"}], ";", 
         "\[IndentingNewLine]", 
         RowBox[{"SetDirectory", "[", 
          RowBox[{"ParentDirectory", "[", "]"}], "]"}]}], 
        "\[IndentingNewLine]", "}"}]}], "\[IndentingNewLine]", " ", "]"}]}], 
    "\[IndentingNewLine]", "}"}]}], "\[IndentingNewLine]", "]"}]], "Input",
 CellChangeTimes->{{3.7334970521694527`*^9, 3.7334970617299023`*^9}, {
   3.7334976480385942`*^9, 3.7334976546676207`*^9}, {3.733498204104307*^9, 
   3.733498213519082*^9}, {3.73350155207281*^9, 3.733501556150497*^9}, {
   3.7335016574958143`*^9, 3.733501785271907*^9}, {3.7335018373079567`*^9, 
   3.7335019080579987`*^9}, {3.7335019498199167`*^9, 3.733501965328895*^9}, {
   3.73350201633457*^9, 3.733502019448845*^9}, {3.7335020692908173`*^9, 
   3.7335020822405376`*^9}, {3.733502113503186*^9, 3.7335021397049932`*^9}, {
   3.733502208217431*^9, 3.733502212573146*^9}, {3.733502333523135*^9, 
   3.7335023369459047`*^9}, {3.733827524563004*^9, 3.7338275302418222`*^9}, 
   3.733829196888096*^9, {3.7338568604334793`*^9, 3.733856860591955*^9}, {
   3.733856899143198*^9, 3.7338569112788267`*^9}, {3.734626322370565*^9, 
   3.734626381455419*^9}, {3.734626483197896*^9, 3.73462649518288*^9}, {
   3.734627808719995*^9, 3.7346278911461353`*^9}, 3.734627927018008*^9, {
   3.734674379659286*^9, 3.734674380087924*^9}, {3.7348515532451572`*^9, 
   3.7348515533151007`*^9}, 3.734855063402506*^9, {3.738847607895138*^9, 
   3.738847624444182*^9}, {3.738851348667437*^9, 3.738851348950672*^9}, 
   3.7388615605070763`*^9, {3.738864651256185*^9, 3.7388646756487017`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"16.598594`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388646938303556`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.869209`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864694067657*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"16.755794`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388647103251762`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.841151`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388647105483007`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"16.913782`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864727241395*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.711411`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864727450699*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"16.647962`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388647435556383`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.864512`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388647437848063`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"16.722702`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864760359438*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.851618`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864760584558*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"16.497239`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388647768948917`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.880244`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388647771302843`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"17.179023`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864793894553*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.8944`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864794133144*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"16.888455`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864810599944*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.928849`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388648108457003`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"16.968418`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864827694887*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.870804`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864827926239*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"16.808911`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388648442794523`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.867243`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.73886484451162*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"17.311024`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864861254653*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.907113`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864861495775*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"15.469377`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864876914337*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.86073`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864877143137*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"17.262178`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864893887034*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.874305`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388648941200933`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"17.021548`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864910532497*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.853957`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864910760824*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"17.35543`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.73886492759542*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.718526`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388649278069143`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"17.017137`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864944421485*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.805256`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864944636354*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"17.191297`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388649614891376`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.734578`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388649617022877`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"16.757912`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388649782948112`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.8772`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388649785294037`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"16.974203`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738864995395095*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.938842`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.7388649956443987`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"16.892678`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738865012148679*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.88797`", ",", "Null"}], "}"}]], "Print",
 CellChangeTimes->{3.738865012385442*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"CalculateMinSpace", "[", 
   RowBox[{
    RowBox[{"-", "1"}], ",", "Nmax", ",", "\"\<fout\>\""}], "]"}], "//", 
  "Timing"}]], "Input",
 CellChangeTimes->{{3.733474960806519*^9, 3.733474974318275*^9}, 
   3.7334765763065233`*^9, {3.7334766200295477`*^9, 3.733476623643774*^9}, {
   3.733476704484713*^9, 3.733476729093017*^9}, 3.7334767645158453`*^9, 
   3.733476834793095*^9, 3.733477084881919*^9, {3.7334771348301277`*^9, 
   3.733477148789997*^9}, 3.733477332331402*^9, {3.7334782107054768`*^9, 
   3.733478211964508*^9}, {3.733499666692904*^9, 3.733499666862095*^9}, 
   3.733499914036415*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"9.450332`", ",", "Null"}], "}"}]], "Output",
 CellChangeTimes->{
  3.733474974943034*^9, {3.733475025738473*^9, 3.733475036282394*^9}, {
   3.733475079128872*^9, 3.73347518298205*^9}, {3.73347521534824*^9, 
   3.733475226059507*^9}, 3.733475256691387*^9, 3.733475352291507*^9, {
   3.733475643891171*^9, 3.7334756622249603`*^9}, {3.733476570839061*^9, 
   3.73347658662114*^9}, {3.733476632651669*^9, 3.733476656522292*^9}, {
   3.733476689947769*^9, 3.733476740235032*^9}, 3.733476776157453*^9, {
   3.733476832805376*^9, 3.733476871730125*^9}, 3.7334769438140783`*^9, {
   3.733477023257536*^9, 3.733477041867982*^9}, {3.73347708313078*^9, 
   3.733477096126808*^9}, {3.733477130928164*^9, 3.733477160080041*^9}, 
   3.733477216862383*^9, {3.733477324781494*^9, 3.733477343584322*^9}, 
   3.733477377240155*^9, 3.733477416485096*^9, {3.7334780138641567`*^9, 
   3.7334780371620417`*^9}, 3.733478140687796*^9, 3.7334782259263697`*^9, 
   3.733481012431415*^9, 3.733482158706955*^9, 3.733496984493058*^9, 
   3.733498499372426*^9, 3.733499677219412*^9, {3.733499911929318*^9, 
   3.733499924264049*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"CalculateMinSpace", "[", 
   RowBox[{"1", ",", "1", ",", 
    RowBox[{"\[Pi]", "/", "2"}], ",", "Nmax", ",", "\"\<fth\>\""}], "]"}], "//",
   "Timing"}]], "Input",
 CellChangeTimes->{{3.733499693844915*^9, 3.733499735476611*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"40.880326`", ",", "Null"}], "}"}]], "Output",
 CellChangeTimes->{{3.7334997038243113`*^9, 3.733499749668066*^9}, 
   3.733499819205304*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"t1", "=", 
   RowBox[{"{", 
    RowBox[{"1", ",", "2", ",", "3", ",", "4", ",", "5"}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"t2", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"2", ",", "3"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"1", ",", "4"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"4", ",", "1"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"3", ",", "9"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"5", ",", "0"}], "}"}]}], "}"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.733472665397895*^9, 3.733472674875557*^9}, {
  3.733473223386043*^9, 3.733473236625614*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"t2", "=", 
  RowBox[{"Sort", "[", 
   RowBox[{"t2", ",", 
    RowBox[{
     RowBox[{
      RowBox[{"#1", "[", 
       RowBox[{"[", "1", "]"}], "]"}], "<", 
      RowBox[{"#2", "[", 
       RowBox[{"[", "1", "]"}], "]"}]}], "&"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.7334732401639023`*^9, 3.733473242212328*^9}, {
  3.733473302663866*^9, 3.7334733200663223`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"1", ",", "4"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"2", ",", "3"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"3", ",", "9"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"4", ",", "1"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"5", ",", "0"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.733473243520998*^9, {3.733473311725566*^9, 3.733473324163012*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"t1", "[", 
  RowBox[{"[", "1", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.733472900942195*^9, 3.7334729030040894`*^9}}],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{3.733472903441875*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"FirstPosition", "[", 
  RowBox[{"t2", ",", "1"}], "]"}]], "Input",
 CellChangeTimes->{{3.7334729162338057`*^9, 3.733472933106374*^9}}],

Cell[BoxData[
 RowBox[{"{", "2", "}"}]], "Output",
 CellChangeTimes->{{3.7334729216033382`*^9, 3.733472933501977*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"FirstPosition", "[", 
     RowBox[{"t2", ",", 
      RowBox[{"t1", "[", 
       RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "5"}], "}"}]}], "]"}], "//", 
  "Flatten"}]], "Input",
 CellChangeTimes->{{3.7334728710267153`*^9, 3.733472891755481*^9}, {
  3.733472936775338*^9, 3.733472946440288*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"2", ",", "1", ",", "4", ",", "3", ",", "5"}], "}"}]], "Output",
 CellChangeTimes->{
  3.733472892359647*^9, {3.7334729392795753`*^9, 3.73347294672099*^9}}]
}, Open  ]],

Cell[BoxData["Which"], "Input",
 CellChangeTimes->{{3.7334727644595957`*^9, 3.733472764824739*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Print", "[", "a", "]"}], ";", " ", 
  RowBox[{"Print", "[", "b", "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.733474063370008*^9, 3.733474069252907*^9}}],

Cell[CellGroupData[{

Cell[BoxData["a"], "Print",
 CellChangeTimes->{3.733474069676841*^9}],

Cell[BoxData["b"], "Print",
 CellChangeTimes->{3.733474069677539*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Length", "[", "H0p", "]"}]], "Input",
 CellChangeTimes->{{3.733481159500161*^9, 3.733481212845375*^9}}],

Cell[BoxData["21"], "Output",
 CellChangeTimes->{{3.733481187423502*^9, 3.733481230199687*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"L", "=", "5"}], ";"}]], "Input",
 CellChangeTimes->{{3.733481963483753*^9, 3.7334819657821207`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Mp", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"MinCosMat", "[", 
      RowBox[{
       RowBox[{"2", "l1"}], ",", 
       RowBox[{"2", "l2"}], ",", "L"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"l1", ",", "0", ",", "Nmax"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"l2", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"M2p", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"zmu", "[", 
       RowBox[{"[", 
        RowBox[{"L", ",", "2"}], "]"}], "]"}], 
      RowBox[{"N", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"Mp", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"i", "+", "1"}], ",", 
            RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "/", 
         SqrtBox[
          RowBox[{
           RowBox[{"Factorial", "[", 
            RowBox[{"2", "i"}], "]"}], 
           RowBox[{"Factorial", "[", 
            RowBox[{"2", "j"}], "]"}]}]]}], ",", "20"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Mo", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"MinCosMat", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"2", "l1"}], "+", "1"}], ",", 
       RowBox[{
        RowBox[{"2", "l2"}], "+", "1"}], ",", "L"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"l1", ",", "0", ",", "Nmax"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"l2", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"M2o", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"zmu", "[", 
       RowBox[{"[", 
        RowBox[{"L", ",", "2"}], "]"}], "]"}], 
      RowBox[{"N", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"Mo", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"i", "+", "1"}], ",", 
            RowBox[{"j", "+", "1"}]}], "]"}], "]"}], "/", 
         SqrtBox[
          RowBox[{
           RowBox[{"Factorial", "[", 
            RowBox[{
             RowBox[{"2", "i"}], "+", "1"}], "]"}], 
           RowBox[{"Factorial", "[", 
            RowBox[{
             RowBox[{"2", "j"}], "+", "1"}], "]"}]}]]}], ",", "20"}], "]"}]}],
      ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "0", ",", "Nmax"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "0", ",", "Nmax"}], "}"}]}], "]"}]}], ";"}]}], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"H", "=", 
   RowBox[{
    RowBox[{"ArrayFlatten", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"H0p", ",", 
         RowBox[{"\[Theta]", 
          SqrtBox[
           FractionBox["L", 
            RowBox[{"8", "\[Pi]"}]]], "ThetaBl"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"\[Theta]", 
          SqrtBox[
           FractionBox["L", 
            RowBox[{"8", "\[Pi]"}]]], 
          RowBox[{"Transpose", "[", "ThetaBl", "]"}]}], ",", "H0o"}], "}"}]}],
       "}"}], "]"}], "/.", 
    RowBox[{"\[Theta]", "\[Rule]", "0"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.733481246780224*^9, 3.733481282065559*^9}, {
  3.733481919796259*^9, 3.7334819316901817`*^9}, {3.7334819717626877`*^9, 
  3.733481973573126*^9}, {3.7334820073546753`*^9, 3.7334820225330343`*^9}}]
},
WindowSize->{1920, 1031},
WindowMargins->{{0, Automatic}, {0, Automatic}},
FrontEndVersion->"11.0 for Linux x86 (64-bit) (September 21, 2016)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 895, 20, 72, "Input"],
Cell[CellGroupData[{
Cell[1478, 44, 620, 16, 50, "Input"],
Cell[2101, 62, 6065, 194, 195, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8203, 261, 633, 10, 34, "Input"],
Cell[8839, 273, 200, 3, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9076, 281, 127, 2, 32, "Input"],
Cell[9206, 285, 105, 1, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9348, 291, 151, 2, 32, "Input"],
Cell[9502, 295, 161, 2, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9700, 302, 309, 5, 32, "Input"],
Cell[10012, 309, 459, 6, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10508, 320, 168, 3, 32, "Input"],
Cell[10679, 325, 180, 3, 32, "Output"]
}, Open  ]],
Cell[10874, 331, 1385, 35, 73, "Input"],
Cell[CellGroupData[{
Cell[12284, 370, 639, 16, 34, "Input"],
Cell[12926, 388, 2042, 45, 225, "Output"]
}, Open  ]],
Cell[14983, 436, 1383, 38, 73, "Input"],
Cell[16369, 476, 1242, 38, 91, "Input"],
Cell[17614, 516, 908, 14, 50, "Input"],
Cell[18525, 532, 613, 17, 73, "Input"],
Cell[19141, 551, 503, 15, 73, "Input"],
Cell[19647, 568, 5639, 149, 355, "Input"],
Cell[25289, 719, 4890, 136, 286, "Input"],
Cell[30182, 857, 227, 3, 32, "Input"],
Cell[30412, 862, 223, 4, 34, "Input"],
Cell[30638, 868, 242, 7, 32, "Input"],
Cell[30883, 877, 10692, 253, 620, "Input"],
Cell[41578, 1132, 818, 24, 65, "Input"],
Cell[CellGroupData[{
Cell[42421, 1160, 135, 2, 32, "Input"],
Cell[42559, 1164, 1355, 40, 214, "Output"]
}, Open  ]],
Cell[43929, 1207, 9378, 214, 435, "Input"],
Cell[53310, 1423, 2018, 64, 149, "Input"],
Cell[CellGroupData[{
Cell[55353, 1491, 193, 3, 32, "Input"],
Cell[55549, 1496, 127, 1, 32, "Output"]
}, Open  ]],
Cell[55691, 1500, 144, 2, 32, "Input"],
Cell[CellGroupData[{
Cell[55860, 1506, 139, 2, 32, "Input"],
Cell[56002, 1510, 99, 1, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[56138, 1516, 2235, 52, 262, "Input"],
Cell[CellGroupData[{
Cell[58398, 1572, 244, 5, 26, "Print"],
Cell[58645, 1579, 244, 5, 26, "Print"],
Cell[58892, 1586, 242, 5, 26, "Print"],
Cell[59137, 1593, 244, 5, 26, "Print"],
Cell[59384, 1600, 243, 5, 26, "Print"],
Cell[59630, 1607, 244, 5, 26, "Print"],
Cell[59877, 1614, 244, 5, 26, "Print"],
Cell[60124, 1621, 244, 5, 26, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[60417, 1632, 127, 2, 32, "Input"],
Cell[60547, 1636, 149, 2, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[60733, 1643, 163, 3, 32, "Input"],
Cell[60899, 1648, 124, 1, 32, "Output"]
}, Open  ]],
Cell[61038, 1652, 185, 4, 32, "Input"],
Cell[CellGroupData[{
Cell[61248, 1660, 601, 15, 123, "Input"],
Cell[CellGroupData[{
Cell[61874, 1679, 151, 3, 26, "Print"],
Cell[62028, 1684, 151, 3, 26, "Print"],
Cell[62182, 1689, 151, 3, 26, "Print"],
Cell[62336, 1694, 151, 3, 26, "Print"],
Cell[62490, 1699, 151, 3, 26, "Print"],
Cell[62644, 1704, 153, 3, 26, "Print"],
Cell[62800, 1709, 150, 3, 26, "Print"],
Cell[62953, 1714, 151, 3, 26, "Print"],
Cell[63107, 1719, 151, 3, 26, "Print"],
Cell[63261, 1724, 150, 3, 26, "Print"],
Cell[63414, 1729, 153, 3, 26, "Print"],
Cell[63570, 1734, 151, 3, 26, "Print"],
Cell[63724, 1739, 151, 3, 26, "Print"],
Cell[63878, 1744, 151, 3, 26, "Print"],
Cell[64032, 1749, 151, 3, 26, "Print"],
Cell[64186, 1754, 153, 3, 26, "Print"],
Cell[64342, 1759, 151, 3, 26, "Print"],
Cell[64496, 1764, 153, 3, 26, "Print"],
Cell[64652, 1769, 151, 3, 26, "Print"],
Cell[64806, 1774, 151, 3, 26, "Print"]
}, Open  ]]
}, Open  ]],
Cell[64984, 1781, 3592, 82, 416, "Input"],
Cell[CellGroupData[{
Cell[68601, 1867, 123, 2, 32, "Input"],
Cell[68727, 1871, 109, 1, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[68873, 1877, 168, 3, 32, "Input"],
Cell[69044, 1882, 122, 1, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[69203, 1888, 251, 4, 34, "Input"],
Cell[69457, 1894, 182, 2, 32, "Output"]
}, Open  ]],
Cell[69654, 1899, 93, 1, 32, "Input"],
Cell[69750, 1902, 361, 8, 34, "Input"],
Cell[CellGroupData[{
Cell[70136, 1914, 3337, 67, 325, "Input"],
Cell[CellGroupData[{
Cell[73498, 1985, 128, 3, 26, "Print"],
Cell[73629, 1990, 125, 3, 26, "Print"],
Cell[73757, 1995, 128, 3, 26, "Print"],
Cell[73888, 2000, 127, 3, 26, "Print"],
Cell[74018, 2005, 126, 3, 26, "Print"],
Cell[74147, 2010, 125, 3, 26, "Print"],
Cell[74275, 2015, 128, 3, 26, "Print"],
Cell[74406, 2020, 127, 3, 26, "Print"],
Cell[74536, 2025, 126, 3, 26, "Print"],
Cell[74665, 2030, 125, 3, 26, "Print"],
Cell[74793, 2035, 128, 3, 26, "Print"],
Cell[74924, 2040, 127, 3, 26, "Print"],
Cell[75054, 2045, 126, 3, 26, "Print"],
Cell[75183, 2050, 123, 3, 26, "Print"],
Cell[75309, 2055, 126, 3, 26, "Print"],
Cell[75438, 2060, 127, 3, 26, "Print"],
Cell[75568, 2065, 126, 3, 26, "Print"],
Cell[75697, 2070, 125, 3, 26, "Print"],
Cell[75825, 2075, 128, 3, 26, "Print"],
Cell[75956, 2080, 124, 3, 26, "Print"],
Cell[76083, 2085, 126, 3, 26, "Print"],
Cell[76212, 2090, 125, 3, 26, "Print"],
Cell[76340, 2095, 126, 3, 26, "Print"],
Cell[76469, 2100, 124, 3, 26, "Print"],
Cell[76596, 2105, 126, 3, 26, "Print"],
Cell[76725, 2110, 127, 3, 26, "Print"],
Cell[76855, 2115, 126, 3, 26, "Print"],
Cell[76984, 2120, 125, 3, 26, "Print"],
Cell[77112, 2125, 124, 3, 26, "Print"],
Cell[77239, 2130, 127, 3, 26, "Print"],
Cell[77369, 2135, 126, 3, 26, "Print"],
Cell[77498, 2140, 125, 3, 26, "Print"],
Cell[77626, 2145, 128, 3, 26, "Print"],
Cell[77757, 2150, 127, 3, 26, "Print"],
Cell[77887, 2155, 128, 3, 26, "Print"],
Cell[78018, 2160, 125, 3, 26, "Print"],
Cell[78146, 2165, 126, 3, 26, "Print"],
Cell[78275, 2170, 127, 3, 26, "Print"],
Cell[78405, 2175, 126, 3, 26, "Print"],
Cell[78534, 2180, 124, 3, 26, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[78707, 2189, 632, 12, 34, "Input"],
Cell[79342, 2203, 1143, 18, 34, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[80522, 2226, 264, 6, 34, "Input"],
Cell[80789, 2234, 179, 4, 34, "Output"]
}, Open  ]],
Cell[80983, 2241, 680, 21, 57, "Input"],
Cell[CellGroupData[{
Cell[81688, 2266, 393, 11, 34, "Input"],
Cell[82084, 2279, 439, 14, 34, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[82560, 2298, 146, 3, 32, "Input"],
Cell[82709, 2303, 70, 1, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[82816, 2309, 158, 3, 34, "Input"],
Cell[82977, 2314, 118, 2, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[83132, 2321, 423, 12, 34, "Input"],
Cell[83558, 2335, 195, 4, 34, "Output"]
}, Open  ]],
Cell[83768, 2342, 99, 1, 32, "Input"],
Cell[CellGroupData[{
Cell[83892, 2347, 189, 4, 34, "Input"],
Cell[CellGroupData[{
Cell[84106, 2355, 69, 1, 24, "Print"],
Cell[84178, 2358, 69, 1, 24, "Print"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[84296, 2365, 127, 2, 32, "Input"],
Cell[84426, 2369, 95, 1, 32, "Output"]
}, Open  ]],
Cell[84536, 2373, 135, 3, 34, "Input"],
Cell[84674, 2378, 2588, 84, 141, "Input"],
Cell[87265, 2464, 861, 24, 76, "Input"]
}
]
*)