(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 10.1' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 158, 7] NotebookDataLength[ 16852, 441] NotebookOptionsPosition[ 16272, 415] NotebookOutlinePosition[ 16616, 430] CellTagsIndexPosition[ 16573, 427] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["\<\ Animation of the particle\[CloseCurlyQuote]s movement and velocity vector (in \ green)\ \>", "Section", CellChangeTimes->{{3.7492883389002237`*^9, 3.7492883998897123`*^9}, { 3.7492884791402445`*^9, 3.7492884843685436`*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"\[Omega]", "=", "0.45"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Rho]", "=", "3"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"r", "[", "t_", "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{"\[Rho]", " ", RowBox[{"Cos", "[", RowBox[{"\[Omega]", " ", "t"}], "]"}]}], ",", RowBox[{"\[Rho]", " ", RowBox[{"Sin", "[", RowBox[{"\[Omega]", " ", "t"}], "]"}]}]}], "}"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"r", "[", "l_", "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{"\[Rho]", " ", RowBox[{"Cos", "[", RowBox[{"\[Omega]", " ", RowBox[{"l", "/", "2"}]}], "]"}]}], ",", RowBox[{"\[Rho]", " ", RowBox[{"Sin", "[", RowBox[{"\[Omega]", " ", RowBox[{"l", "/", "2"}]}], "]"}]}]}], "}"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Manipulate", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"Quiet", "[", "\[IndentingNewLine]", RowBox[{"Show", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Arrow", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"r", "[", "k", "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"Green", ",", RowBox[{"Arrow", "[", RowBox[{"{", RowBox[{ RowBox[{"r", "[", "k", "]"}], ",", RowBox[{ RowBox[{"r", "[", "k", "]"}], "+", RowBox[{ RowBox[{"r", "'"}], "[", "k", "]"}]}]}], "}"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", " ", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"r", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "50"}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"r", "[", "l", "]"}], ",", RowBox[{"{", RowBox[{"l", ",", "0", ",", "100"}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"r", "[", "k", "]"}], "}"}], ",", RowBox[{"PlotMarkers", "\[Rule]", "\[FilledCircle]"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Green"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "}"}]}]}], "\[IndentingNewLine]", "]"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"k", ",", "0", ",", "100", ",", "0.1"}], "}"}]}], "\[IndentingNewLine]", "]"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.7488447940972614`*^9, 3.748845102364893*^9}, { 3.7488451337506886`*^9, 3.7488452594558783`*^9}, {3.7488452905226555`*^9, 3.748845438259105*^9}, 3.74884548117556*^9, {3.748845512835371*^9, 3.748845546085273*^9}, {3.7488456125240726`*^9, 3.7488456595247607`*^9}, { 3.7488456995530505`*^9, 3.7488457045443363`*^9}, {3.748845745857699*^9, 3.748845761457591*^9}, {3.749288194074991*^9, 3.7492882021014504`*^9}, { 3.749288300752042*^9, 3.7492883050412865`*^9}, {3.7492884032619047`*^9, 3.7492884584260597`*^9}}], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`k$$ = 13.3, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`k$$], 0, 100, 0.1}}, Typeset`size$$ = { 360., {178., 182.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`k$273097$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`k$$ = 0}, "ControllerVariables" :> { Hold[$CellContext`k$$, $CellContext`k$273097$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> Quiet[ Show[ Graphics[{ Arrow[{{0, 0}, $CellContext`r[$CellContext`k$$]}], {Green, Arrow[{ $CellContext`r[$CellContext`k$$], \ $CellContext`r[$CellContext`k$$] + Derivative[1][$CellContext`r][$CellContext`k$$]}]}}], ParametricPlot[ $CellContext`r[$CellContext`t], {$CellContext`t, 0, 50}], ParametricPlot[ $CellContext`r[$CellContext`l], {$CellContext`l, 0, 100}], ListPlot[{ $CellContext`r[$CellContext`k$$]}, PlotMarkers -> $CellContext`\[FilledCircle], PlotStyle -> Green], PlotRange -> {{-4, 4}, {-4, 4}}]], "Specifications" :> {{$CellContext`k$$, 0, 100, 0.1}}, "Options" :> {}, "DefaultOptions" :> {}], ImageSizeCache->{411., {236., 243.}}, 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.749288442318139*^9, 3.749288458544067*^9}}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Particle, local basis (in red), velocity vector (in green) Remark: use the \[OpenCurlyDoubleQuote]+\[OpenCurlyDoubleQuote] button \ instead of running the animation, because animation takes long to compute.\ \>", "Section", CellChangeTimes->{{3.7492883389002237`*^9, 3.749288363575635*^9}, { 3.7492885008334856`*^9, 3.749288564270042*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"\[Omega]", "=", "0.75"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Rho]", "=", "3"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"r", "[", "t_", "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{"\[Rho]", " ", RowBox[{"Cos", "[", RowBox[{"\[Omega]", " ", "t"}], "]"}]}], ",", RowBox[{"\[Rho]", " ", RowBox[{"Sin", "[", RowBox[{"\[Omega]", " ", "t"}], "]"}]}]}], "}"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"r", "[", "l_", "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{"\[Rho]", " ", RowBox[{"Cos", "[", RowBox[{"\[Omega]", " ", RowBox[{"l", "/", "2"}]}], "]"}]}], ",", RowBox[{"\[Rho]", " ", RowBox[{"Sin", "[", RowBox[{"\[Omega]", " ", RowBox[{"l", "/", "2"}]}], "]"}]}]}], "}"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Manipulate", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"Quiet", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"T", "[", "k_", "]"}], ":=", RowBox[{ RowBox[{ RowBox[{"r", "'"}], "[", "k", "]"}], "/", SqrtBox[ RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"r", "'"}], "[", "k", "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "^", "2"}], "+", RowBox[{ RowBox[{ RowBox[{ RowBox[{"r", "'"}], "[", "k", "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], "^", "2"}]}]]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"n", "[", "k_", "]"}], ":=", RowBox[{ RowBox[{ RowBox[{"T", "'"}], "[", "k", "]"}], "/", SqrtBox[ RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"T", "'"}], "[", "k", "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], "^", "2"}], "+", RowBox[{ RowBox[{ RowBox[{ RowBox[{"T", "'"}], "[", "k", "]"}], "[", RowBox[{"[", "2", "]"}], "]"}], "^", "2"}]}]]}]}], ";", "\[IndentingNewLine]", RowBox[{"Show", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Arrow", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"r", "[", "k", "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"Green", ",", RowBox[{"Arrow", "[", RowBox[{"{", RowBox[{ RowBox[{"r", "[", "k", "]"}], ",", RowBox[{ RowBox[{"r", "[", "k", "]"}], "+", RowBox[{ RowBox[{"r", "'"}], "[", "k", "]"}]}]}], "}"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"Red", ",", RowBox[{"Arrow", "[", RowBox[{"{", RowBox[{ RowBox[{"r", "[", "k", "]"}], ",", RowBox[{ RowBox[{"r", "[", "k", "]"}], "+", RowBox[{"T", "[", "k", "]"}]}]}], "}"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"Red", ",", RowBox[{"Arrow", "[", RowBox[{"{", RowBox[{ RowBox[{"r", "[", "k", "]"}], ",", RowBox[{ RowBox[{"r", "[", "k", "]"}], "+", RowBox[{"n", "[", "k", "]"}]}]}], "}"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", " ", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"r", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "50"}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"r", "[", "l", "]"}], ",", RowBox[{"{", RowBox[{"l", ",", "0", ",", "100"}], "}"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"r", "[", "k", "]"}], "}"}], ",", RowBox[{"PlotMarkers", "\[Rule]", "\[FilledCircle]"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Green"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "}"}]}]}], "\[IndentingNewLine]", "]"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{"k", ",", "0", ",", "100", ",", "0.1"}], "}"}]}], "\[IndentingNewLine]", "]"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.7488447940972614`*^9, 3.748845102364893*^9}, { 3.7488451337506886`*^9, 3.7488452594558783`*^9}, {3.7488452905226555`*^9, 3.748845438259105*^9}, 3.74884548117556*^9, {3.748845512835371*^9, 3.748845546085273*^9}, {3.7488456125240726`*^9, 3.7488456595247607`*^9}, { 3.7488456995530505`*^9, 3.7488457045443363`*^9}, {3.748845745857699*^9, 3.748845761457591*^9}, {3.749288194074991*^9, 3.7492882021014504`*^9}, { 3.749288300752042*^9, 3.7492883050412865`*^9}}], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`k$$ = 0, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`k$$], 0, 100, 0.1}}, Typeset`size$$ = { 360., {178., 182.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`k$357193$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`k$$ = 0}, "ControllerVariables" :> { Hold[$CellContext`k$$, $CellContext`k$357193$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> Quiet[$CellContext`T[ Pattern[$CellContext`k, Blank[]]] := Derivative[1][$CellContext`r][$CellContext`k]/(Part[ Derivative[1][$CellContext`r][$CellContext`k], 1]^2 + Part[ Derivative[1][$CellContext`r][$CellContext`k], 2]^2)^ Rational[1, 2]; $CellContext`n[ Pattern[$CellContext`k, Blank[]]] := Derivative[1][$CellContext`T][$CellContext`k]/(Part[ Derivative[1][$CellContext`T][$CellContext`k], 1]^2 + Part[ Derivative[1][$CellContext`T][$CellContext`k], 2]^2)^ Rational[1, 2]; Show[ Graphics[{ Arrow[{{0, 0}, $CellContext`r[$CellContext`k$$]}], {Green, Arrow[{ $CellContext`r[$CellContext`k$$], \ $CellContext`r[$CellContext`k$$] + Derivative[1][$CellContext`r][$CellContext`k$$]}]}, {Red, Arrow[{ $CellContext`r[$CellContext`k$$], \ $CellContext`r[$CellContext`k$$] + $CellContext`T[$CellContext`k$$]}]}, {Red, Arrow[{ $CellContext`r[$CellContext`k$$], \ $CellContext`r[$CellContext`k$$] + $CellContext`n[$CellContext`k$$]}]}}], ParametricPlot[ $CellContext`r[$CellContext`t], {$CellContext`t, 0, 50}], ParametricPlot[ $CellContext`r[$CellContext`l], {$CellContext`l, 0, 100}], ListPlot[{ $CellContext`r[$CellContext`k$$]}, PlotMarkers -> $CellContext`\[FilledCircle], PlotStyle -> Green], PlotRange -> {{-4, 4}, {-4, 4}}]], "Specifications" :> {{$CellContext`k$$, 0, 100, 0.1}}, "Options" :> {}, "DefaultOptions" :> {}], ImageSizeCache->{411., {223., 230.}}, 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.748845170237775*^9, 3.7488452600379114`*^9}, { 3.748845349247014*^9, 3.7488453592665873`*^9}, {3.748845398627838*^9, 3.748845409473459*^9}, 3.7488454398851986`*^9, 3.748845482033609*^9, { 3.748845514813484*^9, 3.7488455470143256`*^9}, 3.7488455889557247`*^9, { 3.74884563117614*^9, 3.748845660433813*^9}, 3.7488457061834297`*^9, { 3.7488457365561666`*^9, 3.7488457623486423`*^9}, {3.749288202901496*^9, 3.749288226103823*^9}, 3.7492883061153483`*^9, 3.749288567272214*^9}] }, Open ]] }, Open ]] }, WindowSize->{1536, 781}, WindowMargins->{{17, Automatic}, {Automatic, 30}}, FrontEndVersion->"10.1 for Microsoft Windows (64-bit) (March 23, 2015)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[580, 22, 236, 5, 70, "Section"], Cell[CellGroupData[{ Cell[841, 31, 3624, 94, 352, "Input"], Cell[4468, 127, 2246, 48, 496, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[6763, 181, 352, 6, 111, "Section"], Cell[CellGroupData[{ Cell[7140, 191, 5506, 149, 442, "Input"], Cell[12649, 342, 3595, 69, 470, "Output"] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)