(* 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[ 31995, 944] NotebookOptionsPosition[ 30563, 892] NotebookOutlinePosition[ 30928, 908] CellTagsIndexPosition[ 30885, 905] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[{ RowBox[{ RowBox[{"F", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "1", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"3", " ", SuperscriptBox["w", "2"]}], ",", "0", ",", "0", ",", RowBox[{"2", " ", "w"}]}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{ RowBox[{"-", "2"}], " ", "w"}], ",", "0", ",", "0"}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Q", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"q1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "q2"}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"G", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"P", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"p11", ",", "p12", ",", "p13", ",", "p14"}], "}"}], ",", RowBox[{"{", RowBox[{"p12", ",", "p22", ",", "p23", ",", "p24"}], "}"}], ",", RowBox[{"{", RowBox[{"p13", ",", "p23", ",", "p33", ",", "p34"}], "}"}], ",", RowBox[{"{", RowBox[{"p14", ",", "p24", ",", "p34", ",", "p44"}], "}"}]}], "}"}]}], ";"}]}], "Input", CellChangeTimes->{{3.4855072535210676`*^9, 3.485507297303158*^9}, { 3.485507375164028*^9, 3.485507385007967*^9}, {3.4855074349776764`*^9, 3.485507552042424*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"Part", " ", RowBox[{"(", "a", ")"}]}], ":"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"H", "=", RowBox[{"{", RowBox[{"{", RowBox[{"0", ",", "0", ",", "1", ",", "0"}], "}"}], "}"}]}], ";"}]}]], "Input", CellChangeTimes->{{3.485507690420081*^9, 3.4855077283114333`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"F", ".", "P"}], " ", "+", " ", RowBox[{"P", ".", RowBox[{"Transpose", "[", "F", "]"}]}], "+", RowBox[{"G", " ", ".", " ", "Q", ".", " ", RowBox[{"Transpose", "[", "G", "]"}]}], "-", RowBox[{ FractionBox["1", "r"], RowBox[{"P", ".", RowBox[{"Transpose", "[", "H", "]"}], ".", "H", ".", "P"}]}]}], "//", "MatrixForm"}]], "Input", CellChangeTimes->{{3.485507638153452*^9, 3.4855076761854324`*^9}, { 3.48550777320292*^9, 3.485507774140438*^9}, 3.485509448200675*^9, { 3.48550988622471*^9, 3.4855098895372734`*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{ RowBox[{"2", " ", "p12"}], "-", FractionBox[ SuperscriptBox["p13", "2"], "r"]}], RowBox[{"p22", "-", FractionBox[ RowBox[{"p13", " ", "p23"}], "r"], "+", RowBox[{"2", " ", "p14", " ", "w"}], "+", RowBox[{"3", " ", "p11", " ", SuperscriptBox["w", "2"]}]}], RowBox[{"p14", "+", "p23", "-", FractionBox[ RowBox[{"p13", " ", "p33"}], "r"]}], RowBox[{"p24", "-", FractionBox[ RowBox[{"p13", " ", "p34"}], "r"], "-", RowBox[{"2", " ", "p12", " ", "w"}]}]}, { RowBox[{"p22", "-", FractionBox[ RowBox[{"p13", " ", "p23"}], "r"], "+", RowBox[{"2", " ", "p14", " ", "w"}], "+", RowBox[{"3", " ", "p11", " ", SuperscriptBox["w", "2"]}]}], RowBox[{"q1", "-", FractionBox[ SuperscriptBox["p23", "2"], "r"], "+", RowBox[{"4", " ", "p24", " ", "w"}], "+", RowBox[{"6", " ", "p12", " ", SuperscriptBox["w", "2"]}]}], RowBox[{"p24", "-", FractionBox[ RowBox[{"p23", " ", "p33"}], "r"], "+", RowBox[{"2", " ", "p34", " ", "w"}], "+", RowBox[{"3", " ", "p13", " ", SuperscriptBox["w", "2"]}]}], RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"p23", " ", "p34"}], "r"]}], "-", RowBox[{"2", " ", "p22", " ", "w"}], "+", RowBox[{"2", " ", "p44", " ", "w"}], "+", RowBox[{"3", " ", "p14", " ", SuperscriptBox["w", "2"]}]}]}, { RowBox[{"p14", "+", "p23", "-", FractionBox[ RowBox[{"p13", " ", "p33"}], "r"]}], RowBox[{"p24", "-", FractionBox[ RowBox[{"p23", " ", "p33"}], "r"], "+", RowBox[{"2", " ", "p34", " ", "w"}], "+", RowBox[{"3", " ", "p13", " ", SuperscriptBox["w", "2"]}]}], RowBox[{ RowBox[{"2", " ", "p34"}], "-", FractionBox[ SuperscriptBox["p33", "2"], "r"]}], RowBox[{"p44", "-", FractionBox[ RowBox[{"p33", " ", "p34"}], "r"], "-", RowBox[{"2", " ", "p23", " ", "w"}]}]}, { RowBox[{"p24", "-", FractionBox[ RowBox[{"p13", " ", "p34"}], "r"], "-", RowBox[{"2", " ", "p12", " ", "w"}]}], RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"p23", " ", "p34"}], "r"]}], "-", RowBox[{"2", " ", "p22", " ", "w"}], "+", RowBox[{"2", " ", "p44", " ", "w"}], "+", RowBox[{"3", " ", "p14", " ", SuperscriptBox["w", "2"]}]}], RowBox[{"p44", "-", FractionBox[ RowBox[{"p33", " ", "p34"}], "r"], "-", RowBox[{"2", " ", "p23", " ", "w"}]}], RowBox[{"q2", "-", FractionBox[ SuperscriptBox["p34", "2"], "r"], "-", RowBox[{"4", " ", "p24", " ", "w"}]}]} }, 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.485510366595351*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "solving", " ", "these", " ", "is", " ", "taking", " ", "a", " ", "very", " ", "long", " ", "time"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"Solve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{ RowBox[{"2", " ", "p12"}], "-", FractionBox[ SuperscriptBox["p11", "2"], "r"]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{"p22", "-", FractionBox[ RowBox[{"p11", " ", "p12"}], "r"], "+", RowBox[{"2", " ", "p14", " ", "w"}], "+", RowBox[{"3", " ", "p11", " ", SuperscriptBox["w", "2"]}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{"p14", "+", "p23", "-", FractionBox[ RowBox[{"p11", " ", "p13"}], "r"]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{"p24", "-", FractionBox[ RowBox[{"p11", " ", "p14"}], "r"], "-", RowBox[{"2", " ", "p12", " ", "w"}]}], "\[Equal]", "0"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"q1", "-", FractionBox[ SuperscriptBox["p23", "2"], "r"], "+", RowBox[{"4", " ", "p24", " ", "w"}], "+", RowBox[{"6", " ", "p12", " ", SuperscriptBox["w", "2"]}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{"p24", "-", FractionBox[ RowBox[{"p23", " ", "p33"}], "r"], "+", RowBox[{"2", " ", "p34", " ", "w"}], "+", RowBox[{"3", " ", "p13", " ", SuperscriptBox["w", "2"]}]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"p23", " ", "p34"}], "r"]}], "-", RowBox[{"2", " ", "p22", " ", "w"}], "+", RowBox[{"2", " ", "p44", " ", "w"}], "+", RowBox[{"3", " ", "p14", " ", SuperscriptBox["w", "2"]}]}], "\[Equal]", "0"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"2", " ", "p34"}], "-", FractionBox[ SuperscriptBox["p33", "2"], "r"]}], "\[Equal]", "0"}], ",", RowBox[{ RowBox[{"p44", "-", FractionBox[ RowBox[{"p33", " ", "p34"}], "r"], "-", RowBox[{"2", " ", "p23", " ", "w"}]}], "\[Equal]", "0"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"q2", "-", FractionBox[ SuperscriptBox["p34", "2"], "r"], "-", RowBox[{"4", " ", "p24", " ", "w"}]}], "\[Equal]", "0"}]}], "\[IndentingNewLine]", "}"}], ",", RowBox[{"{", RowBox[{ "p11", ",", "p12", ",", "p13", ",", "p14", ",", "p22", ",", "p23", ",", "p24", ",", "p33", ",", "p34"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.4855099101157937`*^9, 3.4855100159459505`*^9}, { 3.4855106362585044`*^9, 3.4855106838378477`*^9}, {3.4855113209947834`*^9, 3.4855113290887403`*^9}}], Cell[BoxData["$Aborted"], "Output", CellChangeTimes->{3.485510016914719*^9, 3.48551131794783*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "create", " ", "the", " ", "block", " ", "matrix", " ", "in", " ", "the", " ", "solution", " ", "of", " ", "the", " ", "Riccati", " ", "equation"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"bigM", "=", RowBox[{"ArrayFlatten", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"Transpose", "[", "F", "]"}]}], ",", RowBox[{ FractionBox["1", "r"], RowBox[{ RowBox[{"Transpose", "[", "H", "]"}], ".", "H"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"G", ".", "Q", ".", RowBox[{"Transpose", "[", "G", "]"}]}], ",", "F"}], "}"}]}], "}"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", "bigM", "]"}]}]}]], "Input", CellChangeTimes->{{3.4855113498392715`*^9, 3.48551146560786*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", RowBox[{ RowBox[{"-", "3"}], " ", SuperscriptBox["w", "2"]}], "0", "0", "0", "0", "0", "0"}, { RowBox[{"-", "1"}], "0", "0", RowBox[{"2", " ", "w"}], "0", "0", "0", "0"}, {"0", "0", "0", "0", "0", "0", FractionBox["1", "r"], "0"}, {"0", RowBox[{ RowBox[{"-", "2"}], " ", "w"}], RowBox[{"-", "1"}], "0", "0", "0", "0", "0"}, {"0", "0", "0", "0", "0", "1", "0", "0"}, {"0", "q1", "0", "0", RowBox[{"3", " ", SuperscriptBox["w", "2"]}], "0", "0", RowBox[{"2", " ", "w"}]}, {"0", "0", "0", "0", "0", "0", "0", "1"}, {"0", "0", "0", "q2", "0", RowBox[{ RowBox[{"-", "2"}], " ", "w"}], "0", "0"} }, 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.485511369605403*^9, 3.4855113748086605`*^9}, { 3.485511429013173*^9, 3.485511465998495*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"bigPhi", "=", RowBox[{"Exp", "[", RowBox[{"bigM", " ", "t"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", "bigPhi", "]"}]}], "Input", CellChangeTimes->{{3.4855114745612144`*^9, 3.4855115171404295`*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "3"}], " ", "t", " ", SuperscriptBox["w", "2"]}]], "1", "1", "1", "1", "1", "1"}, { SuperscriptBox["\[ExponentialE]", RowBox[{"-", "t"}]], "1", "1", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "t", " ", "w"}]], "1", "1", "1", "1"}, {"1", "1", "1", "1", "1", "1", SuperscriptBox["\[ExponentialE]", FractionBox["t", "r"]], "1"}, {"1", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "2"}], " ", "t", " ", "w"}]], SuperscriptBox["\[ExponentialE]", RowBox[{"-", "t"}]], "1", "1", "1", "1", "1"}, {"1", "1", "1", "1", "1", SuperscriptBox["\[ExponentialE]", "t"], "1", "1"}, {"1", SuperscriptBox["\[ExponentialE]", RowBox[{"q1", " ", "t"}]], "1", "1", SuperscriptBox["\[ExponentialE]", RowBox[{"3", " ", "t", " ", SuperscriptBox["w", "2"]}]], "1", "1", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "t", " ", "w"}]]}, {"1", "1", "1", "1", "1", "1", "1", SuperscriptBox["\[ExponentialE]", "t"]}, {"1", "1", "1", SuperscriptBox["\[ExponentialE]", RowBox[{"q2", " ", "t"}]], "1", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "2"}], " ", "t", " ", "w"}]], "1", "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.48551151754669*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData["bigPhi"], "Input", CellChangeTimes->{{3.48551163826853*^9, 3.485511639315432*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "3"}], " ", "t", " ", SuperscriptBox["w", "2"]}]], ",", "1", ",", "1", ",", "1", ",", "1", ",", "1", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{ SuperscriptBox["\[ExponentialE]", RowBox[{"-", "t"}]], ",", "1", ",", "1", ",", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "t", " ", "w"}]], ",", "1", ",", "1", ",", "1", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "1", ",", "1", ",", "1", ",", "1", ",", "1", ",", SuperscriptBox["\[ExponentialE]", FractionBox["t", "r"]], ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "2"}], " ", "t", " ", "w"}]], ",", SuperscriptBox["\[ExponentialE]", RowBox[{"-", "t"}]], ",", "1", ",", "1", ",", "1", ",", "1", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "1", ",", "1", ",", "1", ",", "1", ",", SuperscriptBox["\[ExponentialE]", "t"], ",", "1", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", SuperscriptBox["\[ExponentialE]", RowBox[{"q1", " ", "t"}]], ",", "1", ",", "1", ",", SuperscriptBox["\[ExponentialE]", RowBox[{"3", " ", "t", " ", SuperscriptBox["w", "2"]}]], ",", "1", ",", "1", ",", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "t", " ", "w"}]]}], "}"}], ",", RowBox[{"{", RowBox[{ "1", ",", "1", ",", "1", ",", "1", ",", "1", ",", "1", ",", "1", ",", SuperscriptBox["\[ExponentialE]", "t"]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "1", ",", "1", ",", SuperscriptBox["\[ExponentialE]", RowBox[{"q2", " ", "t"}]], ",", "1", ",", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "2"}], " ", "t", " ", "w"}]], ",", "1", ",", "1"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.485511639893572*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ "extract", " ", "out", " ", "the", " ", "pieces", " ", "we", " ", "need", " ", "to", " ", "solve", " ", "for", " ", "P", RowBox[{"(", "infty", ")"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"phiLambdaY", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "1", ",", "1", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", SuperscriptBox["E", RowBox[{"q1", " ", "t"}]], ",", "1", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "1", ",", "1", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "1", ",", "1", ",", SuperscriptBox["E", RowBox[{"q2", " ", "t"}]]}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", "phiLambdaY", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"phiLambdaLambda", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", SuperscriptBox["E", "t"], ",", "1", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{ SuperscriptBox["E", RowBox[{"3", " ", "t", " ", SuperscriptBox["w", "2"]}]], ",", "1", ",", "1", ",", SuperscriptBox["E", RowBox[{"2", " ", "t", " ", "w"}]]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", "1", ",", "1", ",", SuperscriptBox["E", "t"]}], "}"}], ",", RowBox[{"{", RowBox[{"1", ",", SuperscriptBox["E", RowBox[{ RowBox[{"-", "2"}], " ", "t", " ", "w"}]], ",", "1", ",", "1"}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", "phiLambdaLambda", "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"phiLambdaY", " ", ".", " ", "phiLambdaLambda"}], "//", "Simplify"}], "//", "MatrixForm"}]}]}]], "Input", CellChangeTimes->{{3.4855115680479827`*^9, 3.485511858492918*^9}, { 3.485511947604574*^9, 3.485511957589205*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "1", "1", "1"}, {"1", SuperscriptBox["\[ExponentialE]", RowBox[{"q1", " ", "t"}]], "1", "1"}, {"1", "1", "1", "1"}, {"1", "1", "1", SuperscriptBox["\[ExponentialE]", RowBox[{"q2", " ", "t"}]]} }, 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.4855117319115524`*^9, 3.485511760381031*^9}, { 3.485511835257948*^9, 3.485511859321064*^9}, {3.4855119522453175`*^9, 3.4855119578548365`*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", SuperscriptBox["\[ExponentialE]", "t"], "1", "1"}, { SuperscriptBox["\[ExponentialE]", RowBox[{"3", " ", "t", " ", SuperscriptBox["w", "2"]}]], "1", "1", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "t", " ", "w"}]]}, {"1", "1", "1", SuperscriptBox["\[ExponentialE]", "t"]}, {"1", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "2"}], " ", "t", " ", "w"}]], "1", "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.4855117319115524`*^9, 3.485511760381031*^9}, { 3.485511835257948*^9, 3.485511859321064*^9}, {3.4855119522453175`*^9, 3.4855119578704615`*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{"3", "+", SuperscriptBox["\[ExponentialE]", RowBox[{"3", " ", "t", " ", SuperscriptBox["w", "2"]}]]}], RowBox[{"2", "+", SuperscriptBox["\[ExponentialE]", "t"], "+", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "2"}], " ", "t", " ", "w"}]]}], "4", RowBox[{"2", "+", SuperscriptBox["\[ExponentialE]", "t"], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "t", " ", "w"}]]}]}, { RowBox[{"3", "+", SuperscriptBox["\[ExponentialE]", RowBox[{"t", " ", RowBox[{"(", RowBox[{"q1", "+", RowBox[{"3", " ", SuperscriptBox["w", "2"]}]}], ")"}]}]]}], RowBox[{"1", "+", SuperscriptBox["\[ExponentialE]", "t"], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"q1", " ", "t"}]], "+", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "2"}], " ", "t", " ", "w"}]]}], RowBox[{"3", "+", SuperscriptBox["\[ExponentialE]", RowBox[{"q1", " ", "t"}]]}], RowBox[{"2", "+", SuperscriptBox["\[ExponentialE]", "t"], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"t", " ", RowBox[{"(", RowBox[{"q1", "+", RowBox[{"2", " ", "w"}]}], ")"}]}]]}]}, { RowBox[{"3", "+", SuperscriptBox["\[ExponentialE]", RowBox[{"3", " ", "t", " ", SuperscriptBox["w", "2"]}]]}], RowBox[{"2", "+", SuperscriptBox["\[ExponentialE]", "t"], "+", SuperscriptBox["\[ExponentialE]", RowBox[{ RowBox[{"-", "2"}], " ", "t", " ", "w"}]]}], "4", RowBox[{"2", "+", SuperscriptBox["\[ExponentialE]", "t"], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "t", " ", "w"}]]}]}, { RowBox[{"2", "+", SuperscriptBox["\[ExponentialE]", RowBox[{"q2", " ", "t"}]], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"3", " ", "t", " ", SuperscriptBox["w", "2"]}]]}], RowBox[{"2", "+", SuperscriptBox["\[ExponentialE]", "t"], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"t", " ", RowBox[{"(", RowBox[{"q2", "-", RowBox[{"2", " ", "w"}]}], ")"}]}]]}], RowBox[{"3", "+", SuperscriptBox["\[ExponentialE]", RowBox[{"q2", " ", "t"}]]}], RowBox[{"1", "+", SuperscriptBox["\[ExponentialE]", "t"], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"q2", " ", "t"}]], "+", SuperscriptBox["\[ExponentialE]", RowBox[{"2", " ", "t", " ", "w"}]]}]} }, 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.4855117319115524`*^9, 3.485511760381031*^9}, { 3.485511835257948*^9, 3.485511859321064*^9}, {3.4855119522453175`*^9, 3.4855119578704615`*^9}}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ "this", " ", "looks", " ", "like", " ", "it", " ", "diverges", " ", "as", " ", "we", " ", "take", " ", "the", " ", "limit", " ", "of", " ", "t", " ", "going", " ", "to", " ", RowBox[{"infinity", ".", " ", "Can"}], " ", "I", " ", "use", " ", "this", " ", RowBox[{"method", "?", " ", "What"}], " ", "are", " ", "the", " ", "eigenvalues", " ", "of", " ", RowBox[{"F", "?"}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.4855120157781944`*^9, 3.4855120464508543`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Eigenvalues", "[", "F", "]"}]], "Input", CellChangeTimes->{{3.485511981636695*^9, 3.4855119980277395`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{ RowBox[{"-", "\[ImaginaryI]"}], " ", "w"}], ",", RowBox[{"\[ImaginaryI]", " ", "w"}]}], "}"}]], "Output", CellChangeTimes->{{3.485511993824507*^9, 3.4855119982933717`*^9}}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ "Since", " ", "they", " ", "do", " ", "not", " ", "have", " ", "negative", " ", "real", " ", "parts", " ", "we", " ", "cannot", " ", "us", " ", "the", " ", "above", " ", "limit", " ", "to", " ", "determine", " ", "P", RowBox[{"(", "infty", ")"}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.485512049825941*^9, 3.485512066795125*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"Part", " ", RowBox[{"(", "b", ")"}]}], ":"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"H", "=", RowBox[{"{", RowBox[{"{", RowBox[{"1", ",", "0", ",", "0", ",", "0"}], "}"}], "}"}]}], ";"}]}]], "Input", CellChangeTimes->{{3.485507690420081*^9, 3.4855077283114333`*^9}, { 3.4855077982190256`*^9, 3.4855078045316467`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"F", ".", "P"}], " ", "+", " ", RowBox[{"P", ".", RowBox[{"Transpose", "[", "F", "]"}]}], "+", RowBox[{"G", " ", ".", " ", "Q", ".", " ", RowBox[{"Transpose", "[", "G", "]"}]}], "-", RowBox[{ FractionBox["1", "r"], RowBox[{"P", ".", RowBox[{"Transpose", "[", "H", "]"}], ".", "H", ".", "P"}]}]}], "//", "MatrixForm"}]], "Input", CellChangeTimes->{{3.485507638153452*^9, 3.4855076761854324`*^9}, { 3.48550777320292*^9, 3.485507774140438*^9}}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { RowBox[{ RowBox[{"2", " ", "p12"}], "-", FractionBox[ SuperscriptBox["p11", "2"], "r"]}], RowBox[{"p22", "-", FractionBox[ RowBox[{"p11", " ", "p12"}], "r"], "+", RowBox[{"2", " ", "p14", " ", "w"}], "+", RowBox[{"3", " ", "p11", " ", SuperscriptBox["w", "2"]}]}], RowBox[{"p14", "+", "p23", "-", FractionBox[ RowBox[{"p11", " ", "p13"}], "r"]}], RowBox[{"p24", "-", FractionBox[ RowBox[{"p11", " ", "p14"}], "r"], "-", RowBox[{"2", " ", "p12", " ", "w"}]}]}, { RowBox[{"p22", "-", FractionBox[ RowBox[{"p11", " ", "p12"}], "r"], "+", RowBox[{"2", " ", "p14", " ", "w"}], "+", RowBox[{"3", " ", "p11", " ", SuperscriptBox["w", "2"]}]}], RowBox[{"q1", "-", FractionBox[ SuperscriptBox["p12", "2"], "r"], "+", RowBox[{"4", " ", "p24", " ", "w"}], "+", RowBox[{"6", " ", "p12", " ", SuperscriptBox["w", "2"]}]}], RowBox[{"p24", "-", FractionBox[ RowBox[{"p12", " ", "p13"}], "r"], "+", RowBox[{"2", " ", "p34", " ", "w"}], "+", RowBox[{"3", " ", "p13", " ", SuperscriptBox["w", "2"]}]}], RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"p12", " ", "p14"}], "r"]}], "-", RowBox[{"2", " ", "p22", " ", "w"}], "+", RowBox[{"2", " ", "p44", " ", "w"}], "+", RowBox[{"3", " ", "p14", " ", SuperscriptBox["w", "2"]}]}]}, { RowBox[{"p14", "+", "p23", "-", FractionBox[ RowBox[{"p11", " ", "p13"}], "r"]}], RowBox[{"p24", "-", FractionBox[ RowBox[{"p12", " ", "p13"}], "r"], "+", RowBox[{"2", " ", "p34", " ", "w"}], "+", RowBox[{"3", " ", "p13", " ", SuperscriptBox["w", "2"]}]}], RowBox[{ RowBox[{"2", " ", "p34"}], "-", FractionBox[ SuperscriptBox["p13", "2"], "r"]}], RowBox[{"p44", "-", FractionBox[ RowBox[{"p13", " ", "p14"}], "r"], "-", RowBox[{"2", " ", "p23", " ", "w"}]}]}, { RowBox[{"p24", "-", FractionBox[ RowBox[{"p11", " ", "p14"}], "r"], "-", RowBox[{"2", " ", "p12", " ", "w"}]}], RowBox[{ RowBox[{"-", FractionBox[ RowBox[{"p12", " ", "p14"}], "r"]}], "-", RowBox[{"2", " ", "p22", " ", "w"}], "+", RowBox[{"2", " ", "p44", " ", "w"}], "+", RowBox[{"3", " ", "p14", " ", SuperscriptBox["w", "2"]}]}], RowBox[{"p44", "-", FractionBox[ RowBox[{"p13", " ", "p14"}], "r"], "-", RowBox[{"2", " ", "p23", " ", "w"}]}], RowBox[{"q2", "-", FractionBox[ SuperscriptBox["p14", "2"], "r"], "-", RowBox[{"4", " ", "p24", " ", "w"}]}]} }, 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.485507644403572*^9, 3.4855076771073256`*^9}, 3.4855077747967005`*^9, 3.48550780704732*^9}] }, Open ]] }, WindowSize->{1105, 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, 1782, 55, 94, "Input"], Cell[2330, 77, 387, 12, 52, "Input"], Cell[CellGroupData[{ Cell[2742, 93, 604, 15, 47, "Input"], Cell[3349, 110, 3469, 100, 142, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6855, 215, 2918, 79, 214, "Input"], Cell[9776, 296, 98, 1, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9911, 302, 980, 27, 88, "Input"], Cell[10894, 331, 1404, 39, 162, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[12335, 375, 272, 6, 52, "Input"], Cell[12610, 383, 2004, 53, 183, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[14651, 441, 97, 1, 31, "Input"], Cell[14751, 444, 2081, 53, 71, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[16869, 502, 2098, 55, 145, "Input"], Cell[18970, 559, 971, 26, 92, "Output"], Cell[19944, 587, 1207, 32, 100, "Output"], Cell[21154, 621, 3464, 94, 116, "Output"] }, Open ]], Cell[24633, 718, 544, 11, 52, "Input"], Cell[CellGroupData[{ Cell[25202, 733, 132, 2, 31, "Input"], Cell[25337, 737, 253, 6, 30, "Output"] }, Open ]], Cell[25605, 746, 397, 7, 31, "Input"], Cell[26005, 755, 440, 13, 52, "Input"], Cell[CellGroupData[{ Cell[26470, 772, 530, 14, 47, "Input"], Cell[27003, 788, 3544, 101, 128, "Output"] }, Open ]] } ] *) (* End of internal cache information *)