(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 13.1' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 42377, 1145] NotebookOptionsPosition[ 37437, 1058] NotebookOutlinePosition[ 37878, 1076] CellTagsIndexPosition[ 37835, 1073] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell["\<\ Christoffel Computation David Hilditch & Anton Khirnov, Feb. 2023.\ \>", "Subsubsection", CellChangeTimes->{{3.886238482476939*^9, 3.8862385051531477`*^9}, { 3.886238556279541*^9, 3.886238570895651*^9}},ExpressionUUID->"373d7a92-84df-4cc8-94b3-\ ee4da185a64c"], Cell[CellGroupData[{ Cell["Define coordinates, metric and inverse metric", "Text", CellChangeTimes->{{3.88623852761099*^9, 3.886238544018025*^9}},ExpressionUUID->"7b3be704-efcd-4b3c-a35c-\ 5fafcf85c005"], Cell[BoxData[ RowBox[{ RowBox[{"X", "=", RowBox[{"{", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}]}], ";"}]], "Input", CellChangeTimes->{{3.886171931510335*^9, 3.886171933623374*^9}, { 3.886171972930284*^9, 3.886171975545806*^9}, {3.88617265745574*^9, 3.886172663671729*^9}, {3.886173418537775*^9, 3.8861734198173847`*^9}}, CellLabel->"In[1]:=",ExpressionUUID->"33ae6975-ec1d-4e40-8c00-336b4482e361"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"g", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"gtt", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gtx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gty", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gtz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"gtx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gxx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gxy", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gxz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"gty", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gxy", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gyy", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gyz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"gtz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gxz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gyz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gzz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}]}], "}"}]}], ";", RowBox[{"g", "//", "MatrixForm"}]}]], "Input", CellChangeTimes->{{3.8861719282262793`*^9, 3.886171929989098*^9}, { 3.886171987992716*^9, 3.88617200715347*^9}, {3.886172316594161*^9, 3.886172333075539*^9}, {3.886172647434791*^9, 3.886172651808694*^9}, { 3.886172682276737*^9, 3.886172909660861*^9}, {3.8861729438025417`*^9, 3.886172948517481*^9}, {3.8861733110110073`*^9, 3.8861734651163197`*^9}, { 3.8861744456232347`*^9, 3.886174447108243*^9}, 3.886174614151684*^9}, CellLabel->"In[2]:=",ExpressionUUID->"bfbd2c04-efbd-49b8-940f-e8e2a48cee58"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"gtt", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gtx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gty", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gtz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}, { RowBox[{"gtx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gxx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gxy", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gxz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}, { RowBox[{"gty", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gxy", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gyy", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gyz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}, { RowBox[{"gtz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gxz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gyz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gzz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]} }, GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{{3.886172749895281*^9, 3.8861727752028303`*^9}, { 3.886172893537417*^9, 3.886172912569895*^9}, 3.886172949054069*^9, 3.886173465656674*^9, 3.886174204828216*^9, 3.886174447409831*^9, 3.8861746149973183`*^9, 3.886238766426577*^9}, CellLabel-> "Out[2]//MatrixForm=",ExpressionUUID->"afb62946-6d22-40db-87c7-\ 06d8aaa5eafb"] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"ginv", "=", RowBox[{ RowBox[{"g", "//", "Inverse"}], "//", "Simplify"}]}], ";"}]], "Input", CellChangeTimes->{{3.886173014702764*^9, 3.88617303331397*^9}}, CellLabel->"In[3]:=",ExpressionUUID->"ff8d5500-b0da-4dea-aed3-005d3f38ce38"] }, Closed]], Cell[CellGroupData[{ Cell["Compute Christoffels (ddd & udd)", "Text", CellChangeTimes->{{3.8862386105171843`*^9, 3.8862386296497107`*^9}},ExpressionUUID->"0e39b00c-5443-4cc9-9420-\ 976408ee456d"], Cell[BoxData[ RowBox[{ RowBox[{"Christoffelsddd", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"1", "/", "2"}], RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"X", "[", RowBox[{"[", "\[Mu]", "]"}], "]"}]], RowBox[{"g", "[", RowBox[{"[", RowBox[{"\[Nu]", ",", "\[Delta]"}], "]"}], "]"}]}], "+", RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"X", "[", RowBox[{"[", "\[Nu]", "]"}], "]"}]], RowBox[{"g", "[", RowBox[{"[", RowBox[{"\[Mu]", ",", "\[Delta]"}], "]"}], "]"}]}], "-", RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"X", "[", RowBox[{"[", "\[Delta]", "]"}], "]"}]], RowBox[{"g", "[", RowBox[{"[", RowBox[{"\[Mu]", ",", "\[Nu]"}], "]"}], "]"}]}]}], ")"}]}], ",", RowBox[{"{", RowBox[{"\[Delta]", ",", "4"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Mu]", ",", "4"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Nu]", ",", "4"}], "}"}]}], "]"}], "//", "Simplify"}]}], ";"}]], "Input", CellChangeTimes->{{3.8861743287685547`*^9, 3.886174365030443*^9}}, CellLabel->"In[4]:=",ExpressionUUID->"ac869ef8-5c04-4ff7-b6ff-063dc324d8f4"], Cell[BoxData[ RowBox[{ RowBox[{"Christoffelsudd", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"1", "/", "2"}], RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"ginv", "[", RowBox[{"[", RowBox[{"\[Alpha]", ",", "\[Delta]"}], "]"}], "]"}], RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"X", "[", RowBox[{"[", "\[Mu]", "]"}], "]"}]], RowBox[{"g", "[", RowBox[{"[", RowBox[{"\[Nu]", ",", "\[Delta]"}], "]"}], "]"}]}], "+", RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"X", "[", RowBox[{"[", "\[Nu]", "]"}], "]"}]], RowBox[{"g", "[", RowBox[{"[", RowBox[{"\[Mu]", ",", "\[Delta]"}], "]"}], "]"}]}], "-", RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"X", "[", RowBox[{"[", "\[Delta]", "]"}], "]"}]], RowBox[{"g", "[", RowBox[{"[", RowBox[{"\[Mu]", ",", "\[Nu]"}], "]"}], "]"}]}]}], ")"}]}], ",", RowBox[{"{", RowBox[{"\[Delta]", ",", "4"}], "}"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"\[Alpha]", ",", "4"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Mu]", ",", "4"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Nu]", ",", "4"}], "}"}]}], "]"}], "//", "Simplify"}]}], ";"}]], "Input", CellChangeTimes->{{3.8861729804852133`*^9, 3.886173003793852*^9}, { 3.886173056335821*^9, 3.886173105859511*^9}, {3.886173175235386*^9, 3.88617327140637*^9}, 3.886173478420867*^9, {3.886173835805682*^9, 3.8861738550901203`*^9}, {3.886174098114808*^9, 3.886174113642325*^9}, { 3.886174235751086*^9, 3.8861742386193647`*^9}}, CellLabel->"In[5]:=",ExpressionUUID->"e7e016e2-c00b-42d2-91be-a38628e97843"] }, Closed]], Cell[CellGroupData[{ Cell["\<\ Rules for twist-free axisymmetry with reflection symmetry, on the symmetry \ axis and equatorial plane\ \>", "Text", CellChangeTimes->{{3.8862386444103394`*^9, 3.88623870518185*^9}},ExpressionUUID->"c11437e5-f609-4db8-bb2d-\ e4ba921912d9"], Cell[BoxData[ RowBox[{ RowBox[{"KillMeOnAxisRules", "=", RowBox[{"{", RowBox[{ RowBox[{"gtx", "->", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], "]"}]}], ",", RowBox[{"gty", "->", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], "]"}]}], ",", RowBox[{"gxy", "->", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], "]"}]}], ",", RowBox[{"gyz", "->", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], "]"}]}], ",", RowBox[{"gxz", "\[Rule]", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], "]"}]}]}], "}"}]}], ";"}]], "Input", CellChangeTimes->{{3.88617348439834*^9, 3.886173509870784*^9}, { 3.886173553075624*^9, 3.8861735746880627`*^9}, {3.886173613403329*^9, 3.886173656011073*^9}, {3.886173711669825*^9, 3.886173730830943*^9}, { 3.886174771825536*^9, 3.886174773975914*^9}}, CellLabel->"In[6]:=",ExpressionUUID->"1f3a5dca-b140-4e51-93fd-93ff73818e14"], Cell[BoxData[ RowBox[{ RowBox[{"KillMeEquatorRules", "=", RowBox[{"{", RowBox[{ RowBox[{"gtz", "->", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], "]"}]}], ",", RowBox[{"gty", "->", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], "]"}]}], ",", RowBox[{"gxy", "->", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], "]"}]}], ",", RowBox[{"gyz", "->", RowBox[{"Function", "[", RowBox[{ RowBox[{"{", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "}"}], ",", "0"}], "]"}]}], ",", RowBox[{ RowBox[{"gxz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "->", "0"}]}], "}"}]}], ";"}]], "Input", CellChangeTimes->{{3.886174790055666*^9, 3.886174820279228*^9}, { 3.886174952284816*^9, 3.886174960675373*^9}}, CellLabel->"In[7]:=",ExpressionUUID->"9461b4d0-c922-4076-a54e-0d59230005dd"] }, Closed]], Cell[CellGroupData[{ Cell["Examine metric on axis and in equator :", "Text", CellChangeTimes->{{3.886238741422525*^9, 3.8862387597451887`*^9}},ExpressionUUID->"9b144fb2-3580-49bd-9bde-\ 8f9e9d66dc8a"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"g", "/.", "KillMeEquatorRules"}], "//", "MatrixForm"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"g", "/.", "KillMeOnAxisRules"}], "//", "MatrixForm"}]}], "Input", CellChangeTimes->{{3.886239099478282*^9, 3.8862391391676826`*^9}}, CellLabel->"In[12]:=",ExpressionUUID->"bd17db4b-99a8-4a6c-9eb7-49accaf01cdd"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"gtt", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gtx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "0", "0"}, { RowBox[{"gtx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gxx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "0", "0"}, {"0", "0", RowBox[{"gyy", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "0"}, {"0", "0", "0", RowBox[{"gzz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]} }, GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{{3.8862391317376204`*^9, 3.8862391395839252`*^9}}, CellLabel-> "Out[12]//MatrixForm=",ExpressionUUID->"bb34dbee-88b7-4abe-9f81-\ 02f49b097c80"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"gtt", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "0", "0", RowBox[{"gtz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}, {"0", RowBox[{"gxx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "0", "0"}, {"0", "0", RowBox[{"gyy", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "0"}, { RowBox[{"gtz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], "0", "0", RowBox[{"gzz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]} }, GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{{3.8862391317376204`*^9, 3.8862391395856543`*^9}}, CellLabel-> "Out[13]//MatrixForm=",ExpressionUUID->"d86ab8bb-5d04-472e-a292-\ 412d5e32b5d3"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Collect the relevant components on Axis:", "Text", CellChangeTimes->{{3.886238862895907*^9, 3.886238894746866*^9}},ExpressionUUID->"ac3861a3-e496-4de1-ad12-\ e033917d7b88"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Christoffelsddd", "[", RowBox[{"[", RowBox[{"1", ",", "1", ",", "1"}], "]"}], "]"}], ",", RowBox[{"Christoffelsddd", "[", RowBox[{"[", RowBox[{"4", ",", "1", ",", "1"}], "]"}], "]"}], ",", RowBox[{"Christoffelsddd", "[", RowBox[{"[", RowBox[{"1", ",", "4", ",", "4"}], "]"}], "]"}], ",", RowBox[{"Christoffelsddd", "[", RowBox[{"[", RowBox[{"4", ",", "4", ",", "4"}], "]"}], "]"}], ",", RowBox[{"Christoffelsddd", "[", RowBox[{"[", RowBox[{"1", ",", "1", ",", "4"}], "]"}], "]"}], ",", RowBox[{"Christoffelsddd", "[", RowBox[{"[", RowBox[{"4", ",", "1", ",", "4"}], "]"}], "]"}]}], "}"}], "/.", "KillMeOnAxisRules"}], "//", "Simplify"}], ";"}]], "Input", CellChangeTimes->{{3.886173999514325*^9, 3.8861740245071373`*^9}, 3.886174252714652*^9, 3.886174383803364*^9, {3.8861749956068773`*^9, 3.886174997474677*^9}, {3.886175279412593*^9, 3.886175386312916*^9}, { 3.886175446572671*^9, 3.886175449889043*^9}, 3.886239866253716*^9}, CellLabel->"In[62]:=",ExpressionUUID->"49dcf80e-4d26-47f8-aad5-d314e4df0e69"], Cell["\<\ Double check that we can get the udd part of the metric by multiplying with \ the inverse of just one block of the inverse:\ \>", "Text", CellChangeTimes->{{3.886238889485704*^9, 3.886238900333441*^9}, { 3.886239153187326*^9, 3.8862391961989613`*^9}},ExpressionUUID->"92d83e45-e318-4d91-be89-\ af13b6883cab"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"g", "/.", "KillMeOnAxisRules"}]], "Input", CellChangeTimes->{{3.8862392389049063`*^9, 3.886239242557054*^9}}, CellLabel->"In[14]:=",ExpressionUUID->"de8bd87c-f1a9-413f-9eb9-806bac2b9685"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"gtt", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", "0", ",", "0", ",", RowBox[{"gtz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{"gxx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{"gyy", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"gtz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", "0", ",", "0", ",", RowBox[{"gzz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.886239242979628*^9}, CellLabel->"Out[14]=",ExpressionUUID->"76943539-2335-4cc0-9e7e-b0a0226fa3ad"] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"XAxi", "=", RowBox[{"{", RowBox[{"t", ",", "z"}], "}"}]}], ";"}]], "Input", CellChangeTimes->{{3.886239392262786*^9, 3.886239399279251*^9}}, CellLabel->"In[21]:=",ExpressionUUID->"01ae6f5e-d085-4a3b-928b-925b844bd0ac"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"gOnAxisBlock", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"gtt", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gtz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"gtz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gzz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}]}], "}"}]}], ";", RowBox[{"gOnAxisBlock", "//", "MatrixForm"}]}]], "Input", CellChangeTimes->{{3.886239203828929*^9, 3.886239300860549*^9}}, CellLabel->"In[19]:=",ExpressionUUID->"8d6186c2-bad4-4e3e-9b47-61d5bbf88a87"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"gtt", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gtz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}, { RowBox[{"gtz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gzz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]} }, GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.886239301125842*^9}, CellLabel-> "Out[19]//MatrixForm=",ExpressionUUID->"2977af70-299b-46e5-9314-\ 010cb4214309"] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"gOnAxisBlockinv", "=", RowBox[{ RowBox[{"gOnAxisBlock", "//", "Inverse"}], "//", "Simplify"}]}], ";"}]], "Input", CellChangeTimes->{{3.886239279053372*^9, 3.8862393025269327`*^9}}, CellLabel->"In[20]:=",ExpressionUUID->"26c014b5-335e-41f3-a343-c5c438670da3"], Cell[BoxData[ RowBox[{ RowBox[{"ChristoffelsAxisBlockddd", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"1", "/", "2"}], RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"XAxi", "[", RowBox[{"[", "\[Mu]", "]"}], "]"}]], RowBox[{"gOnAxisBlock", "[", RowBox[{"[", RowBox[{"\[Nu]", ",", "\[Delta]"}], "]"}], "]"}]}], "+", RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"XAxi", "[", RowBox[{"[", "\[Nu]", "]"}], "]"}]], RowBox[{"gOnAxisBlock", "[", RowBox[{"[", RowBox[{"\[Mu]", ",", "\[Delta]"}], "]"}], "]"}]}], "-", RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"XAxi", "[", RowBox[{"[", "\[Delta]", "]"}], "]"}]], RowBox[{"gOnAxisBlock", "[", RowBox[{"[", RowBox[{"\[Mu]", ",", "\[Nu]"}], "]"}], "]"}]}]}], ")"}]}], ",", RowBox[{"{", RowBox[{"\[Delta]", ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Mu]", ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Nu]", ",", "2"}], "}"}]}], "]"}], "//", "Simplify"}]}], ";"}]], "Input", CellChangeTimes->{{3.8862393767955112`*^9, 3.886239447571889*^9}}, CellLabel->"In[25]:=",ExpressionUUID->"de4e5987-1579-4a27-b9fc-70eb429269ae"], Cell[BoxData[ RowBox[{ RowBox[{"ChristoffelsAxisBlockudd", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"gOnAxisBlockinv", "[", RowBox[{"[", RowBox[{"\[Alpha]", ",", "\[Delta]"}], "]"}], "]"}], RowBox[{"ChristoffelsAxisBlockddd", "[", RowBox[{"[", RowBox[{"\[Delta]", ",", "\[Mu]", ",", "\[Nu]"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"\[Delta]", ",", "2"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"\[Alpha]", ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Mu]", ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Nu]", ",", "2"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.8862394559284563`*^9, 3.886239526238406*^9}, { 3.886239680673327*^9, 3.886239685179586*^9}}, CellLabel->"In[36]:=",ExpressionUUID->"640e7de3-3183-400a-8794-02a2c9d32492"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"ChristoffelsAxisBlockudd", "[", RowBox[{"[", RowBox[{"1", ",", "1", ",", "1"}], "]"}], "]"}], "-", RowBox[{"Christoffelsudd", "[", RowBox[{"[", RowBox[{"1", ",", "1", ",", "1"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{"ChristoffelsAxisBlockudd", "[", RowBox[{"[", RowBox[{"1", ",", "1", ",", "2"}], "]"}], "]"}], "-", RowBox[{"Christoffelsudd", "[", RowBox[{"[", RowBox[{"1", ",", "1", ",", "4"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{"ChristoffelsAxisBlockudd", "[", RowBox[{"[", RowBox[{"1", ",", "2", ",", "2"}], "]"}], "]"}], "-", RowBox[{"Christoffelsudd", "[", RowBox[{"[", RowBox[{"1", ",", "4", ",", "4"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{"ChristoffelsAxisBlockudd", "[", RowBox[{"[", RowBox[{"2", ",", "1", ",", "1"}], "]"}], "]"}], "-", RowBox[{"Christoffelsudd", "[", RowBox[{"[", RowBox[{"4", ",", "1", ",", "1"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{"ChristoffelsAxisBlockudd", "[", RowBox[{"[", RowBox[{"2", ",", "1", ",", "2"}], "]"}], "]"}], "-", RowBox[{"Christoffelsudd", "[", RowBox[{"[", RowBox[{"4", ",", "1", ",", "4"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{"ChristoffelsAxisBlockudd", "[", RowBox[{"[", RowBox[{"2", ",", "2", ",", "2"}], "]"}], "]"}], "-", RowBox[{"Christoffelsudd", "[", RowBox[{"[", RowBox[{"4", ",", "4", ",", "4"}], "]"}], "]"}]}]}], "}"}], "/.", "KillMeOnAxisRules"}], "//", "Simplify"}]], "Input", CellChangeTimes->{{3.886239550225062*^9, 3.88623958206875*^9}, { 3.886239692524212*^9, 3.8862398032153063`*^9}, {3.886240229048933*^9, 3.886240259620287*^9}}, CellLabel->"In[74]:=",ExpressionUUID->"06a60ae0-778e-4912-85dc-d32667a00362"], Cell[BoxData[ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}]], "Output", CellChangeTimes->{{3.886239569587137*^9, 3.8862395822670918`*^9}, { 3.8862396864278193`*^9, 3.886239734674528*^9}, {3.8862397688325167`*^9, 3.8862398035105886`*^9}, {3.886240252178359*^9, 3.886240264435549*^9}}, CellLabel->"Out[74]=",ExpressionUUID->"9359abe8-5143-4c48-81db-e7d805acb053"] }, Open ]], Cell["All is well!", "Text", CellChangeTimes->{{3.886239880987774*^9, 3.886239885684226*^9}},ExpressionUUID->"221e1d75-6a4b-4bc8-9548-\ 1d2a9bac055f"] }, Closed]], Cell[CellGroupData[{ Cell["Collect the relevant components at the equator:", "Text", CellChangeTimes->{{3.886238862895907*^9, 3.886238894746866*^9}, { 3.8862399139244413`*^9, 3.886239918455309*^9}},ExpressionUUID->"0c813465-1111-450c-abee-\ 81cee2d5e354"], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Christoffelsddd", "[", RowBox[{"[", RowBox[{"1", ",", "1", ",", "1"}], "]"}], "]"}], ",", RowBox[{"Christoffelsddd", "[", RowBox[{"[", RowBox[{"2", ",", "1", ",", "1"}], "]"}], "]"}], ",", RowBox[{"Christoffelsddd", "[", RowBox[{"[", RowBox[{"1", ",", "2", ",", "2"}], "]"}], "]"}], ",", RowBox[{"Christoffelsddd", "[", RowBox[{"[", RowBox[{"2", ",", "2", ",", "2"}], "]"}], "]"}], ",", RowBox[{"Christoffelsddd", "[", RowBox[{"[", RowBox[{"1", ",", "1", ",", "2"}], "]"}], "]"}], ",", RowBox[{"Christoffelsddd", "[", RowBox[{"[", RowBox[{"2", ",", "1", ",", "2"}], "]"}], "]"}]}], "}"}], "/.", "KillMeEquatorRules"}], "//", "Simplify"}], ";"}]], "Input", CellChangeTimes->{{3.886173999514325*^9, 3.8861740245071373`*^9}, 3.886174252714652*^9, 3.886174383803364*^9, {3.8861749956068773`*^9, 3.886174997474677*^9}, {3.886175279412593*^9, 3.886175386312916*^9}, { 3.886175446572671*^9, 3.886175449889043*^9}, 3.886239866253716*^9, { 3.8862399323925657`*^9, 3.886239950169238*^9}}, CellLabel->"In[63]:=",ExpressionUUID->"d230e578-d7a3-4696-81a1-55a1e97d698f"], Cell["\<\ Double check that we can get the udd part of the metric by multiplying with \ the inverse of just one block of the inverse:\ \>", "Text", CellChangeTimes->{{3.886238889485704*^9, 3.886238900333441*^9}, { 3.886239153187326*^9, 3.8862391961989613`*^9}},ExpressionUUID->"1b418052-d781-4b3a-886c-\ 17105256f20d"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"g", "/.", "KillMeEquatorRules"}]], "Input", CellChangeTimes->{{3.8862392389049063`*^9, 3.886239242557054*^9}, { 3.886239955096569*^9, 3.8862399559882517`*^9}}, CellLabel->"In[75]:=",ExpressionUUID->"49bce66e-c2ac-4737-a18b-83b9af6efec3"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"gtt", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gtx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"gtx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gxx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{"gyy", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0", ",", RowBox[{"gzz", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.886239242979628*^9, 3.886239956472385*^9, 3.886240293779134*^9}, CellLabel->"Out[75]=",ExpressionUUID->"dcbac780-95b9-4925-9a81-911c83eabaee"] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"XEq", "=", RowBox[{"{", RowBox[{"t", ",", "x"}], "}"}]}], ";"}]], "Input", CellChangeTimes->{{3.886239392262786*^9, 3.886239399279251*^9}, { 3.8862399619870443`*^9, 3.886239964013056*^9}}, CellLabel->"In[76]:=",ExpressionUUID->"2f1d3424-d70f-481b-a66e-b1449c0fcc65"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"gEquatorBlock", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"gtt", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gtx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"gtx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], ",", RowBox[{"gxx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}], "}"}]}], "}"}]}], ";", RowBox[{"gEquatorBlock", "//", "MatrixForm"}]}]], "Input", CellChangeTimes->{{3.886239203828929*^9, 3.886239300860549*^9}, { 3.8862399680307703`*^9, 3.886239991444262*^9}}, CellLabel->"In[77]:=",ExpressionUUID->"4c2fad4b-2113-406e-9140-500b5acb0e14"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"gtt", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gtx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]}, { RowBox[{"gtx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}], RowBox[{"gxx", "[", RowBox[{"t", ",", "x", ",", "y", ",", "z"}], "]"}]} }, GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.886239301125842*^9, 3.886239992311648*^9, 3.8862402203987093`*^9, 3.88624029644741*^9}, CellLabel-> "Out[77]//MatrixForm=",ExpressionUUID->"4debc2a1-1480-4b9c-9db2-\ 294463da8105"] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"gEquatorBlockinv", "=", RowBox[{ RowBox[{"gEquatorBlock", "//", "Inverse"}], "//", "Simplify"}]}], ";"}]], "Input", CellChangeTimes->{{3.886239279053372*^9, 3.8862393025269327`*^9}, { 3.8862400002414837`*^9, 3.886240007499395*^9}}, CellLabel->"In[78]:=",ExpressionUUID->"9f98fd8a-6fa5-4aaa-b85b-5ca08fae24a6"], Cell[BoxData[ RowBox[{ RowBox[{"ChristoffelsEquatorBlockddd", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"1", "/", "2"}], RowBox[{"(", RowBox[{ RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"XEq", "[", RowBox[{"[", "\[Mu]", "]"}], "]"}]], RowBox[{"gEquatorBlock", "[", RowBox[{"[", RowBox[{"\[Nu]", ",", "\[Delta]"}], "]"}], "]"}]}], "+", RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"XEq", "[", RowBox[{"[", "\[Nu]", "]"}], "]"}]], RowBox[{"gEquatorBlock", "[", RowBox[{"[", RowBox[{"\[Mu]", ",", "\[Delta]"}], "]"}], "]"}]}], "-", RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"XEq", "[", RowBox[{"[", "\[Delta]", "]"}], "]"}]], RowBox[{"gEquatorBlock", "[", RowBox[{"[", RowBox[{"\[Mu]", ",", "\[Nu]"}], "]"}], "]"}]}]}], ")"}]}], ",", RowBox[{"{", RowBox[{"\[Delta]", ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Mu]", ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Nu]", ",", "2"}], "}"}]}], "]"}], "//", "Simplify"}]}], ";"}]], "Input", CellChangeTimes->{{3.8862393767955112`*^9, 3.886239447571889*^9}, { 3.8862400174952917`*^9, 3.886240060859981*^9}}, CellLabel->"In[80]:=",ExpressionUUID->"f8b887e2-1f20-49ec-8a4e-6bdd7b719e87"], Cell[BoxData[ RowBox[{ RowBox[{"ChristoffelsEquatorBlockudd", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"gEquatorBlockinv", "[", RowBox[{"[", RowBox[{"\[Alpha]", ",", "\[Delta]"}], "]"}], "]"}], RowBox[{"ChristoffelsEquatorBlockddd", "[", RowBox[{"[", RowBox[{"\[Delta]", ",", "\[Mu]", ",", "\[Nu]"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"\[Delta]", ",", "2"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"\[Alpha]", ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Mu]", ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Nu]", ",", "2"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.8862394559284563`*^9, 3.886239526238406*^9}, { 3.886239680673327*^9, 3.886239685179586*^9}, {3.8862400721140003`*^9, 3.88624009482312*^9}}, CellLabel->"In[81]:=",ExpressionUUID->"c5674cf3-b1f7-4213-ab57-9659993e8b80"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"ChristoffelsEquatorBlockudd", "[", RowBox[{"[", RowBox[{"1", ",", "1", ",", "1"}], "]"}], "]"}], "-", RowBox[{"Christoffelsudd", "[", RowBox[{"[", RowBox[{"1", ",", "1", ",", "1"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{"ChristoffelsEquatorBlockudd", "[", RowBox[{"[", RowBox[{"1", ",", "1", ",", "2"}], "]"}], "]"}], "-", RowBox[{"Christoffelsudd", "[", RowBox[{"[", RowBox[{"1", ",", "1", ",", "2"}], "]"}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"ChristoffelsEquatorBlockudd", "[", RowBox[{"[", RowBox[{"1", ",", "2", ",", "2"}], "]"}], "]"}], "-", RowBox[{"Christoffelsudd", "[", RowBox[{"[", RowBox[{"1", ",", "2", ",", "2"}], "]"}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"ChristoffelsEquatorBlockudd", "[", RowBox[{"[", RowBox[{"2", ",", "1", ",", "1"}], "]"}], "]"}], "-", RowBox[{"Christoffelsudd", "[", RowBox[{"[", RowBox[{"2", ",", "1", ",", "1"}], "]"}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"ChristoffelsEquatorBlockudd", "[", RowBox[{"[", RowBox[{"2", ",", "1", ",", "2"}], "]"}], "]"}], "-", RowBox[{"Christoffelsudd", "[", RowBox[{"[", RowBox[{"2", ",", "1", ",", "2"}], "]"}], "]"}]}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"ChristoffelsEquatorBlockudd", "[", RowBox[{"[", RowBox[{"2", ",", "2", ",", "2"}], "]"}], "]"}], "-", RowBox[{"Christoffelsudd", "[", RowBox[{"[", RowBox[{"2", ",", "2", ",", "2"}], "]"}], "]"}]}]}], "}"}], "/.", "KillMeEquatorRules"}], "//", "Simplify"}]], "Input", CellChangeTimes->{{3.886239550225062*^9, 3.88623958206875*^9}, { 3.886239692524212*^9, 3.8862398032153063`*^9}, {3.886240109340292*^9, 3.886240207302966*^9}}, CellLabel->"In[82]:=",ExpressionUUID->"cd24a2c2-cec1-4e13-8ee2-25af439fceaa"], Cell[BoxData[ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}]], "Output", CellChangeTimes->{{3.886239569587137*^9, 3.8862395822670918`*^9}, { 3.8862396864278193`*^9, 3.886239734674528*^9}, {3.8862397688325167`*^9, 3.8862398035105886`*^9}, 3.886240207733821*^9, 3.886240302195738*^9}, CellLabel->"Out[82]=",ExpressionUUID->"33e9b9d5-5aeb-406f-b752-17986db4beb9"] }, Open ]], Cell["All is well!", "Text", CellChangeTimes->{{3.886239880987774*^9, 3.886239885684226*^9}},ExpressionUUID->"b4df87a0-82a6-4fbd-a126-\ 8cedb791ad67"] }, Closed]] }, CellGrouping->Manual, WindowSize->{1440, 852}, WindowMargins->{{0, Automatic}, {0, Automatic}}, Magnification:>1.5 Inherited, FrontEndVersion->"13.1 for Linux x86 (64-bit) (June 16, 2022)", StyleDefinitions->"Default.nb", ExpressionUUID->"4f361d4c-d8b0-44e2-b48c-8cc05b0570b8" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[558, 20, 276, 7, 107, "Subsubsection",ExpressionUUID->"373d7a92-84df-4cc8-94b3-ee4da185a64c"], Cell[CellGroupData[{ Cell[859, 31, 186, 3, 44, "Text",ExpressionUUID->"7b3be704-efcd-4b3c-a35c-5fafcf85c005"], Cell[1048, 36, 425, 8, 45, "Input",ExpressionUUID->"33ae6975-ec1d-4e40-8c00-336b4482e361"], Cell[CellGroupData[{ Cell[1498, 48, 2332, 53, 173, "Input",ExpressionUUID->"bfbd2c04-efbd-49b8-940f-e8e2a48cee58"], Cell[3833, 103, 2223, 56, 147, "Output",ExpressionUUID->"afb62946-6d22-40db-87c7-06d8aaa5eafb"] }, Open ]], Cell[6071, 162, 277, 6, 45, "Input",ExpressionUUID->"ff8d5500-b0da-4dea-aed3-005d3f38ce38"] }, Closed]], Cell[CellGroupData[{ Cell[6385, 173, 178, 3, 44, "Text",ExpressionUUID->"0e39b00c-5443-4cc9-9420-976408ee456d"], Cell[6566, 178, 1392, 39, 49, "Input",ExpressionUUID->"ac869ef8-5c04-4ff7-b6ff-063dc324d8f4"], Cell[7961, 219, 2017, 52, 80, "Input",ExpressionUUID->"e7e016e2-c00b-42d2-91be-a38628e97843"] }, Closed]], Cell[CellGroupData[{ Cell[10015, 276, 253, 6, 44, "Text",ExpressionUUID->"c11437e5-f609-4db8-bb2d-e4ba921912d9"], Cell[10271, 284, 1420, 39, 77, "Input",ExpressionUUID->"1f3a5dca-b140-4e51-93fd-93ff73818e14"], Cell[11694, 325, 1202, 35, 77, "Input",ExpressionUUID->"9461b4d0-c922-4076-a54e-0d59230005dd"] }, Closed]], Cell[CellGroupData[{ Cell[12933, 365, 183, 3, 44, "Text",ExpressionUUID->"9b144fb2-3580-49bd-9bde-8f9e9d66dc8a"], Cell[CellGroupData[{ Cell[13141, 372, 352, 7, 77, "Input",ExpressionUUID->"bd17db4b-99a8-4a6c-9eb7-49accaf01cdd"], Cell[13496, 381, 1214, 33, 147, "Output",ExpressionUUID->"bb34dbee-88b7-4abe-9f81-02f49b097c80"], Cell[14713, 416, 1214, 33, 147, "Output",ExpressionUUID->"d86ab8bb-5d04-472e-a292-412d5e32b5d3"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[15976, 455, 182, 3, 44, "Text",ExpressionUUID->"ac3861a3-e496-4de1-ad12-e033917d7b88"], Cell[16161, 460, 1255, 29, 77, "Input",ExpressionUUID->"49dcf80e-4d26-47f8-aad5-d314e4df0e69"], Cell[17419, 491, 324, 7, 52, "Text",ExpressionUUID->"92d83e45-e318-4d91-be89-af13b6883cab"], Cell[CellGroupData[{ Cell[17768, 502, 213, 3, 45, "Input",ExpressionUUID->"de8bd87c-f1a9-413f-9eb9-806bac2b9685"], Cell[17984, 507, 995, 29, 83, "Output",ExpressionUUID->"76943539-2335-4cc0-9e7e-b0a0226fa3ad"] }, Open ]], Cell[18994, 539, 264, 6, 45, "Input",ExpressionUUID->"01ae6f5e-d085-4a3b-928b-925b844bd0ac"], Cell[CellGroupData[{ Cell[19283, 549, 768, 20, 45, "Input",ExpressionUUID->"8d6186c2-bad4-4e3e-9b47-61d5bbf88a87"], Cell[20054, 571, 948, 27, 94, "Output",ExpressionUUID->"2977af70-299b-46e5-9314-010cb4214309"] }, Open ]], Cell[21017, 601, 306, 7, 45, "Input",ExpressionUUID->"26c014b5-335e-41f3-a343-c5c438670da3"], Cell[21326, 610, 1444, 39, 115, "Input",ExpressionUUID->"de4e5987-1579-4a27-b9fc-70eb429269ae"], Cell[22773, 651, 959, 25, 77, "Input",ExpressionUUID->"640e7de3-3183-400a-8794-02a2c9d32492"], Cell[CellGroupData[{ Cell[23757, 680, 2000, 51, 205, "Input",ExpressionUUID->"06a60ae0-778e-4912-85dc-d32667a00362"], Cell[25760, 733, 413, 7, 51, "Output",ExpressionUUID->"9359abe8-5143-4c48-81db-e7d805acb053"] }, Open ]], Cell[26188, 743, 154, 3, 52, "Text",ExpressionUUID->"221e1d75-6a4b-4bc8-9548-1d2a9bac055f"] }, Closed]], Cell[CellGroupData[{ Cell[26379, 751, 240, 4, 44, "Text",ExpressionUUID->"0c813465-1111-450c-abee-81cee2d5e354"], Cell[26622, 757, 1308, 30, 77, "Input",ExpressionUUID->"d230e578-d7a3-4696-81a1-55a1e97d698f"], Cell[27933, 789, 324, 7, 52, "Text",ExpressionUUID->"1b418052-d781-4b3a-886c-17105256f20d"], Cell[CellGroupData[{ Cell[28282, 800, 265, 4, 45, "Input",ExpressionUUID->"49bce66e-c2ac-4737-a18b-83b9af6efec3"], Cell[28550, 806, 1036, 29, 83, "Output",ExpressionUUID->"dcbac780-95b9-4925-9a81-911c83eabaee"] }, Open ]], Cell[29601, 838, 314, 7, 45, "Input",ExpressionUUID->"2f1d3424-d70f-481b-a66e-b1449c0fcc65"], Cell[CellGroupData[{ Cell[29940, 849, 821, 21, 45, "Input",ExpressionUUID->"4c2fad4b-2113-406e-9140-500b5acb0e14"], Cell[30764, 872, 1018, 28, 94, "Output",ExpressionUUID->"4debc2a1-1480-4b9c-9db2-294463da8105"] }, Open ]], Cell[31797, 903, 359, 8, 45, "Input",ExpressionUUID->"9f98fd8a-6fa5-4aaa-b85b-5ca08fae24a6"], Cell[32159, 913, 1498, 40, 115, "Input",ExpressionUUID->"f8b887e2-1f20-49ec-8a4e-6bdd7b719e87"], Cell[33660, 955, 1016, 26, 77, "Input",ExpressionUUID->"c5674cf3-b1f7-4213-ab57-9659993e8b80"], Cell[CellGroupData[{ Cell[34701, 985, 2135, 55, 237, "Input",ExpressionUUID->"cd24a2c2-cec1-4e13-8ee2-25af439fceaa"], Cell[36839, 1042, 413, 7, 51, "Output",ExpressionUUID->"33e9b9d5-5aeb-406f-b752-17986db4beb9"] }, Open ]], Cell[37267, 1052, 154, 3, 52, "Text",ExpressionUUID->"b4df87a0-82a6-4fbd-a126-8cedb791ad67"] }, Closed]] } ] *)