(* 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[     35520,       1079]
NotebookOptionsPosition[     32702,        996]
NotebookOutlinePosition[     33107,       1013]
CellTagsIndexPosition[     33064,       1010]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Consider the baker's map", "Subtitle"],

Cell["First the matrix", "Text",
 CellChangeTimes->{{3.5007520748105416`*^9, 3.5007520780097246`*^9}}],

Cell[BoxData[
 RowBox[{"m1", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"2", ",", "0"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"0", ",", "a"}], "}"}]}], "}"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"MatrixForm", "[", "m1", "]"}]], "Input",
 CellChangeTimes->{{3.5007520827979984`*^9, 3.500752087538269*^9}}],

Cell[BoxData[
 RowBox[{"Det", "[", "m1", "]"}]], "Input"],

Cell[BoxData[
 RowBox[{"Eigenvalues", "[", "m1", "]"}]], "Input"],

Cell["e.g., ", "Text",
 CellChangeTimes->{{3.500752172114107*^9, 3.5007521741942263`*^9}}],

Cell[BoxData[
 RowBox[{"a", "=", 
  RowBox[{"1", "/", "2"}]}]], "Input",
 CellChangeTimes->{{3.4693254727125*^9, 3.4693254739275*^9}}],

Cell[BoxData[
 RowBox[{"Eigenvalues", "[", "m1", "]"}]], "Input"],

