Skip to content
Snippets Groups Projects
Commit 55fb376d authored by Nitschke, Ingo's avatar Nitschke, Ingo
Browse files

2023-11-09-18-02-42

parent eec97d1b
No related branches found
No related tags found
No related merge requests found
(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 13.2' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 15704, 421]
NotebookOptionsPosition[ 13328, 379]
NotebookOutlinePosition[ 13728, 395]
CellTagsIndexPosition[ 13685, 392]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"$Assumptions", "=",
RowBox[{"r", ">=", "0"}]}]], "Input",
CellChangeTimes->{{3.908537042470517*^9, 3.908537054965324*^9}},
CellLabel->"In[11]:=",ExpressionUUID->"07d7ef80-8072-4d1f-adc2-6dc778f3a6d4"],
Cell[BoxData[
RowBox[{"r", "\[GreaterEqual]", "0"}]], "Output",
CellChangeTimes->{3.908537056694439*^9},
CellLabel->"Out[11]=",ExpressionUUID->"7e3882c2-5e75-4860-8715-7287c13aef00"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{"X", "=",
RowBox[{"r", "*",
RowBox[{"{",
RowBox[{
RowBox[{"Cos", "[", "\[CurlyPhi]", "]"}], ",",
RowBox[{"Sin", "[", "\[CurlyPhi]", "]"}]}],
"}"}]}]}], "\[IndentingNewLine]",
RowBox[{"DX", "=",
RowBox[{
RowBox[{"Transpose", "@",
RowBox[{"D", "[",
RowBox[{"X", ",",
RowBox[{"{",
RowBox[{"{",
RowBox[{"r", ",", "\[CurlyPhi]"}], "}"}], "}"}]}], "]"}]}], "//",
"Simplify"}]}], "\[IndentingNewLine]",
RowBox[{"g", "=",
RowBox[{
RowBox[{"Outer", "[",
RowBox[{"Dot", ",", "DX", ",", "DX", ",", "1"}], "]"}], "//",
"Simplify"}]}], "\[IndentingNewLine]",
RowBox[{"gInv", "=",
RowBox[{
RowBox[{"Inverse", "[", "g", "]"}], "//",
"Simplify"}]}], "\[IndentingNewLine]",
RowBox[{"mu", "=",
RowBox[{
RowBox[{"Sqrt", "[",
RowBox[{"Det", "[", "g", "]"}], "]"}], "//", "Simplify"}]}]}], "Input",
CellChangeTimes->{{3.908536677795877*^9, 3.908536735890113*^9}, {
3.908536782826435*^9, 3.908536804025585*^9}, 3.908536883968605*^9, {
3.908536934274844*^9, 3.908536953442299*^9}, {3.908536993822443*^9,
3.908537037389576*^9}, 3.9085374203476763`*^9},
CellLabel->"In[31]:=",ExpressionUUID->"749ea6e9-500d-45de-ac40-5a8e6d66bae0"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"r", " ",
RowBox[{"Cos", "[", "\[CurlyPhi]", "]"}]}], ",",
RowBox[{"r", " ",
RowBox[{"Sin", "[", "\[CurlyPhi]", "]"}]}]}], "}"}]], "Output",
CellChangeTimes->{
3.908536804844842*^9, 3.908536954069023*^9, {3.908537037947226*^9,
3.908537060163313*^9}, 3.908537420856242*^9},
CellLabel->"Out[31]=",ExpressionUUID->"0e258462-abb6-42d2-909e-9e3ec13c6c53"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{
RowBox[{"Cos", "[", "\[CurlyPhi]", "]"}], ",",
RowBox[{"Sin", "[", "\[CurlyPhi]", "]"}]}], "}"}], ",",
RowBox[{"{",
RowBox[{
RowBox[{
RowBox[{"-", "r"}], " ",
RowBox[{"Sin", "[", "\[CurlyPhi]", "]"}]}], ",",
RowBox[{"r", " ",
RowBox[{"Cos", "[", "\[CurlyPhi]", "]"}]}]}], "}"}]}], "}"}]], "Output",\
CellChangeTimes->{
3.908536804844842*^9, 3.908536954069023*^9, {3.908537037947226*^9,
3.908537060163313*^9}, 3.9085374208629513`*^9},
CellLabel->"Out[32]=",ExpressionUUID->"c807491c-cb7a-4848-a6ac-23b609c9b524"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"1", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",",
SuperscriptBox["r", "2"]}], "}"}]}], "}"}]], "Output",
CellChangeTimes->{
3.908536804844842*^9, 3.908536954069023*^9, {3.908537037947226*^9,
3.908537060163313*^9}, 3.908537420870977*^9},
CellLabel->"Out[33]=",ExpressionUUID->"3cbbe53f-bc46-402b-a476-5431244fbe7e"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"{",
RowBox[{"1", ",", "0"}], "}"}], ",",
RowBox[{"{",
RowBox[{"0", ",",
FractionBox["1",
SuperscriptBox["r", "2"]]}], "}"}]}], "}"}]], "Output",
CellChangeTimes->{
3.908536804844842*^9, 3.908536954069023*^9, {3.908537037947226*^9,
3.908537060163313*^9}, 3.908537420878971*^9},
CellLabel->"Out[34]=",ExpressionUUID->"f6d3fba2-b471-4115-932c-5f3211d7820e"],
Cell[BoxData["r"], "Output",
CellChangeTimes->{
3.908536804844842*^9, 3.908536954069023*^9, {3.908537037947226*^9,
3.908537060163313*^9}, 3.9085374208865843`*^9},
CellLabel->"Out[35]=",ExpressionUUID->"7d456091-063f-412d-adb3-5d648892e6b6"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{"vecS", "=",
RowBox[{"{",
RowBox[{"0", ",",
RowBox[{
RowBox[{"1", "/", "r"}], "+", "r"}]}], "}"}]}], "\[IndentingNewLine]",
RowBox[{"vecF", "=",
RowBox[{
RowBox[{"vecS", ".", "g"}], "//", "Simplify"}]}], "\[IndentingNewLine]",
RowBox[{"vecR", "=",
RowBox[{
RowBox[{"vecS", ".", "DX"}], "//",
"Simplify"}]}], "\[IndentingNewLine]", "foobar", "\[IndentingNewLine]",
RowBox[{"vecMUS", "=",
RowBox[{"Integrate", "[",
RowBox[{
RowBox[{"vecS", "*", "mu"}], ",",
RowBox[{"{",
RowBox[{"r", ",", "0", ",", "1"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\[CurlyPhi]", ",", "0", ",",
RowBox[{"2", "\[Pi]"}]}], "}"}]}], "]"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"vecMLF", "=",
RowBox[{"Integrate", "[",
RowBox[{
RowBox[{"vecF", "*", "mu"}], ",",
RowBox[{"{",
RowBox[{"r", ",", "0", ",", "1"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\[CurlyPhi]", ",", "0", ",",
RowBox[{"2", "\[Pi]"}]}], "}"}]}], "]"}]}],
";"}], "\[IndentingNewLine]",
RowBox[{"vecMLS", "=",
RowBox[{
RowBox[{"vecMLF", ".", "gInv"}], "//",
"Simplify"}]}], "\[IndentingNewLine]", "foobar", "\[IndentingNewLine]",
RowBox[{"vecMUR", "=",
RowBox[{
RowBox[{"vecMUS", ".", "DX"}], "//", "Simplify"}]}], "\[IndentingNewLine]",
RowBox[{"vecMLR", "=",
RowBox[{
RowBox[{"vecMLS", ".", "DX"}], "//", "Simplify"}]}], "\[IndentingNewLine]",
RowBox[{"vecMMR", "=",
RowBox[{"Integrate", "[",
RowBox[{
RowBox[{"vecR", "*", "mu"}], ",",
RowBox[{"{",
RowBox[{"r", ",", "0", ",", "1"}], "}"}], ",",
RowBox[{"{",
RowBox[{"\[CurlyPhi]", ",", "0", ",",
RowBox[{"2", "\[Pi]"}]}], "}"}]}], "]"}]}]}], "Input",
CellChangeTimes->{{3.9085370677503443`*^9, 3.908537395232176*^9}, {
3.908537436380686*^9, 3.9085374814871607`*^9}, {3.90853753363913*^9,
3.908537667685548*^9}, {3.908537988809556*^9, 3.908537990825014*^9}, {
3.908538062375228*^9, 3.908538083247596*^9}},
CellLabel->
"In[122]:=",ExpressionUUID->"a3153d31-308b-4081-a837-d78b3f1b207e"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"0", ",",
RowBox[{
FractionBox["1", "r"], "+", "r"}]}], "}"}]], "Output",
CellChangeTimes->{
3.908537148247705*^9, 3.908537246724028*^9, 3.908537352752882*^9,
3.908537395821574*^9, 3.908537482306665*^9, {3.90853754579412*^9,
3.908537551164577*^9}, {3.908537634592929*^9, 3.908537668296736*^9},
3.908537996566593*^9, {3.908538063423064*^9, 3.908538084495439*^9}},
CellLabel->
"Out[122]=",ExpressionUUID->"20bbf761-0741-4840-b2ea-ff9d690e8e8d"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"0", ",",
RowBox[{"r", "+",
SuperscriptBox["r", "3"]}]}], "}"}]], "Output",
CellChangeTimes->{
3.908537148247705*^9, 3.908537246724028*^9, 3.908537352752882*^9,
3.908537395821574*^9, 3.908537482306665*^9, {3.90853754579412*^9,
3.908537551164577*^9}, {3.908537634592929*^9, 3.908537668296736*^9},
3.908537996566593*^9, {3.908538063423064*^9, 3.908538084502393*^9}},
CellLabel->
"Out[123]=",ExpressionUUID->"e9b6fc32-7212-43f2-9f5a-8321b630ff86"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"-",
RowBox[{"(",
RowBox[{
RowBox[{"(",
RowBox[{"1", "+",
SuperscriptBox["r", "2"]}], ")"}], " ",
RowBox[{"Sin", "[", "\[CurlyPhi]", "]"}]}], ")"}]}], ",",
RowBox[{
RowBox[{"(",
RowBox[{"1", "+",
SuperscriptBox["r", "2"]}], ")"}], " ",
RowBox[{"Cos", "[", "\[CurlyPhi]", "]"}]}]}], "}"}]], "Output",
CellChangeTimes->{
3.908537148247705*^9, 3.908537246724028*^9, 3.908537352752882*^9,
3.908537395821574*^9, 3.908537482306665*^9, {3.90853754579412*^9,
3.908537551164577*^9}, {3.908537634592929*^9, 3.908537668296736*^9},
3.908537996566593*^9, {3.908538063423064*^9, 3.908538084510021*^9}},
CellLabel->
"Out[124]=",ExpressionUUID->"4e1ebc3c-5fd3-47e4-9fc8-422c4704d6aa"],
Cell[BoxData["foobar"], "Output",
CellChangeTimes->{
3.908537148247705*^9, 3.908537246724028*^9, 3.908537352752882*^9,
3.908537395821574*^9, 3.908537482306665*^9, {3.90853754579412*^9,
3.908537551164577*^9}, {3.908537634592929*^9, 3.908537668296736*^9},
3.908537996566593*^9, {3.908538063423064*^9, 3.9085380845176897`*^9}},
CellLabel->
"Out[125]=",ExpressionUUID->"0eaa44c3-657b-47b3-8085-af77fb4f72f2"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"0", ",",
FractionBox[
RowBox[{"8", " ", "\[Pi]"}], "3"]}], "}"}]], "Output",
CellChangeTimes->{
3.908537148247705*^9, 3.908537246724028*^9, 3.908537352752882*^9,
3.908537395821574*^9, 3.908537482306665*^9, {3.90853754579412*^9,
3.908537551164577*^9}, {3.908537634592929*^9, 3.908537668296736*^9},
3.908537996566593*^9, {3.908538063423064*^9, 3.9085380848485413`*^9}},
CellLabel->
"Out[126]=",ExpressionUUID->"76b4e534-de19-49bf-963b-f09a13633841"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"0", ",",
FractionBox[
RowBox[{"16", " ", "\[Pi]"}],
RowBox[{"15", " ",
SuperscriptBox["r", "2"]}]]}], "}"}]], "Output",
CellChangeTimes->{
3.908537148247705*^9, 3.908537246724028*^9, 3.908537352752882*^9,
3.908537395821574*^9, 3.908537482306665*^9, {3.90853754579412*^9,
3.908537551164577*^9}, {3.908537634592929*^9, 3.908537668296736*^9},
3.908537996566593*^9, {3.908538063423064*^9, 3.9085380850635138`*^9}},
CellLabel->
"Out[128]=",ExpressionUUID->"eab2d716-e921-4595-8825-8f2f64d0f87e"],
Cell[BoxData["foobar"], "Output",
CellChangeTimes->{
3.908537148247705*^9, 3.908537246724028*^9, 3.908537352752882*^9,
3.908537395821574*^9, 3.908537482306665*^9, {3.90853754579412*^9,
3.908537551164577*^9}, {3.908537634592929*^9, 3.908537668296736*^9},
3.908537996566593*^9, {3.908538063423064*^9, 3.908538085065731*^9}},
CellLabel->
"Out[129]=",ExpressionUUID->"56dde299-1bff-4d32-8852-95df0ebc1b94"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{
RowBox[{"-",
FractionBox["8", "3"]}], " ", "\[Pi]", " ", "r", " ",
RowBox[{"Sin", "[", "\[CurlyPhi]", "]"}]}], ",",
RowBox[{
FractionBox["8", "3"], " ", "\[Pi]", " ", "r", " ",
RowBox[{"Cos", "[", "\[CurlyPhi]", "]"}]}]}], "}"}]], "Output",
CellChangeTimes->{
3.908537148247705*^9, 3.908537246724028*^9, 3.908537352752882*^9,
3.908537395821574*^9, 3.908537482306665*^9, {3.90853754579412*^9,
3.908537551164577*^9}, {3.908537634592929*^9, 3.908537668296736*^9},
3.908537996566593*^9, {3.908538063423064*^9, 3.908538085067708*^9}},
CellLabel->
"Out[130]=",ExpressionUUID->"38a4182b-cc74-443b-80bb-949e276c691e"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"-",
FractionBox[
RowBox[{"16", " ", "\[Pi]", " ",
RowBox[{"Sin", "[", "\[CurlyPhi]", "]"}]}],
RowBox[{"15", " ", "r"}]]}], ",",
FractionBox[
RowBox[{"16", " ", "\[Pi]", " ",
RowBox[{"Cos", "[", "\[CurlyPhi]", "]"}]}],
RowBox[{"15", " ", "r"}]]}], "}"}]], "Output",
CellChangeTimes->{
3.908537148247705*^9, 3.908537246724028*^9, 3.908537352752882*^9,
3.908537395821574*^9, 3.908537482306665*^9, {3.90853754579412*^9,
3.908537551164577*^9}, {3.908537634592929*^9, 3.908537668296736*^9},
3.908537996566593*^9, {3.908538063423064*^9, 3.9085380850694923`*^9}},
CellLabel->
"Out[131]=",ExpressionUUID->"d31ad6d6-9665-47ab-b68f-24fb09b0826b"],
Cell[BoxData[
RowBox[{"{",
RowBox[{"0", ",", "0"}], "}"}]], "Output",
CellChangeTimes->{
3.908537148247705*^9, 3.908537246724028*^9, 3.908537352752882*^9,
3.908537395821574*^9, 3.908537482306665*^9, {3.90853754579412*^9,
3.908537551164577*^9}, {3.908537634592929*^9, 3.908537668296736*^9},
3.908537996566593*^9, {3.908538063423064*^9, 3.908538085547634*^9}},
CellLabel->
"Out[132]=",ExpressionUUID->"28f4f3a7-8fdb-4d45-92af-37191d82ef82"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"vecMUR", "+", "vecMLR"}], ")"}], "/", "2"}], "//",
"Simplify"}]], "Input",
CellChangeTimes->{{3.908537816537592*^9, 3.908537833361465*^9}},
CellLabel->
"In[110]:=",ExpressionUUID->"82847a7b-d6f0-4ad8-bed9-e5b078eb174e"],
Cell[BoxData[
RowBox[{"{",
RowBox[{
RowBox[{"-",
FractionBox[
RowBox[{"\[Pi]", " ",
RowBox[{"(",
RowBox[{"1", "+",
RowBox[{"3", " ",
SuperscriptBox["r", "2"]}]}], ")"}], " ",
RowBox[{"Sin", "[", "\[CurlyPhi]", "]"}]}],
RowBox[{"3", " ", "r"}]]}], ",",
RowBox[{"\[Pi]", " ",
RowBox[{"(",
RowBox[{
FractionBox["1",
RowBox[{"3", " ", "r"}]], "+", "r"}], ")"}], " ",
RowBox[{"Cos", "[", "\[CurlyPhi]", "]"}]}]}], "}"}]], "Output",
CellChangeTimes->{3.908537834632143*^9, 3.908538041202775*^9},
CellLabel->
"Out[110]=",ExpressionUUID->"ac5d2414-3db3-4311-931d-ff92e979f3fd"]
}, Open ]]
},
WindowSize->{1440., 773.25},
WindowMargins->{{-1440, Automatic}, {0, Automatic}},
FrontEndVersion->"13.2 for Linux x86 (64-bit) (January 30, 2023)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"49ce43f7-438b-4079-bd44-acfdb88c333a"
]
(* End of Notebook Content *)
(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 229, 4, 22, "Input",ExpressionUUID->"07d7ef80-8072-4d1f-adc2-6dc778f3a6d4"],
Cell[812, 28, 185, 3, 25, "Output",ExpressionUUID->"7e3882c2-5e75-4860-8715-7287c13aef00"]
}, Open ]],
Cell[CellGroupData[{
Cell[1034, 36, 1259, 34, 87, "Input",ExpressionUUID->"749ea6e9-500d-45de-ac40-5a8e6d66bae0"],
Cell[2296, 72, 423, 10, 28, "Output",ExpressionUUID->"0e258462-abb6-42d2-909e-9e3ec13c6c53"],
Cell[2722, 84, 641, 18, 28, "Output",ExpressionUUID->"c807491c-cb7a-4848-a6ac-23b609c9b524"],
Cell[3366, 104, 416, 11, 27, "Output",ExpressionUUID->"3cbbe53f-bc46-402b-a476-5431244fbe7e"],
Cell[3785, 117, 441, 12, 40, "Output",ExpressionUUID->"f6d3fba2-b471-4115-932c-5f3211d7820e"],
Cell[4229, 131, 248, 4, 25, "Output",ExpressionUUID->"7d456091-063f-412d-adb3-5d648892e6b6"]
}, Open ]],
Cell[CellGroupData[{
Cell[4514, 140, 2109, 57, 187, "Input",ExpressionUUID->"a3153d31-308b-4081-a837-d78b3f1b207e"],
Cell[6626, 199, 507, 11, 50, "Output",ExpressionUUID->"20bbf761-0741-4840-b2ea-ff9d690e8e8d"],
Cell[7136, 212, 510, 11, 38, "Output",ExpressionUUID->"e9b6fc32-7212-43f2-9f5a-8321b630ff86"],
Cell[7649, 225, 804, 21, 40, "Output",ExpressionUUID->"4e1ebc3c-5fd3-47e4-9fc8-422c4704d6aa"],
Cell[8456, 248, 422, 7, 36, "Output",ExpressionUUID->"0eaa44c3-657b-47b3-8085-af77fb4f72f2"],
Cell[8881, 257, 513, 11, 53, "Output",ExpressionUUID->"76b4e534-de19-49bf-963b-f09a13633841"],
Cell[9397, 270, 567, 13, 59, "Output",ExpressionUUID->"eab2d716-e921-4595-8825-8f2f64d0f87e"],
Cell[9967, 285, 420, 7, 36, "Output",ExpressionUUID->"56dde299-1bff-4d32-8852-95df0ebc1b94"],
Cell[10390, 294, 707, 16, 50, "Output",ExpressionUUID->"38a4182b-cc74-443b-80bb-949e276c691e"],
Cell[11100, 312, 747, 18, 59, "Output",ExpressionUUID->"d31ad6d6-9665-47ab-b68f-24fb09b0826b"],
Cell[11850, 332, 460, 9, 62, "Output",ExpressionUUID->"28f4f3a7-8fdb-4d45-92af-37191d82ef82"]
}, Open ]],
Cell[CellGroupData[{
Cell[12347, 346, 291, 8, 22, "Input",ExpressionUUID->"82847a7b-d6f0-4ad8-bed9-e5b078eb174e"],
Cell[12641, 356, 671, 20, 60, "Output",ExpressionUUID->"ac5d2414-3db3-4311-931d-ff92e979f3fd"]
}, Open ]]
}
]
*)
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment