(* Content-type: application/mathematica *)

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

(* CreatedBy='Mathematica 7.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     46915,       1201]
NotebookOptionsPosition[     44565,       1131]
NotebookOutlinePosition[     45009,       1149]
CellTagsIndexPosition[     44966,       1146]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Duffing", "Section",
 CellChangeTimes->{{3.468856993903*^9, 3.468857001199*^9}, 
   3.468862832278*^9, {3.4692952581190643`*^9, 3.469295264015751*^9}}],

Cell[CellGroupData[{

Cell["Part a)", "Subsection",
 CellChangeTimes->{{3.4692953123436213`*^9, 3.469295314387182*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Clear", "[", "x", "]"}], ";", 
  RowBox[{"Clear", "[", "x1", "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.469295506972704*^9, 3.4692955302318563`*^9}, {
  3.4692957167886686`*^9, 3.4692957229817495`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"x0", "=", "0"}], ";", 
  RowBox[{"xD0", "=", "0"}], ";", 
  RowBox[{"\[Gamma]", "=", "0.2"}], ";", 
  RowBox[{"f", "=", "10"}], ";", 
  RowBox[{
   SubscriptBox["\[Omega]", "D"], "=", "1"}], ";"}]], "Input",
 CellChangeTimes->{
  3.467841596705*^9, 3.467841751833*^9, {3.4678606728055*^9, 
   3.4678606753824997`*^9}, {3.4678608832904997`*^9, 3.467860883666*^9}, {
   3.467861140986*^9, 3.467861141658*^9}, {3.467906653521*^9, 
   3.467906659006*^9}, {3.4679067469315*^9, 3.4679067673445*^9}, {
   3.4679068282165003`*^9, 3.4679068368570004`*^9}, 3.4679069646324997`*^9, {
   3.4679171801115*^9, 3.46791718571*^9}, {3.4679173367485*^9, 
   3.4679173369055*^9}, 3.4679173798945*^9, {3.4679174573085003`*^9, 
   3.4679174580825*^9}, 3.467917517924*^9, {3.4679180034065*^9, 
   3.4679180046055*^9}, {3.467919220707*^9, 3.4679192313535*^9}, {
   3.467919405164*^9, 3.4679194080345*^9}, {3.4679195475295*^9, 
   3.4679195477135*^9}, {3.467919710981*^9, 3.467919711572*^9}, 
   3.4679203548265*^9, {3.4679209260030003`*^9, 3.4679209262019997`*^9}, {
   3.4679211605179996`*^9, 3.467921182793*^9}, {3.4679217631305*^9, 
   3.4679217643915*^9}, {3.4679245217539997`*^9, 3.4679245505655003`*^9}, {
   3.46802846966*^9, 3.468028511933*^9}, 3.4680285519765*^9, 
   3.4680288580175*^9, {3.468029178826*^9, 3.468029193855*^9}, {
   3.4680292246275*^9, 3.4680292496805*^9}, 3.468029294163*^9, {
   3.468029382799*^9, 3.468029434977*^9}, 3.4680294803780003`*^9, {
   3.468029536988*^9, 3.4680295542825003`*^9}, {3.468029664835*^9, 
   3.4680296658265*^9}, {3.4680301131305*^9, 3.4680301490439997`*^9}, 
   3.4680302945945*^9, {3.4680304101875*^9, 3.4680304104674997`*^9}, {
   3.468030935449*^9, 3.4680309503905*^9}, {3.46803099105*^9, 
   3.4680310040235*^9}, {3.4688527359119997`*^9, 3.468852781068*^9}, {
   3.468852908909*^9, 3.46885292338*^9}, {3.468852969242*^9, 
   3.4688529869849997`*^9}, 3.468853068932*^9, 3.468855115584*^9, {
   3.468857014001*^9, 3.468857018488*^9}, {3.468857072519*^9, 
   3.468857083652*^9}, 3.4688628322790003`*^9, {3.4692952743427525`*^9, 
   3.469295345976575*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"sol1", "=", 
   RowBox[{"NDSolve", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"x", "''"}], "[", "t", "]"}], "+", 
         RowBox[{"\[Gamma]", " ", 
          RowBox[{
           RowBox[{"x", "'"}], "[", "t", "]"}]}], "+", 
         SuperscriptBox[
          RowBox[{"x", "[", "t", "]"}], "3"]}], "\[Equal]", 
        RowBox[{"f", " ", 
         RowBox[{"Cos", "[", 
          RowBox[{
           SubscriptBox["\[Omega]", "D"], "t"}], "]"}]}]}], ",", 
       RowBox[{
        RowBox[{"x", "[", "0", "]"}], "\[Equal]", "x0"}], ",", 
       RowBox[{
        RowBox[{
         RowBox[{"x", "'"}], "[", "0", "]"}], "\[Equal]", "xD0"}]}], "}"}], 
     ",", "x", ",", 
     RowBox[{"{", 
      RowBox[{"t", ",", "0", ",", 
       RowBox[{"100", " ", "\[Pi]"}]}], "}"}], ",", 
     RowBox[{"MaxSteps", "\[Rule]", "500000"}]}], "]"}]}], ";", 
  RowBox[{
   RowBox[{"x1", "[", "t_", "]"}], ":=", 
   RowBox[{
    RowBox[{"x", "[", "t", "]"}], "/.", 
    RowBox[{"sol1", "[", 
     RowBox[{"[", "1", "]"}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.4680285262365*^9, 3.4680285394585*^9}, {
   3.4680285703585*^9, 3.468028573014*^9}, 3.4680288969895*^9, {
   3.4680290470745*^9, 3.4680290478665*^9}, {3.4680293288145*^9, 
   3.46802933095*^9}, {3.4680296779625*^9, 3.468029679322*^9}, {
   3.4680297339630003`*^9, 3.468029793699*^9}, {3.468029832959*^9, 
   3.468030002321*^9}, {3.4680309622795*^9, 3.4680309640550003`*^9}, {
   3.468031146472*^9, 3.468031148959*^9}, {3.4688527931210003`*^9, 
   3.4688528048310003`*^9}, {3.468852872124*^9, 3.468852894954*^9}, {
   3.46885293219*^9, 3.46885295179*^9}, {3.468853007074*^9, 
   3.468853039137*^9}, {3.468855137976*^9, 3.4688551490950003`*^9}, 
   3.4688628322799997`*^9, {3.469295340407482*^9, 3.469295448341806*^9}, {
   3.469323022951*^9, 3.469323024593*^9}}],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"x1", "[", "t", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"t", ",", "0", ",", " ", 
     RowBox[{"30", "\[Pi]"}]}], "}"}], ",", 
   RowBox[{"Ticks", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{"5", "\[Pi]"}], ",", 
       RowBox[{"10", "\[Pi]"}], ",", " ", 
       RowBox[{"15", "\[Pi]"}], ",", 
       RowBox[{"20", "\[Pi]"}], ",", 
       RowBox[{"25", "\[Pi]"}], ",", " ", 
       RowBox[{"30", "\[Pi]"}]}], "}"}], "}"}]}], ",", 
   RowBox[{"BaseStyle", "->", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
      RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<t\>\"", ",", "\"\<x\>\""}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.4692957800922513`*^9, 3.4692957973611193`*^9}, {
  3.46929620942362*^9, 3.4692962559887247`*^9}, {3.469297547020994*^9, 
  3.4692975512797127`*^9}, {3.469323035863*^9, 3.4693230613564997`*^9}}],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"x1", "[", "t", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"t", ",", 
     RowBox[{"30", "\[Pi]"}], ",", " ", 
     RowBox[{"60", "\[Pi]"}]}], "}"}], ",", 
   RowBox[{"Ticks", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"30", "\[Pi]"}], ",", 
       RowBox[{"35", "\[Pi]"}], ",", 
       RowBox[{"40", "\[Pi]"}], ",", " ", 
       RowBox[{"45", "\[Pi]"}], ",", 
       RowBox[{"50", "\[Pi]"}], ",", 
       RowBox[{"55", "\[Pi]"}], ",", " ", 
       RowBox[{"60", "\[Pi]"}]}], "}"}], "}"}]}], ",", 
   RowBox[{"BaseStyle", "->", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
      RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<t\>\"", ",", "\"\<x\>\""}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.4692957800922513`*^9, 3.4692957973611193`*^9}, {
  3.46929620942362*^9, 3.4692962559887247`*^9}, {3.469297547020994*^9, 
  3.4692975512797127`*^9}, {3.469323035863*^9, 3.4693231111254997`*^9}}],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"x1", "[", "t", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"t", ",", 
     RowBox[{"60", "\[Pi]"}], ",", " ", 
     RowBox[{"90", "\[Pi]"}]}], "}"}], ",", 
   RowBox[{"Ticks", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"60", "\[Pi]"}], ",", 
       RowBox[{"65", "\[Pi]"}], ",", 
       RowBox[{"70", "\[Pi]"}], ",", " ", 
       RowBox[{"75", "\[Pi]"}], ",", 
       RowBox[{"80", "\[Pi]"}], ",", 
       RowBox[{"85", "\[Pi]"}], ",", " ", 
       RowBox[{"90", "\[Pi]"}]}], "}"}], "}"}]}], ",", 
   RowBox[{"BaseStyle", "->", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
      RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<t\>\"", ",", "\"\<x\>\""}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.4692957800922513`*^9, 3.4692957973611193`*^9}, {
  3.46929620942362*^9, 3.4692962559887247`*^9}, {3.469297547020994*^9, 
  3.4692975512797127`*^9}, {3.469323035863*^9, 3.4693231111254997`*^9}, {
  3.4693231785575*^9, 3.4693232011540003`*^9}}],

Cell["\<\
Can see that things have settled down to periodic at the end (really when \
past about 20 \[Pi]) and the last period is \
\>", "Text",
 CellChangeTimes->{{3.4692976369532647`*^9, 3.469297683799164*^9}, {
  3.4693232193164997`*^9, 3.4693232296115*^9}, {3.500751099442754*^9, 
  3.500751103994014*^9}}],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"x1", "[", "t", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"t", ",", 
     RowBox[{"84", " ", "\[Pi]"}], ",", " ", 
     RowBox[{"90", "\[Pi]"}]}], "}"}], ",", 
   RowBox[{"Ticks", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"85", "\[Pi]"}], ",", " ", 
       RowBox[{"90", "\[Pi]"}]}], "}"}], "}"}]}], ",", 
   RowBox[{"BaseStyle", "->", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
      RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<t\>\"", ",", "\"\<x\>\""}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.4692957800922513`*^9, 3.4692957973611193`*^9}, {
  3.46929620942362*^9, 3.4692962559887247`*^9}, {3.469297547020994*^9, 
  3.4692975512797127`*^9}, {3.4692976938609705`*^9, 3.46929771677693*^9}, {
  3.5007511832205453`*^9, 3.5007512125962257`*^9}}],