Cell[CellGroupData[{

Cell["Showing the regions using Mma", "Subsection",
 CellChangeTimes->{{3.4693748738272934`*^9, 3.4693748779456406`*^9}, {
  3.5007522451382837`*^9, 3.500752247313408*^9}}],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{
   RowBox[{"Graphics", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"GrayLevel", "[", ".9", "]"}], ",", 
      RowBox[{"Rectangle", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.005", ",", "0.001"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"0.5", ",", "0.5"}], "}"}]}], "]"}], ",", 
      RowBox[{"GrayLevel", "[", ".4", "]"}], ",", 
      RowBox[{"Rectangle", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.5", ",", "0"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"1.0", ",", "0.5"}], "}"}]}], "]"}]}], "}"}], "]"}], ",", 
   RowBox[{"Axes", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"0", ",", "1"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", 
   RowBox[{"AspectRatio", "\[Rule]", 
    RowBox[{"1", "/", "1"}]}]}], "]"}]], "Input"],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{
   RowBox[{"Graphics", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"GrayLevel", "[", ".9", "]"}], ",", 
      RowBox[{"Rectangle", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.005", ",", "0.001"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"1.0", ",", "0.25"}], "}"}]}], "]"}], ",", 
      RowBox[{"GrayLevel", "[", ".4", "]"}], ",", 
      RowBox[{"Rectangle", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "0.5"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"1.0", ",", "0.75"}], "}"}]}], "]"}]}], "}"}], "]"}], ",", 
   RowBox[{"Axes", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"0", ",", "1"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", 
   RowBox[{"AspectRatio", "\[Rule]", 
    RowBox[{"1", "/", "1"}]}]}], "]"}]], "Input"],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{
   RowBox[{"Graphics", "[", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"GrayLevel", "[", ".9", "]"}], ",", 
      RowBox[{"Rectangle", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.005", ",", "0.001"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"1.0", ",", "0.125"}], "}"}]}], "]"}], ",", 
      RowBox[{"GrayLevel", "[", ".4", "]"}], ",", 
      RowBox[{"Rectangle", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "0.25"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"1.0", ",", "0.375"}], "}"}]}], "]"}], ",", 
      RowBox[{"GrayLevel", "[", ".9", "]"}], ",", 
      RowBox[{"Rectangle", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0.005", ",", "0.5"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"1.0", ",", "0.625"}], "}"}]}], "]"}], ",", 
      RowBox[{"GrayLevel", "[", ".4", "]"}], ",", 
      RowBox[{"Rectangle", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"0", ",", "0.75"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"1.0", ",", "0.875"}], "}"}]}], "]"}]}], "}"}], "]"}], ",", 
   RowBox[{"Axes", "\[Rule]", "True"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"0", ",", "1"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", 
   RowBox[{"AspectRatio", "\[Rule]", 
    RowBox[{"1", "/", "1"}]}]}], "]"}]], "Input"]
}, Closed]],

Cell[CellGroupData[{

Cell["Specific initial points and the 2 exponents", "Subsection",
 CellChangeTimes->{{3.469374424287857*^9, 3.469374441416437*^9}}],

Cell["\<\
Express the map using complex notation (handy for 2-D case) - initial point \
(0.1, 0) \
\>", "Text",
 CellChangeTimes->{{3.4693263000295*^9, 3.4693263130115*^9}, {
  3.469372489996502*^9, 3.4693724912598915`*^9}, {3.469373575847619*^9, 
  3.469373591993309*^9}, {3.5007526472822847`*^9, 3.500752653929665*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"A", "=", 
   RowBox[{"NestList", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"Mod", "[", 
        RowBox[{
         RowBox[{"2", "*", 
          RowBox[{"Re", "[", "#", "]"}]}], ",", "1"}], "]"}], "+", 
       RowBox[{"\[ImaginaryI]", "*", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Re", "[", "#", "]"}], "<", "0.5"}], ",", 
          RowBox[{"a", "*", 
           RowBox[{"Im", "[", "#", "]"}]}], ",", 
          RowBox[{
           RowBox[{"a", "*", 
            RowBox[{"Im", "[", "#", "]"}]}], "+", "0.5"}]}], "]"}]}]}], "&"}],
      ",", "0.1", ",", "40"}], "]"}]}], ";", 
  RowBox[{"LA", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Re", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A", ",", "n"}], "]"}], "]"}], ",", 
       RowBox[{"Im", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A", ",", "n"}], "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", "1", ",", "41"}], "}"}]}], "]"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{"plot1", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"LA", ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"Joined", "\[Rule]", "True"}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"Thickness", "[", "0.02`", "]"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}]}], "]"}]}]], "Input"],

Cell["Now displace the initial condition in y, (0.1, 0.35)", "Text",
 CellChangeTimes->{{3.46937250512598*^9, 3.4693725097896028`*^9}, {
   3.4693726478849573`*^9, 3.469372664794165*^9}, {3.4693735988259773`*^9, 
   3.4693736228495154`*^9}, 3.469536526173563*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"A2", "=", 
   RowBox[{"NestList", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"Mod", "[", 
        RowBox[{
         RowBox[{"2", "*", 
          RowBox[{"Re", "[", "#", "]"}]}], ",", "1"}], "]"}], "+", 
       RowBox[{"\[ImaginaryI]", "*", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Re", "[", "#", "]"}], "<", "0.5"}], ",", 
          RowBox[{"a", "*", 
           RowBox[{"Im", "[", "#", "]"}]}], ",", 
          RowBox[{
           RowBox[{"a", "*", 
            RowBox[{"Im", "[", "#", "]"}]}], "+", "0.5"}]}], "]"}]}]}], "&"}],
      ",", 
     RowBox[{"0.1", "+", 
      RowBox[{"\[ImaginaryI]", "*", "0.35"}]}], ",", "40"}], "]"}]}], ";", 
  RowBox[{"LA2", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Re", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A2", ",", "n"}], "]"}], "]"}], ",", 
       RowBox[{"Im", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A2", ",", "n"}], "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", "1", ",", "41"}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{3.469325889317*^9}],

Cell[BoxData[
 RowBox[{"plot2", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"LA2", ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"Joined", "\[Rule]", "True"}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Thickness", "[", "0.02`", "]"}], ",", 
       RowBox[{"Hue", "[", "0.2`", "]"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}]}], "]"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"plot1", ",", "plot2"}], "]"}]], "Input"],

Cell["Clearly flowing to the same limit cycle.", "Text",
 CellChangeTimes->{{3.4693726852755213`*^9, 3.469372700624839*^9}}],

Cell["\<\
The log of the 2-D differences between the trajectories at each iteration is \
given by\
\>", "Text",
 CellChangeTimes->{{3.4693263289040003`*^9, 3.469326378947*^9}, {
  3.469372717409255*^9, 3.469372728921243*^9}, {3.5007523887704988`*^9, 
  3.5007523926817226`*^9}}],

Cell[BoxData[
 RowBox[{"diff12", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"n", ",", 
      RowBox[{"Log", "[", 
       RowBox[{"Sqrt", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"LA", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "1"}], "]"}], "]"}], "-", 
            RowBox[{"LA2", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"LA", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "2"}], "]"}], "]"}], "-", 
            RowBox[{"LA2", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], 
        "]"}], "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "41"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.469326383873*^9, 3.469326546041*^9}, {3.469326602354*^9,
   3.469326609249*^9}}],

