(* 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[ 32237, 647] NotebookOptionsPosition[ 15630, 344] NotebookOutlinePosition[ 32035, 638] CellTagsIndexPosition[ 31992, 635] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[TextData[{ "Damped pendulum swinging in a fluid\n", StyleBox["by Olaf Ellers, Ph.D., Research Associate, Biology and Math \ Departments, Bowdoin College", "Text", FontFamily->"Arial", FontSize->14], " \n", StyleBox["This model is of a pendulum swinging in a fluid in earth's gravity \ and is based on Neill et al., 2007. See below for things to try, assumptions \ and range of variables.\n", FontSize->14] }], "Title", Evaluatable->False, CellChangeTimes->{{3.449845031651*^9, 3.449845036105*^9}, { 3.4500062297759447`*^9, 3.4500062341149445`*^9}, 3.450015257863945*^9, { 3.450016247826945*^9, 3.450016329838945*^9}, {3.450183332638337*^9, 3.4501833638683367`*^9}, {3.450183609898337*^9, 3.4501836254133368`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"\[Rho]o", ">", "\[Rho]f"}], ",", RowBox[{"Plot", "[", RowBox[{ RowBox[{"l", " ", FractionBox[ RowBox[{"2", " ", "\[Pi]"}], "360"], "\[Alpha]", " ", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "t"}], " ", "\[Lambda]"}]], " ", RowBox[{"Cos", "[", RowBox[{ RowBox[{ FractionBox[ RowBox[{"2", " ", "\[Pi]"}], "360"], "\[Delta]"}], "-", " ", RowBox[{"t", " ", SqrtBox[ RowBox[{ RowBox[{"-", SuperscriptBox["\[Lambda]", "2"]}], "+", FractionBox[ RowBox[{"9.8", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "\[Rho]f"}], "+", "\[Rho]o"}], ")"}]}], RowBox[{"l", " ", RowBox[{"(", RowBox[{ RowBox[{"\[Kappa]", " ", "\[Rho]f"}], "+", "\[Rho]o"}], ")"}]}]]}]]}]}], "]"}]}], " ", ",", RowBox[{"{", RowBox[{"t", ",", "0", ",", "100"}], "}"}], ",", RowBox[{"PlotStyle", "\[Rule]", "Red"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}], ",", RowBox[{"AxesOrigin", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], ",", RowBox[{"Frame", "\[Rule]", "True"}], ",", RowBox[{"FrameLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\