Cell["Compare", "Text",
 CellChangeTimes->{{3.5007512322493496`*^9, 3.500751233793438*^9}}],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"x1", "[", "t", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"t", ",", 
     RowBox[{"24", " ", "\[Pi]"}], ",", " ", 
     RowBox[{"30", "\[Pi]"}]}], "}"}], ",", 
   RowBox[{"Ticks", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"25", "\[Pi]"}], ",", " ", 
       RowBox[{"30", "\[Pi]"}]}], "}"}], "}"}]}], ",", 
   RowBox[{"BaseStyle", "->", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
      RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<t\>\"", ",", "\"\<x\>\""}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.4692957800922513`*^9, 3.4692957973611193`*^9}, {
  3.46929620942362*^9, 3.4692962559887247`*^9}, {3.469297547020994*^9, 
  3.4692975512797127`*^9}, {3.4692976938609705`*^9, 3.46929771677693*^9}}],

Cell["Nearly identical!", "Text",
 CellChangeTimes->{{3.500751240177803*^9, 3.5007512462811527`*^9}}],

Cell["\<\
The period is 6 \[Pi] or 3 times the driving period of 2 \[Pi], i.e., this \
illustrates period tripling.  We can look at this directly in phase space\
\>", "Text",
 CellChangeTimes->{{3.469297733343811*^9, 3.469297788785145*^9}, {
  3.469297839140977*^9, 3.46929785562986*^9}}],

Cell[BoxData[
 RowBox[{"pp1", "=", 
  RowBox[{"ParametricPlot", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"x1", "[", "t", "]"}], ",", 
      RowBox[{
       RowBox[{"x1", "'"}], "[", "t", "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"t", ",", "0", ",", 
      RowBox[{"30", " ", "\[Pi]"}]}], "}"}], ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "7"}], ",", "7"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", "Blue", "}"}]}], ",", 
    RowBox[{"Ticks", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "\[Pi]"}], ",", 
        RowBox[{
         RowBox[{"-", "\[Pi]"}], "/", "2"}], ",", "0", ",", 
        RowBox[{"\[Pi]", "/", "2"}], ",", "\[Pi]"}], "}"}], "}"}]}], ",", 
    RowBox[{"BaseStyle", "->", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<x\>\"", ",", "\"\<\!\(\*OverscriptBox[\(x\), \(.\)]\)\>\""}], 
      "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4680286449795*^9, 3.4680286850055*^9}, {
   3.4680287633094997`*^9, 3.468028768587*^9}, 3.4680289031945*^9, {
   3.4680290752205*^9, 3.468029075908*^9}, {3.468029338283*^9, 
   3.468029338875*^9}, {3.468030174417*^9, 3.4680302386730003`*^9}, {
   3.4680306396565*^9, 3.4680306624455*^9}, {3.468031030751*^9, 
   3.4680310547320004`*^9}, 3.468031125557*^9, {3.4680312544744997`*^9, 
   3.4680312743035*^9}, {3.4688530827609997`*^9, 3.468853111349*^9}, {
   3.468853193159*^9, 3.468853215678*^9}, {3.468853248796*^9, 
   3.4688532831870003`*^9}, {3.468853387758*^9, 3.4688534372019997`*^9}, {
   3.468853468929*^9, 3.4688534793599997`*^9}, {3.468853562748*^9, 
   3.468853619392*^9}, {3.468853950534*^9, 3.468853962612*^9}, 
   3.468862832283*^9, {3.4688774688120003`*^9, 3.4688774809630003`*^9}, {
   3.469296433404113*^9, 3.469296454427534*^9}, {3.469296534485194*^9, 
   3.469296590051326*^9}, {3.469307402037056*^9, 3.469307417636756*^9}, {
   3.4693232466225*^9, 3.469323254117*^9}, 3.46932364966*^9}],

Cell["At late times", "Text",
 CellChangeTimes->{{3.5007512650502257`*^9, 3.500751267801383*^9}}],

