(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 6.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 14604, 402] NotebookOptionsPosition[ 13489, 361] NotebookOutlinePosition[ 13830, 376] CellTagsIndexPosition[ 13787, 373] WindowFrame->Normal ContainsDynamic->True *) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"ClearAll", "[", RowBox[{"\"\\"", ",", " ", "K"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"Clear", "[", "K", "]"}], ";"}], "*)"}]}], "\n", RowBox[{ RowBox[{"Remove", "[", "\"\\"", "]"}], ";"}]}], "Input", CellChangeTimes->{{3.442232427734375*^9, 3.442232430375*^9}, { 3.442232475140625*^9, 3.442232496640625*^9}, 3.44224014571875*^9, { 3.4422402084375*^9, 3.44224024675*^9}, {3.4422471935625*^9, 3.442247218953125*^9}, {3.442247264015625*^9, 3.442247288015625*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"a", "=", "5"}], ";"}]], "Input", CellChangeTimes->{{3.442240158640625*^9, 3.442240161125*^9}}], Cell[CellGroupData[{ Cell[BoxData["a"], "Input", CellChangeTimes->{3.442240167046875*^9}], Cell[BoxData["5"], "Output", CellChangeTimes->{ 3.442240167640625*^9, 3.4422403534375*^9, 3.4422471179375*^9, 3.4422471909375*^9, {3.442247267921875*^9, 3.442247291296875*^9}, 3.442457982515625*^9, 3.44245806396875*^9}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"u1", "[", "t_", "]"}], ":=", RowBox[{".1", "t", " ", RowBox[{"Cos", "[", "t", "]"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"u2", "[", "t_", "]"}], ":=", RowBox[{".1", "t", " ", RowBox[{"Sin", "[", "t", "]"}]}]}]}], "Input", CellChangeTimes->{{3.441830967109375*^9, 3.441831042625*^9}, { 3.4418310731875*^9, 3.441831087265625*^9}, {3.4418311854375*^9, 3.441831238390625*^9}, {3.441831287546875*^9, 3.441831287609375*^9}, { 3.441831454984375*^9, 3.44183146028125*^9}, {3.441832003859375*^9, 3.4418320218125*^9}, 3.4418321413125*^9, {3.44183733746875*^9, 3.441837369703125*^9}, 3.441837852046875*^9, 3.441837902265625*^9, { 3.441841096828125*^9, 3.441841105140625*^9}, {3.441841243046875*^9, 3.44184127271875*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"u1", "[", "2", "]"}]], "Input", CellChangeTimes->{{3.4422323836875*^9, 3.44223238828125*^9}}], Cell[BoxData[ RowBox[{"-", "0.08322936730942848`"}]], "Output", CellChangeTimes->{{3.442232384265625*^9, 3.442232388796875*^9}, 3.44224035346875*^9, 3.442247117953125*^9, 3.442247190953125*^9, { 3.4422472679375*^9, 3.4422472913125*^9}, 3.442457982546875*^9, 3.442458063984375*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"\[Alpha]", "[", "t_", "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{"x", "[", "t", "]"}], ",", " ", RowBox[{"y", "[", "t", "]"}]}], "}"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Alpha]", "[", "t_", "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{"u1", "[", "t", "]"}], ",", " ", RowBox[{"u2", "[", "t", "]"}]}], "}"}]}]}]], "Input", CellChangeTimes->{{3.441834332125*^9, 3.4418343909375*^9}, { 3.44183731315625*^9, 3.441837332046875*^9}, {3.44183819834375*^9, 3.441838205625*^9}, {3.441838250875*^9, 3.441838258609375*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"T", "[", "t_", "]"}], ":=", " ", RowBox[{"N", "[", RowBox[{ RowBox[{ RowBox[{"\[Alpha]", "'"}], "[", "t", "]"}], "/", RowBox[{"Norm", "[", RowBox[{ RowBox[{"\[Alpha]", "'"}], "[", "t", "]"}], "]"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.44183439409375*^9, 3.441834440875*^9}, { 3.441837686296875*^9, 3.441837695171875*^9}, {3.44184023040625*^9, 3.44184023259375*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"K", "[", "t_", "]"}], ":=", RowBox[{"Abs", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[Alpha]", "'"}], "[", "t", "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], RowBox[{ RowBox[{ RowBox[{"\[Alpha]", "''"}], "[", "t", "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}], "-", RowBox[{ RowBox[{ RowBox[{ RowBox[{"\[Alpha]", "''"}], "[", "t", "]"}], "[", RowBox[{"[", "1", "]"}], "]"}], RowBox[{ RowBox[{ RowBox[{"\[Alpha]", "'"}], "[", "t", "]"}], "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ")"}], "/", SuperscriptBox[ RowBox[{"Norm", "[", RowBox[{ RowBox[{"\[Alpha]", "'"}], "[", "t", "]"}], "]"}], "3"]}], "]"}]}]], "Input", CellChangeTimes->{{3.441837079734375*^9, 3.441837119515625*^9}, { 3.44183716771875*^9, 3.44183717225*^9}, {3.44183723078125*^9, 3.4418372835625*^9}, 3.441837701703125*^9, {3.441837892484375*^9, 3.441837893609375*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"Nvec", "[", "t_", "]"}], ":=", RowBox[{ RowBox[{"N", "[", RowBox[{"FunctionExpand", "[", RowBox[{"Normalize", "[", RowBox[{ RowBox[{ RowBox[{"T", "'"}], "[", "t", "]"}], "/", RowBox[{"Norm", "[", RowBox[{ RowBox[{"\[Alpha]", "'"}], "[", "t", "]"}], "]"}]}], " ", "]"}], "]"}], "]"}], "/.", RowBox[{"{", RowBox[{"t0", "\[Rule]", "t"}], "}"}]}]}]], "Input", CellChangeTimes->{{3.4418375784375*^9, 3.4418376678125*^9}, { 3.441837738484375*^9, 3.441837751765625*^9}, {3.441838279546875*^9, 3.44183842453125*^9}, {3.441838564828125*^9, 3.4418386005*^9}, { 3.441839121828125*^9, 3.441839149453125*^9}, {3.4418395068125*^9, 3.441839510734375*^9}, {3.441840121671875*^9, 3.4418401509375*^9}, { 3.441840204671875*^9, 3.441840209078125*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Norm", "[", RowBox[{"Nvec", "[", "1", "]"}], "]"}]], "Input", CellChangeTimes->{{3.441839593515625*^9, 3.4418395994375*^9}, 3.441839758015625*^9, {3.441839822078125*^9, 3.441839839984375*^9}, { 3.441839950046875*^9, 3.44183999503125*^9}, {3.441840054453125*^9, 3.441840105015625*^9}, {3.4418402233125*^9, 3.44184024084375*^9}, { 3.441840361703125*^9, 3.441840381671875*^9}}], Cell[BoxData["1.`"], "Output", CellChangeTimes->{{3.441839825875*^9, 3.441839840828125*^9}, { 3.441839958859375*^9, 3.44183999584375*^9}, {3.4418400604375*^9, 3.441840105875*^9}, {3.4418402260625*^9, 3.441840241640625*^9}, { 3.44184036809375*^9, 3.441840382109375*^9}, 3.441929420765625*^9, 3.442232308578125*^9, 3.442240353515625*^9, 3.442247118171875*^9, 3.442247191015625*^9, {3.442247268*^9, 3.442247291359375*^9}, 3.442457983140625*^9, 3.442458064046875*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"OscCenter", " ", "[", "t_", "]"}], ":=", RowBox[{ RowBox[{"\[Alpha]", "[", "t", "]"}], "+", RowBox[{ RowBox[{"Nvec", "[", "t", "]"}], "/", RowBox[{"K", "[", "t", "]"}]}]}]}]], "Input", CellChangeTimes->{{3.441840392671875*^9, 3.441840484296875*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{"{", RowBox[{ RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"\[Alpha]", "[", "t", "]"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", " ", "0", ",", " ", RowBox[{"4", "\[Pi]"}]}], "}"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"Directive", "[", RowBox[{"Red", ",", " ", "Thick"}], "]"}], "}"}]}]}], "]"}], ",", "\[IndentingNewLine]", "\t\t\t\t", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Blue", ",", " ", "Thick", ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"OscCenter", "[", "t0", "]"}], ",", " ", RowBox[{"1", "/", RowBox[{"K", "[", "t0", "]"}]}]}], "]"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", "\t\t\t\t", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", "Large", "]"}], ",", "Black", ",", RowBox[{"Point", "[", RowBox[{"\[Alpha]", "[", "t0", "]"}], "]"}]}], "}"}], "]"}]}], "}"}], "\[IndentingNewLine]", "\t\t\t ", "]"}], ",", " ", RowBox[{"{", RowBox[{"t0", ",", " ", "0", ",", " ", RowBox[{"4", "\[Pi]"}]}], "}"}]}], "\[IndentingNewLine]", " ", "]"}]], "Input", CellChangeTimes->{{3.44184049265625*^9, 3.441840505234375*^9}, { 3.44184054478125*^9, 3.441840682921875*^9}, {3.441840713828125*^9, 3.441840738625*^9}, {3.44184081521875*^9, 3.441840919109375*^9}, { 3.441840951671875*^9, 3.44184099109375*^9}, {3.441841132640625*^9, 3.441841184859375*^9}, {3.441841295484375*^9, 3.441841308171875*^9}, { 3.442458056671875*^9, 3.4424581778125*^9}}], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`t0$$ = 5.730265000147783, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{ Hold[$CellContext`t0$$], 0, 4 Pi}}, Typeset`size$$ = { 360., {153., 157.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`t0$4907$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`t0$$ = 0}, "ControllerVariables" :> { Hold[$CellContext`t0$$, $CellContext`t0$4907$$, 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" :> Show[{ ParametricPlot[ $CellContext`\[Alpha][$CellContext`t], {$CellContext`t, 0, 4 Pi}, PlotStyle -> { Directive[Red, Thick]}], Graphics[{Blue, Thick, Circle[ $CellContext`OscCenter[$CellContext`t0$$], 1/ K[$CellContext`t0$$]]}], Graphics[{ PointSize[Large], Black, Point[ $CellContext`\[Alpha][$CellContext`t0$$]]}]}], "Specifications" :> {{$CellContext`t0$$, 0, 4 Pi}}, "Options" :> {}, "DefaultOptions" :> {}], ImageSizeCache->{406., {203., 208.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Manipulate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellChangeTimes->{ 3.441840506296875*^9, 3.4418405503125*^9, 3.441840601375*^9, { 3.441840670484375*^9, 3.441840683828125*^9}, 3.4418407396875*^9, { 3.441840898375*^9, 3.44184091996875*^9}, 3.441840952515625*^9, 3.441840993140625*^9, 3.441841111953125*^9, 3.441841145671875*^9, 3.4418411858125*^9, {3.44184127590625*^9, 3.441841309125*^9}, 3.441929421546875*^9, {3.44223231021875*^9, 3.4422323341875*^9}, 3.442240353546875*^9, 3.44224711821875*^9, 3.4422471910625*^9, { 3.442247268046875*^9, 3.44224729140625*^9}, 3.442457984140625*^9, { 3.442458064078125*^9, 3.442458119109375*^9}, {3.44245817359375*^9, 3.44245817859375*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Export", "[", RowBox[{"\"\\"", ",", RowBox[{"Show", "[", RowBox[{"{", RowBox[{ RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"\[Alpha]", "[", "t", "]"}], ",", " ", RowBox[{"{", RowBox[{"t", ",", " ", "0", ",", " ", RowBox[{"4", "\[Pi]"}]}], "}"}], ",", " ", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"Directive", "[", RowBox[{"Red", ",", " ", "Thick"}], "]"}], "}"}]}], ",", RowBox[{"Axes", "\[Rule]", "False"}]}], "]"}], ",", "\[IndentingNewLine]", "\t\t\t\t", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Blue", ",", " ", "Thick", ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"OscCenter", "[", "5.7", "]"}], ",", " ", RowBox[{"1", "/", RowBox[{"K", "[", "5.7", "]"}]}]}], "]"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", "\t\t\t\t", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", "Large", "]"}], ",", "Black", ",", RowBox[{"Point", "[", RowBox[{"\[Alpha]", "[", "5.7", "]"}], "]"}]}], "}"}], "]"}]}], "}"}], "\[IndentingNewLine]", "\t\t\t ", "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.442458208109375*^9, 3.442458257703125*^9}, { 3.442458319140625*^9, 3.442458422421875*^9}, {3.442458496109375*^9, 3.442458567609375*^9}}], Cell[BoxData["\<\"Osculating.jpg\"\>"], "Output", CellChangeTimes->{ 3.44245823246875*^9, {3.442458302375*^9, 3.44245842509375*^9}, 3.442458497609375*^9, 3.442458535703125*^9, 3.44245857078125*^9}], Cell[BoxData["\<\"C:\\\\Documents and Settings\\\\Jon\\\\My Documents\"\>"], \ "Output", CellChangeTimes->{3.442458609265625*^9}] }, Open ]] }, WindowSize->{763, 730}, WindowMargins->{{64, Automatic}, {Automatic, 34}}, FrontEndVersion->"6.0 for Microsoft Windows (32-bit) (June 19, 2007)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[568, 21, 597, 14, 72, "Input"], Cell[1168, 37, 130, 3, 31, "Input"], Cell[CellGroupData[{ Cell[1323, 44, 69, 1, 31, "Input"], Cell[1395, 47, 230, 4, 30, "Output"] }, Open ]], Cell[1640, 54, 797, 16, 52, "Input"], Cell[CellGroupData[{ Cell[2462, 74, 118, 2, 31, "Input"], Cell[2583, 78, 293, 5, 30, "Output"] }, Open ]], Cell[2891, 86, 654, 18, 52, "Input"], Cell[3548, 106, 445, 12, 31, "Input"], Cell[3996, 120, 1098, 33, 33, "Input"], Cell[5097, 155, 859, 21, 31, "Input"], Cell[CellGroupData[{ Cell[5981, 180, 416, 7, 31, "Input"], Cell[6400, 189, 486, 7, 30, "Output"] }, Open ]], Cell[6901, 199, 304, 8, 31, "Input"], Cell[CellGroupData[{ Cell[7230, 211, 1809, 42, 112, "Input"], Cell[9042, 255, 2559, 53, 428, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11638, 313, 1495, 36, 132, "Input"], Cell[13136, 351, 204, 3, 30, "Output"], Cell[13343, 356, 130, 2, 30, "Output"] }, Open ]] } ] *) (* End of internal cache information *)