(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 7.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 28021, 793] NotebookOptionsPosition[ 25765, 714] NotebookOutlinePosition[ 26162, 731] CellTagsIndexPosition[ 26119, 728] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["2. Klassische Mechanik und Computeralgebra", "Subtitle", CellChangeTimes->{{3.494908211331071*^9, 3.4949082312250957`*^9}, { 3.4978460915860558`*^9, 3.497846092295192*^9}}], Cell[CellGroupData[{ Cell["Kepler - Problem", "Section", CellChangeTimes->{{3.494908197098586*^9, 3.494908200099038*^9}, { 3.4978460944861717`*^9, 3.497846095978907*^9}, {3.5314631278878803`*^9, 3.531463128914563*^9}}], Cell["\<\ Betrachte Zwei - K\[ODoubleDot]rperproblem im Zentralfeld Hier : Planet mit Masse m bewegt sich um Sonne mit Masse M >> m\ \>", "Text", CellChangeTimes->{{3.494908259561584*^9, 3.494908328459826*^9}}], Cell[CellGroupData[{ Cell["Bewegungsgleichung", "Subsection", CellChangeTimes->{{3.49490833098764*^9, 3.494908334576943*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ SuperscriptBox["d", "2"], RowBox[{ OverscriptBox["r", "\[RightVector]"], " ", "/", " ", "d"}], " ", SuperscriptBox["t", "2"]}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"-", " ", RowBox[{"(", RowBox[{ RowBox[{"G", "*", RowBox[{"M", "/"}]}], "|", OverscriptBox["r", "\[RightVector]"], " ", "|", RowBox[{"^", "2"}]}], " ", ")"}]}], " ", RowBox[{ OverscriptBox["r", "\[RightVector]"], " ", "/"}]}], "|", OverscriptBox["r", "\[RightVector]"], " ", "|"}]}]], "Text", CellChangeTimes->{ 3.494908360236998*^9, {3.494908442027588*^9, 3.4949086984665937`*^9}}], Cell["\<\ Im Folgenden: G=1, M=1 Aufgrund der Drehimpulserhaltung findet die Bewegung in der Ebene statt \[DoubleLongRightArrow] 2 gekoppelte Differentialgleichungen 2. Ordnung\ \>", "Text", CellChangeTimes->{{3.494908721053933*^9, 3.494908802199401*^9}, { 3.494909097966279*^9, 3.4949091071560793`*^9}, {3.499483144375312*^9, 3.499483169075562*^9}}], Cell[BoxData[{ RowBox[{ OverscriptBox["x", "\[DoubleDot]"], " ", "=", " ", RowBox[{ RowBox[{"-", "x"}], "/", RowBox[{"r", "^", "3"}]}]}], "\[IndentingNewLine]", RowBox[{ OverscriptBox["y", "\[DoubleDot]"], " ", "=", " ", RowBox[{ RowBox[{"-", "y"}], "/", RowBox[{"r", "^", "3"}]}]}], "\[IndentingNewLine]", RowBox[{"r", " ", "=", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"x", "^", "2"}], "+", RowBox[{"y", "^", "2"}]}], "]"}]}]}], "Text", CellChangeTimes->{{3.4949089680084257`*^9, 3.4949090120730963`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["Erhaltungsgr\[ODoubleDot]\[SZ]en", "Subsection", CellChangeTimes->{{3.494909032825118*^9, 3.494909048317815*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Energie", ":", " ", "E"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"dx", "/", "dt"}], ")"}], "^", "2"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{"dy", "/", "dt"}], ")"}], "^", "2"}]}], ")"}], "/", "2"}], " ", "-", " ", RowBox[{"1", "/", "r"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"z", "-", RowBox[{"Komponente", " ", "des", " ", RowBox[{"Drehimpuls", ":", "Lz"}]}]}], " ", "=", " ", RowBox[{"m", RowBox[{"(", RowBox[{ RowBox[{"x", " ", RowBox[{"dy", "/", "dt"}]}], "-", RowBox[{"y", " ", RowBox[{"dx", "/", "dt"}]}]}], ")"}]}]}]}], "Text", CellChangeTimes->{{3.4949090526172457`*^9, 3.49490915696227*^9}, { 3.499570793231776*^9, 3.499570918690576*^9}, 3.5314648653026733`*^9}] }, Open ]], Cell[CellGroupData[{ Cell["3. Keplersches Gesetz", "Subsection", CellChangeTimes->{{3.494909198245583*^9, 3.494909204066551*^9}, { 3.499483241803776*^9, 3.4994832420965*^9}}], Cell["\<\ Die Quadrate der Umlaufzeiten verhalten sich wie die Kuben der gro\[SZ]en \ Halbache: T^2 ~ a^3\ \>", "Text", CellChangeTimes->{{3.49490922200858*^9, 3.4949092595430593`*^9}, { 3.4994348648953114`*^9, 3.4994348652279453`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["Implementierung", "Subsection", CellChangeTimes->{{3.4949095524767323`*^9, 3.494909555239992*^9}}], Cell[CellGroupData[{ Cell["Anfangsbedingungen", "Subsubsection"], Cell[BoxData[ RowBox[{"\n", RowBox[{ RowBox[{ RowBox[{"x0", " ", "=", " ", "0.5"}], ";", " ", RowBox[{"vx0", " ", "=", " ", "0"}], ";", " ", RowBox[{"y0", " ", "=", " ", "0"}], ";", " ", RowBox[{"vy0", " ", "=", " ", "1.2"}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"**", RowBox[{"*", "\[IndentingNewLine]", "x0"}]}], " ", "=", " ", "0.5"}], ";", " ", RowBox[{"vx0", " ", "=", " ", "0"}], ";", " ", RowBox[{"y0", " ", "=", " ", "0"}], ";", " ", RowBox[{"vy0", " ", "=", " ", "0.8"}], ";", "\[IndentingNewLine]", " ", RowBox[{"x0", " ", "=", " ", "0.5"}], ";", " ", RowBox[{"vx0", " ", "=", " ", "0"}], ";", " ", RowBox[{"y0", " ", "=", " ", "0"}], ";", " ", RowBox[{"vy0", " ", "=", " ", "1.2"}], ";"}], "\[IndentingNewLine]", "**"}], "*"}], " ", "*)"}]}]}]], "Input", CellChangeTimes->{ 3.49490958906889*^9, 3.494909619198526*^9, {3.494909652011972*^9, 3.494909786288755*^9}, {3.4949120172502003`*^9, 3.494912021759781*^9}, { 3.494912087403051*^9, 3.494912087478599*^9}, {3.494912139507712*^9, 3.494912139817025*^9}, {3.499435118127083*^9, 3.4994351210449266`*^9}, { 3.499435175121997*^9, 3.499435176230954*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["Numerische Loesung der beiden gekoppelten DGLen :", "Subsubsection"], Cell[BoxData[ RowBox[{ RowBox[{"erg", " ", "=", " ", "\n", " ", RowBox[{"NDSolve", "[", RowBox[{ RowBox[{"{", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"x", "[", "t", "]"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", " ", "2"}], "}"}]}], "]"}], " ", "==", " ", RowBox[{ RowBox[{"-", RowBox[{"x", "[", "t", "]"}]}], "/", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"x", "[", "t", "]"}], "^", "2"}], " ", "+", " ", RowBox[{ RowBox[{"y", "[", "t", "]"}], "^", "2"}]}], ")"}], "^", RowBox[{"(", RowBox[{"3", "/", "2"}], ")"}]}]}]}], ",", "\n", " ", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{"y", "[", "t", "]"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", " ", "2"}], "}"}]}], "]"}], " ", "==", " ", RowBox[{ RowBox[{"-", RowBox[{"y", "[", "t", "]"}]}], "/", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"x", "[", "t", "]"}], "^", "2"}], " ", "+", " ", RowBox[{ RowBox[{"y", "[", "t", "]"}], "^", "2"}]}], ")"}], "^", RowBox[{"(", RowBox[{"3", "/", "2"}], ")"}]}]}]}], ",", "\n", " ", RowBox[{ RowBox[{"x", "[", "0", "]"}], " ", "==", " ", "x0"}], ",", " ", RowBox[{ RowBox[{"y", "[", "0", "]"}], " ", "==", " ", "y0"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"x", "'"}], "[", "0", "]"}], " ", "==", " ", "vx0"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"y", "'"}], "[", "0", "]"}], " ", "==", " ", "vy0"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{ RowBox[{"x", "[", "t", "]"}], ",", " ", RowBox[{"y", "[", "t", "]"}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", " ", "0", ",", " ", "10"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{ 3.49490958906889*^9, 3.494909619198526*^9, {3.494909652011972*^9, 3.4949096932857523`*^9}, {3.494909751683422*^9, 3.49490975454599*^9}, { 3.494909810359324*^9, 3.494909861509883*^9}, {3.4994833165163517`*^9, 3.499483317153082*^9}}], Cell[BoxData["erg"], "Input", CellChangeTimes->{ 3.49490958906889*^9, 3.494909619198526*^9, {3.494909652011972*^9, 3.4949096932857523`*^9}, {3.494909751683422*^9, 3.49490975454599*^9}, { 3.494909810359324*^9, 3.4949098668739*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["Graphische Darstellung von x(t) und y(t)", "Subsubsection", CellChangeTimes->{{3.494909891503182*^9, 3.494909912809421*^9}}], Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"x", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}], ",", " ", RowBox[{ RowBox[{"y", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}]}], "}"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", " ", "0", ",", " ", "3"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{ 3.49490958906889*^9, 3.494909619198526*^9, {3.494909652011972*^9, 3.4949096932857523`*^9}, {3.494909751683422*^9, 3.49490975454599*^9}, { 3.494909810359324*^9, 3.4949098824803247`*^9}, {3.494909918144837*^9, 3.494909922642071*^9}, 3.4949099530651617`*^9}], Cell[BoxData[ RowBox[{"parpl", " ", "=", " ", RowBox[{"ParametricPlot", "[", "\n", " ", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"x", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}], ",", " ", RowBox[{ RowBox[{"y", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}]}], "}"}], "\[IndentingNewLine]", ",", RowBox[{"{", RowBox[{"t", ",", " ", "0", ",", " ", "5"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{ 3.49490958906889*^9, 3.494909619198526*^9, {3.494909652011972*^9, 3.4949096932857523`*^9}, {3.494909751683422*^9, 3.49490975454599*^9}, { 3.494909810359324*^9, 3.4949098824803247`*^9}, {3.494909918144837*^9, 3.494909922642071*^9}, {3.4949099530651617`*^9, 3.494910037177412*^9}, { 3.4994834369967413`*^9, 3.49948343835669*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["Definition des \"Planeten\" als Graphikobjekt :", "Subsubsection", CellChangeTimes->{{3.494910046957212*^9, 3.4949100489044437`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"planet", "[", RowBox[{"x_", ",", " ", "y_"}], "]"}], " ", ":=", " ", RowBox[{"Graphics", "[", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{"x", ",", " ", "y"}], "}"}], ",", " ", ".01"}], "]"}], "]"}]}], ";"}]], "Input", CellChangeTimes->{ 3.49490958906889*^9, 3.494909619198526*^9, {3.494909652011972*^9, 3.4949096932857523`*^9}, {3.494909751683422*^9, 3.49490975454599*^9}, { 3.494909810359324*^9, 3.4949098824803247`*^9}, {3.494909918144837*^9, 3.494909922642071*^9}, {3.4949099530651617`*^9, 3.494910060534709*^9}, 3.494910132964657*^9}] }, Open ]], Cell[CellGroupData[{ Cell["Definiere Norm - Koordinatensystem :", "Subsubsection", CellChangeTimes->{ 3.49490958906889*^9, 3.494909619198526*^9, {3.494909652011972*^9, 3.4949096932857523`*^9}, {3.494909751683422*^9, 3.49490975454599*^9}, { 3.494909810359324*^9, 3.4949098824803247`*^9}, {3.494909918144837*^9, 3.494909922642071*^9}, {3.4949099530651617`*^9, 3.494910060534709*^9}, { 3.494910132964657*^9, 3.494910142084738*^9}}], Cell[BoxData[{ RowBox[{"pl", " ", ":=", " ", RowBox[{"Plot", "[", RowBox[{"0", ",", " ", RowBox[{"{", RowBox[{"x", ",", " ", RowBox[{"-", ".5"}], ",", " ", ".5"}], "}"}], ",", RowBox[{"PlotRange", " ", "->", " ", RowBox[{"{", RowBox[{ RowBox[{"-", ".5"}], ",", " ", ".5"}], "}"}]}]}], "]"}]}], "\n", RowBox[{"Show", "[", RowBox[{"pl", ",", "parpl", ",", " ", RowBox[{"planet", "[", RowBox[{"x0", ",", " ", "y0"}], "]"}]}], "]"}]}], "Input", CellChangeTimes->{ 3.49490958906889*^9, 3.494909619198526*^9, {3.494909652011972*^9, 3.4949096932857523`*^9}, {3.494909751683422*^9, 3.49490975454599*^9}, { 3.494909810359324*^9, 3.4949098824803247`*^9}, {3.494909918144837*^9, 3.494909922642071*^9}, {3.4949099530651617`*^9, 3.494910060534709*^9}, { 3.494910132964657*^9, 3.494910161028717*^9}, {3.499434904486491*^9, 3.4994349051857567`*^9}, {3.499435014065721*^9, 3.4994350176082487`*^9}, { 3.4994350600096684`*^9, 3.499435064335682*^9}, {3.49948419035607*^9, 3.499484195284699*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["Definiere Energie und z - Komponente des Drehimpulses", "Subsubsection", CellChangeTimes->{3.4994352344405737`*^9}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"energie", "[", "t_", "]"}], "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}], ",", " ", "t"}], "]"}], "^", "2"}], " ", "+", " ", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{ RowBox[{"y", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}], ",", " ", "t"}], "]"}], "^", "2"}]}], ")"}], "/", "2"}], " ", "-", "\[IndentingNewLine]", RowBox[{"1", "/", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}], ")"}], "^", "2"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}]}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"lz", "[", "t_", "]"}], " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}], ")"}], "*", RowBox[{"D", "[", RowBox[{ RowBox[{ RowBox[{"y", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}], ",", " ", "t"}], "]"}]}], " ", "-", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}], ")"}], "*", " ", RowBox[{"D", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}], ",", " ", "t"}], "]"}]}]}]}], " ", ";"}]}], "Input", CellChangeTimes->{ 3.49490958906889*^9, 3.494909619198526*^9, {3.494909652011972*^9, 3.4949096932857523`*^9}, {3.494909751683422*^9, 3.49490975454599*^9}, { 3.494909810359324*^9, 3.4949098824803247`*^9}, {3.494909918144837*^9, 3.494909922642071*^9}, {3.4949099530651617`*^9, 3.494910019911171*^9}, { 3.494910171079239*^9, 3.494910243663876*^9}, {3.494910278183054*^9, 3.494910320808909*^9}, {3.494910361893346*^9, 3.494910370064024*^9}, { 3.494910541678844*^9, 3.4949105486835012`*^9}, {3.494910860169392*^9, 3.494910870733488*^9}, {3.4949109464482927`*^9, 3.494910976504025*^9}, { 3.494911406623733*^9, 3.4949114467356663`*^9}}], Cell[BoxData[{ RowBox[{"energie", "[", ".3", "]"}], "\[IndentingNewLine]", RowBox[{"lz", "[", ".3", "]"}]}], "Input", CellChangeTimes->{{3.494910917108083*^9, 3.494910919516596*^9}, { 3.494910962658945*^9, 3.494910990781973*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Animierte Darstellung der L\[ODoubleDot]sung\ \>", "Subsubsection", CellChangeTimes->{{3.494911122985497*^9, 3.494911130901196*^9}, { 3.494911639254579*^9, 3.4949116447935133`*^9}, {3.499435070725974*^9, 3.4994350721853733`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"graph", "[", "tt_", "]"}], " ", ":=", " ", "\n", " ", RowBox[{"Show", "[", RowBox[{"pl", ",", "parpl", ",", " ", "\n", " ", RowBox[{ RowBox[{"planet", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}], ",", " ", RowBox[{ RowBox[{"y", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}]}], "]"}], " ", "/.", " ", RowBox[{"{", RowBox[{"t", " ", "->", " ", "tt"}], "}"}]}], ",", " ", "\n", " ", RowBox[{"PlotLabel", " ", "->", " ", "\n", " ", RowBox[{"\"\\"", "<>", RowBox[{"ToString", "[", RowBox[{"PaddedForm", "[", RowBox[{ RowBox[{"N", "[", "tt", "]"}], ",", " ", RowBox[{"{", RowBox[{"7", ",", " ", "5"}], "}"}]}], "]"}], "]"}], " ", "<>", " ", "\"\< -- \>\"", " ", "<>", " ", "\n", " ", "\"\\"", " ", "<>", RowBox[{"ToString", "[", RowBox[{"PaddedForm", "[", RowBox[{ RowBox[{"energie", "[", "tt", "]"}], ",", " ", RowBox[{"{", RowBox[{"10", ",", " ", "7"}], "}"}]}], "]"}], "]"}], " ", "<>", " ", "\"\< -- \>\"", " ", "<>", " ", "\n", " ", "\"\\"", "<>", RowBox[{"ToString", "[", RowBox[{"PaddedForm", "[", RowBox[{ RowBox[{"lz", "[", "tt", "]"}], ",", " ", RowBox[{"{", RowBox[{"10", ",", " ", "7"}], "}"}]}], "]"}], "]"}]}]}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{ 3.49490958906889*^9, 3.494909619198526*^9, {3.494909652011972*^9, 3.4949096932857523`*^9}, {3.494909751683422*^9, 3.49490975454599*^9}, { 3.494909810359324*^9, 3.4949098824803247`*^9}, {3.494909918144837*^9, 3.494909922642071*^9}, {3.4949099530651617`*^9, 3.494910019911171*^9}, { 3.494910171079239*^9, 3.494910243663876*^9}, 3.4949104449058933`*^9, { 3.4949111381139383`*^9, 3.494911140454699*^9}, {3.49491119937873*^9, 3.494911284585759*^9}, {3.499434926966313*^9, 3.499434927057315*^9}, { 3.499434978442395*^9, 3.499434988200617*^9}, 3.499435053348344*^9, { 3.499435090149701*^9, 3.499435091833229*^9}}], Cell[BoxData[ RowBox[{"graph", "[", ".3", "]"}]], "Input", CellChangeTimes->{{3.4949111733636103`*^9, 3.494911180281497*^9}}], Cell[BoxData[ RowBox[{"\n", RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"graph", "[", "t", "]"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", " ", "0", ",", " ", "6", ",", " ", "\n", " ", RowBox[{"AppearanceElements", " ", "->", " ", RowBox[{"{", RowBox[{ "\"\\"", ",", " ", "\"\\"", ",", " ", "\n", " ", "\"\\"", ",", " ", "\"\\""}], "}"}]}]}], "}"}], ",", " ", "\[IndentingNewLine]", RowBox[{"ControlType", " ", "->", " ", "Animator"}]}], "]"}]}]], "Input", CellChangeTimes->{ 3.49490958906889*^9, 3.494909619198526*^9, {3.494909652011972*^9, 3.4949096932857523`*^9}, {3.494909751683422*^9, 3.49490975454599*^9}, { 3.494909810359324*^9, 3.4949098824803247`*^9}, {3.494909918144837*^9, 3.494909922642071*^9}, {3.4949099530651617`*^9, 3.494910019911171*^9}, { 3.494910171079239*^9, 3.494910243663876*^9}, {3.4949104449058933`*^9, 3.4949104471795807`*^9}, {3.4949113418795023`*^9, 3.494911347203251*^9}, { 3.4949115286022*^9, 3.494911530689999*^9}, {3.494920518804624*^9, 3.494920521371744*^9}, {3.49943527015025*^9, 3.4994352717490683`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Bestimme gro\[SZ]e Halbachse: Suche Minimum der Geschwindigkeit\ \>", "Subsubsection", CellChangeTimes->{{3.494911572512436*^9, 3.494911583199009*^9}, { 3.4949116613671722`*^9, 3.4949116800271397`*^9}}], Cell[BoxData[ RowBox[{"\n", RowBox[{ RowBox[{ RowBox[{ RowBox[{"v", "[", "t_", "]"}], " ", "=", " ", "\n", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{ RowBox[{"x", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}], ",", " ", "t"}], "]"}], "^", "2"}], " ", "+", " ", RowBox[{ RowBox[{"D", "[", RowBox[{ RowBox[{ RowBox[{"y", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}], ",", " ", "t"}], "]"}], "^", "2"}]}], "]"}]}], " ", ";"}], "\n", RowBox[{"v", "[", "0", "]"}], "\n", RowBox[{"Plot", "[", RowBox[{ RowBox[{"v", "[", "t1", "]"}], ",", " ", RowBox[{"{", RowBox[{"t1", ",", " ", "0", ",", " ", "2"}], "}"}]}], "]"}]}]}]], "Input", CellChangeTimes->{ 3.49490958906889*^9, 3.494909619198526*^9, {3.494909652011972*^9, 3.4949096932857523`*^9}, {3.494909751683422*^9, 3.49490975454599*^9}, { 3.494909810359324*^9, 3.4949098824803247`*^9}, {3.494909918144837*^9, 3.494909922642071*^9}, {3.4949099530651617`*^9, 3.494910019911171*^9}, { 3.494910171079239*^9, 3.494910243663876*^9}, {3.4949104449058933`*^9, 3.4949104471795807`*^9}, {3.4949113418795023`*^9, 3.494911347203251*^9}, { 3.494911670477119*^9, 3.494911702344118*^9}}], Cell["Ableitung der Geschwindigkeit; bestimme Periode", "Text", CellChangeTimes->{{3.494911741843398*^9, 3.494911779843813*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"dv", "[", "t_", "]"}], " ", "=", " ", RowBox[{"D", "[", RowBox[{ RowBox[{"v", "[", "t", "]"}], ",", " ", "t"}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"t1", " ", "=", " ", RowBox[{"t", " ", "/.", " ", RowBox[{ RowBox[{"FindRoot", "[", RowBox[{ RowBox[{ RowBox[{"dv", "[", "t", "]"}], " ", "==", " ", "0"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", " ", "0"}], "}"}]}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\n", RowBox[{ RowBox[{"t2", " ", "=", " ", RowBox[{"t", " ", "/.", " ", RowBox[{ RowBox[{"FindRoot", "[", RowBox[{ RowBox[{ RowBox[{"dv", "[", "t", "]"}], " ", "==", " ", "0"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", " ", "1."}], "}"}]}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\n", RowBox[{"tPeriode", " ", "=", " ", RowBox[{"t2", " ", "-", " ", "t1"}]}]}], "Input", CellChangeTimes->{ 3.49490958906889*^9, 3.494909619198526*^9, {3.494909652011972*^9, 3.4949096932857523`*^9}, {3.494909751683422*^9, 3.49490975454599*^9}, { 3.494909810359324*^9, 3.4949098824803247`*^9}, {3.494909918144837*^9, 3.494909922642071*^9}, {3.4949099530651617`*^9, 3.494910019911171*^9}, { 3.494910171079239*^9, 3.494910243663876*^9}, {3.4949104449058933`*^9, 3.4949104471795807`*^9}, {3.4949113418795023`*^9, 3.494911347203251*^9}, { 3.494911670477119*^9, 3.494911733993441*^9}, {3.494911783979638*^9, 3.494911790678986*^9}}], Cell["Gro\[SZ]e Halbachse :", "Text", CellChangeTimes->{{3.499435346140337*^9, 3.499435358603764*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"r", "[", "t_", "]"}], " ", "=", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}], ")"}], "^", "2"}], " ", "+", " ", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y", "[", "t", "]"}], " ", "/.", " ", RowBox[{"erg", "[", RowBox[{"[", "1", "]"}], "]"}]}], " ", ")"}], "^", "2"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"a", " ", "=", " ", RowBox[{ RowBox[{"1", "/", "2"}], "*", RowBox[{"(", RowBox[{ RowBox[{"r", "[", "t1", "]"}], " ", "+", " ", RowBox[{"r", "[", RowBox[{"t1", " ", "+", " ", RowBox[{"tPeriode", "/", "2"}]}], "]"}]}], ")"}]}]}], ";"}], "\n", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", " ", "tPeriode"}], "]"}], ";"}], "\n", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", " ", "a"}], "]"}], ";"}], "\n", RowBox[{ RowBox[{"Print", "[", RowBox[{"\"\\"", ",", " ", RowBox[{ RowBox[{"tPeriode", "^", "2"}], "/", RowBox[{"a", "^", "3"}]}]}], "]"}], ";"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{ 3.49490958906889*^9, 3.494909619198526*^9, {3.494909652011972*^9, 3.4949096932857523`*^9}, {3.494909751683422*^9, 3.49490975454599*^9}, { 3.494909810359324*^9, 3.4949098824803247`*^9}, {3.494909918144837*^9, 3.494909922642071*^9}, {3.4949099530651617`*^9, 3.494910019911171*^9}, { 3.494910171079239*^9, 3.494910243663876*^9}, {3.4949104449058933`*^9, 3.4949104471795807`*^9}, {3.4949113418795023`*^9, 3.494911347203251*^9}, { 3.494911670477119*^9, 3.494911733993441*^9}, {3.494911783979638*^9, 3.494911892095868*^9}, {3.494911927257243*^9, 3.494911931238872*^9}, 3.4949119811118813`*^9}] }, Open ]] }, Open ]] }, Open ]] }, Open ]] }, ScreenStyleEnvironment->"Presentation", WindowSize->{1154, 907}, WindowMargins->{{24, Automatic}, {Automatic, 0}}, ShowSelection->True, FrontEndVersion->"8.0 for Linux x86 (64-bit) (October 10, 2011)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[567, 22, 181, 2, 54, "Subtitle"], Cell[CellGroupData[{ Cell[773, 28, 203, 3, 74, "Section"], Cell[979, 33, 211, 4, 49, "Text"], Cell[CellGroupData[{ Cell[1215, 41, 105, 1, 39, "Subsection"], Cell[1323, 44, 676, 20, 88, "Text"], Cell[2002, 66, 357, 7, 185, "Text"], Cell[2362, 75, 556, 16, 124, "Text"] }, Open ]], Cell[CellGroupData[{ Cell[2955, 96, 120, 1, 67, "Subsection"], Cell[3078, 99, 894, 27, 87, "Text"] }, Open ]], Cell[CellGroupData[{ Cell[4009, 131, 156, 2, 67, "Subsection"], Cell[4168, 135, 239, 6, 117, "Text"] }, Open ]], Cell[CellGroupData[{ Cell[4444, 146, 105, 1, 67, "Subsection"], Cell[CellGroupData[{ Cell[4574, 151, 43, 0, 48, "Subsubsection"], Cell[4620, 153, 1335, 30, 238, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[5992, 188, 74, 0, 48, "Subsubsection"], Cell[6069, 190, 2403, 66, 238, "Input"], Cell[8475, 258, 240, 4, 52, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[8752, 267, 131, 1, 48, "Subsubsection"], Cell[8886, 270, 760, 19, 52, "Input"], Cell[9649, 291, 925, 21, 127, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[10611, 317, 140, 1, 48, "Subsubsection"], Cell[10754, 320, 648, 16, 52, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[11439, 341, 424, 6, 48, "Subsubsection"], Cell[11866, 349, 1069, 23, 90, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[12972, 377, 122, 1, 48, "Subsubsection"], Cell[13097, 380, 2852, 79, 423, "Input"], Cell[15952, 461, 235, 4, 50, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[16224, 470, 246, 5, 28, "Subsubsection"], Cell[16473, 477, 2347, 52, 145, "Input"], Cell[18823, 531, 127, 2, 30, "Input"], Cell[18953, 535, 1232, 23, 107, "Input"] }, Open ]], Cell[CellGroupData[{ Cell[20222, 563, 215, 4, 28, "Subsubsection"], Cell[20440, 569, 1485, 38, 107, "Input"], Cell[21928, 609, 129, 1, 30, "Text"], Cell[22060, 612, 1569, 39, 88, "Input"], Cell[23632, 653, 103, 1, 30, "Text"], Cell[23738, 656, 1975, 52, 126, "Input"] }, Open ]] }, Open ]] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)