(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 8.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 157, 7] NotebookDataLength[ 61501, 1849] NotebookOptionsPosition[ 58811, 1764] NotebookOutlinePosition[ 59166, 1780] CellTagsIndexPosition[ 59123, 1777] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[TextData[{ "Some Basics of Plotting in ", StyleBox["Mathematica", FontSlant->"Italic"], "\n", StyleBox["SIAM Computers Group \nAlanna Hoyer-Letizel", FontSize->16], "\n", StyleBox["April 11, 2012", FontSize->16] }], "Title", CellChangeTimes->{{3.543091816089097*^9, 3.5430918667457247`*^9}, { 3.543091972775078*^9, 3.5430919737516127`*^9}, {3.5430922145246067`*^9, 3.5430922187689342`*^9}}, TextAlignment->Center, Background->RGBColor[0.88, 1, 0.88]], Cell[TextData[{ "The biggest thing about using ", StyleBox["Mathematica", FontSlant->"Italic"], " is using the Documentation Center. I wouldn\[CloseCurlyQuote]t be able to \ do anything without it." }], "Text", CellChangeTimes->{{3.543092016246893*^9, 3.543092042698243*^9}, { 3.5480731816833344`*^9, 3.5480731827364883`*^9}}], Cell[CellGroupData[{ Cell["\<\ 2D Plots\ \>", "Section", CellChangeTimes->{{3.543090931476693*^9, 3.5430909467888527`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"Plotting", " ", "lines", " ", "in", " ", "2", RowBox[{"D", ".", " ", "The"}], " ", "Plot", " ", "Command", " ", "plots", " ", "graphs"}], "*)"}], "\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{ RowBox[{"Cos", "[", "x", "]"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "Pi"}], ",", "Pi"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.54298573440154*^9, 3.542985812869171*^9}, { 3.542985898853923*^9, 3.54298592665984*^9}, {3.5430894409274473`*^9, 3.543089450876*^9}, {3.543089488970437*^9, 3.543089492875338*^9}, { 3.5430901023233547`*^9, 3.543090140401409*^9}, {3.543092255131518*^9, 3.5430922627128363`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"Changing", " ", "the", " ", "PlotStyle"}], "*)"}], "\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{ RowBox[{"Cos", "[", "x", "]"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "Pi"}], ",", "Pi"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"Thick", ",", "Red"}], "}"}]}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"ContourPlot", " ", "in", " ", "2", "D"}], " ", "-", " ", RowBox[{ "plotting", " ", "curves", " ", "that", " ", "are", " ", "not", " ", "graphs"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"a", "=", RowBox[{"ContourPlot", "[", RowBox[{ RowBox[{ RowBox[{"x", "^", "2"}], "+", RowBox[{"y", "^", "2"}]}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", RowBox[{"{", RowBox[{"y", ",", RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", RowBox[{"ColorFunction", "\[Rule]", "\"\\""}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"ColorFunction", "\[Rule]", "\"\\""}], "*)"}]}]}]], "Input", CellChangeTimes->{{3.542985825923691*^9, 3.542985888848565*^9}, { 3.54298596880299*^9, 3.542985994382778*^9}, {3.542986026195225*^9, 3.542986054073496*^9}, {3.542986175707121*^9, 3.542986176645179*^9}, { 3.543065466706132*^9, 3.543065470603066*^9}, {3.543089463334189*^9, 3.543089524866724*^9}, {3.543091747727336*^9, 3.543091748206562*^9}, { 3.543092380277893*^9, 3.543092418857643*^9}, {3.543166203592572*^9, 3.543166204727873*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"One", " ", "Contour"}], "*)"}], RowBox[{"b", "=", RowBox[{"ContourPlot", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"x", "^", "2"}], "+", RowBox[{"y", "^", "2"}]}], "==", "1"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", RowBox[{"{", RowBox[{"y", ",", RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"ContourStyle", "\[Rule]", RowBox[{"{", RowBox[{"Red", ",", "Dashed", ",", "Thick"}], "}"}]}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.543092424714575*^9, 3.543092428514104*^9}, { 3.543166218535244*^9, 3.5431662225729856`*^9}, {3.543166286244215*^9, 3.5431662894412947`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"Putting", " ", "two", " ", "plots", " ", "together"}], "*)"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"a", ",", "b"}], "]"}]}]], "Input", CellChangeTimes->{{3.543089602453998*^9, 3.54308961813947*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", "ParametricPlot", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ "Plotting", " ", "lots", " ", "of", " ", "curves", " ", "in", " ", "2", " ", "different", " ", "ways"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"first", " ", "way"}], " ", "-", " ", RowBox[{"using", " ", "Show"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"Define", " ", "funtions", " ", "for", " ", "each", " ", RowBox[{"line", ".", " ", "t_"}], " ", "is", " ", "important"}], " ", "-", " ", RowBox[{ "tells", " ", "Mathematica", " ", "what", " ", "the", " ", "variable", " ", "of", " ", "the", " ", "function", " ", "is"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"f1", "[", "t_", "]"}], "=", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", "1"}], "}"}], "+", RowBox[{"t", "*", RowBox[{"(", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "2"}], "}"}], "-", RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", "1"}], "}"}]}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"a1", "=", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"f1", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "1"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"f2", "[", "t_", "]"}], "=", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "2"}], "}"}], "+", RowBox[{"t", "*", RowBox[{"(", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], "-", RowBox[{"{", RowBox[{"0", ",", "2"}], "}"}]}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"a2", "=", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"f2", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "1"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"f3", "[", "t_", "]"}], "=", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}], "+", RowBox[{"t", "*", RowBox[{"(", RowBox[{ RowBox[{"{", RowBox[{"1", ",", RowBox[{"-", "1"}]}], "}"}], "-", RowBox[{"{", RowBox[{"2", ",", "1"}], "}"}]}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"a3", "=", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"f3", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "1"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"f4", "[", "t_", "]"}], "=", RowBox[{ RowBox[{"{", RowBox[{"1", ",", RowBox[{"-", "1"}]}], "}"}], "+", RowBox[{"t", "*", RowBox[{"(", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", RowBox[{"-", "1"}]}], "}"}], "-", RowBox[{"{", RowBox[{"1", ",", RowBox[{"-", "1"}]}], "}"}]}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"a4", "=", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"f4", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "1"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"f5", "[", "t_", "]"}], "=", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", RowBox[{"-", "1"}]}], "}"}], "+", RowBox[{"t", "*", RowBox[{"(", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", "1"}], "}"}], "-", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", RowBox[{"-", "1"}]}], "}"}]}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"a5", "=", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"f5", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "1"}], "}"}]}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.543089674069655*^9, 3.543089868687202*^9}, 3.543089935348817*^9, {3.543090039829364*^9, 3.543090044399179*^9}, { 3.543090158383133*^9, 3.5430901652300663`*^9}, 3.543091753892076*^9, { 3.543092494572003*^9, 3.543092519057927*^9}}], Cell[BoxData[{ RowBox[{"Show", "[", RowBox[{"a1", ",", "a2", ",", "a3", ",", "a4", ",", "a5", ",", RowBox[{"PlotRange", "\[Rule]", "All"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}]}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{"Where", " ", "are", " ", "all", " ", "the", " ", "lines", RowBox[{"??", " ", RowBox[{ RowBox[{ RowBox[{"Add", " ", "PlotRange"}], "\[Rule]", "All"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}]}]}]}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5430898446168203`*^9, 3.543089853079336*^9}, { 3.5430898889907303`*^9, 3.543089923100087*^9}, {3.543090055885289*^9, 3.543090057237306*^9}, {3.5430925378903933`*^9, 3.5430925388901167`*^9}, { 3.5430925704289083`*^9, 3.543092587745303*^9}, {3.5431663451991568`*^9, 3.543166380077364*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"Second", " ", "way"}], " ", "-", " ", RowBox[{ "listing", " ", "them", " ", "all", " ", "in", " ", "ParametricPlot"}]}], "*)"}], RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"f1", "[", "t", "]"}], ",", RowBox[{"f2", "[", "t", "]"}], ",", RowBox[{"f3", "[", "t", "]"}], ",", RowBox[{"f4", "[", "t", "]"}], ",", RowBox[{"f5", "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "1"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.543089807552866*^9, 3.543089807834818*^9}, { 3.54309007790355*^9, 3.543090088791411*^9}, 3.5430917587832108`*^9, { 3.54316640500539*^9, 3.543166412850724*^9}}] }, Closed]], Cell[CellGroupData[{ Cell["\<\ 3D Plots\ \>", "Section", CellChangeTimes->{{3.543090963668992*^9, 3.543090969074733*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"Plotting", " ", "in", " ", "3", "D"}], " ", "-", " ", RowBox[{"Plot3D", " ", "is", " ", "good", " ", "for", " ", "graphs"}]}], ",", " ", RowBox[{ "ContourPlot3D", " ", "is", " ", "good", " ", "for", " ", "functions", " ", "that", " ", "are", " ", "not", " ", "graphs"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"surf1", "=", RowBox[{"Plot3D", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "^", "2"}], "-", RowBox[{"y", "^", "2"}]}], ")"}], "/", "10"}], "+", "3"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "3"}], ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"y", ",", RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"cylinder", "=", RowBox[{"ContourPlot3D", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"x", "^", "2"}], "+", RowBox[{"y", "^", "2"}]}], "==", "9"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "3"}], ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"y", ",", RowBox[{"-", "3"}], ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"z", ",", "0", ",", "9"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"surf1", ",", "cylinder"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.5429860584479237`*^9, 3.54298615109634*^9}, { 3.5429861866281633`*^9, 3.542986288767783*^9}, {3.5430901873999777`*^9, 3.5430902512130413`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "Use", " ", "ParametricPlot3D", " ", "to", " ", "get", " ", "a", " ", "curve", " ", "in", " ", "the", " ", "3", "D", " ", "by", " ", "only", " ", "having", " ", "one", " ", "parameter"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"intersection", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"3", RowBox[{"Cos", "[", "t", "]"}]}], ",", RowBox[{"3", RowBox[{"Sin", "[", "t", "]"}]}], ",", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"9", "/", "10"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Cos", "[", "t", "]"}], "^", "2"}], "-", RowBox[{ RowBox[{"Sin", "[", "t", "]"}], "^", "2"}]}], ")"}]}], "+", "3"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", RowBox[{"2", "Pi"}]}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"Thick", ",", "Green"}], "}"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"surf1", ",", "cylinder", ",", "intersection"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.542986296911333*^9, 3.542986382893647*^9}, { 3.543090253891777*^9, 3.5430902739712057`*^9}}] }, Closed]], Cell[CellGroupData[{ Cell["\<\ Animate command\ \>", "Section", CellChangeTimes->{{3.543090980421053*^9, 3.543090989793755*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"Animating", " ", "a", " ", "curve"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"Want", " ", "to", " ", "animate", " ", "a", " ", RowBox[{"point", "/", "particle"}], " ", "moving", " ", "along", " ", "this", " ", "half", " ", "circle"}], "*)"}], "\[IndentingNewLine]", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", "t", "]"}], ",", RowBox[{"Sin", "[", "t", "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "Pi"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.542986044024562*^9, 3.5429860457903633`*^9}, { 3.543090321894269*^9, 3.543090357885377*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"Plot", " ", "a", " ", "point", " ", "using", " ", RowBox[{"Graphics", "/", "Graphics3D"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", "Large", "]"}], ",", RowBox[{"Point", "[", RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", "0", "]"}], ",", RowBox[{"Sin", "[", "0", "]"}], ",", "0"}], "}"}], "]"}]}], "}"}], "]"}]}]], "Input", CellChangeTimes->{{3.543090404680251*^9, 3.5430904459674788`*^9}}], Cell[BoxData[ RowBox[{"Show", "[", RowBox[{ RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", "t", "]"}], ",", RowBox[{"Sin", "[", "t", "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "Pi"}], "}"}]}], "]"}], ",", RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", "Large", "]"}], ",", RowBox[{"Point", "[", RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", "0", "]"}], ",", RowBox[{"Sin", "[", "0", "]"}], ",", "0"}], "}"}], "]"}]}], "}"}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.5430927954342127`*^9, 3.543092832636244*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"Animate", " ", "the", " ", "Point", " ", "moving"}], "*)"}], "\[IndentingNewLine]", RowBox[{"Animate", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", "Large", "]"}], ",", RowBox[{"Point", "[", RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", "a", "]"}], ",", RowBox[{"Sin", "[", "a", "]"}], ",", "0"}], "}"}], "]"}]}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "}"}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"a", ",", "0", ",", "Pi"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.54309287141733*^9, 3.543092905998776*^9}, { 3.54316666103086*^9, 3.543166694157305*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"Animate", " ", "command"}], " ", "-", " ", RowBox[{"Have", " ", "the", " ", "line", " ", "plotted"}]}], ",", " ", RowBox[{ "and", " ", "have", " ", "the", " ", "point", " ", "moving", " ", "along", " ", "the", " ", "curve"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"Syntax", " ", "-", " ", RowBox[{ "need", " ", "Animate", " ", "outside", " ", "of", " ", "Show"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"Animate", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{ RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", "t", "]"}], ",", RowBox[{"Sin", "[", "t", "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "Pi"}], "}"}]}], "]"}], ",", RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", "Large", "]"}], ",", RowBox[{"Point", "[", RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", "a", "]"}], ",", RowBox[{"Sin", "[", "a", "]"}], ",", "0"}], "}"}], "]"}]}], "}"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"a", ",", "0", ",", "Pi"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.542980531943995*^9, 3.5429805334284573`*^9}, { 3.54298076576464*^9, 3.542980813272526*^9}, {3.542981079790057*^9, 3.5429811349844007`*^9}, {3.542981214220645*^9, 3.5429813374432917`*^9}, { 3.542981409269133*^9, 3.542981439677223*^9}, {3.5430902856261044`*^9, 3.5430903184406424`*^9}, {3.5430903643653393`*^9, 3.543090425981266*^9}, { 3.543090458200186*^9, 3.543090467233033*^9}, 3.543091765445591*^9}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"The", " ", "point", " ", "draws", " ", "the", " ", "line"}], " ", "-", " ", RowBox[{ "change", " ", "the", " ", "range", " ", "of", " ", "the", " ", "parametric", " ", "plot", " ", "to", " ", "have", " ", "the", " ", "animate", " ", "parameter"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"PlotRange", " ", "issue"}], ",", " ", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "}"}]}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"Animate", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{ RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", "t", "]"}], ",", RowBox[{"Sin", "[", "t", "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "a"}], "}"}]}], "]"}], ",", RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", "Large", "]"}], ",", RowBox[{"Point", "[", RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", "a", "]"}], ",", RowBox[{"Sin", "[", "a", "]"}], ",", "0"}], "}"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", "1"}], "}"}]}], "}"}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"a", ",", "0", ",", "Pi"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.543090472663525*^9, 3.543090613834465*^9}, { 3.54316674899695*^9, 3.543166749372306*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "Another", " ", "curve", " ", "with", " ", "the", " ", "same", " ", "animation", " ", "idea"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"p", "[", "t_", "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"-", "t"}], ",", RowBox[{ RowBox[{"Exp", "[", RowBox[{ RowBox[{"-", "1"}], "/", RowBox[{"(", RowBox[{"1", "-", RowBox[{"t", "^", "2"}]}], ")"}]}], "]"}], RowBox[{"Cos", "[", RowBox[{"20", " ", "t"}], "]"}]}], ",", RowBox[{ RowBox[{"Exp", "[", RowBox[{ RowBox[{"-", "1"}], "/", RowBox[{"(", RowBox[{"1", "-", RowBox[{"t", "^", "2"}]}], ")"}]}], "]"}], RowBox[{"Sin", "[", RowBox[{"20", " ", "t"}], "]"}]}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Animate", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{ RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"p", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{"-", ".99"}], ",", ".99"}], "}"}]}], "]"}], ",", RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", "Large", "]"}], ",", RowBox[{"Point", "[", RowBox[{"p", "[", "a", "]"}], "]"}]}], "}"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"a", ",", RowBox[{"-", ".99"}], ",", ".99"}], "}"}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.542982579908327*^9, 3.542982660912483*^9}, { 3.543090619057645*^9, 3.5430906500802794`*^9}}] }, Closed]], Cell[CellGroupData[{ Cell["\<\ Manipulate Command\ \>", "Section", CellChangeTimes->{{3.543091015282033*^9, 3.543091019991971*^9}}], Cell[CellGroupData[{ Cell["Kepler Problem", "Subsubsection", CellChangeTimes->{{3.543091032680092*^9, 3.543091034343277*^9}}], Cell[TextData[{ "Idea: As an astronaut circles the earth, she throws the lens cap of her \ camera at the earth. How does the lens cap move?\n\nThe Kepler Problem is \ given by the second order ODE ", Cell[BoxData[ FormBox[ RowBox[{" ", RowBox[{ RowBox[{ SuperscriptBox["d", "2"], RowBox[{"q", "/", SuperscriptBox[ StyleBox[ RowBox[{"d", "t"}]], "2"]}]}], "=", RowBox[{ RowBox[{ RowBox[{"-", "q"}], "/"}], "|", "q", "|", "3"}]}]}], TraditionalForm]], FormatType->"TraditionalForm"], ". One solution is a circle of radius 1, travelled at constant angular \ speed of 1, given by ", StyleBox["f(t)", FontSlant->"Italic"], " below. ", StyleBox["g(t)", FontSlant->"Italic"], " is the solution of the linear variational equation for this orbit, i.e. it \ gives the linear approximation of the motion of the lens cap thrown with \ inital velocity ", StyleBox["b. ", FontSlant->"Italic"], "Since this is the linearization, the results would only be accurate for a \ small about of time." }], "Text", CellChangeTimes->{{3.5430907243854027`*^9, 3.543090769269052*^9}, { 3.543152816510016*^9, 3.543152829055748*^9}, {3.5431528659960623`*^9, 3.543153093336309*^9}, {3.5431531729501762`*^9, 3.543153295917983*^9}, { 3.5431533711561327`*^9, 3.543153418064148*^9}, {3.543153460008553*^9, 3.543153477821621*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{ "solution", " ", "to", " ", "linear", " ", "variational", " ", "equation", " ", "of", " ", "a", " ", "circle", " ", "orbit", " ", "for", " ", "the", " ", "kepler", " ", "problem"}], ",", " ", RowBox[{ RowBox[{"i", ".", "e", ".", " ", "an"}], " ", "astronaut", " ", "throwing", " ", "a", " ", "lens", " ", "cap"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"b", "=", ".5"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"f", "[", "t_", "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", "t", "]"}], ",", RowBox[{"Sin", "[", "t", "]"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"g", "[", "t_", "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"b", " ", RowBox[{"Sin", "[", "t", "]"}]}]}], ")"}], RowBox[{"Cos", "[", RowBox[{"t", "-", RowBox[{"2", "b", " ", RowBox[{"Cos", "[", "t", "]"}]}], "+", RowBox[{"2", "b"}]}], "]"}]}], ",", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"b", " ", RowBox[{"Sin", "[", "t", "]"}]}]}], ")"}], RowBox[{"Sin", "[", RowBox[{"t", "-", RowBox[{"2", "b", " ", RowBox[{"Cos", "[", "t", "]"}]}], "+", RowBox[{"2", "b"}]}], "]"}]}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"f", "[", "t", "]"}], ",", RowBox[{"g", "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", RowBox[{"2", "Pi"}]}], "}"}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.5429872895228033`*^9, 3.542987301165572*^9}, { 3.54298765511399*^9, 3.542987656305252*^9}, {3.543153109750434*^9, 3.543153127570046*^9}, {3.543166986136785*^9, 3.543166988707027*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{ "Check", " ", "different", " ", "values", " ", "for", " ", "b", " ", "without", " ", "having", " ", "to", " ", "change", " ", "and", " ", "evaluate", " ", "the", " ", "cell", " ", "each", " ", RowBox[{"time", ".", " ", "Using"}], " ", "a"}], "=", "b"}], "*)"}], RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"a", " ", RowBox[{"Sin", "[", "t", "]"}]}]}], ")"}], RowBox[{"Cos", "[", RowBox[{"t", "-", RowBox[{"2", "a", " ", RowBox[{"Cos", "[", "t", "]"}]}], "+", RowBox[{"2", "a"}]}], "]"}]}], ",", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{"a", " ", RowBox[{"Sin", "[", "t", "]"}]}]}], ")"}], RowBox[{"Sin", "[", RowBox[{"t", "-", RowBox[{"2", "a", " ", RowBox[{"Cos", "[", "t", "]"}]}], "+", RowBox[{"2", "a"}]}], "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", "t", "]"}], ",", RowBox[{"Sin", "[", "t", "]"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", RowBox[{"2", "Pi"}]}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}], ",", RowBox[{"{", RowBox[{"a", ",", "0", ",", "2"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.543153132326734*^9, 3.54315316887643*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "pick", " ", "a", " ", "value", " ", "for", " ", "b", " ", "and", " ", "Animate", " ", "the", " ", "motion"}], "*)"}], RowBox[{ RowBox[{"Animate", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{ RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"f", "[", "t", "]"}], ",", RowBox[{"g", "[", "t", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", RowBox[{"2", "Pi"}]}], "}"}]}], "]"}], ",", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", "Large", "]"}], ",", RowBox[{"Point", "[", RowBox[{"f", "[", "a", "]"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", "Large", "]"}], ",", RowBox[{"Point", "[", RowBox[{"g", "[", "a", "]"}], "]"}]}], "}"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"a", ",", "0", ",", RowBox[{"2", "Pi"}]}], "}"}]}], "]"}], "\[IndentingNewLine]"}]}]], "Input", CellChangeTimes->{{3.542987319748034*^9, 3.542987359604966*^9}, { 3.54298739819851*^9, 3.542987534418174*^9}, {3.543153335291047*^9, 3.543153350185092*^9}}] }, Closed]], Cell[CellGroupData[{ Cell["Tangent/Normal vectors", "Subsubsection", CellChangeTimes->{{3.543091056320116*^9, 3.5430910676614733`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "moving", " ", "around", " ", "tangent", " ", "vectors", " ", "on", " ", "a", " ", "surface"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"f", "[", RowBox[{"x_", ",", "y_"}], "]"}], " ", "=", RowBox[{ RowBox[{"Cos", "[", "x", "]"}], RowBox[{"Sin", "[", "y", "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"surf", " ", "=", " ", RowBox[{"Plot3D", "[", RowBox[{ RowBox[{"f", "[", RowBox[{"x", ",", "y"}], "]"}], ",", " ", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "Pi"}], ",", "Pi"}], "}"}], ",", RowBox[{"{", RowBox[{"y", ",", RowBox[{"-", "Pi"}], ",", "Pi"}], "}"}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\"", ",", "\"\\""}], "}"}]}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"Opacity", "[", "0.7", "]"}]}]}], "]"}]}]}]}]], "Input", CellChangeTimes->{ 3.5429815318878727`*^9, {3.54298159335992*^9, 3.542981595861565*^9}, { 3.5429816406354923`*^9, 3.542981742758642*^9}, 3.542981774382258*^9, { 3.542981831456635*^9, 3.542981835869978*^9}, {3.542987632345478*^9, 3.542987640712852*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"fx", "[", RowBox[{"x_", ",", "y_"}], "]"}], "=", RowBox[{"D", "[", RowBox[{ RowBox[{"f", "[", RowBox[{"x", ",", "y"}], "]"}], ",", "x"}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"fy", "[", RowBox[{"x_", ",", "y_"}], "]"}], "=", RowBox[{"D", "[", RowBox[{ RowBox[{"f", "[", RowBox[{"x", ",", "y"}], "]"}], ",", "y"}], "]"}]}]}], "Input", CellChangeTimes->{{3.542981845167671*^9, 3.542982009749681*^9}, { 3.542982109113448*^9, 3.542982161907308*^9}, {3.542982210576249*^9, 3.542982214457073*^9}, {3.542982251372095*^9, 3.542982285620368*^9}, { 3.543167789622699*^9, 3.543167792716675*^9}}], Cell[BoxData[ RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{"surf", ",", RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{"Red", ",", RowBox[{"Thickness", "[", ".02", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", RowBox[{"f", "[", RowBox[{"a", ",", "b"}], "]"}]}], "}"}], ",", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", RowBox[{"f", "[", RowBox[{"a", ",", "b"}], "]"}]}], "}"}], "+", RowBox[{"{", RowBox[{"1", ",", "0", ",", RowBox[{"fx", "[", RowBox[{"a", ",", "b"}], "]"}]}], "}"}]}]}], "}"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{"Green", ",", RowBox[{"Thickness", "[", ".02", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", RowBox[{"f", "[", RowBox[{"a", ",", "b"}], "]"}]}], "}"}], ",", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "b", ",", RowBox[{"f", "[", RowBox[{"a", ",", "b"}], "]"}]}], "}"}], "+", RowBox[{"{", RowBox[{"0", ",", "1", ",", RowBox[{"fy", "[", RowBox[{"a", ",", "b"}], "]"}]}], "}"}]}]}], "}"}], "]"}]}], "}"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"a", ",", RowBox[{"-", "Pi"}], ",", "Pi"}], "}"}], ",", RowBox[{"{", RowBox[{"b", ",", RowBox[{"-", "Pi"}], ",", "Pi"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.5429822173726273`*^9, 3.542982217839198*^9}, { 3.542982254969059*^9, 3.5429822892397213`*^9}, {3.5429823727156763`*^9, 3.542982418912768*^9}}], Cell[BoxData[ RowBox[{"(*", RowBox[{ RowBox[{ "plot", " ", "a", " ", "normal", " ", "vector", " ", "on", " ", "a", " ", "surface"}], " ", "-", " ", RowBox[{"cylinder", " ", "and", " ", "ellipsoid"}]}], "*)"}]], "Input", CellChangeTimes->{{3.542988261115035*^9, 3.542988310419178*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"Parameterize", " ", "a", " ", "cylinder"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"cyl", "[", RowBox[{"t_", ",", "z_"}], "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", "t", "]"}], ",", RowBox[{"Sin", "[", "t", "]"}], ",", "z"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"cylPlot", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"cyl", "[", RowBox[{"t", ",", "z"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", RowBox[{"2", "Pi"}]}], "}"}], ",", RowBox[{"{", RowBox[{"z", ",", "0", ",", "3"}], "}"}], ",", RowBox[{"PlotStyle", "->", RowBox[{"Opacity", "[", "0.3", "]"}]}]}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.543071790477228*^9, 3.5430717985451107`*^9}, { 3.54307298970611*^9, 3.543073077029077*^9}, {3.543073835097546*^9, 3.5430738388288507`*^9}, {3.543073882547213*^9, 3.543073885649953*^9}, { 3.543074034523075*^9, 3.54307404569663*^9}, 3.543090851456023*^9}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"Normal", " ", "vector", " ", "expression"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Dcdt", "[", RowBox[{"t_", ",", "z_"}], "]"}], "=", RowBox[{"D", "[", RowBox[{ RowBox[{"cyl", "[", RowBox[{"t", ",", "z"}], "]"}], ",", "t"}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Dcdz", "[", RowBox[{"t_", ",", "z_"}], "]"}], "=", RowBox[{"D", "[", RowBox[{ RowBox[{"cyl", "[", RowBox[{"t", ",", "z"}], "]"}], ",", "z"}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"normal1", "[", RowBox[{"t_", ",", "z_"}], "]"}], "=", RowBox[{ RowBox[{"Cross", "[", RowBox[{ RowBox[{"Dcdt", "[", RowBox[{"t", ",", "z"}], "]"}], ",", RowBox[{"Dcdz", "[", RowBox[{"t", ",", "z"}], "]"}]}], "]"}], "/", RowBox[{"Norm", "[", RowBox[{"Cross", "[", RowBox[{ RowBox[{"Dcdt", "[", RowBox[{"t", ",", "z"}], "]"}], ",", RowBox[{"Dcdz", "[", RowBox[{"t", ",", "z"}], "]"}]}], "]"}], "]"}]}]}]}]}]], "Input", CellChangeTimes->{{3.5430730864685717`*^9, 3.543073128608835*^9}, { 3.5430731725332947`*^9, 3.5430731999092903`*^9}, {3.5430736816058397`*^9, 3.5430736943968887`*^9}, 3.5430740990729322`*^9, 3.543074130375354*^9}], Cell[BoxData[ RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{"cylPlot", ",", RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{"Thick", ",", "Red", ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"cyl", "[", RowBox[{"a", ",", "b"}], "]"}], ",", RowBox[{ RowBox[{"cyl", "[", RowBox[{"a", ",", "b"}], "]"}], "+", RowBox[{".5", RowBox[{"normal1", "[", RowBox[{"a", ",", "b"}], "]"}]}]}]}], "}"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "3"}], "}"}]}], "}"}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"a", ",", "0", ",", RowBox[{"2", "Pi"}]}], "}"}], ",", RowBox[{"{", RowBox[{"b", ",", "0", ",", "3"}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.543073213749926*^9, 3.543073261353788*^9}, { 3.543073355510126*^9, 3.543073522573155*^9}, {3.543073628817011*^9, 3.543073631751145*^9}, {3.543073764985094*^9, 3.543073765487995*^9}, { 3.543073800646401*^9, 3.543073818420796*^9}, {3.543073956911873*^9, 3.543074015891156*^9}, {3.543074052659634*^9, 3.5430740760553617`*^9}, { 3.543074118151785*^9, 3.5430741188625937`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"Parameterize", " ", "a", " ", "ellipsoid"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"ell", "[", RowBox[{"s_", ",", "t_"}], "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"3", RowBox[{"Cos", "[", "t", "]"}], RowBox[{"Sin", "[", "s", "]"}]}], ",", RowBox[{"2", RowBox[{"Sin", "[", "t", "]"}], RowBox[{"Sin", "[", "s", "]"}]}], ",", RowBox[{"Cos", "[", "s", "]"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"ellPlot", "=", RowBox[{"ParametricPlot3D", "[", RowBox[{ RowBox[{"ell", "[", RowBox[{"s", ",", "t"}], "]"}], ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", RowBox[{"2", "Pi"}]}], "}"}], ",", RowBox[{"{", RowBox[{"s", ",", "0", ",", "Pi"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"Opacity", "[", "0.3", "]"}]}]}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.543090883702993*^9, 3.5430908977238483`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"Normal", " ", "Vector", " ", "expression"}], " ", "-", " ", RowBox[{"same", " ", "idea", " ", "as", " ", "before"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Deds", "[", RowBox[{"s_", ",", "t_"}], "]"}], "=", RowBox[{"D", "[", RowBox[{ RowBox[{"ell", "[", RowBox[{"s", ",", "t"}], "]"}], ",", "s"}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Dedt", "[", RowBox[{"s_", ",", "t_"}], "]"}], "=", RowBox[{"D", "[", RowBox[{ RowBox[{"ell", "[", RowBox[{"s", ",", "t"}], "]"}], ",", "t"}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"normal2", "[", RowBox[{"s_", ",", "t_"}], "]"}], "=", RowBox[{"Cross", "[", RowBox[{ RowBox[{"Deds", "[", RowBox[{"s", ",", "t"}], "]"}], ",", RowBox[{"Dedt", "[", RowBox[{"s", ",", "t"}], "]"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{"Simplify", "[", "%", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"normal22", "[", RowBox[{"s_", ",", "t_"}], "]"}], "=", RowBox[{ RowBox[{"normal2", "[", RowBox[{"s", ",", "t"}], "]"}], "/", RowBox[{"Norm", "[", RowBox[{"normal2", "[", RowBox[{"s", ",", "t"}], "]"}], "]"}]}]}], ";"}]}]}]], "Input", CellChangeTimes->{{3.543090904634996*^9, 3.5430909083611298`*^9}}], Cell[BoxData[ RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{"ellPlot", ",", RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{"Thick", ",", "Red", ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"ell", "[", RowBox[{"a", ",", "b"}], "]"}], ",", RowBox[{ RowBox[{"ell", "[", RowBox[{"a", ",", "b"}], "]"}], "+", RowBox[{"normal22", "[", RowBox[{"a", ",", "b"}], "]"}]}]}], "}"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}]}], "]"}], ",", RowBox[{"{", RowBox[{"a", ",", ".001", ",", RowBox[{"Pi", "-", ".001"}]}], "}"}], ",", RowBox[{"{", RowBox[{"b", ",", "0", ",", RowBox[{"2", "Pi"}]}], "}"}]}], "]"}]], "Input"] }, Closed]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Cycloid Animation\ \>", "Section", CellChangeTimes->{{3.54309109192589*^9, 3.54309109595606*^9}}], Cell[CellGroupData[{ Cell["Cycloid", "Subsubsection", CellChangeTimes->{{3.543091124332377*^9, 3.543091138481805*^9}, { 3.54309118227136*^9, 3.543091233084634*^9}, {3.54315396764279*^9, 3.5431539705540943`*^9}, 3.543154326040881*^9, {3.543155061473106*^9, 3.543155065307933*^9}}], Cell[TextData[StyleBox["A cycloid is a way of writing out the time \ parameterizations of radius term in Kepler\[CloseCurlyQuote]s problem. \ Newtown actually found this without any formulas. ", FontWeight->"Plain"]], "Text", CellChangeTimes->{3.543155071480131*^9}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "Can", " ", "we", " ", "plot", " ", "the", " ", "circle", " ", "that", " ", "generates", " ", "the", " ", "cycloid", RowBox[{"??", "?"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"cycloid", "[", RowBox[{"a_", ",", "b_"}], "]"}], "[", "t_", "]"}], ":=", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"a", " ", "t"}], "-", RowBox[{"b", " ", RowBox[{"Sin", "[", "t", "]"}]}]}], ",", RowBox[{"a", "-", RowBox[{"b", " ", RowBox[{"Cos", "[", "t", "]"}]}]}]}], "}"}]}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"at", ",", "a"}], "}"}], " ", "is", " ", "the", " ", "center", " ", "of", " ", "the", " ", "circle"}], ",", " ", RowBox[{"b", " ", "is", " ", "length", " ", "of", " ", RowBox[{"line", "."}]}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"the", " ", ":=", " ", RowBox[{ RowBox[{ "is", " ", "a", " ", "way", " ", "to", " ", "define", " ", "variables", " ", "in", " ", "a", " ", RowBox[{"more", " ", "'"}], "flexibile", " ", RowBox[{"way", "'"}], " ", "to", " ", "define", " ", "a", " ", "function"}], " ", "-", " ", RowBox[{ "look", " ", "at", " ", "documentation", " ", "center", " ", "here"}]}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"\[Alpha]", ":=", RowBox[{"cycloid", "[", RowBox[{"1", ",", "1"}], "]"}]}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"Plot", " ", "of", " ", "the", " ", "cycloid"}], "*)"}], "\[IndentingNewLine]", RowBox[{"aaa", "=", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"\[Alpha]", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{ RowBox[{"-", "3"}], "Pi"}], ",", RowBox[{"3", "Pi"}]}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "3"}], "Pi"}], ",", RowBox[{"3", "Pi"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", "4"}], "}"}]}], "}"}]}]}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.542982936922441*^9, 3.542983078745019*^9}, { 3.542984360342574*^9, 3.542984361125194*^9}, {3.542986590229257*^9, 3.542986615049432*^9}, {3.542986683476039*^9, 3.542986707218588*^9}, { 3.542987673912116*^9, 3.542987674960492*^9}, {3.54309126163953*^9, 3.543091276444491*^9}, {3.543091328579352*^9, 3.543091370256719*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "My", " ", "Goal", " ", "is", " ", "to", " ", "get", " ", "a", " ", "picture", " ", "like", " ", "the", " ", "one", " ", "on", " ", "the", " ", "wikipedia", " ", "page", " ", "for", " ", "cycloid"}], "*)"}], RowBox[{"(*", RowBox[{ "Plotting", " ", "a", " ", "circle", " ", "with", " ", "Graphics"}], "*)"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"aaa", ",", RowBox[{"Graphics", "[", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}], ",", "1"}], "]"}], "]"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.542984283335067*^9, 3.5429843653547783`*^9}, { 3.543091376472145*^9, 3.543091455045549*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"Animating", " ", "the", " ", "moving", " ", "circle"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ "Cheating", " ", "a", " ", "little", " ", "because", " ", "the", " ", "circle", " ", RowBox[{"isn", "'"}], "t", " ", "spinning"}], ",", " ", RowBox[{ RowBox[{"it", "'"}], "s", " ", "just", " ", "sliding", " ", "across"}]}], "*)"}], RowBox[{"Animate", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{"aaa", ",", RowBox[{"Graphics", "[", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "1"}], "}"}], ",", "1"}], "]"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "Pi"}], ",", RowBox[{"2", "Pi"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", "4"}], "}"}]}], "}"}]}], ",", RowBox[{"Axes", "\[Rule]", "True"}]}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"a", ",", "0", ",", RowBox[{"2", "Pi"}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.542984511086886*^9, 3.542984594399383*^9}, { 3.542984628480743*^9, 3.5429846415302057`*^9}, {3.542984681433155*^9, 3.542984682655438*^9}, {3.542984726494389*^9, 3.5429847684558983`*^9}, { 3.542984890663534*^9, 3.5429849049883413`*^9}, {3.543091390881584*^9, 3.5430914005124474`*^9}, {3.543091458493557*^9, 3.543091458877694*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "Adding", " ", "a", " ", "point", " ", "on", " ", "the", " ", "circle"}], "*)"}], RowBox[{"Animate", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{"aaa", ",", RowBox[{"Graphics", "[", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "1"}], "}"}], ",", "1"}], "]"}], "]"}], ",", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", "Large", "]"}], ",", RowBox[{"Point", "[", RowBox[{"{", RowBox[{ RowBox[{"a", "-", RowBox[{"Sin", "[", "a", "]"}]}], ",", RowBox[{"1", "-", RowBox[{"Cos", "[", "a", "]"}]}]}], "}"}], "]"}]}], "}"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"a", ",", "0", ",", RowBox[{"2", "Pi"}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.542984776776166*^9, 3.542984862857094*^9}, { 3.542984946577058*^9, 3.542985162346569*^9}, {3.542985217844492*^9, 3.5429852180655193`*^9}, {3.542985249530497*^9, 3.542985267725499*^9}, { 3.543091463494931*^9, 3.543091472198501*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "Adding", " ", "a", " ", "line", " ", "on", " ", "the", " ", "radius"}], "*)"}], RowBox[{"Animate", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{"aaa", ",", RowBox[{"Graphics", "[", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "1"}], "}"}], ",", "1"}], "]"}], "]"}], ",", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", "Medium", "]"}], ",", RowBox[{"Point", "[", RowBox[{"{", RowBox[{ RowBox[{"a", "-", RowBox[{"Sin", "[", "a", "]"}]}], ",", RowBox[{"1", "-", RowBox[{"Cos", "[", "a", "]"}]}]}], "}"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"Graphics", "[", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "1"}], "}"}], ",", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "1"}], "}"}], "+", RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"Sin", "[", "a", "]"}]}], ",", RowBox[{"-", RowBox[{"Cos", "[", "a", "]"}]}]}], "}"}]}]}], "}"}], "]"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"a", ",", RowBox[{ RowBox[{"-", "3"}], "Pi"}], ",", RowBox[{"3", "Pi"}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.542985306055427*^9, 3.5429853800243807`*^9}, { 3.5429854131293707`*^9, 3.542985447264612*^9}, {3.5430914742459497`*^9, 3.543091483060041*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"Having", " ", "the", " ", RowBox[{"point", " ", "'"}], RowBox[{"draw", "'"}], " ", "the", " ", "cycloid"}], "*)"}], RowBox[{"Animate", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{ RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"\[Alpha]", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{ RowBox[{"-", "3"}], "Pi"}], ",", "a"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "3"}], "Pi"}], ",", RowBox[{"3", "Pi"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", "4"}], "}"}]}], "}"}]}]}], "]"}], ",", RowBox[{"Graphics", "[", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "1"}], "}"}], ",", "1"}], "]"}], "]"}], ",", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", "Medium", "]"}], ",", RowBox[{"Point", "[", RowBox[{"{", RowBox[{ RowBox[{"a", "-", RowBox[{"Sin", "[", "a", "]"}]}], ",", RowBox[{"1", "-", RowBox[{"Cos", "[", "a", "]"}]}]}], "}"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"Graphics", "[", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "1"}], "}"}], ",", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "1"}], "}"}], "+", RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"Sin", "[", "a", "]"}]}], ",", RowBox[{"-", RowBox[{"Cos", "[", "a", "]"}]}]}], "}"}]}]}], "}"}], "]"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"a", ",", RowBox[{ RowBox[{"-", "3"}], "Pi"}], ",", RowBox[{"3", "Pi"}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.542985513246793*^9, 3.542985613612606*^9}, { 3.543091486045701*^9, 3.543091499817669*^9}, 3.54315509952092*^9}] }, Closed]], Cell[CellGroupData[{ Cell["Prolate Cycloid", "Subsubsection", CellChangeTimes->{{3.5430911077967987`*^9, 3.543091111517544*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"A", " ", "Prolate", " ", "Cycloid"}], " ", "-", " ", RowBox[{ "point", " ", "is", " ", "outside", " ", "of", " ", "the", " ", "circle"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Beta]", ":=", RowBox[{"cycloid", "[", RowBox[{"1", ",", "3"}], "]"}]}], "\[IndentingNewLine]", RowBox[{"\[Beta]", "[", "t", "]"}], "\[IndentingNewLine]", RowBox[{"bbb", "=", RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"\[Beta]", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{ RowBox[{"-", "3"}], "Pi"}], ",", RowBox[{"3", "Pi"}]}], "}"}]}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.5429833272203617`*^9, 3.542983375799705*^9}, { 3.542984349028757*^9, 3.542984350028599*^9}, {3.542986757801737*^9, 3.5429867625393047`*^9}, {3.542986811430359*^9, 3.542986818598353*^9}, { 3.5429877337617083`*^9, 3.542987739133597*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "Prolate", " ", "cycloid", " ", "is", " ", "when", " ", "the", " ", "point", " ", "is", " ", "a", " ", "fixed", " ", "distance", " ", "outside", " ", "the", " ", "circle"}], "*)"}], RowBox[{"Show", "[", RowBox[{"bbb", ",", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}], ",", "1"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", "Medium", "]"}], ",", "Red", ",", RowBox[{"Point", "[", RowBox[{"{", RowBox[{"0", ",", RowBox[{"1", "-", "3"}]}], "}"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"Thick", ",", "Red", ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}], ",", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}], "+", RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", "3"}]}], "}"}]}]}], "}"}], "]"}]}], "}"}]}], "}"}], "]"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.5429877484329767`*^9, 3.542987896601365*^9}, { 3.543091509842142*^9, 3.543091610212529*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "Circle", " ", "and", " ", "line", " ", "drawing", " ", "out", " ", "the", " ", "cycloid"}], "*)"}], RowBox[{"Animate", "[", RowBox[{ RowBox[{"Show", "[", RowBox[{ RowBox[{"ParametricPlot", "[", RowBox[{ RowBox[{"\[Beta]", "[", "t", "]"}], ",", RowBox[{"{", RowBox[{"t", ",", RowBox[{ RowBox[{"-", "3"}], "Pi"}], ",", "a"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "3"}], "Pi"}], ",", RowBox[{"3", "Pi"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "4"}], "}"}]}], "}"}]}]}], "]"}], ",", RowBox[{"Graphics", "[", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "1"}], "}"}], ",", "1"}], "]"}], "]"}], ",", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"PointSize", "[", "Medium", "]"}], ",", "Red", ",", RowBox[{"Point", "[", RowBox[{"{", RowBox[{ RowBox[{"a", "-", RowBox[{"3", RowBox[{"Sin", "[", "a", "]"}]}]}], ",", RowBox[{"1", "-", RowBox[{"3", RowBox[{"Cos", "[", "a", "]"}]}]}]}], "}"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Thick", ",", "Red", ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "1"}], "}"}], ",", RowBox[{ RowBox[{"{", RowBox[{"a", ",", "1"}], "}"}], "+", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "3"}], RowBox[{"Sin", "[", "a", "]"}]}], ",", RowBox[{ RowBox[{"-", "3"}], RowBox[{"Cos", "[", "a", "]"}]}]}], "}"}]}]}], "}"}], "]"}]}], "}"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"a", ",", RowBox[{ RowBox[{"-", "3"}], "Pi"}], ",", RowBox[{"3", "Pi"}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.54298672531723*^9, 3.542986729972666*^9}, { 3.542986776742844*^9, 3.542986779974863*^9}, {3.543091617873549*^9, 3.543091641889699*^9}, {3.5430916972717876`*^9, 3.543091707665392*^9}}] }, Closed]] }, Open ]] }, Open ]] }, WindowSize->{800, 503}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, FrontEndVersion->"8.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (February 23, \ 2011)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[579, 22, 475, 15, 173, "Title"], Cell[1057, 39, 334, 8, 26, "Text"], Cell[CellGroupData[{ Cell[1416, 51, 104, 4, 89, "Section"], Cell[1523, 57, 729, 16, 43, "Input"], Cell[2255, 75, 426, 13, 43, "Input"], Cell[2684, 90, 1291, 32, 58, "Input"], Cell[3978, 124, 782, 23, 43, "Input"], Cell[4763, 149, 278, 7, 43, "Input"], Cell[5044, 158, 4458, 140, 238, "Input"], Cell[9505, 300, 949, 21, 43, "Input"], Cell[10457, 323, 784, 21, 43, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[11278, 349, 102, 4, 59, "Section"], Cell[11383, 355, 1743, 51, 88, "Input"], Cell[13129, 408, 1448, 41, 73, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[14614, 454, 109, 4, 59, "Section"], Cell[14726, 460, 744, 18, 58, "Input"], Cell[15473, 480, 574, 16, 43, "Input"], Cell[16050, 498, 729, 21, 43, "Input"], Cell[16782, 521, 1094, 32, 58, "Input"], Cell[17879, 555, 1806, 44, 73, "Input"], Cell[19688, 601, 2116, 61, 103, "Input"], Cell[21807, 664, 1741, 53, 73, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[23585, 722, 112, 4, 59, "Section"], Cell[CellGroupData[{ Cell[23722, 730, 105, 1, 24, "Subsubsection"], Cell[23830, 733, 1399, 38, 107, "Text"], Cell[25232, 773, 2084, 60, 103, "Input"], Cell[27319, 835, 1732, 47, 88, "Input"], Cell[29054, 884, 1380, 39, 88, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[30471, 928, 115, 1, 18, "Subsubsection"], Cell[30589, 931, 1305, 35, 73, "Input"], Cell[31897, 968, 690, 19, 43, "Input"], Cell[32590, 989, 1992, 55, 88, "Input"], Cell[34585, 1046, 304, 7, 27, "Input"], Cell[34892, 1055, 1132, 30, 58, "Input"], Cell[36027, 1087, 1390, 41, 73, "Input"], Cell[37420, 1130, 1553, 41, 43, "Input"], Cell[38976, 1173, 1086, 32, 58, "Input"], Cell[40065, 1207, 1463, 44, 103, "Input"], Cell[41531, 1253, 1174, 36, 43, "Input"] }, Closed]] }, Open ]], Cell[CellGroupData[{ Cell[42754, 1295, 109, 4, 89, "Section"], Cell[CellGroupData[{ Cell[42888, 1303, 270, 4, 24, "Subsubsection"], Cell[43161, 1309, 268, 4, 41, "Text"], Cell[43432, 1315, 2718, 73, 133, "Input"], Cell[46153, 1390, 762, 20, 58, "Input"], Cell[46918, 1412, 1602, 42, 73, "Input"], Cell[48523, 1456, 1200, 33, 58, "Input"], Cell[49726, 1491, 1672, 50, 73, "Input"], Cell[51401, 1543, 2275, 67, 73, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[53713, 1615, 108, 1, 18, "Subsubsection"], Cell[53824, 1618, 1003, 25, 73, "Input"], Cell[54830, 1645, 1368, 38, 73, "Input"], Cell[56201, 1685, 2570, 74, 88, "Input"] }, Closed]] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)