Cell[BoxData[
 RowBox[{"pp1a", "=", 
  RowBox[{"ParametricPlot", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"x1", "[", "t", "]"}], ",", 
      RowBox[{
       RowBox[{"x1", "'"}], "[", "t", "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"t", ",", 
      RowBox[{"22", "\[Pi]"}], ",", 
      RowBox[{"30", " ", "\[Pi]"}]}], "}"}], ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "7"}], ",", "7"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", "Blue", "}"}]}], ",", 
    RowBox[{"Ticks", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "\[Pi]"}], ",", 
        RowBox[{
         RowBox[{"-", "\[Pi]"}], "/", "2"}], ",", "0", ",", 
        RowBox[{"\[Pi]", "/", "2"}], ",", "\[Pi]"}], "}"}], "}"}]}], ",", 
    RowBox[{"BaseStyle", "->", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<x\>\"", ",", "\"\<\!\(\*OverscriptBox[\(x\), \(.\)]\)\>\""}], 
      "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4680286449795*^9, 3.4680286850055*^9}, {
   3.4680287633094997`*^9, 3.468028768587*^9}, 3.4680289031945*^9, {
   3.4680290752205*^9, 3.468029075908*^9}, {3.468029338283*^9, 
   3.468029338875*^9}, {3.468030174417*^9, 3.4680302386730003`*^9}, {
   3.4680306396565*^9, 3.4680306624455*^9}, {3.468031030751*^9, 
   3.4680310547320004`*^9}, 3.468031125557*^9, {3.4680312544744997`*^9, 
   3.4680312743035*^9}, {3.4688530827609997`*^9, 3.468853111349*^9}, {
   3.468853193159*^9, 3.468853215678*^9}, {3.468853248796*^9, 
   3.4688532831870003`*^9}, {3.468853387758*^9, 3.4688534372019997`*^9}, {
   3.468853468929*^9, 3.4688534793599997`*^9}, {3.468853562748*^9, 
   3.468853619392*^9}, {3.468853950534*^9, 3.468853962612*^9}, 
   3.468862832283*^9, {3.4688774688120003`*^9, 3.4688774809630003`*^9}, {
   3.469296433404113*^9, 3.469296454427534*^9}, {3.469296534485194*^9, 
   3.469296590051326*^9}, {3.469296625150651*^9, 3.46929667201215*^9}, 
   3.4693075964405174`*^9, {3.4693232782345*^9, 3.4693232955205*^9}}],

Cell["Go to later times", "Text",
 CellChangeTimes->{{3.469323312729*^9, 3.469323315664*^9}}],

Cell[BoxData[
 RowBox[{"pp1b", "=", 
  RowBox[{"ParametricPlot", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"x1", "[", "t", "]"}], ",", 
      RowBox[{
       RowBox[{"x1", "'"}], "[", "t", "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"t", ",", 
      RowBox[{"80", "\[Pi]"}], ",", 
      RowBox[{"100", " ", "\[Pi]"}]}], "}"}], ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "7"}], ",", "7"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", "Blue", "}"}]}], ",", 
    RowBox[{"Ticks", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "\[Pi]"}], ",", 
        RowBox[{
         RowBox[{"-", "\[Pi]"}], "/", "2"}], ",", "0", ",", 
        RowBox[{"\[Pi]", "/", "2"}], ",", "\[Pi]"}], "}"}], "}"}]}], ",", 
    RowBox[{"BaseStyle", "->", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<x\>\"", ",", "\"\<\!\(\*OverscriptBox[\(x\), \(.\)]\)\>\""}], 
      "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4680286449795*^9, 3.4680286850055*^9}, {
   3.4680287633094997`*^9, 3.468028768587*^9}, 3.4680289031945*^9, {
   3.4680290752205*^9, 3.468029075908*^9}, {3.468029338283*^9, 
   3.468029338875*^9}, {3.468030174417*^9, 3.4680302386730003`*^9}, {
   3.4680306396565*^9, 3.4680306624455*^9}, {3.468031030751*^9, 
   3.4680310547320004`*^9}, 3.468031125557*^9, {3.4680312544744997`*^9, 
   3.4680312743035*^9}, {3.4688530827609997`*^9, 3.468853111349*^9}, {
   3.468853193159*^9, 3.468853215678*^9}, {3.468853248796*^9, 
   3.4688532831870003`*^9}, {3.468853387758*^9, 3.4688534372019997`*^9}, {
   3.468853468929*^9, 3.4688534793599997`*^9}, {3.468853562748*^9, 
   3.468853619392*^9}, {3.468853950534*^9, 3.468853962612*^9}, 
   3.468862832283*^9, {3.4688774688120003`*^9, 3.4688774809630003`*^9}, {
   3.469296433404113*^9, 3.469296454427534*^9}, {3.469296534485194*^9, 
   3.469296590051326*^9}, {3.469296625150651*^9, 3.46929667201215*^9}, 
   3.4693075964405174`*^9, {3.4693232782345*^9, 3.4693233255804996`*^9}, {
   3.4693234736825*^9, 3.46932347429*^9}}],

Cell["\<\
A bit more stable - make a table of points - the stobed sections\
\>", "Text",
 CellChangeTimes->{{3.469323341533*^9, 3.469323354171*^9}, {
  3.500751283186263*^9, 3.500751303433421*^9}}],

Cell[BoxData[
 RowBox[{"tab1", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"x1", "[", 
       RowBox[{"n", " ", "2", 
        RowBox[{"\[Pi]", "/", 
         SubscriptBox["\[Omega]", "D"]}]}], "]"}], ",", 
      RowBox[{
       RowBox[{"x1", "'"}], "[", 
       RowBox[{"n", " ", "2", 
        RowBox[{"\[Pi]", "/", 
         SubscriptBox["\[Omega]", "D"]}]}], "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "40", ",", "50"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.468857249308*^9, 3.468857251201*^9}, {3.46885728616*^9, 
   3.468857348133*^9}, {3.468857380652*^9, 3.4688574111940002`*^9}, {
   3.468857445472*^9, 3.468857486389*^9}, {3.468857912574*^9, 
   3.468857924229*^9}, {3.468858375179*^9, 3.468858397817*^9}, {
   3.468859646458*^9, 3.468859646602*^9}, 3.4688600821140003`*^9, {
   3.468860361554*^9, 3.468860361825*^9}, 3.468862832336*^9, {
   3.469296747655095*^9, 3.469296784470387*^9}, {3.4693233715150003`*^9, 
   3.4693234007200003`*^9}}],

Cell["\<\
Note these results are not completely stable (i.e., identical), but pretty \
good for numerical integration!\
\>", "Text",
 CellChangeTimes->{{3.4693234944495*^9, 3.4693235172305*^9}, {
  3.5007513185382853`*^9, 3.5007513264257364`*^9}}],

Cell[BoxData[
 RowBox[{"phplot1a", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"tab1", ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "7"}], ",", "7"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<x\>\"", ",", "\"\<\!\(\*OverscriptBox[\(x\), \(.\)]\)\>\""}], 
      "}"}]}], ",", 
    RowBox[{"BaseStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "\[Rule]", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"Red", ",", 
       RowBox[{"PointSize", "[", "0.02`", "]"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4692972454943933`*^9, 3.469297280921312*^9}, {
   3.469297318173395*^9, 3.4692973483900146`*^9}, 3.4693076066115217`*^9}],

