(* 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[ 19460, 640] NotebookOptionsPosition[ 18166, 590] NotebookOutlinePosition[ 18530, 606] CellTagsIndexPosition[ 18487, 603] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "the", " ", "iterations", " ", "we", " ", "want", " ", "to", " ", "perform"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"f", "[", "M_", "]"}], ":=", RowBox[{"Simplify", "[", RowBox[{"Inverse", "[", " ", RowBox[{ RowBox[{"Inverse", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "T"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], " ", "}"}], ".", " ", "M", " ", ".", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"T", ",", "1"}], "}"}]}], "}"}]}], "]"}], " ", "+", " ", RowBox[{ FractionBox["1", "r"], RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], "}"}]}]}], "]"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.48223500595582*^9, 3.4822351187854586`*^9}, 3.482235182146456*^9, {3.4822356320798483`*^9, 3.48223563486117*^9}, { 3.4822356811436043`*^9, 3.4822356835030403`*^9}, {3.48464872992892*^9, 3.4846487392260327`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"iterate", " ", "the", " ", "first", " ", "time"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"P1", "=", RowBox[{"f", "[", RowBox[{ FractionBox["1", "eps"], RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", "P1", "]"}]}]}]], "Input", CellChangeTimes->{{3.4822350543775597`*^9, 3.482235057393262*^9}, { 3.482235101410014*^9, 3.4822351053944907`*^9}, {3.482235175380657*^9, 3.482235178005725*^9}, {3.4822356482521377`*^9, 3.48223569501896*^9}, { 3.482236946721062*^9, 3.4822369802998323`*^9}, {3.48464875238262*^9, 3.4846487547108045`*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { FractionBox[ RowBox[{"r", " ", RowBox[{"(", RowBox[{"1", "+", SuperscriptBox["T", "2"]}], ")"}]}], RowBox[{"1", "+", RowBox[{"eps", " ", "r"}], "+", SuperscriptBox["T", "2"]}]], FractionBox[ RowBox[{"r", " ", "T"}], RowBox[{"1", "+", RowBox[{"eps", " ", "r"}], "+", SuperscriptBox["T", "2"]}]]}, { FractionBox[ RowBox[{"r", " ", "T"}], RowBox[{"1", "+", RowBox[{"eps", " ", "r"}], "+", SuperscriptBox["T", "2"]}]], FractionBox[ RowBox[{"1", "+", RowBox[{"eps", " ", "r"}]}], RowBox[{"eps", "+", RowBox[{ SuperscriptBox["eps", "2"], " ", "r"}], "+", RowBox[{"eps", " ", SuperscriptBox["T", "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.4822351059257545`*^9, {3.482235165880414*^9, 3.482235183755872*^9}, { 3.4822356651275697`*^9, 3.4822356953158426`*^9}, 3.4822369482523413`*^9, 3.482236980799842*^9, 3.4846487583046465`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"P1", " ", RowBox[{"(", " ", RowBox[{"1", " ", "+", " ", RowBox[{"eps", " ", "r"}], " ", "+", " ", SuperscriptBox["T", "2"]}], ")"}]}], "//", "Simplify"}], "//", "MatrixForm"}]], "Input", CellChangeTimes->{{3.4822369498461223`*^9, 3.4822369901750216`*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"r", " ", RowBox[{"(", RowBox[{"1", "+", SuperscriptBox["T", "2"]}], ")"}]}], RowBox[{"r", " ", "T"}]}, { RowBox[{"r", " ", "T"}], RowBox[{ FractionBox["1", "eps"], "+", "r"}]} }, 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.4822369630807514`*^9, 3.4822369903781505`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"iterate", " ", "a", " ", "second", " ", "time"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"P2", "=", RowBox[{"f", "[", RowBox[{"f", "[", RowBox[{ FractionBox["1", "eps"], RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], "]"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.4822359458378806`*^9, 3.4822360071988263`*^9}, { 3.4822370904894476`*^9, 3.4822370927707415`*^9}, {3.4846487839303026`*^9, 3.484648788571046*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ FractionBox[ RowBox[{"r", " ", RowBox[{"(", RowBox[{ SuperscriptBox["T", "2"], "+", RowBox[{"eps", " ", RowBox[{"(", RowBox[{"r", "+", RowBox[{"4", " ", "r", " ", SuperscriptBox["T", "2"]}]}], ")"}]}]}], ")"}]}], RowBox[{ RowBox[{ SuperscriptBox["eps", "2"], " ", SuperscriptBox["r", "2"]}], "+", SuperscriptBox["T", "2"], "+", RowBox[{"eps", " ", "r", " ", RowBox[{"(", RowBox[{"2", "+", RowBox[{"5", " ", SuperscriptBox["T", "2"]}]}], ")"}]}]}]], ",", FractionBox[ RowBox[{"r", " ", RowBox[{"(", RowBox[{"T", "+", RowBox[{"2", " ", "eps", " ", "r", " ", "T"}]}], ")"}]}], RowBox[{ RowBox[{ SuperscriptBox["eps", "2"], " ", SuperscriptBox["r", "2"]}], "+", SuperscriptBox["T", "2"], "+", RowBox[{"eps", " ", "r", " ", RowBox[{"(", RowBox[{"2", "+", RowBox[{"5", " ", SuperscriptBox["T", "2"]}]}], ")"}]}]}]]}], "}"}], ",", RowBox[{"{", RowBox[{ FractionBox[ RowBox[{"r", " ", RowBox[{"(", RowBox[{"T", "+", RowBox[{"2", " ", "eps", " ", "r", " ", "T"}]}], ")"}]}], RowBox[{ RowBox[{ SuperscriptBox["eps", "2"], " ", SuperscriptBox["r", "2"]}], "+", SuperscriptBox["T", "2"], "+", RowBox[{"eps", " ", "r", " ", RowBox[{"(", RowBox[{"2", "+", RowBox[{"5", " ", SuperscriptBox["T", "2"]}]}], ")"}]}]}]], ",", FractionBox[ RowBox[{"r", " ", RowBox[{"(", RowBox[{"2", "+", RowBox[{"eps", " ", "r"}]}], ")"}]}], RowBox[{ RowBox[{ SuperscriptBox["eps", "2"], " ", SuperscriptBox["r", "2"]}], "+", SuperscriptBox["T", "2"], "+", RowBox[{"eps", " ", "r", " ", RowBox[{"(", RowBox[{"2", "+", RowBox[{"5", " ", SuperscriptBox["T", "2"]}]}], ")"}]}]}]]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.4822359537599583`*^9, {3.482235995136018*^9, 3.48223600752696*^9}, 3.482237094130142*^9, 3.484648790586723*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "lets", " ", "take", " ", "eps", " ", "equal", " ", "to", " ", "zero", " ", "now"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"P2p", " ", "=", " ", RowBox[{"P2", " ", "/.", RowBox[{"{", RowBox[{"eps", "\[Rule]", "0"}], "}"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", "P2p", "]"}]}]}]], "Input", CellChangeTimes->{{3.482237097864589*^9, 3.4822371151149206`*^9}, { 3.482237176600476*^9, 3.482237188960088*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"r", FractionBox["r", "T"]}, { FractionBox["r", "T"], FractionBox[ RowBox[{"2", " ", "r"}], SuperscriptBox["T", "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.4822371089273014`*^9, 3.4822371154274263`*^9}, 3.48223718931947*^9, 3.4846487997588325`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"iterate", " ", "a", " ", "third", " ", "time"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"f", "[", "P2p", "]"}], "//", "MatrixForm"}]}]], "Input", CellChangeTimes->{{3.482237211085513*^9, 3.4822372183825283`*^9}, { 3.4846488119778953`*^9, 3.484648816821769*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { FractionBox[ RowBox[{"5", " ", "r"}], "6"], FractionBox["r", RowBox[{"2", " ", "T"}]]}, { FractionBox["r", RowBox[{"2", " ", "T"}]], FractionBox["r", RowBox[{"2", " ", SuperscriptBox["T", "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.482237213382432*^9, 3.482237218648158*^9}, 3.4846488172280293`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"a", " ", "fourth", " ", "time"}], " ", "..."}], " ", "we", " ", "begin", " ", "to", " ", "see", " ", "the", " ", "pattern"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"f", "[", RowBox[{"f", "[", "P2p", "]"}], "]"}], "//", "MatrixForm"}]}]], "Input", CellChangeTimes->{{3.482237211085513*^9, 3.4822372183825283`*^9}, { 3.4822372785711837`*^9, 3.482237280883728*^9}, {3.484648823775072*^9, 3.4846488290252066`*^9}, {3.484649021498884*^9, 3.484649027952174*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { FractionBox[ RowBox[{"7", " ", "r"}], "10"], FractionBox[ RowBox[{"3", " ", "r"}], RowBox[{"10", " ", "T"}]]}, { FractionBox[ RowBox[{"3", " ", "r"}], RowBox[{"10", " ", "T"}]], FractionBox["r", RowBox[{"5", " ", SuperscriptBox["T", "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.482237213382432*^9, 3.482237218648158*^9}, 3.482237281211859*^9, 3.4846488293689656`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"Lets", " ", "do", " ", "the", " ", "general", " ", "case"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Clear", "[", "Pkp1", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Pkp1", "[", "k_", "]"}], ":=", RowBox[{"2", " ", FractionBox["r", RowBox[{ RowBox[{"(", RowBox[{"k", "+", "1"}], ")"}], RowBox[{"(", RowBox[{"k", "+", "2"}], ")"}]}]], RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"2", "k"}], "+", "1"}], ",", FractionBox["3", "T"]}], "}"}], ",", RowBox[{"{", RowBox[{ FractionBox["3", "T"], ",", FractionBox["6", RowBox[{"k", " ", SuperscriptBox["T", "2"]}]]}], "}"}]}], "}"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", RowBox[{"Pkp1", "[", "k", "]"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.4846488444631023`*^9, 3.484648904667768*^9}, { 3.484648944450036*^9, 3.4846489618254814`*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{"1", "+", RowBox[{"2", " ", "k"}]}], ")"}], " ", "r"}], RowBox[{ RowBox[{"(", RowBox[{"1", "+", "k"}], ")"}], " ", RowBox[{"(", RowBox[{"2", "+", "k"}], ")"}]}]], FractionBox[ RowBox[{"6", " ", "r"}], RowBox[{ RowBox[{"(", RowBox[{"1", "+", "k"}], ")"}], " ", RowBox[{"(", RowBox[{"2", "+", "k"}], ")"}], " ", "T"}]]}, { FractionBox[ RowBox[{"6", " ", "r"}], RowBox[{ RowBox[{"(", RowBox[{"1", "+", "k"}], ")"}], " ", RowBox[{"(", RowBox[{"2", "+", "k"}], ")"}], " ", "T"}]], FractionBox[ RowBox[{"12", " ", "r"}], RowBox[{"k", " ", RowBox[{"(", RowBox[{"1", "+", "k"}], ")"}], " ", RowBox[{"(", RowBox[{"2", "+", "k"}], ")"}], " ", SuperscriptBox["T", "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.4846488914018035`*^9, 3.48464890571467*^9}, { 3.484648948590768*^9, 3.484648962153615*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"value", "=", RowBox[{"f", "[", RowBox[{"Pkp1", "[", "k", "]"}], "]"}]}]], "Input", CellChangeTimes->{{3.4846489107773*^9, 3.4846489335122566`*^9}, { 3.4846489689506636`*^9, 3.484648969435051*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{"3", "+", RowBox[{"2", " ", "k"}]}], ")"}], " ", "r"}], RowBox[{"6", "+", RowBox[{"5", " ", "k"}], "+", SuperscriptBox["k", "2"]}]], ",", FractionBox[ RowBox[{"6", " ", "r"}], RowBox[{ RowBox[{"(", RowBox[{"6", "+", RowBox[{"5", " ", "k"}], "+", SuperscriptBox["k", "2"]}], ")"}], " ", "T"}]]}], "}"}], ",", RowBox[{"{", RowBox[{ FractionBox[ RowBox[{"6", " ", "r"}], RowBox[{ RowBox[{"(", RowBox[{"6", "+", RowBox[{"5", " ", "k"}], "+", SuperscriptBox["k", "2"]}], ")"}], " ", "T"}]], ",", FractionBox[ RowBox[{"12", " ", "r"}], RowBox[{ RowBox[{"(", RowBox[{"6", "+", RowBox[{"11", " ", "k"}], "+", RowBox[{"6", " ", SuperscriptBox["k", "2"]}], "+", SuperscriptBox["k", "3"]}], ")"}], " ", SuperscriptBox["T", "2"]}]]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{{3.4846489137930017`*^9, 3.484648935059171*^9}, 3.4846489701538196`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"Does", " ", "this", " ", "equal", " ", RowBox[{ RowBox[{"Pkp1", "[", RowBox[{"k", "+", "1"}], "]"}], "?", " ", "Yes"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"value", " ", "-", " ", RowBox[{"Pkp1", "[", RowBox[{"k", "+", "1"}], "]"}]}], "//", "Simplify"}]}]], "Input", CellChangeTimes->{{3.484648936215451*^9, 3.484648938559261*^9}, { 3.4846489731226454`*^9, 3.4846490036703024`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{{3.484648982747892*^9, 3.4846489964044914`*^9}}] }, Open ]] }, WindowSize->{773, 750}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, ShowSelection->True, FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (February 18, 2009)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[545, 20, 1316, 38, 67, "Input"], Cell[CellGroupData[{ Cell[1886, 62, 869, 23, 90, "Input"], Cell[2758, 87, 1601, 47, 89, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4396, 139, 331, 9, 33, "Input"], Cell[4730, 150, 879, 26, 60, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5646, 181, 641, 18, 69, "Input"], Cell[6290, 201, 2362, 75, 100, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8689, 281, 557, 15, 72, "Input"], Cell[9249, 298, 836, 24, 75, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[10122, 327, 345, 8, 52, "Input"], Cell[10470, 337, 921, 28, 77, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11428, 370, 581, 13, 52, "Input"], Cell[12012, 385, 1002, 30, 79, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[13051, 420, 1155, 35, 110, "Input"], Cell[14209, 457, 1668, 51, 81, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[15914, 513, 230, 5, 31, "Input"], Cell[16147, 520, 1230, 40, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[17414, 565, 494, 13, 52, "Input"], Cell[17911, 580, 239, 7, 30, "Output"] }, Open ]] } ] *) (* End of internal cache information *)