Cell["which clearly shrinks exponentially", "Text",
 CellChangeTimes->{{3.4693727340376825`*^9, 3.4693727446917315`*^9}}],

Cell[BoxData[
 RowBox[{"pd1", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"diff12", ",", 
    RowBox[{"BaseStyle", "->", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<n\>\"", ",", 
       "\"\<ln(\!\(\*SubscriptBox[\(\[Rho]\), \(n\)]\))\>\""}], "}"}]}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.4693265589605*^9, 3.4693265858255*^9}, {
  3.4693266207415*^9, 3.4693266214455*^9}, {3.4693267745220003`*^9, 
  3.4693267958645*^9}}],

Cell["\<\
in a fashion consistent with an exponent of ln[1/2].  To demonstrate this \
last point consider\
\>", "Text",
 CellChangeTimes->{{3.469372758403164*^9, 3.4693727883686514`*^9}, {
  3.500752414401965*^9, 3.500752431569947*^9}}],

Cell[BoxData[
 RowBox[{"pd2", "=", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"(", "n", ")"}], " ", 
      RowBox[{"Log", "[", 
       RowBox[{"1", "/", "2"}], "]"}]}], "-", "0.2"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "41"}], "}"}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.469326684568*^9, 3.4693267457980003`*^9}, {
  3.469326804631*^9, 3.4693268247720003`*^9}, {3.46932686096*^9, 
  3.4693268783335*^9}, {3.469372801222145*^9, 3.469372806853348*^9}}],

Cell["and superimpose", "Text",
 CellChangeTimes->{{3.500752439769416*^9, 3.5007524473698506`*^9}}],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"pd1", ",", "pd2"}], "]"}]], "Input",
 CellChangeTimes->{{3.469326807318*^9, 3.469326813653*^9}}],

Cell["Pretty good agreement!!", "Text",
 CellChangeTimes->{{3.5007524591945267`*^9, 3.500752464161811*^9}}],