Cell["Put the points on the attractor", "Text",
 CellChangeTimes->{{3.5007513567294693`*^9, 3.500751370409252*^9}}],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"pp1b", ",", "phplot1a"}], "]"}]], "Input",
 CellChangeTimes->{{3.4692973289059887`*^9, 3.4692973536783123`*^9}, {
  3.4693234327145*^9, 3.469323433354*^9}}]
}, Closed]],

Cell[CellGroupData[{

Cell["Part b)", "Subsection",
 CellChangeTimes->{{3.4692953123436213`*^9, 3.469295314387182*^9}, {
  3.469323992674*^9, 3.469323993218*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Clear", "[", "x", "]"}], ";", 
  RowBox[{"Clear", "[", "x1", "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.469295506972704*^9, 3.4692955302318563`*^9}, {
  3.4692957167886686`*^9, 3.4692957229817495`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"x0", "=", "0"}], ";", 
  RowBox[{"xD0", "=", "0"}], ";", 
  RowBox[{"\[Gamma]", "=", "0.25"}], ";", 
  RowBox[{"f", "=", "11"}], ";", 
  RowBox[{
   SubscriptBox["\[Omega]", "D"], "=", "1"}], ";"}]], "Input",
 CellChangeTimes->{
  3.467841596705*^9, 3.467841751833*^9, {3.4678606728055*^9, 
   3.4678606753824997`*^9}, {3.4678608832904997`*^9, 3.467860883666*^9}, {
   3.467861140986*^9, 3.467861141658*^9}, {3.467906653521*^9, 
   3.467906659006*^9}, {3.4679067469315*^9, 3.4679067673445*^9}, {
   3.4679068282165003`*^9, 3.4679068368570004`*^9}, 3.4679069646324997`*^9, {
   3.4679171801115*^9, 3.46791718571*^9}, {3.4679173367485*^9, 
   3.4679173369055*^9}, 3.4679173798945*^9, {3.4679174573085003`*^9, 
   3.4679174580825*^9}, 3.467917517924*^9, {3.4679180034065*^9, 
   3.4679180046055*^9}, {3.467919220707*^9, 3.4679192313535*^9}, {
   3.467919405164*^9, 3.4679194080345*^9}, {3.4679195475295*^9, 
   3.4679195477135*^9}, {3.467919710981*^9, 3.467919711572*^9}, 
   3.4679203548265*^9, {3.4679209260030003`*^9, 3.4679209262019997`*^9}, {
   3.4679211605179996`*^9, 3.467921182793*^9}, {3.4679217631305*^9, 
   3.4679217643915*^9}, {3.4679245217539997`*^9, 3.4679245505655003`*^9}, {
   3.46802846966*^9, 3.468028511933*^9}, 3.4680285519765*^9, 
   3.4680288580175*^9, {3.468029178826*^9, 3.468029193855*^9}, {
   3.4680292246275*^9, 3.4680292496805*^9}, 3.468029294163*^9, {
   3.468029382799*^9, 3.468029434977*^9}, 3.4680294803780003`*^9, {
   3.468029536988*^9, 3.4680295542825003`*^9}, {3.468029664835*^9, 
   3.4680296658265*^9}, {3.4680301131305*^9, 3.4680301490439997`*^9}, 
   3.4680302945945*^9, {3.4680304101875*^9, 3.4680304104674997`*^9}, {
   3.468030935449*^9, 3.4680309503905*^9}, {3.46803099105*^9, 
   3.4680310040235*^9}, {3.4688527359119997`*^9, 3.468852781068*^9}, {
   3.468852908909*^9, 3.46885292338*^9}, {3.468852969242*^9, 
   3.4688529869849997`*^9}, 3.468853068932*^9, 3.468855115584*^9, {
   3.468857014001*^9, 3.468857018488*^9}, {3.468857072519*^9, 
   3.468857083652*^9}, 3.4688628322790003`*^9, {3.4692952743427525`*^9, 
   3.469295345976575*^9}, {3.4693240013785*^9, 3.46932400341*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"sol2", "=", 
   RowBox[{"NDSolve", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"x", "''"}], "[", "t", "]"}], "+", 
         RowBox[{"\[Gamma]", " ", 
          RowBox[{
           RowBox[{"x", "'"}], "[", "t", "]"}]}], "+", 
         SuperscriptBox[
          RowBox[{"x", "[", "t", "]"}], "3"]}], "\[Equal]", 
        RowBox[{"f", " ", 
         RowBox[{"Cos", "[", 
          RowBox[{
           SubscriptBox["\[Omega]", "D"], "t"}], "]"}]}]}], ",", 
       RowBox[{
        RowBox[{"x", "[", "0", "]"}], "\[Equal]", "x0"}], ",", 
       RowBox[{
        RowBox[{
         RowBox[{"x", "'"}], "[", "0", "]"}], "\[Equal]", "xD0"}]}], "}"}], 
     ",", "x", ",", 
     RowBox[{"{", 
      RowBox[{"t", ",", "0", ",", 
       RowBox[{"100", " ", "\[Pi]"}]}], "}"}], ",", 
     RowBox[{"MaxSteps", "\[Rule]", "500000"}]}], "]"}]}], ";", 
  RowBox[{
   RowBox[{"x2", "[", "t_", "]"}], ":=", 
   RowBox[{
    RowBox[{"x", "[", "t", "]"}], "/.", 
    RowBox[{"sol2", "[", 
     RowBox[{"[", "1", "]"}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.4680285262365*^9, 3.4680285394585*^9}, {
   3.4680285703585*^9, 3.468028573014*^9}, 3.4680288969895*^9, {
   3.4680290470745*^9, 3.4680290478665*^9}, {3.4680293288145*^9, 
   3.46802933095*^9}, {3.4680296779625*^9, 3.468029679322*^9}, {
   3.4680297339630003`*^9, 3.468029793699*^9}, {3.468029832959*^9, 
   3.468030002321*^9}, {3.4680309622795*^9, 3.4680309640550003`*^9}, {
   3.468031146472*^9, 3.468031148959*^9}, {3.4688527931210003`*^9, 
   3.4688528048310003`*^9}, {3.468852872124*^9, 3.468852894954*^9}, {
   3.46885293219*^9, 3.46885295179*^9}, {3.468853007074*^9, 
   3.468853039137*^9}, {3.468855137976*^9, 3.4688551490950003`*^9}, 
   3.4688628322799997`*^9, {3.469295340407482*^9, 3.469295448341806*^9}, {
   3.469323022951*^9, 3.469323024593*^9}, {3.469324016532*^9, 
   3.469324022163*^9}}],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"x2", "[", "t", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"t", ",", "0", ",", " ", 
     RowBox[{"60", "\[Pi]"}]}], "}"}], ",", 
   RowBox[{"Ticks", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{"5", "\[Pi]"}], ",", 
       RowBox[{"10", "\[Pi]"}], ",", " ", 
       RowBox[{"15", "\[Pi]"}], ",", 
       RowBox[{"20", "\[Pi]"}], ",", 
       RowBox[{"25", "\[Pi]"}], ",", " ", 
       RowBox[{"30", "\[Pi]"}], ",", 
       RowBox[{"35", "\[Pi]"}], ",", 
       RowBox[{"40", "\[Pi]"}], ",", " ", 
       RowBox[{"45", "\[Pi]"}], ",", 
       RowBox[{"50", "\[Pi]"}], ",", 
       RowBox[{"55", "\[Pi]"}], ",", " ", 
       RowBox[{"60", "\[Pi]"}]}], "}"}], "}"}]}], ",", 
   RowBox[{"BaseStyle", "->", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
      RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<t\>\"", ",", "\"\<x\>\""}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.4692957800922513`*^9, 3.4692957973611193`*^9}, {
  3.46929620942362*^9, 3.4692962559887247`*^9}, {3.469297547020994*^9, 
  3.4692975512797127`*^9}, {3.469323035863*^9, 3.4693230613564997`*^9}, {
  3.4693240287609997`*^9, 3.4693240775715*^9}}],

Cell["Look at various time segments", "Text",
 CellChangeTimes->{{3.500751400489972*^9, 3.50075141163361*^9}}],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"x2", "[", "t", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"t", ",", 
     RowBox[{"30", "\[Pi]"}], ",", " ", 
     RowBox[{"60", "\[Pi]"}]}], "}"}], ",", 
   RowBox[{"Ticks", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"30", "\[Pi]"}], ",", 
       RowBox[{"35", "\[Pi]"}], ",", 
       RowBox[{"40", "\[Pi]"}], ",", " ", 
       RowBox[{"45", "\[Pi]"}], ",", 
       RowBox[{"50", "\[Pi]"}], ",", 
       RowBox[{"55", "\[Pi]"}], ",", " ", 
       RowBox[{"60", "\[Pi]"}]}], "}"}], "}"}]}], ",", 
   RowBox[{"BaseStyle", "->", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
      RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<t\>\"", ",", "\"\<x\>\""}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.4692957800922513`*^9, 3.4692957973611193`*^9}, {
  3.46929620942362*^9, 3.4692962559887247`*^9}, {3.469297547020994*^9, 
  3.4692975512797127`*^9}, {3.469323035863*^9, 3.4693231111254997`*^9}, {
  3.4693241043754997`*^9, 3.4693241045039997`*^9}}],

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"x2", "[", "t", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"t", ",", 
     RowBox[{"60", "\[Pi]"}], ",", " ", 
     RowBox[{"90", "\[Pi]"}]}], "}"}], ",", 
   RowBox[{"Ticks", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"60", "\[Pi]"}], ",", 
       RowBox[{"65", "\[Pi]"}], ",", 
       RowBox[{"70", "\[Pi]"}], ",", " ", 
       RowBox[{"75", "\[Pi]"}], ",", 
       RowBox[{"80", "\[Pi]"}], ",", 
       RowBox[{"85", "\[Pi]"}], ",", " ", 
       RowBox[{"90", "\[Pi]"}]}], "}"}], "}"}]}], ",", 
   RowBox[{"BaseStyle", "->", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
      RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<t\>\"", ",", "\"\<x\>\""}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.4692957800922513`*^9, 3.4692957973611193`*^9}, {
  3.46929620942362*^9, 3.4692962559887247`*^9}, {3.469297547020994*^9, 
  3.4692975512797127`*^9}, {3.469323035863*^9, 3.4693231111254997`*^9}, {
  3.4693231785575*^9, 3.4693232011540003`*^9}, {3.4693241135345*^9, 
  3.4693241138145*^9}}],

