(* 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[ 13581, 496] NotebookOptionsPosition[ 12348, 448] NotebookOutlinePosition[ 12753, 465] CellTagsIndexPosition[ 12710, 462] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["Program 04-02:Transition matrix: Limiting Power", "Section", FontColor->RGBColor[1, 0, 0]], Cell[TextData[{ "This program works through Example 4.3 in the text, and finally determines \ lim ", Cell[BoxData[ FormBox[ SuperscriptBox["T", "n"], TraditionalForm]]], " as n\[Rule]\[Infinity]. Obviously other matrices can be substituted for T." }], "Text", FontSlant->"Italic", FontColor->RGBColor[0, 0, 1]], Cell[BoxData[ RowBox[{ RowBox[{"t", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"1", "/", "4"}], ",", RowBox[{"1", "/", "2"}], ",", RowBox[{"1", "/", "4"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"1", "/", "2"}], ",", RowBox[{"1", "/", "4"}], ",", RowBox[{"1", "/", "4"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"1", "/", "4"}], ",", RowBox[{"1", "/", "4"}], ",", RowBox[{"1", "/", "2"}]}], "}"}]}], "}"}]}], ";"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"MatrixForm", "[", "t", "]"}]], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { FractionBox["1", "4"], FractionBox["1", "2"], FractionBox["1", "4"]}, { FractionBox["1", "2"], FractionBox["1", "4"], FractionBox["1", "4"]}, { FractionBox["1", "4"], FractionBox["1", "4"], FractionBox["1", "2"]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.442564612002*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Eigenvalues", "[", "t", "]"}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{"1", ",", RowBox[{"-", FractionBox["1", "4"]}], ",", FractionBox["1", "4"]}], "}"}]], "Output", CellChangeTimes->{3.442564612038*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"c", "=", RowBox[{"Transpose", "[", RowBox[{"Eigenvectors", "[", "t", "]"}], "]"}]}], ";"}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"MatrixForm", "[", "c", "]"}]], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", RowBox[{"-", "1"}], RowBox[{"-", FractionBox["1", "2"]}]}, {"1", "1", RowBox[{"-", FractionBox["1", "2"]}]}, {"1", "0", "1"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.44256461209*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d", "=", RowBox[{"DiagonalMatrix", "[", RowBox[{"Eigenvalues", "[", "t", "]"}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{"-", FractionBox["1", "4"]}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", FractionBox["1", "4"]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.44256461213*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"tt", "[", "n_", "]"}], "=", RowBox[{"c", ".", RowBox[{"MatrixPower", "[", RowBox[{"d", ",", "n"}], "]"}], ".", RowBox[{"Inverse", "[", "c", "]"}]}]}]], "Input", CellChangeTimes->{{3.4407481176219997`*^9, 3.4407481581664*^9}, { 3.4407481886956*^9, 3.4407481972911997`*^9}, {3.4407482462752*^9, 3.4407482507524*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ FractionBox["1", "3"], "+", RowBox[{ FractionBox["1", "3"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "n"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}]}], ",", RowBox[{ FractionBox["1", "3"], "+", RowBox[{ FractionBox["1", "3"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}], "-", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "n"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}]}], ",", RowBox[{ FractionBox["1", "3"], "-", FractionBox[ SuperscriptBox["2", RowBox[{ RowBox[{"-", "2"}], " ", "n"}]], "3"]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ FractionBox["1", "3"], "+", RowBox[{ FractionBox["1", "3"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}], "-", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "n"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}]}], ",", RowBox[{ FractionBox["1", "3"], "+", RowBox[{ FractionBox["1", "3"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "n"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}]}], ",", RowBox[{ FractionBox["1", "3"], "-", FractionBox[ SuperscriptBox["2", RowBox[{ RowBox[{"-", "2"}], " ", "n"}]], "3"]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ FractionBox["1", "3"], "-", FractionBox[ SuperscriptBox["4", RowBox[{"-", "n"}]], "3"]}], ",", RowBox[{ FractionBox["1", "3"], "-", FractionBox[ SuperscriptBox["4", RowBox[{"-", "n"}]], "3"]}], ",", RowBox[{ FractionBox["1", "3"], "+", RowBox[{ FractionBox["1", "3"], " ", SuperscriptBox["2", RowBox[{"1", "-", RowBox[{"2", " ", "n"}]}]]}]}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.4425646121949997`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"MatrixForm", "[", "%", "]"}]], "Input", CellChangeTimes->{{3.4407482770072002`*^9, 3.4407482834344*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{ FractionBox["1", "3"], "+", RowBox[{ FractionBox["1", "3"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "n"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}]}], RowBox[{ FractionBox["1", "3"], "+", RowBox[{ FractionBox["1", "3"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}], "-", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "n"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}]}], RowBox[{ FractionBox["1", "3"], "-", FractionBox[ SuperscriptBox["2", RowBox[{ RowBox[{"-", "2"}], " ", "n"}]], "3"]}]}, { RowBox[{ FractionBox["1", "3"], "+", RowBox[{ FractionBox["1", "3"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}], "-", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "n"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}]}], RowBox[{ FractionBox["1", "3"], "+", RowBox[{ FractionBox["1", "3"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"-", "1"}], ")"}], "n"], " ", SuperscriptBox["2", RowBox[{ RowBox[{"-", "1"}], "-", RowBox[{"2", " ", "n"}]}]]}]}], RowBox[{ FractionBox["1", "3"], "-", FractionBox[ SuperscriptBox["2", RowBox[{ RowBox[{"-", "2"}], " ", "n"}]], "3"]}]}, { RowBox[{ FractionBox["1", "3"], "-", FractionBox[ SuperscriptBox["4", RowBox[{"-", "n"}]], "3"]}], RowBox[{ FractionBox["1", "3"], "-", FractionBox[ SuperscriptBox["4", RowBox[{"-", "n"}]], "3"]}], RowBox[{ FractionBox["1", "3"], "+", RowBox[{ FractionBox["1", "3"], " ", SuperscriptBox["2", RowBox[{"1", "-", RowBox[{"2", " ", "n"}]}]]}]}]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.442564612259*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"MatrixForm", "[", RowBox[{"Limit", "[", RowBox[{ RowBox[{"tt", "[", "n", "]"}], ",", RowBox[{"n", "\[Rule]", "Infinity"}]}], "]"}], "]"}]], "Input", CellChangeTimes->{{3.4407482871940002`*^9, 3.44074831309*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { FractionBox["1", "3"], FractionBox["1", "3"], FractionBox["1", "3"]}, { FractionBox["1", "3"], FractionBox["1", "3"], FractionBox["1", "3"]}, { FractionBox["1", "3"], FractionBox["1", "3"], FractionBox["1", "3"]} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.4425646126*^9}] }, Open ]] }, Open ]] }, WindowSize->{796, 599}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic}, FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (November 10, 2008)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[567, 22, 97, 1, 71, "Section"], Cell[667, 25, 320, 9, 29, "Text"], Cell[990, 36, 567, 19, 31, "Input"], Cell[CellGroupData[{ Cell[1582, 59, 63, 1, 31, "Input"], Cell[1648, 62, 881, 28, 101, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[2566, 95, 64, 1, 31, "Input"], Cell[2633, 98, 186, 6, 47, "Output"] }, Open ]], Cell[2834, 107, 145, 4, 31, "Input"], Cell[CellGroupData[{ Cell[3004, 115, 63, 1, 31, "Input"], Cell[3070, 118, 773, 24, 91, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3880, 147, 129, 3, 31, "Input"], Cell[4012, 152, 368, 12, 47, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4417, 169, 371, 9, 31, "Input"], Cell[4791, 180, 2872, 102, 89, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7700, 287, 129, 2, 24, "Input"], Cell[7832, 291, 3317, 112, 69, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11186, 408, 252, 6, 20, "Input"], Cell[11441, 416, 879, 28, 68, "Output"] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)