(* 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[     16074,        470]
NotebookOptionsPosition[     14602,        424]
NotebookOutlinePosition[     14949,        439]
CellTagsIndexPosition[     14906,        436]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Consider the cat map", "Subtitle"],

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

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

Cell[BoxData[
 RowBox[{"v1", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"1", "/", "3"}], ",", 
    RowBox[{"2", "/", "3"}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.468526978538*^9, 3.4685269948719997`*^9}}],

Cell[BoxData[
 RowBox[{"m1", ".", "v1"}]], "Input",
 CellChangeTimes->{{3.468526999824*^9, 3.4685270045179996`*^9}}],

Cell["Consider the properties of this matrix", "Text",
 CellChangeTimes->{{3.4685281398929996`*^9, 3.46852814754*^9}}],

Cell[BoxData[
 RowBox[{"mu", "=", 
  RowBox[{"Eigenvalues", "[", "m1", "]"}]}]], "Input",
 CellChangeTimes->{{3.468719541454*^9, 3.46871955059*^9}}],

Cell[BoxData[
 RowBox[{"\[Lambda]1", "=", 
  RowBox[{"N", "[", 
   RowBox[{"Log", "[", 
    RowBox[{"mu", "[", 
     RowBox[{"[", "1", "]"}], "]"}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.468719555046*^9, 3.468719574309*^9}, {
  3.4687196256809998`*^9, 3.468719627514*^9}}],

Cell[BoxData[
 RowBox[{"\[Lambda]2", "=", 
  RowBox[{"N", "[", 
   RowBox[{"Log", "[", 
    RowBox[{"mu", "[", 
     RowBox[{"[", "2", "]"}], "]"}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.468719577701*^9, 3.468719584892*^9}, {
  3.4687196392320004`*^9, 3.4687196416099997`*^9}}],

Cell[BoxData[
 RowBox[{"m1", ".", "m1", ".", "m1"}]], "Input",
 CellChangeTimes->{{3.4685272625039997`*^9, 3.468527264584*^9}, {
  3.46852730867*^9, 3.4685273282200003`*^9}}],

Cell[BoxData[
 RowBox[{"MatrixPower", "[", 
  RowBox[{"m1", ",", "1"}], "]"}]], "Input",
 CellChangeTimes->{{3.468527430078*^9, 3.468527435566*^9}, {
  3.4685280353570004`*^9, 3.468528036124*^9}}],

Cell[BoxData[
 RowBox[{"MatrixPower", "[", 
  RowBox[{"m1", ",", "2"}], "]"}]], "Input",
 CellChangeTimes->{{3.468527430078*^9, 3.468527435566*^9}, {
  3.4685280353570004`*^9, 3.468528047683*^9}}],

Cell[BoxData[
 RowBox[{"MatrixPower", "[", 
  RowBox[{"m1", ",", "3"}], "]"}]], "Input",
 CellChangeTimes->{{3.468527370109*^9, 3.468527385593*^9}, {3.468527922892*^9,
   3.468527933281*^9}}],

Cell["Check the relation - compare to above", "Text",
 CellChangeTimes->{{3.4685281573640003`*^9, 3.468528158803*^9}, {
  3.500758749833331*^9, 3.500758780834104*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"m11", "[", "n_", "]"}], ":=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Fibonacci", "[", 
       RowBox[{
        RowBox[{"2", "n"}], "+", "1"}], "]"}], ",", 
      RowBox[{"Fibonacci", "[", 
       RowBox[{"2", "n"}], "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Fibonacci", "[", 
       RowBox[{"2", "n"}], "]"}], ",", 
      RowBox[{"Fibonacci", "[", 
       RowBox[{
        RowBox[{"2", "n"}], "-", "1"}], "]"}]}], "}"}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.468527939568*^9, 3.468528016663*^9}}],

Cell[BoxData[
 RowBox[{"m11", "[", "1", "]"}]], "Input",
 CellChangeTimes->{{3.468528021835*^9, 3.46852802366*^9}}],

Cell[BoxData[
 RowBox[{"m11", "[", "2", "]"}]], "Input",
 CellChangeTimes->{{3.468528057115*^9, 3.468528058403*^9}}],

Cell[BoxData[
 RowBox[{"m11", "[", "3", "]"}]], "Input",
 CellChangeTimes->{{3.468528062084*^9, 3.468528063626*^9}}],

Cell["Now evaluate several steps", "Text",
 CellChangeTimes->{{3.4685281760620003`*^9, 3.468528184075*^9}}],

Cell[BoxData[
 RowBox[{"A3", "=", 
  RowBox[{"NestList", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"Mod", "[", 
      RowBox[{
       RowBox[{"m1", ".", "#"}], ",", "1"}], "]"}], "&"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"1", "/", "3"}], ",", 
      RowBox[{"2", "/", "3"}]}], "}"}], ",", "10"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.468526898791*^9, 3.468526953826*^9}, {3.468527096954*^9,
   3.4685271156730003`*^9}, {3.468527146479*^9, 3.4685271484309998`*^9}}],

Cell[BoxData[
 RowBox[{"plot22B", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"A3", ",", 
    RowBox[{"Joined", "\[Rule]", "True"}], ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Thickness", "[", "0.015`", "]"}], ",", "Blue"}], "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}], ",", 
    RowBox[{"BaseStyle", "->", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.46851794723*^9, 3.468517949334*^9}, {3.468518476712*^9, 
  3.468518493678*^9}, {3.4685187604230003`*^9, 3.4685187733190002`*^9}, {
  3.468518809501*^9, 3.468518813508*^9}, {3.468518846363*^9, 
  3.468518906984*^9}, {3.468519068775*^9, 3.46851906927*^9}, {
  3.468519321625*^9, 3.468519395212*^9}, {3.468519981332*^9, 
  3.4685199873459997`*^9}, {3.4685203203459997`*^9, 3.4685203975810003`*^9}, {
  3.468520432035*^9, 3.468520524841*^9}, {3.468522142377*^9, 
  3.4685221744230003`*^9}, {3.468524327068*^9, 3.46852435642*^9}, {
  3.468527202903*^9, 3.468527227038*^9}}],

Cell[BoxData[
 RowBox[{"A4", "=", 
  RowBox[{"NestList", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"Mod", "[", 
      RowBox[{
       RowBox[{"m1", ".", "#"}], ",", "1"}], "]"}], "&"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"11", "/", "30"}], ",", 
      RowBox[{"2", "/", "3"}]}], "}"}], ",", "60"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.468526898791*^9, 3.468526953826*^9}, {3.468527096954*^9,
   3.4685271156730003`*^9}, {3.468527146479*^9, 3.4685271484309998`*^9}, {
  3.468718046132*^9, 3.4687180550360003`*^9}, {3.468718118946*^9, 
  3.4687181197200003`*^9}}],

Cell[BoxData[
 RowBox[{"plot44B", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"A4", ",", 
    RowBox[{"Joined", "\[Rule]", "True"}], ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Thickness", "[", "0.015`", "]"}], ",", "Green"}], "}"}]}], 
    ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}], ",", 
    RowBox[{"BaseStyle", "->", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.46851794723*^9, 3.468517949334*^9}, {3.468518476712*^9, 
  3.468518493678*^9}, {3.4685187604230003`*^9, 3.4685187733190002`*^9}, {
  3.468518809501*^9, 3.468518813508*^9}, {3.468518846363*^9, 
  3.468518906984*^9}, {3.468519068775*^9, 3.46851906927*^9}, {
  3.468519321625*^9, 3.468519395212*^9}, {3.468519981332*^9, 
  3.4685199873459997`*^9}, {3.4685203203459997`*^9, 3.4685203975810003`*^9}, {
  3.468520432035*^9, 3.468520524841*^9}, {3.468522142377*^9, 
  3.4685221744230003`*^9}, {3.468524327068*^9, 3.46852435642*^9}, {
  3.468527202903*^9, 3.468527227038*^9}, {3.4687181309300003`*^9, 
  3.468718147273*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"A5", "=", 
   RowBox[{"NestList", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"Mod", "[", 
       RowBox[{
        RowBox[{"m1", ".", "#"}], ",", "1"}], "]"}], "&"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SqrtBox["0.0002"], ",", "0"}], "}"}], ",", "30"}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.468526898791*^9, 3.468526953826*^9}, {3.468527096954*^9,
   3.4685271156730003`*^9}, {3.468527146479*^9, 3.4685271484309998`*^9}, {
  3.468718315925*^9, 3.468718352244*^9}}],

Cell[BoxData[
 RowBox[{"plot55B", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"A5", ",", 
    RowBox[{"Joined", "\[Rule]", "True"}], ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Thickness", "[", "0.015`", "]"}], ",", "Black"}], "}"}]}], 
    ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}], ",", 
    RowBox[{"BaseStyle", "->", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.46851794723*^9, 3.468517949334*^9}, {3.468518476712*^9, 
   3.468518493678*^9}, {3.4685187604230003`*^9, 3.4685187733190002`*^9}, {
   3.468518809501*^9, 3.468518813508*^9}, {3.468518846363*^9, 
   3.468518906984*^9}, {3.468519068775*^9, 3.46851906927*^9}, {
   3.468519321625*^9, 3.468519395212*^9}, {3.468519981332*^9, 
   3.4685199873459997`*^9}, {3.4685203203459997`*^9, 
   3.4685203975810003`*^9}, {3.468520432035*^9, 3.468520524841*^9}, {
   3.468522142377*^9, 3.4685221744230003`*^9}, {3.468524327068*^9, 
   3.46852435642*^9}, {3.468527202903*^9, 3.468527227038*^9}, {
   3.468718358245*^9, 3.468718369708*^9}, 3.468718490222*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"A6", "=", 
   RowBox[{"NestList", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"Mod", "[", 
       RowBox[{
        RowBox[{"m1", ".", "#"}], ",", "1"}], "]"}], "&"}], ",", 
     RowBox[{"{", 
      RowBox[{
       SqrtBox["0.00021"], ",", "0"}], "}"}], ",", "30"}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.468526898791*^9, 3.468526953826*^9}, {3.468527096954*^9,
   3.4685271156730003`*^9}, {3.468527146479*^9, 3.4685271484309998`*^9}, {
  3.468718315925*^9, 3.468718352244*^9}, {3.4687184947320004`*^9, 
  3.468718499675*^9}}],

Cell[BoxData[
 RowBox[{"plot66B", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"A6", ",", 
    RowBox[{"Joined", "\[Rule]", "True"}], ",", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ",", 
    RowBox[{"PlotStyle", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Thickness", "[", "0.015`", "]"}], ",", "Dashed", ",", "Red"}],
       "}"}]}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", 
     RowBox[{"1", " ", 
      FractionBox["1", "1"]}]}], ",", 
    RowBox[{"BaseStyle", "->", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.46851794723*^9, 3.468517949334*^9}, {3.468518476712*^9, 
  3.468518493678*^9}, {3.4685187604230003`*^9, 3.4685187733190002`*^9}, {
  3.468518809501*^9, 3.468518813508*^9}, {3.468518846363*^9, 
  3.468518906984*^9}, {3.468519068775*^9, 3.46851906927*^9}, {
  3.468519321625*^9, 3.468519395212*^9}, {3.468519981332*^9, 
  3.4685199873459997`*^9}, {3.4685203203459997`*^9, 3.4685203975810003`*^9}, {
  3.468520432035*^9, 3.468520524841*^9}, {3.468522142377*^9, 
  3.4685221744230003`*^9}, {3.468524327068*^9, 3.46852435642*^9}, {
  3.468527202903*^9, 3.468527227038*^9}, {3.468718358245*^9, 
  3.468718369708*^9}, {3.468718507932*^9, 3.46871854333*^9}}],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"plot55B", ",", "plot66B"}], "]"}]], "Input",
 CellChangeTimes->{{3.468718653697*^9, 3.46871866852*^9}}],