Cell["\<\
Can see that things do NOT settle down to periodic (at least for these times \
or probably ever).\
\>", "Text",
 CellChangeTimes->{{3.4692976369532647`*^9, 3.469297683799164*^9}, {
  3.4693232193164997`*^9, 3.4693232296115*^9}, {3.4693241244175*^9, 
  3.4693241582054996`*^9}, {3.5007514295456343`*^9, 3.500751450953859*^9}}],

Cell["\<\
We can look at this directly in phase space\
\>", "Text",
 CellChangeTimes->{{3.469297733343811*^9, 3.469297788785145*^9}, {
  3.469297839140977*^9, 3.46929785562986*^9}, {3.4693241723635*^9, 
  3.4693241770905*^9}}],

Cell[BoxData[
 RowBox[{"pp2", "=", 
  RowBox[{"ParametricPlot", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"x2", "[", "t", "]"}], ",", 
      RowBox[{
       RowBox[{"x2", "'"}], "[", "t", "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"t", ",", "0", ",", 
      RowBox[{"100", " ", "\[Pi]"}]}], "}"}], ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "7.5"}], ",", "7.5"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", "Blue", "}"}]}], ",", 
    RowBox[{"Ticks", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "\[Pi]"}], ",", 
        RowBox[{
         RowBox[{"-", "\[Pi]"}], "/", "2"}], ",", "0", ",", 
        RowBox[{"\[Pi]", "/", "2"}], ",", "\[Pi]"}], "}"}], "}"}]}], ",", 
    RowBox[{"BaseStyle", "->", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<x\>\"", ",", "\"\<\!\(\*OverscriptBox[\(x\), \(.\)]\)\>\""}], 
      "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4680286449795*^9, 3.4680286850055*^9}, {
   3.4680287633094997`*^9, 3.468028768587*^9}, 3.4680289031945*^9, {
   3.4680290752205*^9, 3.468029075908*^9}, {3.468029338283*^9, 
   3.468029338875*^9}, {3.468030174417*^9, 3.4680302386730003`*^9}, {
   3.4680306396565*^9, 3.4680306624455*^9}, {3.468031030751*^9, 
   3.4680310547320004`*^9}, 3.468031125557*^9, {3.4680312544744997`*^9, 
   3.4680312743035*^9}, {3.4688530827609997`*^9, 3.468853111349*^9}, {
   3.468853193159*^9, 3.468853215678*^9}, {3.468853248796*^9, 
   3.4688532831870003`*^9}, {3.468853387758*^9, 3.4688534372019997`*^9}, {
   3.468853468929*^9, 3.4688534793599997`*^9}, {3.468853562748*^9, 
   3.468853619392*^9}, {3.468853950534*^9, 3.468853962612*^9}, 
   3.468862832283*^9, {3.4688774688120003`*^9, 3.4688774809630003`*^9}, {
   3.469296433404113*^9, 3.469296454427534*^9}, {3.469296534485194*^9, 
   3.469296590051326*^9}, {3.469307402037056*^9, 3.469307417636756*^9}, {
   3.4693232466225*^9, 3.469323254117*^9}, 3.46932364966*^9, {
   3.4693241807279997`*^9, 3.469324209788*^9}}],

Cell[BoxData[
 RowBox[{"pp2a", "=", 
  RowBox[{"ParametricPlot", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"x2", "[", "t", "]"}], ",", 
      RowBox[{
       RowBox[{"x2", "'"}], "[", "t", "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"t", ",", 
      RowBox[{"50", "\[Pi]"}], ",", 
      RowBox[{"60", " ", "\[Pi]"}]}], "}"}], ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "7.5"}], ",", "7.5"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", "Blue", "}"}]}], ",", 
    RowBox[{"Ticks", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "\[Pi]"}], ",", 
        RowBox[{
         RowBox[{"-", "\[Pi]"}], "/", "2"}], ",", "0", ",", 
        RowBox[{"\[Pi]", "/", "2"}], ",", "\[Pi]"}], "}"}], "}"}]}], ",", 
    RowBox[{"BaseStyle", "->", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<x\>\"", ",", "\"\<\!\(\*OverscriptBox[\(x\), \(.\)]\)\>\""}], 
      "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4680286449795*^9, 3.4680286850055*^9}, {
   3.4680287633094997`*^9, 3.468028768587*^9}, 3.4680289031945*^9, {
   3.4680290752205*^9, 3.468029075908*^9}, {3.468029338283*^9, 
   3.468029338875*^9}, {3.468030174417*^9, 3.4680302386730003`*^9}, {
   3.4680306396565*^9, 3.4680306624455*^9}, {3.468031030751*^9, 
   3.4680310547320004`*^9}, 3.468031125557*^9, {3.4680312544744997`*^9, 
   3.4680312743035*^9}, {3.4688530827609997`*^9, 3.468853111349*^9}, {
   3.468853193159*^9, 3.468853215678*^9}, {3.468853248796*^9, 
   3.4688532831870003`*^9}, {3.468853387758*^9, 3.4688534372019997`*^9}, {
   3.468853468929*^9, 3.4688534793599997`*^9}, {3.468853562748*^9, 
   3.468853619392*^9}, {3.468853950534*^9, 3.468853962612*^9}, 
   3.468862832283*^9, {3.4688774688120003`*^9, 3.4688774809630003`*^9}, {
   3.469296433404113*^9, 3.469296454427534*^9}, {3.469296534485194*^9, 
   3.469296590051326*^9}, {3.469296625150651*^9, 3.46929667201215*^9}, 
   3.4693075964405174`*^9, {3.4693232782345*^9, 3.4693232955205*^9}, {
   3.469324217948*^9, 3.4693242644135*^9}}],