Cell["\<\
So to see the positive exponent we consider an initial condition that is \
slightly displaced in x, (0.11, 0.35).\
\>", "Text",
 CellChangeTimes->{{3.4693268904185*^9, 3.4693268975369997`*^9}, {
  3.469372840250592*^9, 3.4693728754885073`*^9}, {3.4693735497560964`*^9, 
  3.469373560513114*^9}, {3.4695365321482477`*^9, 3.469536533848615*^9}, {
  3.5007526778020306`*^9, 3.500752679562131*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"A3", "=", 
   RowBox[{"NestList", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"Mod", "[", 
        RowBox[{
         RowBox[{"2", "*", 
          RowBox[{"Re", "[", "#", "]"}]}], ",", "1"}], "]"}], "+", 
       RowBox[{"\[ImaginaryI]", "*", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Re", "[", "#", "]"}], "<", "0.5"}], ",", 
          RowBox[{"a", "*", 
           RowBox[{"Im", "[", "#", "]"}]}], ",", 
          RowBox[{
           RowBox[{"a", "*", 
            RowBox[{"Im", "[", "#", "]"}]}], "+", "0.5"}]}], "]"}]}]}], "&"}],
      ",", 
     RowBox[{"0.11", "+", 
      RowBox[{"\[ImaginaryI]", "*", "0.35"}]}], ",", "40"}], "]"}]}], ";", 
  RowBox[{"LA3", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Re", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A3", ",", "n"}], "]"}], "]"}], ",", 
       RowBox[{"Im", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A3", ",", "n"}], "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", "1", ",", "41"}], "}"}]}], "]"}]}]}]], "Input"],

Cell[BoxData[
 RowBox[{"plot3", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"LA3", ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"Joined", "\[Rule]", "True"}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Thickness", "[", "0.02`", "]"}], ",", 
       RowBox[{"Hue", "[", "0.4`", "]"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}]}], "]"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"plot1", ",", "plot2", ",", "plot3"}], "]"}]], "Input"],

Cell["\<\
Clearly a different final limit cycle.  Consider the (log) 2-D separation \
between the nearby 2nd and 3rd sets of initial conditions\
\>", "Text",
 CellChangeTimes->{{3.469372917231164*^9, 3.4693729838646455`*^9}, 
   3.4695365358141775`*^9, {3.500752497465716*^9, 3.500752500073865*^9}}],

Cell[BoxData[
 RowBox[{"diff23", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"n", ",", 
      RowBox[{"Log", "[", 
       RowBox[{"Sqrt", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"LA3", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "1"}], "]"}], "]"}], "-", 
            RowBox[{"LA2", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"LA3", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "2"}], "]"}], "]"}], "-", 
            RowBox[{"LA2", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], 
        "]"}], "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "41"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.469326383873*^9, 3.469326546041*^9}, {3.469326602354*^9,
   3.469326609249*^9}, {3.469373018511135*^9, 3.469373033673849*^9}}],

Cell["\<\
which clearly grows exponentially, at least initially and before finite size \
effects set in.\
\>", "Text",
 CellChangeTimes->{{3.4693727340376825`*^9, 3.4693727446917315`*^9}, {
   3.4693730458258595`*^9, 3.4693730469490232`*^9}, 3.4695365417264633`*^9, {
   3.5007525164578023`*^9, 3.500752554097955*^9}}],

Cell[BoxData[
 RowBox[{"pd3", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"diff23", ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "20"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "5"}], ",", "0"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"BaseStyle", "->", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<n\>\"", ",", 
       "\"\<ln(\!\(\*SubscriptBox[\(\[Rho]\), \(n\)]\))\>\""}], "}"}]}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.4693265589605*^9, 3.4693265858255*^9}, {
  3.4693266207415*^9, 3.4693266214455*^9}, {3.4693267745220003`*^9, 
  3.4693267958645*^9}, {3.4693730793023868`*^9, 3.4693731350394*^9}, {
  3.469373165676818*^9, 3.4693731658484125`*^9}}],

Cell["in a fashion consistent with an exponent of ln[2]", "Text",
 CellChangeTimes->{{3.469372758403164*^9, 3.4693727883686514`*^9}, {
  3.4693731449294834`*^9, 3.4693731455222645`*^9}}],

Cell[BoxData[
 RowBox[{"pd4", "=", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"(", "n", ")"}], " ", 
      RowBox[{"Log", "[", "2", "]"}]}], "-", "5.3"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "0", ",", "7"}], "}"}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.469326684568*^9, 3.4693267457980003`*^9}, {
  3.469326804631*^9, 3.4693268247720003`*^9}, {3.46932686096*^9, 
  3.4693268783335*^9}, {3.469372801222145*^9, 3.469372806853348*^9}, {
  3.4693731492817435`*^9, 3.469373242083213*^9}, {3.4693732789297047`*^9, 
  3.4693732794132953`*^9}}],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"pd3", ",", "pd4"}], "]"}]], "Input",
 CellChangeTimes->{{3.469326807318*^9, 3.469326813653*^9}, {
  3.469373178343612*^9, 3.4693731808551316`*^9}}],

