From dfda358a7100f5491a4cc537c28d0fc718e5a544 Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Mon, 21 Dec 2020 17:37:37 +0100 Subject: Started moving work to script file. --- new_schofield.nb | 578 ++++++++++++++++++++++++++++++------------------------- schofield.wl | 54 ++++++ 2 files changed, 367 insertions(+), 265 deletions(-) create mode 100644 schofield.wl diff --git a/new_schofield.nb b/new_schofield.nb index 84ee0f3..7aafe40 100644 --- a/new_schofield.nb +++ b/new_schofield.nb @@ -10,36 +10,28 @@ NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] -NotebookDataLength[ 590390, 11609] -NotebookOptionsPosition[ 576572, 11379] -NotebookOutlinePosition[ 576986, 11396] -CellTagsIndexPosition[ 576943, 11393] +NotebookDataLength[ 591593, 11657] +NotebookOptionsPosition[ 577508, 11420] +NotebookOutlinePosition[ 577901, 11436] +CellTagsIndexPosition[ 577858, 11433] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ +Cell[BoxData[ + RowBox[{"Get", "[", + "\"\<~/doc/research/first_order_singularities/mma/schofield.wl\>\"", + "]"}]], "Input", + CellChangeTimes->{{3.817548726786456*^9, 3.817548838429791*^9}}, + CellLabel->"In[3]:=",ExpressionUUID->"e595a667-7f25-4f98-9f12-a321a236d113"], Cell[CellGroupData[{ + Cell["Definitions", "Section", CellChangeTimes->{{3.817119832436002*^9, 3.817119834596055*^9}},ExpressionUUID->"bb531fce-8a9d-4a5f-8d6a-\ ec66b9fbf698"], -Cell[BoxData[ - RowBox[{ - RowBox[{"\[Beta]", "=", - RowBox[{"1", "/", "8"}]}], ";", - RowBox[{"\[Delta]", "=", "15"}], ";", - RowBox[{"\[Alpha]", "=", "0"}], ";", - RowBox[{"\[CapitalDelta]", "=", - RowBox[{"\[Beta]", " ", "\[Delta]"}]}], ";"}]], "Input", - CellChangeTimes->{{3.8171203487765493`*^9, 3.817120359509359*^9}, { - 3.8171253080634003`*^9, 3.817125374665786*^9}, {3.81712545121809*^9, - 3.8171254551141233`*^9}, {3.817209512993226*^9, 3.81720951605719*^9}, { - 3.817216134198064*^9, 3.817216147934895*^9}}, - CellLabel-> - "(WOPR) In[1]:=",ExpressionUUID->"fafcf491-a81f-4f91-b9fe-4d27770e7577"], - Cell[BoxData[ RowBox[{ RowBox[{"$Assumptions", "=", @@ -70,144 +62,193 @@ Cell[BoxData[ CellLabel-> "(WOPR) In[446]:=",ExpressionUUID->"395191f4-75f5-4aa3-87f6-ccf9cbedf73f"], -Cell[BoxData[ - RowBox[{ - RowBox[{ - RowBox[{"f", "[", "n_", "]"}], "[", "\[Theta]_", "]"}], ":=", - RowBox[{ - SuperscriptBox[ - RowBox[{"(", - RowBox[{"\[Theta]", "/", "\[Theta]i"}], ")"}], "2"], "-", - "1"}]}]], "Input", - CellChangeTimes->{{3.8171202124436197`*^9, 3.817120269635832*^9}, { - 3.817120312676557*^9, 3.8171203389089193`*^9}, {3.8171204067661743`*^9, - 3.8171204839196587`*^9}, {3.8171205561371202`*^9, 3.8171205568649817`*^9}, { - 3.817204484948892*^9, 3.817204488556822*^9}, {3.817204991326084*^9, - 3.81720499357379*^9}, {3.817209015952724*^9, 3.8172090687475147`*^9}, { - 3.8172134845351257`*^9, 3.817213486502965*^9}, {3.817269303799518*^9, - 3.817269304543895*^9}, {3.81727116429113*^9, 3.8172712020249577`*^9}}, - CellLabel-> - "(WOPR) In[266]:=",ExpressionUUID->"3e69b8b0-cd07-4524-a460-17b8352b5902"], +Cell[CellGroupData[{ Cell[BoxData[ - RowBox[{ - RowBox[{ - RowBox[{"g", "[", "n_", "]"}], "[", "\[Theta]_", "]"}], ":=", + RowBox[{"Manipulate", "[", "\[IndentingNewLine]", RowBox[{ - RowBox[{"(", - RowBox[{"1", "-", - SuperscriptBox[ - RowBox[{"(", - FractionBox["\[Theta]", "\[Theta]c"], ")"}], "2"]}], ")"}], - RowBox[{"Sum", "[", + RowBox[{"Plot", "[", RowBox[{ RowBox[{ - RowBox[{"gC", "[", "i", "]"}], - SuperscriptBox["\[Theta]", - RowBox[{ - RowBox[{"2", "i"}], "+", "1"}]]}], ",", + RowBox[{ + RowBox[{"R\[ScriptCapitalM]", "[", "2", "]"}], "[", + RowBox[{"\[Gamma]", ",", "B", ",", "\[Theta]c", ",", "1"}], "]"}], "[", + "\[Theta]", "]"}], ",", "\[IndentingNewLine]", RowBox[{"{", - RowBox[{"i", ",", "0", ",", "n"}], "}"}]}], "]"}]}]}]], "Input", - CellChangeTimes->{{3.817119845565673*^9, 3.817119882052971*^9}, { - 3.817119939573969*^9, 3.817120069536271*^9}, {3.817209086426073*^9, - 3.8172090977061653`*^9}, 3.8172693107527933`*^9}, - CellLabel-> - "(WOPR) In[137]:=",ExpressionUUID->"daded195-b1a7-448b-87f6-d258f2dcf5e8"], + RowBox[{"\[Theta]", ",", + RowBox[{"-", "20"}], ",", "20"}], "}"}], ",", "\[IndentingNewLine]", + RowBox[{"Epilog", "\[Rule]", + RowBox[{"{", + RowBox[{"Line", "[", + RowBox[{"{", + RowBox[{ + RowBox[{"{", + RowBox[{"\[Theta]c", ",", + RowBox[{"-", "1000"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"\[Theta]c", ",", "1000"}], "}"}]}], "}"}], "]"}], "}"}]}], + ",", + RowBox[{"WorkingPrecision", "\[Rule]", "30"}]}], "]"}], + "\[IndentingNewLine]", ",", + RowBox[{"{", + RowBox[{"\[Gamma]", ",", + RowBox[{"1", "/", "10"}], ",", "10", ",", + RowBox[{"1", "/", "10"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"B", ",", + RowBox[{"1", "/", "10"}], ",", "10", ",", + RowBox[{"1", "/", "10"}]}], "}"}], ",", + RowBox[{"{", + RowBox[{"\[Theta]c", ",", + RowBox[{"1", "+", + RowBox[{"1", "/", "10"}]}], ",", "2", ",", + RowBox[{"1", "/", "100"}]}], "}"}]}], "]"}]], "Input", + CellChangeTimes->{{3.8175496131650867`*^9, 3.817549762607664*^9}, { + 3.817550982837791*^9, 3.817551030438055*^9}, {3.8175529357604218`*^9, + 3.81755293778432*^9}, {3.8175530988827267`*^9, 3.8175531070583982`*^9}, { + 3.817553402424604*^9, 3.8175534118159523`*^9}, {3.8175536922933187`*^9, + 3.81755370940517*^9}, {3.81755419838216*^9, 3.817554207237907*^9}, { + 3.8175542657350683`*^9, 3.817554403122216*^9}, {3.817555063285273*^9, + 3.8175551707746563`*^9}, {3.817557024056016*^9, 3.817557024895731*^9}, { + 3.8175573533737717`*^9, 3.817557355637789*^9}}, + CellLabel->"In[6]:=",ExpressionUUID->"389cbfb1-8cc7-458c-9f8c-1fd631f97f2c"], Cell[BoxData[ - RowBox[{ - RowBox[{"\[ScriptCapitalM]f", "[", "y_", "]"}], ":=", + TagBox[ + StyleBox[ + DynamicModuleBox[{Schofield`B$$ = Rational[9, 10], Schofield`\[Gamma]$$ = + Rational[22, 5], Schofield`\[Theta]c$$ = Rational[179, 100], + Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, + Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = + 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ + Hold[Schofield`\[Gamma]$$], + Rational[1, 10], 10, + Rational[1, 10]}, { + Hold[Schofield`B$$], + Rational[1, 10], 10, + Rational[1, 10]}, { + Hold[Schofield`\[Theta]c$$], + Rational[11, 10], 2, + Rational[1, 100]}}, Typeset`size$$ = {270., {82., 85.33243899140032}}, + Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True}, + + DynamicBox[Manipulate`ManipulateBoxes[ + 1, StandardForm, + "Variables" :> { + Schofield`B$$ = Rational[1, 10], Schofield`\[Gamma]$$ = + Rational[1, 10], Schofield`\[Theta]c$$ = Rational[11, 10]}, + "ControllerVariables" :> {}, + "OtherVariables" :> { + Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, + Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, + Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, + Typeset`skipInitDone$$}, "Body" :> Plot[ + Schofield`R\[ScriptCapitalM][2][ + Schofield`\[Gamma]$$, Schofield`B$$, Schofield`\[Theta]c$$, 1][ + Schofield`\[Theta]], {Schofield`\[Theta], -20, 20}, Epilog -> { + + Line[{{Schofield`\[Theta]c$$, -1000}, { + Schofield`\[Theta]c$$, 1000}}]}, WorkingPrecision -> 30], + "Specifications" :> {{Schofield`\[Gamma]$$, + Rational[1, 10], 10, + Rational[1, 10]}, {Schofield`B$$, + Rational[1, 10], 10, + Rational[1, 10]}, {Schofield`\[Theta]c$$, + Rational[11, 10], 2, + Rational[1, 100]}}, "Options" :> {}, "DefaultOptions" :> {}], + ImageSizeCache->{307., {134.849609375, 139.150390625}}, + SingleEvaluation->True], + Deinitialization:>None, + DynamicModuleValues:>{}, + SynchronousInitialization->True, + UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$}, + UnsavedVariables:>{Typeset`initDone$$}, + UntrackedVariables:>{Typeset`size$$}], "Manipulate", + Deployed->True, + StripOnInput->False], + Manipulate`InterpretManipulate[1]]], "Output", + CellChangeTimes->{{3.8175543100594673`*^9, 3.817554403402668*^9}, { + 3.817555144679367*^9, 3.817555170954733*^9}, 3.817557025045121*^9, { + 3.817557353599874*^9, 3.81755735592977*^9}}, + CellLabel->"Out[6]=",ExpressionUUID->"2b3d333d-6106-45b7-8d31-f398af6dca84"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ + RowBox[{"Limit", "[", RowBox[{ - RowBox[{"-", - FractionBox["1", "\[Pi]"]}], - FractionBox["1", "y"], - SuperscriptBox["\[ExponentialE]", - FractionBox["1", "y"]], " ", - RowBox[{"(", + RowBox[{ RowBox[{ - RowBox[{"-", "1"}], "+", "y"}], ")"}], " ", - RowBox[{"ExpIntegralEi", "[", - RowBox[{"-", - FractionBox["1", "y"]}], "]"}]}]}]], "Input", - CellChangeTimes->{{3.81712187820217*^9, 3.817121900537793*^9}, - 3.817122299113345*^9, {3.817198483403965*^9, 3.817198509154025*^9}, { - 3.817199440148189*^9, 3.817199440731868*^9}, {3.817199845562872*^9, - 3.817199848402961*^9}, 3.8172001271039267`*^9, {3.8172001950013113`*^9, - 3.8172001956329947`*^9}, {3.817209116226256*^9, 3.817209123026651*^9}, { - 3.817209186739526*^9, 3.817209188013309*^9}, {3.817210278358514*^9, - 3.817210288080538*^9}, {3.817269682319373*^9, 3.817269725334901*^9}, { - 3.817269813956863*^9, 3.817269814867185*^9}, {3.8172765008393383`*^9, - 3.817276503991358*^9}, {3.8172765584000072`*^9, 3.8172765730241413`*^9}, { - 3.817279978610174*^9, 3.8172799852438726`*^9}}, - CellLabel-> - "(WOPR) In[580]:=",ExpressionUUID->"f27ba651-c22d-460a-a64d-186282c0ce66"], + RowBox[{"R\[ScriptCapitalM]", "[", "2", "]"}], "[", + RowBox[{"\[Gamma]", ",", "B", ",", "\[Theta]c", ",", "M0"}], "]"}], "[", + "\[Theta]", "]"}], ",", + RowBox[{"\[Theta]", "\[Rule]", "\[Infinity]"}]}], "]"}]], "Input", + CellChangeTimes->{{3.817557366814589*^9, 3.817557383669715*^9}}, + CellLabel->"In[7]:=",ExpressionUUID->"377f0dcd-6050-40b3-be70-41379e9fda93"], -Cell[BoxData[{ +Cell[BoxData["0"], "Output", + CellChangeTimes->{3.8175573858936357`*^9}, + CellLabel->"Out[7]=",ExpressionUUID->"02602b05-0f8e-4f3a-8e08-cb1c3aa44a3f"] +}, Open ]], + +Cell[CellGroupData[{ + +Cell[BoxData[ RowBox[{ - RowBox[{"R\[ScriptCapitalM]", "[", "\[Theta]_", "]"}], ":=", RowBox[{ - RowBox[{"\[ScriptCapitalM]0", - RowBox[{"(", - RowBox[{ - RowBox[{"\[ScriptCapitalM]f", "[", - RowBox[{"B", - RowBox[{"(", - RowBox[{"\[Theta]c", "-", "\[Theta]"}], ")"}]}], "]"}], "-", - RowBox[{"\[ScriptCapitalM]f", "[", - RowBox[{"B", - RowBox[{"(", - RowBox[{"\[Theta]c", "+", "\[Theta]"}], ")"}]}], "]"}]}], ")"}]}], - "+", - RowBox[{"\[ScriptCapitalM]1", " ", "\[Theta]"}]}]}], "\[IndentingNewLine]", + RowBox[{"eqLow", "[", "2", "]"}], "[", + RowBox[{ + RowBox[{ + RowBox[{"f", "[", "\[Theta]i", "]"}], "[", "0", "]"}], ",", + RowBox[{ + RowBox[{"g", "[", + RowBox[{"gC", ",", "\[Theta]c"}], "]"}], "[", "0", "]"}]}], "]"}], "[", + "0", "]"}]], "Input", + CellChangeTimes->{{3.8175547641765738`*^9, 3.817554807184411*^9}, { + 3.81755721589187*^9, 3.817557223395288*^9}}, + CellLabel->"In[4]:=",ExpressionUUID->"8c3eb7fb-65d2-4a56-a3e3-0dc6e0518fe6"], + +Cell[BoxData[ RowBox[{ - RowBox[{"I\[ScriptCapitalF]", "[", "\[Theta]_", "]"}], ":=", - RowBox[{"\[ScriptCapitalF]1", + RowBox[{"M0", " ", RowBox[{"(", RowBox[{ - RowBox[{"HeavisideTheta", "[", - RowBox[{"\[Theta]", "-", "\[Theta]c"}], "]"}], "+", - RowBox[{"HeavisideTheta", "[", + FractionBox["1", "\[Pi]"], "+", + FractionBox[ RowBox[{ - RowBox[{"-", "\[Theta]"}], "-", "\[Theta]c"}], "]"}]}], ")"}], "B", - " ", + SuperscriptBox["\[ExponentialE]", + FractionBox["1", + RowBox[{"2", " ", "B", " ", "\[Theta]c"}]]], " ", + RowBox[{"(", + RowBox[{"1", "-", + RowBox[{"2", " ", "B", " ", "\[Theta]c"}], "-", + RowBox[{"2", " ", "B", " ", "\[Gamma]", " ", "\[Theta]c"}]}], ")"}], + " ", + RowBox[{"ExpIntegralEi", "[", + RowBox[{"-", + FractionBox["1", + RowBox[{"2", " ", "B", " ", "\[Theta]c"}]]}], "]"}]}], + RowBox[{"2", " ", "B", " ", "\[Pi]", " ", "\[Theta]c"}]]}], ")"}]}], + "+", + FractionBox[ RowBox[{ - RowBox[{"RealAbs", "[", - RowBox[{"\[Theta]", "-", "\[Theta]c"}], "]"}], "/", - RowBox[{"(", - RowBox[{"1", "+", - RowBox[{"\[Gamma]", " ", "B", " ", - RowBox[{"RealAbs", "[", - RowBox[{"\[Theta]", "-", "\[Theta]c"}], "]"}]}]}], ")"}]}], - RowBox[{"Exp", "[", - RowBox[{"-", - SuperscriptBox[ - RowBox[{"(", - RowBox[{"B", " ", - RowBox[{"RealAbs", "[", - RowBox[{"\[Theta]", "-", "\[Theta]c"}], "]"}]}], ")"}], - RowBox[{"-", "1"}]]}], "]"}]}]}]}], "Input", - CellChangeTimes->{{3.817121798544737*^9, 3.817121811575968*^9}, - 3.817121863281023*^9, {3.817121906833729*^9, 3.81712192235432*^9}, { - 3.817121959355006*^9, 3.8171219616748743`*^9}, {3.8171228987400923`*^9, - 3.8171229707415123`*^9}, {3.817198243983283*^9, 3.817198246927249*^9}, { - 3.817198496051778*^9, 3.817198523876305*^9}, {3.817198579461516*^9, - 3.817198579533231*^9}, {3.817198943419979*^9, 3.8171989466357327`*^9}, { - 3.8171993746668243`*^9, 3.817199374738553*^9}, {3.817200189041057*^9, - 3.817200208305159*^9}, {3.817209205684279*^9, 3.817209388039248*^9}, { - 3.81720944395175*^9, 3.817209446047728*^9}, {3.817209930464348*^9, - 3.8172099380644417`*^9}, {3.817210264070428*^9, 3.817210267542055*^9}, { - 3.817210329919582*^9, 3.817210334487361*^9}, {3.817210372639934*^9, - 3.817210377592119*^9}, 3.817212316154291*^9, 3.817212645208345*^9, { - 3.817216916290226*^9, 3.8172169167762327`*^9}, {3.817219212992364*^9, - 3.8172192132224283`*^9}, {3.817266621409737*^9, 3.8172666282959213`*^9}, { - 3.8172683850976753`*^9, 3.817268385343678*^9}, 3.817269739852894*^9, { - 3.81726981683014*^9, 3.817269821291363*^9}, {3.817271032349269*^9, - 3.8172710326758747`*^9}, {3.817271110024907*^9, 3.8172711106148663`*^9}, { - 3.8172719744968967`*^9, 3.817271998377635*^9}, {3.817272634473014*^9, - 3.81727265440335*^9}, {3.817276578733289*^9, 3.81727661335675*^9}, { - 3.817280005679316*^9, 3.817280009539542*^9}}, - CellLabel-> - "(WOPR) In[581]:=",ExpressionUUID->"e32b7309-d4a7-469a-ab41-1800db7276a7"], + SuperscriptBox[ + RowBox[{"(", + RowBox[{ + SuperscriptBox["\[Theta]c", "2"], "-", + SuperscriptBox["\[Theta]i", "2"]}], ")"}], + RowBox[{"1", "/", "8"}]], " ", + RowBox[{ + SuperscriptBox["Gl", "\[Prime]", + MultilineFunction->None], "[", "0", "]"}]}], + SuperscriptBox["\[Theta]i", + RowBox[{"1", "/", "4"}]]]}]], "Output", + CellChangeTimes->{ + 3.817554807406887*^9, {3.817557206327607*^9, 3.81755724215594*^9}, + 3.8175573305621147`*^9}, + CellLabel->"Out[4]=",ExpressionUUID->"e44b7486-110d-476c-b94b-28e4b600a40c"] +}, Open ]], Cell[BoxData[{ RowBox[{ @@ -11377,10 +11418,9 @@ Cell[BoxData[ }, Open ]] }, Open ]] }, -Evaluator->"WOPR", -WindowSize->{1437, 792}, -WindowMargins->{{1.5, Automatic}, {1.5, Automatic}}, -FrontEndVersion->"12.1 for Linux ARM (32-bit) (August 18, 2020)", +WindowSize->{957, 1062}, +WindowMargins->{{Automatic, 1.5}, {1.5, Automatic}}, +FrontEndVersion->"12.1 for Linux x86 (64-bit) (June 19, 2020)", StyleDefinitions->"Default.nb", ExpressionUUID->"6f6c570c-dfc0-4320-bdc7-0f1de62bb24c" ] @@ -11395,220 +11435,228 @@ CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ +Cell[558, 20, 270, 5, 24, "Input",ExpressionUUID->"e595a667-7f25-4f98-9f12-a321a236d113"], +Cell[CellGroupData[{ +Cell[853, 29, 156, 3, 50, "Section",ExpressionUUID->"bb531fce-8a9d-4a5f-8d6a-ec66b9fbf698"], +Cell[1012, 34, 1380, 28, 24, "Input",ExpressionUUID->"395191f4-75f5-4aa3-87f6-ccf9cbedf73f"], Cell[CellGroupData[{ -Cell[580, 22, 156, 3, 68, "Section",ExpressionUUID->"bb531fce-8a9d-4a5f-8d6a-ec66b9fbf698"], -Cell[739, 27, 613, 13, 42, "Input",ExpressionUUID->"fafcf491-a81f-4f91-b9fe-4d27770e7577"], -Cell[1355, 42, 1380, 28, 42, "Input",ExpressionUUID->"395191f4-75f5-4aa3-87f6-ccf9cbedf73f"], -Cell[2738, 72, 838, 17, 45, "Input",ExpressionUUID->"3e69b8b0-cd07-4524-a460-17b8352b5902"], -Cell[3579, 91, 820, 23, 58, "Input",ExpressionUUID->"daded195-b1a7-448b-87f6-d258f2dcf5e8"], -Cell[4402, 116, 1225, 26, 62, "Input",ExpressionUUID->"f27ba651-c22d-460a-a64d-186282c0ce66"], -Cell[5630, 144, 3022, 65, 66, "Input",ExpressionUUID->"e32b7309-d4a7-469a-ab41-1800db7276a7"], -Cell[8655, 211, 8663, 212, 436, "Input",ExpressionUUID->"77a88da0-5b31-458a-9ed8-0b34545e511d"], -Cell[17321, 425, 4921, 136, 171, "Input",ExpressionUUID->"1f309ab1-20a0-46cc-93fe-35aed54ccb3d"], -Cell[22245, 563, 233, 5, 42, "Input",ExpressionUUID->"e7fef682-ee6e-4d26-962a-6946f28d5999"], -Cell[22481, 570, 2312, 53, 64, "Input",ExpressionUUID->"2f5b927d-f164-4051-b8ca-bad9a5ba3b97"], +Cell[2417, 66, 2007, 47, 92, "Input",ExpressionUUID->"389cbfb1-8cc7-458c-9f8c-1fd631f97f2c"], +Cell[4427, 115, 2664, 57, 288, "Output",ExpressionUUID->"2b3d333d-6106-45b7-8d31-f398af6dca84"] +}, Open ]], +Cell[CellGroupData[{ +Cell[7128, 177, 453, 10, 22, "Input",ExpressionUUID->"377f0dcd-6050-40b3-be70-41379e9fda93"], +Cell[7584, 189, 150, 2, 25, "Output",ExpressionUUID->"02602b05-0f8e-4f3a-8e08-cb1c3aa44a3f"] +}, Open ]], +Cell[CellGroupData[{ +Cell[7771, 196, 502, 13, 24, "Input",ExpressionUUID->"8c3eb7fb-65d2-4a56-a3e3-0dc6e0518fe6"], +Cell[8276, 211, 1300, 38, 62, "Output",ExpressionUUID->"e44b7486-110d-476c-b94b-28e4b600a40c"] +}, Open ]], +Cell[9591, 252, 8663, 212, 332, "Input",ExpressionUUID->"77a88da0-5b31-458a-9ed8-0b34545e511d"], +Cell[18257, 466, 4921, 136, 144, "Input",ExpressionUUID->"1f309ab1-20a0-46cc-93fe-35aed54ccb3d"], +Cell[23181, 604, 233, 5, 22, "Input",ExpressionUUID->"e7fef682-ee6e-4d26-962a-6946f28d5999"], +Cell[23417, 611, 2312, 53, 56, "Input",ExpressionUUID->"2f5b927d-f164-4051-b8ca-bad9a5ba3b97"], Cell[CellGroupData[{ -Cell[24818, 627, 4764, 78, 42, "Input",ExpressionUUID->"630fc9c9-229c-4372-8a16-eac1f2d98690"], -Cell[29585, 707, 632, 12, 38, "Message",ExpressionUUID->"20e41565-1936-405a-9504-d530fdcc8152"], -Cell[30220, 721, 1499, 29, 49, "Output",ExpressionUUID->"319eb113-9321-48fc-ad65-5700453ef613"] +Cell[25754, 668, 4764, 78, 41, "Input",ExpressionUUID->"630fc9c9-229c-4372-8a16-eac1f2d98690"], +Cell[30521, 748, 632, 12, 31, "Message",ExpressionUUID->"20e41565-1936-405a-9504-d530fdcc8152"], +Cell[31156, 762, 1499, 29, 44, "Output",ExpressionUUID->"319eb113-9321-48fc-ad65-5700453ef613"] }, Open ]], -Cell[31734, 753, 461, 14, 42, "Input",ExpressionUUID->"ea36fa2d-8c46-4383-be61-eff682dedb14"], +Cell[32670, 794, 461, 14, 22, "Input",ExpressionUUID->"ea36fa2d-8c46-4383-be61-eff682dedb14"], Cell[CellGroupData[{ -Cell[32220, 771, 452, 12, 42, "Input",ExpressionUUID->"5823df11-0739-494d-bf1b-44bf2a14cb2b"], -Cell[32675, 785, 278, 4, 46, "Output",ExpressionUUID->"cf69a3b6-a533-4dea-b33c-8d65e371fe08"] +Cell[33156, 812, 452, 12, 22, "Input",ExpressionUUID->"5823df11-0739-494d-bf1b-44bf2a14cb2b"], +Cell[33611, 826, 278, 4, 25, "Output",ExpressionUUID->"cf69a3b6-a533-4dea-b33c-8d65e371fe08"] }, Open ]], Cell[CellGroupData[{ -Cell[32990, 794, 1740, 46, 105, "Input",ExpressionUUID->"2a7b9f00-b942-4f1e-89df-1d3081b1fc6a"], -Cell[34733, 842, 21020, 380, 218, "Output",ExpressionUUID->"155c39dd-ec1d-4900-9cbf-4fc4b876ac75"] +Cell[33926, 835, 1740, 46, 71, "Input",ExpressionUUID->"2a7b9f00-b942-4f1e-89df-1d3081b1fc6a"], +Cell[35669, 883, 21020, 380, 153, "Output",ExpressionUUID->"155c39dd-ec1d-4900-9cbf-4fc4b876ac75"] }, Open ]], Cell[CellGroupData[{ -Cell[55790, 1227, 966, 27, 45, "Input",ExpressionUUID->"da645597-b6a8-4dfc-83c5-3cec31d750f9"], -Cell[56759, 1256, 469, 10, 46, "Output",ExpressionUUID->"f27b02e6-f3db-42e2-b534-43e5409fe36c"] +Cell[56726, 1268, 966, 27, 24, "Input",ExpressionUUID->"da645597-b6a8-4dfc-83c5-3cec31d750f9"], +Cell[57695, 1297, 469, 10, 25, "Output",ExpressionUUID->"f27b02e6-f3db-42e2-b534-43e5409fe36c"] }, Open ]], Cell[CellGroupData[{ -Cell[57265, 1271, 1847, 47, 109, "Input",ExpressionUUID->"b03f723e-24d4-4097-8ca2-9b1145403bfb"], -Cell[59115, 1320, 519, 11, 22, "Message",ExpressionUUID->"1ff6b3ff-f0d9-48ea-95df-10584520bc02"], -Cell[59637, 1333, 519, 11, 22, "Message",ExpressionUUID->"82682d18-f80c-4d23-8675-d31ca1576f0e"], -Cell[60159, 1346, 521, 11, 22, "Message",ExpressionUUID->"164319c1-a2e0-427a-9f70-94e763f48d8e"], -Cell[60683, 1359, 486, 10, 22, "Message",ExpressionUUID->"5d8fc41c-1820-4591-bdb1-29d65ed68611"], -Cell[61172, 1371, 25289, 433, 243, "Output",ExpressionUUID->"0f099d03-656a-4602-85f6-26b4b2405b3b"] +Cell[58201, 1312, 1847, 47, 73, "Input",ExpressionUUID->"b03f723e-24d4-4097-8ca2-9b1145403bfb"], +Cell[60051, 1361, 519, 11, 19, "Message",ExpressionUUID->"1ff6b3ff-f0d9-48ea-95df-10584520bc02"], +Cell[60573, 1374, 519, 11, 19, "Message",ExpressionUUID->"82682d18-f80c-4d23-8675-d31ca1576f0e"], +Cell[61095, 1387, 521, 11, 19, "Message",ExpressionUUID->"164319c1-a2e0-427a-9f70-94e763f48d8e"], +Cell[61619, 1400, 486, 10, 19, "Message",ExpressionUUID->"5d8fc41c-1820-4591-bdb1-29d65ed68611"], +Cell[62108, 1412, 25289, 433, 172, "Output",ExpressionUUID->"0f099d03-656a-4602-85f6-26b4b2405b3b"] }, Open ]], -Cell[86476, 1807, 1836, 43, 85, "Input",ExpressionUUID->"868d07c1-b3cf-44ff-a070-78c2a21c388b"], +Cell[87412, 1848, 1836, 43, 57, "Input",ExpressionUUID->"868d07c1-b3cf-44ff-a070-78c2a21c388b"], Cell[CellGroupData[{ -Cell[88337, 1854, 4656, 86, 64, "Input",ExpressionUUID->"7c02fde8-7f1a-472a-9e7d-0e95772516eb"], -Cell[92996, 1942, 899, 22, 46, "Output",ExpressionUUID->"6b22737f-c85f-4c36-9b44-d9e694bee115"] +Cell[89273, 1895, 4656, 86, 73, "Input",ExpressionUUID->"7c02fde8-7f1a-472a-9e7d-0e95772516eb"], +Cell[93932, 1983, 899, 22, 25, "Output",ExpressionUUID->"6b22737f-c85f-4c36-9b44-d9e694bee115"] }, Open ]], Cell[CellGroupData[{ -Cell[93932, 1969, 1989, 49, 105, "Input",ExpressionUUID->"0e89b49c-ebf1-4e9a-800f-62ec63f935a1"], -Cell[95924, 2020, 21952, 393, 218, "Output",ExpressionUUID->"7cf0a89e-6d18-442b-b56d-94b97702709f"] +Cell[94868, 2010, 1989, 49, 75, "Input",ExpressionUUID->"0e89b49c-ebf1-4e9a-800f-62ec63f935a1"], +Cell[96860, 2061, 21952, 393, 154, "Output",ExpressionUUID->"7cf0a89e-6d18-442b-b56d-94b97702709f"] }, Open ]], Cell[CellGroupData[{ -Cell[117913, 2418, 1382, 36, 45, "Input",ExpressionUUID->"6c30dc25-997d-4790-bfcd-6e15ebcf7cbc"], -Cell[119298, 2456, 520, 12, 40, "Message",ExpressionUUID->"56028de4-0ddc-4ea7-bf21-a25534c5c68d"], -Cell[119821, 2470, 521, 12, 40, "Message",ExpressionUUID->"3d57804e-381a-4ba2-8dce-4428206df060"], -Cell[120345, 2484, 3161, 69, 259, "Output",ExpressionUUID->"56916602-5965-4962-8188-b6662a5aa173"] +Cell[118849, 2459, 1382, 36, 24, "Input",ExpressionUUID->"6c30dc25-997d-4790-bfcd-6e15ebcf7cbc"], +Cell[120234, 2497, 520, 12, 33, "Message",ExpressionUUID->"56028de4-0ddc-4ea7-bf21-a25534c5c68d"], +Cell[120757, 2511, 521, 12, 40, "Message",ExpressionUUID->"3d57804e-381a-4ba2-8dce-4428206df060"], +Cell[121281, 2525, 3161, 69, 259, "Output",ExpressionUUID->"56916602-5965-4962-8188-b6662a5aa173"] }, Open ]], -Cell[123521, 2556, 1713, 40, 64, "Input",ExpressionUUID->"b31259b6-cc89-48a8-920d-e9dbcee664fb"], +Cell[124457, 2597, 1713, 40, 64, "Input",ExpressionUUID->"b31259b6-cc89-48a8-920d-e9dbcee664fb"], Cell[CellGroupData[{ -Cell[125259, 2600, 3913, 66, 42, "Input",ExpressionUUID->"30d0f855-c86b-4121-855c-813e9ee04f96"], -Cell[129175, 2668, 1064, 19, 28, "Message",ExpressionUUID->"f1cd282b-fb73-4166-8899-1f1b71289512"], -Cell[130242, 2689, 3809, 57, 49, "Output",ExpressionUUID->"3b7bac7f-12db-4bb2-8d83-0f02d8f19997"] +Cell[126195, 2641, 3913, 66, 42, "Input",ExpressionUUID->"30d0f855-c86b-4121-855c-813e9ee04f96"], +Cell[130111, 2709, 1064, 19, 28, "Message",ExpressionUUID->"f1cd282b-fb73-4166-8899-1f1b71289512"], +Cell[131178, 2730, 3809, 57, 49, "Output",ExpressionUUID->"3b7bac7f-12db-4bb2-8d83-0f02d8f19997"] }, Open ]], Cell[CellGroupData[{ -Cell[134088, 2751, 1587, 43, 105, "Input",ExpressionUUID->"e1346f99-b8a8-49f2-9bab-d7d94b52e5d3"], -Cell[135678, 2796, 61160, 1040, 218, 39593, 685, "CachedBoxData", "BoxData", "Output",ExpressionUUID->"24eee5ea-c8f9-4063-9d36-9e38f4d71b1d"] +Cell[135024, 2792, 1587, 43, 105, "Input",ExpressionUUID->"e1346f99-b8a8-49f2-9bab-d7d94b52e5d3"], +Cell[136614, 2837, 61160, 1040, 218, 39593, 685, "CachedBoxData", "BoxData", "Output",ExpressionUUID->"24eee5ea-c8f9-4063-9d36-9e38f4d71b1d"] }, Open ]], Cell[CellGroupData[{ -Cell[196875, 3841, 1714, 46, 105, "Input",ExpressionUUID->"8b99e721-7bc5-4c63-b150-8946bb7d2007"], -Cell[198592, 3889, 30207, 531, 218, "Output",ExpressionUUID->"d79294e8-cb2f-4112-ad68-da965b8a242a"] +Cell[197811, 3882, 1714, 46, 105, "Input",ExpressionUUID->"8b99e721-7bc5-4c63-b150-8946bb7d2007"], +Cell[199528, 3930, 30207, 531, 218, "Output",ExpressionUUID->"d79294e8-cb2f-4112-ad68-da965b8a242a"] }, Open ]], Cell[CellGroupData[{ -Cell[228836, 4425, 491, 13, 29, "Input",ExpressionUUID->"8dd8c18b-504b-469d-8072-a9fb0cc417df"], -Cell[229330, 4440, 539, 12, 22, "Message",ExpressionUUID->"a55da929-28cd-4032-aea8-5daec85d020a"], -Cell[229872, 4454, 539, 12, 22, "Message",ExpressionUUID->"d9f18df2-5058-4626-9d48-5486b1743e96"], -Cell[230414, 4468, 2243, 56, 79, "Output",ExpressionUUID->"7ce119a6-a1c0-48cc-9d86-6e7b6a1eafe3"] +Cell[229772, 4466, 491, 13, 29, "Input",ExpressionUUID->"8dd8c18b-504b-469d-8072-a9fb0cc417df"], +Cell[230266, 4481, 539, 12, 22, "Message",ExpressionUUID->"a55da929-28cd-4032-aea8-5daec85d020a"], +Cell[230808, 4495, 539, 12, 22, "Message",ExpressionUUID->"d9f18df2-5058-4626-9d48-5486b1743e96"], +Cell[231350, 4509, 2243, 56, 79, "Output",ExpressionUUID->"7ce119a6-a1c0-48cc-9d86-6e7b6a1eafe3"] }, Open ]], Cell[CellGroupData[{ -Cell[232694, 4529, 561, 16, 29, "Input",ExpressionUUID->"b0a27cf9-de55-4cf4-9b19-ff7b3add442a"], -Cell[233258, 4547, 913, 27, 75, "Output",ExpressionUUID->"d3fcf7af-82f3-49c4-81c0-3bf29938e0d1"] +Cell[233630, 4570, 561, 16, 29, "Input",ExpressionUUID->"b0a27cf9-de55-4cf4-9b19-ff7b3add442a"], +Cell[234194, 4588, 913, 27, 75, "Output",ExpressionUUID->"d3fcf7af-82f3-49c4-81c0-3bf29938e0d1"] }, Open ]], Cell[CellGroupData[{ -Cell[234208, 4579, 681, 19, 55, "Input",ExpressionUUID->"f49c10e4-6ea3-40fb-85e8-436bcc73196f"], -Cell[234892, 4600, 4071, 86, 235, "Output",ExpressionUUID->"0ede64a7-b62f-453c-b5e6-06e4c6ab27af"] +Cell[235144, 4620, 681, 19, 55, "Input",ExpressionUUID->"f49c10e4-6ea3-40fb-85e8-436bcc73196f"], +Cell[235828, 4641, 4071, 86, 235, "Output",ExpressionUUID->"0ede64a7-b62f-453c-b5e6-06e4c6ab27af"] }, Open ]], Cell[CellGroupData[{ -Cell[239000, 4691, 486, 13, 29, "Input",ExpressionUUID->"67bd5207-9130-427b-afb3-21bda9c7e82e"], -Cell[239489, 4706, 200, 4, 33, "Output",ExpressionUUID->"39e5dba9-9c25-42e4-8876-75a530e46ee5"] +Cell[239936, 4732, 486, 13, 29, "Input",ExpressionUUID->"67bd5207-9130-427b-afb3-21bda9c7e82e"], +Cell[240425, 4747, 200, 4, 33, "Output",ExpressionUUID->"39e5dba9-9c25-42e4-8876-75a530e46ee5"] }, Open ]], Cell[CellGroupData[{ -Cell[239726, 4715, 783, 22, 29, "Input",ExpressionUUID->"6db9bdd5-3bc4-4993-85d9-2bb9c2d0db38"], -Cell[240512, 4739, 3496, 77, 249, "Output",ExpressionUUID->"028520eb-aee9-483d-9152-2416a8a91977"] +Cell[240662, 4756, 783, 22, 29, "Input",ExpressionUUID->"6db9bdd5-3bc4-4993-85d9-2bb9c2d0db38"], +Cell[241448, 4780, 3496, 77, 249, "Output",ExpressionUUID->"028520eb-aee9-483d-9152-2416a8a91977"] }, Open ]], Cell[CellGroupData[{ -Cell[244045, 4821, 1223, 33, 24, "Input",ExpressionUUID->"aab225ac-8e68-4d46-a4da-600fe6364133"], -Cell[245271, 4856, 10880, 231, 176, "Output",ExpressionUUID->"6e880d59-ffe1-440b-af8b-754c17269f6e"] +Cell[244981, 4862, 1223, 33, 24, "Input",ExpressionUUID->"aab225ac-8e68-4d46-a4da-600fe6364133"], +Cell[246207, 4897, 10880, 231, 176, "Output",ExpressionUUID->"6e880d59-ffe1-440b-af8b-754c17269f6e"] }, Open ]], Cell[CellGroupData[{ -Cell[256188, 5092, 2115, 52, 66, "Input",ExpressionUUID->"f01a72f1-91c7-4454-9040-b54a61c9cd07"], -Cell[258306, 5146, 615, 13, 19, "Message",ExpressionUUID->"46cbf58a-9f1b-42bc-8a11-616c84d0e1d7"], -Cell[258924, 5161, 613, 13, 19, "Message",ExpressionUUID->"c25c6435-fc6d-4b46-b41c-a3d203aa09da"], -Cell[259540, 5176, 613, 13, 19, "Message",ExpressionUUID->"3a09025b-666f-454b-9a32-9f324e354788"], -Cell[260156, 5191, 580, 12, 19, "Message",ExpressionUUID->"5b1f7808-f0a6-4257-b7ed-4acdc8c4ac95"] +Cell[257124, 5133, 2115, 52, 66, "Input",ExpressionUUID->"f01a72f1-91c7-4454-9040-b54a61c9cd07"], +Cell[259242, 5187, 615, 13, 19, "Message",ExpressionUUID->"46cbf58a-9f1b-42bc-8a11-616c84d0e1d7"], +Cell[259860, 5202, 613, 13, 19, "Message",ExpressionUUID->"c25c6435-fc6d-4b46-b41c-a3d203aa09da"], +Cell[260476, 5217, 613, 13, 19, "Message",ExpressionUUID->"3a09025b-666f-454b-9a32-9f324e354788"], +Cell[261092, 5232, 580, 12, 19, "Message",ExpressionUUID->"5b1f7808-f0a6-4257-b7ed-4acdc8c4ac95"] }, Open ]], Cell[CellGroupData[{ -Cell[260773, 5208, 536, 10, 24, "Input",ExpressionUUID->"1dfbee5f-7e33-4446-bc10-991ba8838497"], -Cell[261312, 5220, 93071, 1628, 177, "Output",ExpressionUUID->"e57fa90a-27ff-4013-8fad-ca7e3c98d348"] +Cell[261709, 5249, 536, 10, 24, "Input",ExpressionUUID->"1dfbee5f-7e33-4446-bc10-991ba8838497"], +Cell[262248, 5261, 93071, 1628, 177, "Output",ExpressionUUID->"e57fa90a-27ff-4013-8fad-ca7e3c98d348"] }, Open ]], Cell[CellGroupData[{ -Cell[354420, 6853, 314, 8, 22, "Input",ExpressionUUID->"8499af31-a50c-4059-afac-b628095fc761"], -Cell[354737, 6863, 368, 8, 25, "Output",ExpressionUUID->"9aaf347d-0fae-4cb5-8920-f15fee437a9e"] +Cell[355356, 6894, 314, 8, 22, "Input",ExpressionUUID->"8499af31-a50c-4059-afac-b628095fc761"], +Cell[355673, 6904, 368, 8, 25, "Output",ExpressionUUID->"9aaf347d-0fae-4cb5-8920-f15fee437a9e"] }, Open ]], Cell[CellGroupData[{ -Cell[355142, 6876, 577, 17, 24, "Input",ExpressionUUID->"fecf4a67-c7e8-40b5-8e36-582a0bab9fe6"], -Cell[355722, 6895, 753, 23, 48, "Output",ExpressionUUID->"742a0341-4af8-4ba4-a799-e15a1e847a9c"] +Cell[356078, 6917, 577, 17, 24, "Input",ExpressionUUID->"fecf4a67-c7e8-40b5-8e36-582a0bab9fe6"], +Cell[356658, 6936, 753, 23, 48, "Output",ExpressionUUID->"742a0341-4af8-4ba4-a799-e15a1e847a9c"] }, Open ]], Cell[CellGroupData[{ -Cell[356512, 6923, 438, 11, 24, "Input",ExpressionUUID->"9884f3ff-ff57-4cf7-b122-2aebcb709801"], -Cell[356953, 6936, 947, 26, 45, "Output",ExpressionUUID->"5f69419d-b0c1-4cd5-9616-460f429c823a"] +Cell[357448, 6964, 438, 11, 24, "Input",ExpressionUUID->"9884f3ff-ff57-4cf7-b122-2aebcb709801"], +Cell[357889, 6977, 947, 26, 45, "Output",ExpressionUUID->"5f69419d-b0c1-4cd5-9616-460f429c823a"] }, Open ]], Cell[CellGroupData[{ -Cell[357937, 6967, 549, 13, 24, "Input",ExpressionUUID->"c18d621f-7e5f-48db-8ae7-0fda0221c79c"], -Cell[358489, 6982, 599, 17, 45, "Output",ExpressionUUID->"b45884b0-ccdc-4e99-ac15-83498ea86eff"] +Cell[358873, 7008, 549, 13, 24, "Input",ExpressionUUID->"c18d621f-7e5f-48db-8ae7-0fda0221c79c"], +Cell[359425, 7023, 599, 17, 45, "Output",ExpressionUUID->"b45884b0-ccdc-4e99-ac15-83498ea86eff"] }, Open ]], Cell[CellGroupData[{ -Cell[359125, 7004, 338, 7, 22, "Input",ExpressionUUID->"8e16b08d-33f5-4a51-9850-75f98a0e6a08"], -Cell[359466, 7013, 2269, 65, 61, "Output",ExpressionUUID->"21b799e5-0965-4692-b4ba-6083041c0cdd"] +Cell[360061, 7045, 338, 7, 22, "Input",ExpressionUUID->"8e16b08d-33f5-4a51-9850-75f98a0e6a08"], +Cell[360402, 7054, 2269, 65, 61, "Output",ExpressionUUID->"21b799e5-0965-4692-b4ba-6083041c0cdd"] }, Open ]], Cell[CellGroupData[{ -Cell[361772, 7083, 237, 5, 22, "Input",ExpressionUUID->"4c8cb259-b56a-423e-957f-9a65ba700a16"], -Cell[362012, 7090, 2063, 59, 109, "Output",ExpressionUUID->"312f09f3-55f8-4f9d-a42c-05ca1c4aea27"] +Cell[362708, 7124, 237, 5, 22, "Input",ExpressionUUID->"4c8cb259-b56a-423e-957f-9a65ba700a16"], +Cell[362948, 7131, 2063, 59, 109, "Output",ExpressionUUID->"312f09f3-55f8-4f9d-a42c-05ca1c4aea27"] }, Open ]], Cell[CellGroupData[{ -Cell[364112, 7154, 564, 14, 24, "Input",ExpressionUUID->"069edadd-116a-4ccf-831d-0fe6e07540fe"], -Cell[364679, 7170, 1113, 29, 52, "Output",ExpressionUUID->"e6426e3c-3d4d-49f1-becd-dfef45704088"] +Cell[365048, 7195, 564, 14, 24, "Input",ExpressionUUID->"069edadd-116a-4ccf-831d-0fe6e07540fe"], +Cell[365615, 7211, 1113, 29, 52, "Output",ExpressionUUID->"e6426e3c-3d4d-49f1-becd-dfef45704088"] }, Open ]], Cell[CellGroupData[{ -Cell[365829, 7204, 522, 13, 24, "Input",ExpressionUUID->"3cd180b6-0e22-4832-af2b-e9ecec7e0a49"], -Cell[366354, 7219, 5586, 138, 258, "Output",ExpressionUUID->"45f19748-f63a-49e1-8ced-a4f5c5dda272"] +Cell[366765, 7245, 522, 13, 24, "Input",ExpressionUUID->"3cd180b6-0e22-4832-af2b-e9ecec7e0a49"], +Cell[367290, 7260, 5586, 138, 258, "Output",ExpressionUUID->"45f19748-f63a-49e1-8ced-a4f5c5dda272"] }, Open ]], -Cell[371955, 7360, 366, 10, 22, "Input",ExpressionUUID->"3ffc56f5-fc40-4c69-bdce-f91a4aab1c67"], +Cell[372891, 7401, 366, 10, 22, "Input",ExpressionUUID->"3ffc56f5-fc40-4c69-bdce-f91a4aab1c67"], Cell[CellGroupData[{ -Cell[372346, 7374, 838, 23, 42, "Input",ExpressionUUID->"e131ab87-4ef6-4852-9ddd-e4499b39559a"], -Cell[373187, 7399, 597, 16, 51, "Output",ExpressionUUID->"86a22957-07a1-48ab-90d1-d142a27afe22"] +Cell[373282, 7415, 838, 23, 42, "Input",ExpressionUUID->"e131ab87-4ef6-4852-9ddd-e4499b39559a"], +Cell[374123, 7440, 597, 16, 51, "Output",ExpressionUUID->"86a22957-07a1-48ab-90d1-d142a27afe22"] }, Open ]], Cell[CellGroupData[{ -Cell[373821, 7420, 921, 23, 42, "Input",ExpressionUUID->"22a30e99-2bd5-43a5-9353-3a3be14232a6"], -Cell[374745, 7445, 2812, 74, 80, "Output",ExpressionUUID->"470e97c2-b944-4611-8c77-75e5a0915c04"] +Cell[374757, 7461, 921, 23, 42, "Input",ExpressionUUID->"22a30e99-2bd5-43a5-9353-3a3be14232a6"], +Cell[375681, 7486, 2812, 74, 80, "Output",ExpressionUUID->"470e97c2-b944-4611-8c77-75e5a0915c04"] }, Open ]], Cell[CellGroupData[{ -Cell[377594, 7524, 824, 23, 42, "Input",ExpressionUUID->"3752e4d8-bc7c-475c-b88c-a0d99559d99c"], -Cell[378421, 7549, 2580, 68, 78, "Output",ExpressionUUID->"db6b4f65-90c7-4c37-9945-5690d8472368"] +Cell[378530, 7565, 824, 23, 42, "Input",ExpressionUUID->"3752e4d8-bc7c-475c-b88c-a0d99559d99c"], +Cell[379357, 7590, 2580, 68, 78, "Output",ExpressionUUID->"db6b4f65-90c7-4c37-9945-5690d8472368"] }, Open ]], Cell[CellGroupData[{ -Cell[381038, 7622, 463, 12, 40, "Input",ExpressionUUID->"6a34969c-5bbf-4acb-a54f-a9bfbf2c9f87"], -Cell[381504, 7636, 8623, 160, 171, "Output",ExpressionUUID->"70cff46d-010d-4a9a-8cf5-3901b2a30edf"] +Cell[381974, 7663, 463, 12, 40, "Input",ExpressionUUID->"6a34969c-5bbf-4acb-a54f-a9bfbf2c9f87"], +Cell[382440, 7677, 8623, 160, 171, "Output",ExpressionUUID->"70cff46d-010d-4a9a-8cf5-3901b2a30edf"] }, Open ]], Cell[CellGroupData[{ -Cell[390164, 7801, 595, 17, 24, "Input",ExpressionUUID->"6afc45fb-65bc-47fe-b120-408f15ede40b"], -Cell[390762, 7820, 470, 13, 50, "Output",ExpressionUUID->"4cdf78f8-6821-4e67-9385-8a091cc9cf9d"] +Cell[391100, 7842, 595, 17, 24, "Input",ExpressionUUID->"6afc45fb-65bc-47fe-b120-408f15ede40b"], +Cell[391698, 7861, 470, 13, 50, "Output",ExpressionUUID->"4cdf78f8-6821-4e67-9385-8a091cc9cf9d"] }, Open ]], Cell[CellGroupData[{ -Cell[391269, 7838, 236, 6, 22, "Input",ExpressionUUID->"b66a0012-a40b-4041-8ea9-a9cad40db85f"], -Cell[391508, 7846, 301, 8, 39, "Output",ExpressionUUID->"9744135b-7d1f-4023-b02d-1bdf3b29d60f"] +Cell[392205, 7879, 236, 6, 22, "Input",ExpressionUUID->"b66a0012-a40b-4041-8ea9-a9cad40db85f"], +Cell[392444, 7887, 301, 8, 39, "Output",ExpressionUUID->"9744135b-7d1f-4023-b02d-1bdf3b29d60f"] }, Open ]], Cell[CellGroupData[{ -Cell[391846, 7859, 419, 12, 24, "Input",ExpressionUUID->"7ceb4bac-ce7b-4adc-aeb9-95a215777744"], -Cell[392268, 7873, 720, 22, 40, "Output",ExpressionUUID->"9c7df971-6bfc-41f8-a68e-fb9921c0b017"] +Cell[392782, 7900, 419, 12, 24, "Input",ExpressionUUID->"7ceb4bac-ce7b-4adc-aeb9-95a215777744"], +Cell[393204, 7914, 720, 22, 40, "Output",ExpressionUUID->"9c7df971-6bfc-41f8-a68e-fb9921c0b017"] }, Open ]], Cell[CellGroupData[{ -Cell[393025, 7900, 492, 13, 24, "Input",ExpressionUUID->"07932986-b2e1-477a-9d99-1b9701e3af0c"], -Cell[393520, 7915, 510, 11, 19, "Message",ExpressionUUID->"f0c4d2f9-9d6d-48c0-b78e-b8f84a7f6bb2"], -Cell[394033, 7928, 5864, 115, 182, "Output",ExpressionUUID->"82bf5b22-f54c-4da8-bea6-61b7c42c3450"] +Cell[393961, 7941, 492, 13, 24, "Input",ExpressionUUID->"07932986-b2e1-477a-9d99-1b9701e3af0c"], +Cell[394456, 7956, 510, 11, 19, "Message",ExpressionUUID->"f0c4d2f9-9d6d-48c0-b78e-b8f84a7f6bb2"], +Cell[394969, 7969, 5864, 115, 182, "Output",ExpressionUUID->"82bf5b22-f54c-4da8-bea6-61b7c42c3450"] }, Open ]], Cell[CellGroupData[{ -Cell[399934, 8048, 804, 22, 54, "Input",ExpressionUUID->"13e1541b-4ad9-4a40-89f2-be40a0e4d6fd"], -Cell[400741, 8072, 479, 13, 51, "Output",ExpressionUUID->"dff49a57-ddf5-4e9d-aa0a-55006221f3c2"] +Cell[400870, 8089, 804, 22, 54, "Input",ExpressionUUID->"13e1541b-4ad9-4a40-89f2-be40a0e4d6fd"], +Cell[401677, 8113, 479, 13, 51, "Output",ExpressionUUID->"dff49a57-ddf5-4e9d-aa0a-55006221f3c2"] }, Open ]], Cell[CellGroupData[{ -Cell[401257, 8090, 650, 19, 24, "Input",ExpressionUUID->"6def2df0-3761-4f37-b66f-ba9071b93625"], -Cell[401910, 8111, 534, 15, 42, "Output",ExpressionUUID->"175a2c7c-fa54-45c3-bf4f-c4ed89dad110"] +Cell[402193, 8131, 650, 19, 24, "Input",ExpressionUUID->"6def2df0-3761-4f37-b66f-ba9071b93625"], +Cell[402846, 8152, 534, 15, 42, "Output",ExpressionUUID->"175a2c7c-fa54-45c3-bf4f-c4ed89dad110"] }, Open ]], Cell[CellGroupData[{ -Cell[402481, 8131, 302, 7, 24, "Input",ExpressionUUID->"657e079b-66fb-4614-816e-af41b7fabf10"], -Cell[402786, 8140, 454, 11, 19, "Message",ExpressionUUID->"f6532f4e-d02f-402b-985e-f8964b91f90e"], -Cell[403243, 8153, 1129, 29, 32, "Output",ExpressionUUID->"5f2bc73c-3ef7-4874-86b3-dd3ac3e99858"] +Cell[403417, 8172, 302, 7, 24, "Input",ExpressionUUID->"657e079b-66fb-4614-816e-af41b7fabf10"], +Cell[403722, 8181, 454, 11, 19, "Message",ExpressionUUID->"f6532f4e-d02f-402b-985e-f8964b91f90e"], +Cell[404179, 8194, 1129, 29, 32, "Output",ExpressionUUID->"5f2bc73c-3ef7-4874-86b3-dd3ac3e99858"] }, Open ]], Cell[CellGroupData[{ -Cell[404409, 8187, 325, 7, 24, "Input",ExpressionUUID->"49f74264-7484-4d5a-b52f-d0892384e631"], -Cell[404737, 8196, 4933, 99, 175, "Output",ExpressionUUID->"26e2807e-1336-4751-8e74-d4414f081679"] +Cell[405345, 8228, 325, 7, 24, "Input",ExpressionUUID->"49f74264-7484-4d5a-b52f-d0892384e631"], +Cell[405673, 8237, 4933, 99, 175, "Output",ExpressionUUID->"26e2807e-1336-4751-8e74-d4414f081679"] }, Open ]], Cell[CellGroupData[{ -Cell[409707, 8300, 1059, 26, 48, "Input",ExpressionUUID->"3bd71288-3233-4171-a393-57251fb3dc52"], -Cell[410769, 8328, 6890, 131, 180, "Output",ExpressionUUID->"203873a1-0dd1-4567-8d12-97eeeed127f7"] +Cell[410643, 8341, 1059, 26, 48, "Input",ExpressionUUID->"3bd71288-3233-4171-a393-57251fb3dc52"], +Cell[411705, 8369, 6890, 131, 180, "Output",ExpressionUUID->"203873a1-0dd1-4567-8d12-97eeeed127f7"] }, Open ]], Cell[CellGroupData[{ -Cell[417696, 8464, 1358, 35, 48, "Input",ExpressionUUID->"39644cf7-65fa-4952-a674-714ff8619dd6"], -Cell[419057, 8501, 15118, 266, 179, "Output",ExpressionUUID->"16b2f1f1-4bda-4b9c-9a2d-c4721af32a8e"] +Cell[418632, 8505, 1358, 35, 48, "Input",ExpressionUUID->"39644cf7-65fa-4952-a674-714ff8619dd6"], +Cell[419993, 8542, 15118, 266, 179, "Output",ExpressionUUID->"16b2f1f1-4bda-4b9c-9a2d-c4721af32a8e"] }, Open ]], Cell[CellGroupData[{ -Cell[434212, 8772, 1607, 46, 65, "Input",ExpressionUUID->"72bac293-cf66-40a9-809d-4d0874dae037"], -Cell[435822, 8820, 59065, 1048, 175, "Output",ExpressionUUID->"d758ff21-6281-4951-879d-9a8a224caf38"] +Cell[435148, 8813, 1607, 46, 65, "Input",ExpressionUUID->"72bac293-cf66-40a9-809d-4d0874dae037"], +Cell[436758, 8861, 59065, 1048, 175, "Output",ExpressionUUID->"d758ff21-6281-4951-879d-9a8a224caf38"] }, Open ]], Cell[CellGroupData[{ -Cell[494924, 9873, 1002, 31, 42, "Input",ExpressionUUID->"41ea287c-9ee3-4b18-b17d-5dfc6a4e056d"], -Cell[495929, 9906, 3476, 103, 125, "Output",ExpressionUUID->"7b8e54ff-34ea-444e-8e78-ad376945c6d8"] +Cell[495860, 9914, 1002, 31, 42, "Input",ExpressionUUID->"41ea287c-9ee3-4b18-b17d-5dfc6a4e056d"], +Cell[496865, 9947, 3476, 103, 125, "Output",ExpressionUUID->"7b8e54ff-34ea-444e-8e78-ad376945c6d8"] }, Open ]], Cell[CellGroupData[{ -Cell[499442, 10014, 753, 23, 24, "Input",ExpressionUUID->"77ced036-93d3-4558-8eb2-1aeeb8345ed7"], -Cell[500198, 10039, 840, 23, 40, "Output",ExpressionUUID->"49c3b3ed-cf01-4c58-acac-6b3b7b9cb822"] +Cell[500378, 10055, 753, 23, 24, "Input",ExpressionUUID->"77ced036-93d3-4558-8eb2-1aeeb8345ed7"], +Cell[501134, 10080, 840, 23, 40, "Output",ExpressionUUID->"49c3b3ed-cf01-4c58-acac-6b3b7b9cb822"] }, Open ]], Cell[CellGroupData[{ -Cell[501075, 10067, 1262, 35, 41, "Input",ExpressionUUID->"1dc20dd2-2639-4b38-8eb3-e39481c6a15e"], -Cell[502340, 10104, 71734, 1205, 306, "Output",ExpressionUUID->"ae6c01c8-9abd-4358-a6f9-fb8f033a8ba4"] +Cell[502011, 10108, 1262, 35, 41, "Input",ExpressionUUID->"1dc20dd2-2639-4b38-8eb3-e39481c6a15e"], +Cell[503276, 10145, 71734, 1205, 306, "Output",ExpressionUUID->"ae6c01c8-9abd-4358-a6f9-fb8f033a8ba4"] }, Open ]], Cell[CellGroupData[{ -Cell[574111, 11314, 1045, 25, 40, "Input",ExpressionUUID->"25d7aa57-46b2-4829-8aad-b00b9ab2ecff"], -Cell[575159, 11341, 1385, 34, 41, "Output",ExpressionUUID->"5f0c650e-d3ec-4a10-aad2-c391b45db45f"] +Cell[575047, 11355, 1045, 25, 40, "Input",ExpressionUUID->"25d7aa57-46b2-4829-8aad-b00b9ab2ecff"], +Cell[576095, 11382, 1385, 34, 41, "Output",ExpressionUUID->"5f0c650e-d3ec-4a10-aad2-c391b45db45f"] }, Open ]] }, Open ]] } diff --git a/schofield.wl b/schofield.wl new file mode 100644 index 0000000..fdb76e3 --- /dev/null +++ b/schofield.wl @@ -0,0 +1,54 @@ + +BeginPackage["Schofield`"] + +β::usage = "Gives the exponent for the magnetization." +δ::usage = "Exponent." +α::usage = "Exponent." +Δ::usage = "Exponent." + +$Assumptions = {θc > 0, θi > 0, θc > θi, gC[_] ∈ Reals, B > 0, γ > 0} + +β[D_:2] := Piecewise[ + { + {1/8, D == 2}, + {0.326419, D == 3}, + {1/2, D == 4}, + {β, True} + } +] + +δ[D_:2] := Piecewise[ + { + {15, D == 2}, + {4.78984, D == 3}, + {3, D == 4}, + {δ, True} + } +] + +α[D_:2] := Piecewise[ + { + {0, D == 2}, + {0.11008, D == 3}, + {0, D == 4}, + {α, True} + } +] + +Δ[D_:2] := β[D] δ[D] + +f[θi_:1][n_][θ_] := (θ / θi)^2 - 1 +g[gC_:gC, θc_:θc][n_][θ_] := (1 - (θ/θc)^2) Sum[gC[i] θ^(2i+1), {i, 0, n}] + +I\[ScriptCapitalM]f[γ_][y_] := (1 + (1 + γ x) / x) Exp[-1/x] +R\[ScriptCapitalM]f[γ_][y_] := (1 - y - γ y) Exp[1/y] ExpIntegralEi[-1/y] / (π y) + +R\[ScriptCapitalM][2][γ_, B_, θc_, M0_][θ_] := - M0 (R\[ScriptCapitalM]f[γ][B(θc - θ)] - R\[ScriptCapitalM]f[γ][B(θc + θ)]) + +eqLow[D_:2][f_, g_][m_] := SeriesCoefficient[ + R\[ScriptCapitalM][D][γ, B, θc, M0][θ] + f[θ]^β[D] Gl'[g[θ] f[θ]^(-Δ[D])], + {θ, θc, m}, + Assumptions -> Join[$Assumptions, {θ < θc, θ > θi}] +] +EndPackage[] + -- cgit v1.2.3-70-g09d2