Cell[BoxData[
 RowBox[{"tabdif", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"n", ",", 
      RowBox[{"Log", "[", 
       RowBox[{"Sqrt", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"A5", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "1"}], "]"}], "]"}], "-", 
            RowBox[{"A6", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"A5", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "2"}], "]"}], "]"}], "-", 
            RowBox[{"A6", "[", 
             RowBox[{"[", 
              RowBox[{"n", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], 
        "]"}], "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "30"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.468719168615*^9, 3.46871929375*^9}, {3.468719327797*^9, 
  3.468719332748*^9}}],

Cell[BoxData[
 RowBox[{"pd", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"tabdif", ",", 
    RowBox[{"BaseStyle", "->", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"FontFamily", "->", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "}"}]}], ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
      "\"\<n\>\"", ",", 
       "\"\<ln(\!\(\*SubscriptBox[\(\[Rho]\), \(n\)]\))\>\""}], "}"}]}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.468719276698*^9, 3.468719309523*^9}, {3.468719370224*^9,
   3.468719371144*^9}, {3.46871997995*^9, 3.468720057921*^9}, {
  3.4687202805769997`*^9, 3.468720282594*^9}}],

Cell[BoxData[
 RowBox[{"pl", "=", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "7.95975"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"n", "-", "1"}], ")"}], " ", "\[Lambda]1"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"n", ",", "1", ",", "10"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4687196030360003`*^9, 3.468719611354*^9}, {
  3.468719668487*^9, 3.4687197648900003`*^9}}],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"pd", ",", "pl"}], "]"}]], "Input",
 CellChangeTimes->{{3.4687197780290003`*^9, 3.4687197840480003`*^9}}]
}, Open  ]]
},
WindowSize->{1195, 684},
WindowMargins->{{61, Automatic}, {Automatic, 114}},
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, 40, 0, 49, "Subtitle"],
Cell[610, 24, 199, 7, 31, "Input"],
Cell[812, 33, 57, 1, 31, "Input"],
Cell[872, 36, 212, 6, 31, "Input"],
Cell[1087, 44, 116, 2, 31, "Input"],
Cell[1206, 48, 118, 1, 29, "Text"],
Cell[1327, 51, 148, 3, 31, "Input"],
Cell[1478, 56, 282, 7, 31, "Input"],
Cell[1763, 65, 287, 7, 31, "Input"],
Cell[2053, 74, 174, 3, 31, "Input"],
Cell[2230, 79, 196, 4, 31, "Input"],
Cell[2429, 85, 196, 4, 31, "Input"],
Cell[2628, 91, 191, 4, 31, "Input"],
Cell[2822, 97, 167, 2, 29, "Text"],
Cell[2992, 101, 600, 19, 31, "Input"],
Cell[3595, 122, 115, 2, 31, "Input"],
Cell[3713, 126, 116, 2, 31, "Input"],
Cell[3832, 130, 116, 2, 31, "Input"],
Cell[3951, 134, 107, 1, 29, "Text"],
Cell[4061, 137, 487, 13, 31, "Input"],
Cell[4551, 152, 1403, 32, 74, "Input"],
Cell[5957, 186, 585, 15, 31, "Input"],
Cell[6545, 203, 1457, 34, 74, "Input"],
Cell[8005, 239, 527, 15, 40, "Input"],
Cell[8535, 256, 1480, 34, 74, "Input"],
Cell[10018, 292, 576, 16, 40, "Input"],
Cell[10597, 310, 1506, 34, 74, "Input"],
Cell[12106, 346, 151, 3, 31, "Input"],
Cell[12260, 351, 1051, 31, 31, "Input"],
Cell[13314, 384, 662, 17, 31, "Input"],
Cell[13979, 403, 452, 13, 31, "Input"],
Cell[14434, 418, 152, 3, 31, "Input"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)