(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.2' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[      5852,        154]
NotebookOptionsPosition[      5552,        141]
NotebookOutlinePosition[      5898,        156]
CellTagsIndexPosition[      5855,        153]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{"DynamicModule", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"\[Alpha]pr", ",", " ", "Tmax", ",", " ", "\[Mu]pr"}], "}"}], ",",
    "\[IndentingNewLine]", 
   RowBox[{"Manipulate", "[", "\[IndentingNewLine]", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{"Module", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"u", ",", " ", "t", ",", " ", "\[Phi]"}], "}"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"sol", " ", "=", " ", 
          RowBox[{"NDSolve", "[", "\[IndentingNewLine]", 
           RowBox[{
            RowBox[{"{", "\[IndentingNewLine]", 
             RowBox[{
              RowBox[{
               RowBox[{
                RowBox[{
                 RowBox[{"u", "''"}], "[", "t", "]"}], " ", "-", " ", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{"1", " ", "-", " ", "\[Mu]pr"}], ")"}], " ", 
                 RowBox[{"u", "[", "t", "]"}], " ", 
                 SuperscriptBox[
                  RowBox[{
                   RowBox[{"\[Phi]", "'"}], "[", "t", "]"}], "2"]}], "+", " ",
                 "\[Mu]pr"}], " ", "\[Equal]", " ", "0"}], ",", " ", 
              RowBox[{
               RowBox[{
                RowBox[{
                 RowBox[{"u", "[", "t", "]"}], " ", 
                 RowBox[{
                  RowBox[{"\[Phi]", "''"}], "[", "t", "]"}]}], " ", "+", " ", 
                RowBox[{"2", " ", 
                 RowBox[{
                  RowBox[{"u", "'"}], "[", "t", "]"}], " ", 
                 RowBox[{
                  RowBox[{"\[Phi]", "'"}], "[", "t", "]"}]}]}], " ", 
               "\[Equal]", " ", "0"}], ",", " ", "\[IndentingNewLine]", 
              RowBox[{
               RowBox[{"u", "[", "0", "]"}], " ", "\[Equal]", " ", "1"}], ",",
               " ", "\[IndentingNewLine]", 
              RowBox[{
               RowBox[{
                RowBox[{"u", "'"}], "[", "0", "]"}], " ", "\[Equal]", " ", 
               "0"}], ",", "\[IndentingNewLine]", 
              RowBox[{
               RowBox[{"\[Phi]", "[", "0", "]"}], " ", "\[Equal]", "0"}], ",",
               " ", "\[IndentingNewLine]", 
              RowBox[{
               RowBox[{
                RowBox[{"\[Phi]", "'"}], "[", "0", "]"}], " ", "\[Equal]", 
               " ", "\[Alpha]pr"}]}], "\[IndentingNewLine]", "}"}], ",", " ", 
            "\[IndentingNewLine]", 
            RowBox[{"{", 
             RowBox[{"u", ",", " ", "\[Phi]"}], "}"}], ",", " ", 
            "\[IndentingNewLine]", 
            RowBox[{"{", 
             RowBox[{"t", ",", " ", "0", ",", " ", 
              RowBox[{"tmax", " ", "=", " ", "100"}]}], "}"}]}], 
           "\[IndentingNewLine]", "]"}]}], ";"}]}], "\[IndentingNewLine]", 
       "]"}], ";", "\n", 
      RowBox[{"usol", " ", "=", " ", 
       RowBox[{
        RowBox[{
         RowBox[{"sol", "[", 
          RowBox[{"[", "1", "]"}], "]"}], "[", 
         RowBox[{"[", "1", "]"}], "]"}], "[", 
        RowBox[{"[", "2", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"\[Phi]sol", " ", "=", " ", 
       RowBox[{
        RowBox[{
         RowBox[{"sol", "[", 
          RowBox[{"[", "1", "]"}], "]"}], "[", 
         RowBox[{"[", "2", "]"}], "]"}], "[", 
        RowBox[{"[", "2", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
      RowBox[{"ParametricPlot", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"usol", "[", "t", "]"}], " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Cos", "[", 
            RowBox[{"\[Phi]sol", "[", "t", "]"}], "]"}], ",", " ", 
           RowBox[{"Sin", "[", 
            RowBox[{"\[Phi]sol", "[", "t", "]"}], "]"}]}], "}"}]}], ",", " ", 
        RowBox[{"{", 
         RowBox[{"t", ",", " ", "0", ",", " ", "Tmax"}], "}"}]}], "]"}]}], 
     ",", " ", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"\[Alpha]pr", ",", " ", "0.429"}], "}"}], ",", " ", "0.001", 
       ",", " ", "10", ",", " ", "0.001", ",", " ", 
       RowBox[{"Appearance", "\[Rule]", "\"\<Open\>\""}]}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"Tmax", ",", " ", "21.511"}], "}"}], ",", " ", "0.001", ",", 
       " ", "tmax", ",", " ", "0.01", ",", " ", 
       RowBox[{"Appearance", "\[Rule]", "\"\<Open\>\""}]}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"\[Mu]pr", ",", " ", "0.5"}], "}"}], ",", " ", "0.001", ",", 
       " ", "1", ",", " ", "0.001", ",", " ", 
       RowBox[{"Appearance", "\[Rule]", "\"\<Open\>\""}]}], "}"}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"TrackedSymbols", "\[RuleDelayed]", " ", 
      RowBox[{"{", 
       RowBox[{"\[Alpha]pr", ",", " ", "Tmax", ",", " ", "\[Mu]pr"}], 
       "}"}]}]}], "\[IndentingNewLine]", "]"}]}], "\[IndentingNewLine]", 
  "]"}]], "Input",ExpressionUUID->"47d1a10d-09e6-4b41-b659-221db4de48ec"]
},
WindowSize->{1904, 997},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
FrontEndVersion->"11.2 for Microsoft Windows (64-bit) (September 10, 2017)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 4990, 119, 528, "Input",ExpressionUUID->"47d1a10d-09e6-4b41-b659-221db4de48ec"]
}
]
*)