Cell["\<\
But note that there are other (distant) initial conditions, e.g., (0.11, 0), \
that flow to this new limit cycle\
\>", "Text",
 CellChangeTimes->{{3.4693733061667805`*^9, 3.4693733339186473`*^9}, {
   3.469373666918668*^9, 3.469373674515722*^9}, {3.469373942566033*^9, 
   3.4693739484315577`*^9}, 3.4695365468587646`*^9, {3.500752588713935*^9, 
   3.500752595250309*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"A4", "=", 
   RowBox[{"NestList", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"Mod", "[", 
        RowBox[{
         RowBox[{"2", "*", 
          RowBox[{"Re", "[", "#", "]"}]}], ",", "1"}], "]"}], "+", 
       RowBox[{"\[ImaginaryI]", "*", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Re", "[", "#", "]"}], "<", "0.5"}], ",", 
          RowBox[{"a", "*", 
           RowBox[{"Im", "[", "#", "]"}]}], ",", 
          RowBox[{
           RowBox[{"a", "*", 
            RowBox[{"Im", "[", "#", "]"}]}], "+", "0.5"}]}], "]"}]}]}], "&"}],
      ",", "0.11", ",", "40"}], "]"}]}], ";", 
  RowBox[{"LA4", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Re", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A4", ",", "n"}], "]"}], "]"}], ",", 
       RowBox[{"Im", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A4", ",", "n"}], "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", "1", ",", "41"}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.469373908979664*^9, 3.46937391546918*^9}}],

Cell[BoxData[
 RowBox[{"plot4", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"LA4", ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"Joined", "\[Rule]", "True"}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Thickness", "[", "0.02`", "]"}], ",", 
       RowBox[{"Hue", "[", "0.6`", "]"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}]}], "]"}]}]], "Input"],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"plot4", ",", "plot3"}], "]"}]], "Input"],

Cell["\<\
The 2-D differences between the trajectories at each iteration is given by\
\>", "Text",
 CellChangeTimes->{{3.4693263289040003`*^9, 3.469326378947*^9}, {
  3.469372717409255*^9, 3.469372728921243*^9}}],

Cell[BoxData[
 RowBox[{"diff34", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"n", ",", 
      RowBox[{"Log", "[", 
       RowBox[{"Sqrt", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"LA3", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "1"}], "]"}], "]"}], "-", 
            RowBox[{"LA4", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"LA3", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "2"}], "]"}], "]"}], "-", 
            RowBox[{"LA4", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], 
        "]"}], "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "41"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.469326383873*^9, 3.469326546041*^9}, {3.469326602354*^9,
   3.469326609249*^9}, {3.4693734217761574`*^9, 3.4693734394194183`*^9}}],

Cell["which clearly shrinks exponentially", "Text",
 CellChangeTimes->{{3.4693727340376825`*^9, 3.4693727446917315`*^9}}],

Cell[BoxData[
 RowBox[{"pd5", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"diff34", ",", 
    RowBox[{"BaseStyle", "->", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<n\>\"", ",", 
       "\"\<ln(\!\(\*SubscriptBox[\(\[Rho]\), \(n\)]\))\>\""}], "}"}]}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.4693265589605*^9, 3.4693265858255*^9}, {
  3.4693266207415*^9, 3.4693266214455*^9}, {3.4693267745220003`*^9, 
  3.4693267958645*^9}, {3.4693734467200775`*^9, 3.469373452148773*^9}, {
  3.469374007835596*^9, 3.4693740083347898`*^9}}],

Cell["in a fashion consistent with an exponent of ln[1/2]", "Text",
 CellChangeTimes->{{3.469372758403164*^9, 3.4693727883686514`*^9}}],

Cell[BoxData[
 RowBox[{"pd6", "=", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{
     RowBox[{
      RowBox[{"(", "n", ")"}], " ", 
      RowBox[{"Log", "[", 
       RowBox[{"1", "/", "2"}], "]"}]}], "-", "0.2"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "41"}], "}"}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.469326684568*^9, 3.4693267457980003`*^9}, {
  3.469326804631*^9, 3.4693268247720003`*^9}, {3.46932686096*^9, 
  3.4693268783335*^9}, {3.469372801222145*^9, 3.469372806853348*^9}, {
  3.4693739927661896`*^9, 3.4693740142159147`*^9}}],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"pd5", ",", "pd6"}], "]"}]], "Input",
 CellChangeTimes->{{3.469326807318*^9, 3.469326813653*^9}, {
  3.4693740197070436`*^9, 3.469374023575794*^9}}],

Cell["Other initial conditions", "Text",
 CellChangeTimes->{{3.5007527508822107`*^9, 3.5007527532253447`*^9}, {
  3.500752804569281*^9, 3.5007528162739506`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"A5", "=", 
   RowBox[{"NestList", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"Mod", "[", 
        RowBox[{
         RowBox[{"2", "*", 
          RowBox[{"Re", "[", "#", "]"}]}], ",", "1"}], "]"}], "+", 
       RowBox[{"\[ImaginaryI]", "*", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Re", "[", "#", "]"}], "<", "0.5"}], ",", 
          RowBox[{"a", "*", 
           RowBox[{"Im", "[", "#", "]"}]}], ",", 
          RowBox[{
           RowBox[{"a", "*", 
            RowBox[{"Im", "[", "#", "]"}]}], "+", "0.5"}]}], "]"}]}]}], "&"}],
      ",", 
     RowBox[{"0.65", "+", 
      RowBox[{"\[ImaginaryI]", " ", "0.1"}]}], ",", "40"}], "]"}]}], ";", 
  RowBox[{"LA5", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Re", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A5", ",", "n"}], "]"}], "]"}], ",", 
       RowBox[{"Im", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A5", ",", "n"}], "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", "1", ",", "31"}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.469374134849168*^9, 3.4693741498561754`*^9}}],

Cell[BoxData[
 RowBox[{"plot5", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"LA5", ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"Joined", "\[Rule]", "True"}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Thickness", "[", "0.02`", "]"}], ",", 
       RowBox[{"Hue", "[", "0.6`", "]"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}]}], "]"}]}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"A6", "=", 
   RowBox[{"NestList", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"Mod", "[", 
        RowBox[{
         RowBox[{"2", "*", 
          RowBox[{"Re", "[", "#", "]"}]}], ",", "1"}], "]"}], "+", 
       RowBox[{"\[ImaginaryI]", "*", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Re", "[", "#", "]"}], "<", "0.5"}], ",", 
          RowBox[{"a", "*", 
           RowBox[{"Im", "[", "#", "]"}]}], ",", 
          RowBox[{
           RowBox[{"a", "*", 
            RowBox[{"Im", "[", "#", "]"}]}], "+", "0.5"}]}], "]"}]}]}], "&"}],
      ",", 
     RowBox[{"0.66", "+", 
      RowBox[{"\[ImaginaryI]", "*", "0.1"}]}], ",", "40"}], "]"}]}], ";", 
  RowBox[{"LA6", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Re", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A6", ",", "n"}], "]"}], "]"}], ",", 
       RowBox[{"Im", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A6", ",", "n"}], "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", "1", ",", "41"}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.469374184113336*^9, 3.469374190212858*^9}}],

Cell[BoxData[
 RowBox[{"plot6", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"LA6", ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"Joined", "\[Rule]", "True"}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Thickness", "[", "0.02`", "]"}], ",", 
       RowBox[{"Hue", "[", "0.6`", "]"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5007527721094246`*^9, 3.500752777820751*^9}}],

Cell["Now initial conditions that map to the origin.", "Text",
 CellChangeTimes->{{3.500752841306382*^9, 3.5007528575063095`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"A7", "=", 
   RowBox[{"NestList", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"Mod", "[", 
        RowBox[{
         RowBox[{"2", "*", 
          RowBox[{"Re", "[", "#", "]"}]}], ",", "1"}], "]"}], "+", 
       RowBox[{"\[ImaginaryI]", "*", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Re", "[", "#", "]"}], "<", "0.5"}], ",", 
          RowBox[{"a", "*", 
           RowBox[{"Im", "[", "#", "]"}]}], ",", 
          RowBox[{
           RowBox[{"a", "*", 
            RowBox[{"Im", "[", "#", "]"}]}], "+", "0.5"}]}], "]"}]}]}], "&"}],
      ",", 
     RowBox[{"0.5", "+", 
      RowBox[{"\[ImaginaryI]", "*", "0.5"}]}], ",", "40"}], "]"}]}], ";", 
  RowBox[{"LA7", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Re", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A7", ",", "n"}], "]"}], "]"}], ",", 
       RowBox[{"Im", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A7", ",", "n"}], "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", "1", ",", "40"}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.4693742243296204`*^9, 3.4693742388062353`*^9}, {
  3.469374276791748*^9, 3.46937427818013*^9}}],

Cell[BoxData[
 RowBox[{"plot7", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"LA7", ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"Joined", "\[Rule]", "True"}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Thickness", "[", "0.02`", "]"}], ",", 
       RowBox[{"Hue", "[", "0.6`", "]"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4693742976486807`*^9, 3.4693743032022095`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"A8", "=", 
   RowBox[{"NestList", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"Mod", "[", 
        RowBox[{
         RowBox[{"2", "*", 
          RowBox[{"Re", "[", "#", "]"}]}], ",", "1"}], "]"}], "+", 
       RowBox[{"\[ImaginaryI]", "*", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Re", "[", "#", "]"}], "<", "0.5"}], ",", 
          RowBox[{"a", "*", 
           RowBox[{"Im", "[", "#", "]"}]}], ",", 
          RowBox[{
           RowBox[{"a", "*", 
            RowBox[{"Im", "[", "#", "]"}]}], "+", "0.5"}]}], "]"}]}]}], "&"}],
      ",", 
     RowBox[{"0.25", "+", 
      RowBox[{"\[ImaginaryI]", "*", "0.5"}]}], ",", "40"}], "]"}]}], ";", 
  RowBox[{"LA8", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Re", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A8", ",", "n"}], "]"}], "]"}], ",", 
       RowBox[{"Im", "[", 
        RowBox[{"Part", "[", 
         RowBox[{"A8", ",", "n"}], "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"n", ",", "1", ",", "40"}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.4693742243296204`*^9, 3.4693742388062353`*^9}, {
  3.469374276791748*^9, 3.46937427818013*^9}, {3.469374320767584*^9, 
  3.4693743366013813`*^9}, {3.4693743719973273`*^9, 3.4693743758036785`*^9}}],

Cell[BoxData[
 RowBox[{"plot8", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"LA8", ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"Joined", "\[Rule]", "True"}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Thickness", "[", "0.02`", "]"}], ",", 
       RowBox[{"Hue", "[", "0.6`", "]"}]}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4693742976486807`*^9, 3.4693743312038507`*^9}}]
}, Closed]]
}, Open  ]]
},
WindowSize->{848, 757},
WindowMargins->{{3, Automatic}, {Automatic, 5}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
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, 44, 0, 49, "Subtitle"],
Cell[614, 24, 102, 1, 29, "Text"],
Cell[719, 27, 199, 7, 31, "Input"],
Cell[921, 36, 132, 2, 31, "Input"],
Cell[1056, 40, 57, 1, 31, "Input"],
Cell[1116, 43, 65, 1, 31, "Input"],
Cell[1184, 46, 90, 1, 29, "Text"],
Cell[1277, 49, 134, 3, 31, "Input"],
Cell[1414, 54, 65, 1, 31, "Input"],
Cell[CellGroupData[{
Cell[1504, 59, 172, 2, 36, "Subsection"],
Cell[1679, 63, 979, 29, 52, "Input"],
Cell[2661, 94, 981, 29, 52, "Input"],
Cell[3645, 125, 1499, 43, 92, "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[5181, 173, 131, 1, 28, "Subsection"],
Cell[5315, 176, 321, 6, 29, "Text"],
Cell[5639, 184, 1091, 34, 52, "Input"],
Cell[6733, 220, 562, 16, 74, "Input"],
Cell[7298, 238, 263, 3, 29, "Text"],
Cell[7564, 243, 1202, 37, 52, "Input"],
Cell[8769, 282, 656, 19, 74, "Input"],
Cell[9428, 303, 88, 2, 31, "Input"],
Cell[9519, 307, 124, 1, 29, "Text"],
Cell[9646, 310, 278, 6, 29, "Text"],
Cell[9927, 318, 1054, 31, 31, "Input"],
Cell[10984, 351, 121, 1, 29, "Text"],
Cell[11108, 354, 629, 17, 31, "Input"],
Cell[11740, 373, 236, 5, 29, "Text"],
Cell[11979, 380, 563, 14, 31, "Input"],
Cell[12545, 396, 99, 1, 29, "Text"],
Cell[12647, 399, 144, 3, 31, "Input"],
Cell[12794, 404, 107, 1, 29, "Text"],
Cell[12904, 407, 404, 7, 29, "Text"],
Cell[13311, 416, 1164, 36, 52, "Input"],
Cell[14478, 454, 656, 19, 74, "Input"],
Cell[15137, 475, 102, 2, 31, "Input"],
Cell[15242, 479, 299, 5, 29, "Text"],
Cell[15544, 486, 1102, 31, 31, "Input"],
Cell[16649, 519, 318, 6, 29, "Text"],
Cell[16970, 527, 968, 26, 52, "Input"],
Cell[17941, 555, 186, 2, 29, "Text"],
Cell[18130, 559, 638, 15, 31, "Input"],
Cell[18771, 576, 195, 4, 31, "Input"],
Cell[18969, 582, 382, 7, 29, "Text"],
Cell[19354, 591, 1161, 35, 52, "Input"],
Cell[20518, 628, 656, 19, 74, "Input"],
Cell[21177, 649, 88, 2, 31, "Input"],
Cell[21268, 653, 212, 4, 29, "Text"],
Cell[21483, 659, 1106, 31, 31, "Input"],
Cell[22592, 692, 121, 1, 29, "Text"],
Cell[22716, 695, 728, 18, 31, "Input"],
Cell[23447, 715, 135, 1, 29, "Text"],
Cell[23585, 718, 616, 15, 31, "Input"],
Cell[24204, 735, 195, 4, 31, "Input"],
Cell[24402, 741, 161, 2, 29, "Text"],
Cell[24566, 745, 1231, 37, 52, "Input"],
Cell[25800, 784, 656, 19, 74, "Input"],
Cell[26459, 805, 1229, 37, 52, "Input"],
Cell[27691, 844, 724, 20, 74, "Input"],
Cell[28418, 866, 130, 1, 29, "Text"],
Cell[28551, 869, 1280, 38, 52, "Input"],
Cell[29834, 909, 726, 20, 74, "Input"],
Cell[30563, 931, 1382, 39, 52, "Input"],
Cell[31948, 972, 726, 20, 74, "Input"]
}, Closed]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)