Cell["Go to later times", "Text",
 CellChangeTimes->{{3.469323312729*^9, 3.469323315664*^9}}],

Cell[BoxData[
 RowBox[{"pp2b", "=", 
  RowBox[{"ParametricPlot", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"x2", "[", "t", "]"}], ",", 
      RowBox[{
       RowBox[{"x2", "'"}], "[", "t", "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"t", ",", 
      RowBox[{"90", "\[Pi]"}], ",", 
      RowBox[{"100", " ", "\[Pi]"}]}], "}"}], ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "7.5"}], ",", "7.5"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", "Blue", "}"}]}], ",", 
    RowBox[{"Ticks", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "\[Pi]"}], ",", 
        RowBox[{
         RowBox[{"-", "\[Pi]"}], "/", "2"}], ",", "0", ",", 
        RowBox[{"\[Pi]", "/", "2"}], ",", "\[Pi]"}], "}"}], "}"}]}], ",", 
    RowBox[{"BaseStyle", "->", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<x\>\"", ",", "\"\<\!\(\*OverscriptBox[\(x\), \(.\)]\)\>\""}], 
      "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4680286449795*^9, 3.4680286850055*^9}, {
   3.4680287633094997`*^9, 3.468028768587*^9}, 3.4680289031945*^9, {
   3.4680290752205*^9, 3.468029075908*^9}, {3.468029338283*^9, 
   3.468029338875*^9}, {3.468030174417*^9, 3.4680302386730003`*^9}, {
   3.4680306396565*^9, 3.4680306624455*^9}, {3.468031030751*^9, 
   3.4680310547320004`*^9}, 3.468031125557*^9, {3.4680312544744997`*^9, 
   3.4680312743035*^9}, {3.4688530827609997`*^9, 3.468853111349*^9}, {
   3.468853193159*^9, 3.468853215678*^9}, {3.468853248796*^9, 
   3.4688532831870003`*^9}, {3.468853387758*^9, 3.4688534372019997`*^9}, {
   3.468853468929*^9, 3.4688534793599997`*^9}, {3.468853562748*^9, 
   3.468853619392*^9}, {3.468853950534*^9, 3.468853962612*^9}, 
   3.468862832283*^9, {3.4688774688120003`*^9, 3.4688774809630003`*^9}, {
   3.469296433404113*^9, 3.469296454427534*^9}, {3.469296534485194*^9, 
   3.469296590051326*^9}, {3.469296625150651*^9, 3.46929667201215*^9}, 
   3.4693075964405174`*^9, {3.4693232782345*^9, 3.4693233255804996`*^9}, {
   3.4693234736825*^9, 3.46932347429*^9}, {3.469324272365*^9, 
   3.4693242910745*^9}}],

Cell["Still changing!", "Text",
 CellChangeTimes->{{3.469323341533*^9, 3.469323354171*^9}, {3.469324300771*^9,
   3.4693243044585*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"tab2", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"x2", "[", 
        RowBox[{"n", " ", "2", 
         RowBox[{"\[Pi]", "/", 
          SubscriptBox["\[Omega]", "D"]}]}], "]"}], ",", 
       RowBox[{
        RowBox[{"x2", "'"}], "[", 
        RowBox[{"n", " ", "2", 
         RowBox[{"\[Pi]", "/", 
          SubscriptBox["\[Omega]", "D"]}]}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", "0", ",", "50"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.468857249308*^9, 3.468857251201*^9}, {3.46885728616*^9, 
   3.468857348133*^9}, {3.468857380652*^9, 3.4688574111940002`*^9}, {
   3.468857445472*^9, 3.468857486389*^9}, {3.468857912574*^9, 
   3.468857924229*^9}, {3.468858375179*^9, 3.468858397817*^9}, {
   3.468859646458*^9, 3.468859646602*^9}, 3.4688600821140003`*^9, {
   3.468860361554*^9, 3.468860361825*^9}, 3.468862832336*^9, {
   3.469296747655095*^9, 3.469296784470387*^9}, {3.4693233715150003`*^9, 
   3.4693234007200003`*^9}, {3.469324310324*^9, 3.4693243295695*^9}, 
   3.4693244098715*^9}],

Cell["\<\
Note these results (as before) are not completely stable, but pretty good for \
numerical integration!\
\>", "Text",
 CellChangeTimes->{{3.4693234944495*^9, 3.4693235172305*^9}, {
  3.5007514866579013`*^9, 3.5007514917701936`*^9}}],

Cell[BoxData[
 RowBox[{"phplot2a", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"tab2", ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "7.5"}], ",", "7.5"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<x\>\"", ",", "\"\<\!\(\*OverscriptBox[\(x\), \(.\)]\)\>\""}], 
      "}"}]}], ",", 
    RowBox[{"BaseStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "\[Rule]", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"Red", ",", 
       RowBox[{"PointSize", "[", "0.02`", "]"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4692972454943933`*^9, 3.469297280921312*^9}, {
   3.469297318173395*^9, 3.4692973483900146`*^9}, 3.4693076066115217`*^9, {
   3.4693243378529997`*^9, 3.469324353002*^9}}],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"pp2", ",", "phplot2a"}], "]"}]], "Input",
 CellChangeTimes->{{3.4692973289059887`*^9, 3.4692973536783123`*^9}, {
  3.4693234327145*^9, 3.469323433354*^9}, {3.4693243706835003`*^9, 
  3.469324396408*^9}}],

Cell["Probe at twice the frequency", "Text",
 CellChangeTimes->{{3.4693244631025*^9, 3.46932447966*^9}, {
  3.4693245545395*^9, 3.4693245556105003`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"tab3", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"x2", "[", 
        RowBox[{"n", " ", 
         RowBox[{"\[Pi]", "/", 
          SubscriptBox["\[Omega]", "D"]}]}], "]"}], ",", 
       RowBox[{
        RowBox[{"x2", "'"}], "[", 
        RowBox[{"n", " ", 
         RowBox[{"\[Pi]", "/", 
          SubscriptBox["\[Omega]", "D"]}]}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", "0", ",", "100"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.468857249308*^9, 3.468857251201*^9}, {3.46885728616*^9, 
   3.468857348133*^9}, {3.468857380652*^9, 3.4688574111940002`*^9}, {
   3.468857445472*^9, 3.468857486389*^9}, {3.468857912574*^9, 
   3.468857924229*^9}, {3.468858375179*^9, 3.468858397817*^9}, {
   3.468859646458*^9, 3.468859646602*^9}, 3.4688600821140003`*^9, {
   3.468860361554*^9, 3.468860361825*^9}, 3.468862832336*^9, {
   3.469296747655095*^9, 3.469296784470387*^9}, {3.4693233715150003`*^9, 
   3.4693234007200003`*^9}, {3.469324310324*^9, 3.4693243295695*^9}, 
   3.4693244098715*^9, {3.469324522641*^9, 3.4693245622279997`*^9}}],

Cell[BoxData[
 RowBox[{"phplot3a", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"tab3", ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "7.5"}], ",", "7.5"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<x\>\"", ",", "\"\<\!\(\*OverscriptBox[\(x\), \(.\)]\)\>\""}], 
      "}"}]}], ",", 
    RowBox[{"BaseStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "\[Rule]", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"Red", ",", 
       RowBox[{"PointSize", "[", "0.02`", "]"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4692972454943933`*^9, 3.469297280921312*^9}, {
   3.469297318173395*^9, 3.4692973483900146`*^9}, 3.4693076066115217`*^9, {
   3.4693243378529997`*^9, 3.469324353002*^9}, {3.469324565782*^9, 
   3.4693245695810003`*^9}}],

Cell["Here see both lobes in the phase portrait", "Text",
 CellChangeTimes->{{3.4693245917105*^9, 3.469324600901*^9}, {
  3.5007515209778643`*^9, 3.500751526186162*^9}}],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"pp2", ",", "phplot3a"}], "]"}]], "Input",
 CellChangeTimes->{{3.4692973289059887`*^9, 3.4692973536783123`*^9}, {
  3.4693234327145*^9, 3.469323433354*^9}, {3.4693243706835003`*^9, 
  3.469324396408*^9}, {3.4693245772735*^9, 3.469324577873*^9}}]
}, Closed]]
}, Open  ]]
},
WindowSize->{1231, 702},
WindowMargins->{{39, Automatic}, {Automatic, 31}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
Magnification->1.1499998569488525`,
FrontEndVersion->"7.0 for Microsoft Windows (64-bit) (February 18, 2009)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[567, 22, 157, 2, 81, "Section"],
Cell[CellGroupData[{
Cell[749, 28, 97, 1, 41, "Subsection"],
Cell[849, 31, 240, 5, 33, "Input"],
Cell[1092, 38, 2120, 36, 33, "Input"],
Cell[3215, 76, 1917, 47, 57, "Input"],
Cell[5135, 125, 1064, 27, 54, "Input"],
Cell[6202, 154, 1119, 29, 54, "Input"],
Cell[7324, 185, 1168, 30, 62, "Input"],
Cell[8495, 217, 310, 6, 36, "Text"],
Cell[8808, 225, 984, 25, 37, "Input"],
Cell[9795, 252, 91, 1, 36, "Text"],
Cell[9889, 255, 931, 24, 37, "Input"],
Cell[10823, 281, 101, 1, 36, "Text"],
Cell[10927, 284, 288, 5, 36, "Text"],
Cell[11218, 291, 2456, 58, 93, "Input"],
Cell[13677, 351, 97, 1, 36, "Text"],
Cell[13777, 354, 2490, 59, 93, "Input"],
Cell[16270, 415, 93, 1, 36, "Text"],
Cell[16366, 418, 2539, 60, 93, "Input"],
Cell[18908, 480, 197, 4, 36, "Text"],
Cell[19108, 486, 1026, 24, 37, "Input"],
Cell[20137, 512, 247, 5, 36, "Text"],
Cell[20387, 519, 1134, 31, 93, "Input"],
Cell[21524, 552, 115, 1, 36, "Text"],
Cell[21642, 555, 204, 4, 37, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[21883, 564, 140, 2, 32, "Subsection"],
Cell[22026, 568, 240, 5, 33, "Input"],
Cell[22269, 575, 2161, 36, 33, "Input"],
Cell[24433, 613, 1961, 48, 57, "Input"],
Cell[26397, 663, 1351, 34, 54, "Input"],
Cell[27751, 699, 110, 1, 32, "Text"],
Cell[27864, 702, 1172, 30, 54, "Input"],
Cell[29039, 734, 1213, 31, 62, "Input"],
Cell[30255, 767, 335, 6, 36, "Text"],
Cell[30593, 775, 226, 5, 36, "Text"],
Cell[30822, 782, 2510, 59, 93, "Input"],
Cell[33335, 843, 2539, 60, 93, "Input"],
Cell[35877, 905, 93, 1, 36, "Text"],
Cell[35973, 908, 2588, 61, 93, "Input"],
Cell[38564, 971, 135, 2, 36, "Text"],
Cell[38702, 975, 1123, 26, 37, "Input"],
Cell[39828, 1003, 241, 5, 36, "Text"],
Cell[40072, 1010, 1187, 32, 93, "Input"],
Cell[41262, 1044, 251, 5, 37, "Input"],
Cell[41516, 1051, 153, 2, 36, "Text"],
Cell[41672, 1055, 1159, 26, 37, "Input"],
Cell[42834, 1083, 1236, 33, 93, "Input"],
Cell[44073, 1118, 169, 2, 36, "Text"],
Cell[44245, 1122, 292, 5, 37, "Input"]
}, Closed]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)