(* 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[     13187,        445]
NotebookOptionsPosition[     12212,        407]
NotebookOutlinePosition[     12638,        425]
CellTagsIndexPosition[     12595,        422]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"F", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"0", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"1", ",", 
       RowBox[{"-", "b"}]}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"GQGT", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"q", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0"}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"P", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"p11", ",", "p12"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"p12", ",", "p22"}], "}"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"H", "=", 
   RowBox[{"{", 
    RowBox[{"{", 
     RowBox[{"0", ",", "a"}], "}"}], "}"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.484049856211033*^9, 3.484049901227234*^9}, {
  3.484049931540122*^9, 3.48404995439979*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"pDot", "=", 
   RowBox[{
    RowBox[{"F", " ", ".", " ", "P"}], "+", 
    RowBox[{"P", ".", 
     RowBox[{"Transpose", "[", "F", "]"}]}], "+", "GQGT", "-", 
    RowBox[{
     FractionBox["1", "r"], 
     RowBox[{"P", ".", 
      RowBox[{"Transpose", "[", "H", "]"}], ".", "H", ".", "P"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "pDot", "]"}]}], "Input",
 CellChangeTimes->{{3.484049903633515*^9, 3.484049928024452*^9}, {
  3.484049973790663*^9, 3.484049992400276*^9}, {3.484050076510728*^9, 
  3.484050090745285*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      RowBox[{"q", "-", 
       FractionBox[
        RowBox[{
         SuperscriptBox["a", "2"], " ", 
         SuperscriptBox["p12", "2"]}], "r"]}], 
      RowBox[{"p11", "-", 
       RowBox[{"b", " ", "p12"}], "-", 
       FractionBox[
        RowBox[{
         SuperscriptBox["a", "2"], " ", "p12", " ", "p22"}], "r"]}]},
     {
      RowBox[{"p11", "-", 
       RowBox[{"b", " ", "p12"}], "-", 
       FractionBox[
        RowBox[{
         SuperscriptBox["a", "2"], " ", "p12", " ", "p22"}], "r"]}], 
      RowBox[{
       RowBox[{"2", " ", "p12"}], "-", 
       RowBox[{"2", " ", "b", " ", "p22"}], "-", 
       FractionBox[
        RowBox[{
         SuperscriptBox["a", "2"], " ", 
         SuperscriptBox["p22", "2"]}], "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.484049905524164*^9, 3.484049928274455*^9}, {
   3.4840499936971674`*^9, 3.484050023088169*^9}, {3.4840500776826177`*^9, 
   3.4840500913077917`*^9}, 3.484064028838688*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"p12Sol", "=", 
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"pDot", "[", 
      RowBox[{"[", 
       RowBox[{"1", ",", "1"}], "]"}], "]"}], "\[Equal]", "0"}], ",", "p12"}],
    "]"}]}]], "Input",
 CellChangeTimes->{{3.484050059916765*^9, 3.4840501369490013`*^9}, {
  3.4840501812620687`*^9, 3.4840501843871083`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"p12", "\[Rule]", 
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        SqrtBox["q"], " ", 
        SqrtBox["r"]}], "a"]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"p12", "\[Rule]", 
     FractionBox[
      RowBox[{
       SqrtBox["q"], " ", 
       SqrtBox["r"]}], "a"]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.484050103432947*^9, 3.484050137261505*^9, 
  3.4840501848089886`*^9, 3.4840640309168663`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"pDot22", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"pDot", "[", 
     RowBox[{"[", 
      RowBox[{"2", ",", "2"}], "]"}], "]"}], "/.", "p12Sol"}], "//", 
   "Simplify"}]}]], "Input",
 CellChangeTimes->{{3.484050170308803*^9, 3.484050192480962*^9}, {
  3.4840640820119247`*^9, 3.4840640852151313`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"-", "2"}], " ", "b", " ", "p22"}], "-", 
    FractionBox[
     RowBox[{
      SuperscriptBox["a", "2"], " ", 
      SuperscriptBox["p22", "2"]}], "r"], "-", 
    FractionBox[
     RowBox[{"2", " ", 
      SqrtBox["q"], " ", 
      SqrtBox["r"]}], "a"]}], ",", 
   RowBox[{
    RowBox[{
     RowBox[{"-", "2"}], " ", "b", " ", "p22"}], "-", 
    FractionBox[
     RowBox[{
      SuperscriptBox["a", "2"], " ", 
      SuperscriptBox["p22", "2"]}], "r"], "+", 
    FractionBox[
     RowBox[{"2", " ", 
      SqrtBox["q"], " ", 
      SqrtBox["r"]}], "a"]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.484050174152602*^9, 3.4840501927622156`*^9}, 
   3.4840640332606764`*^9, 3.484064087777697*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "we", " ", "have", " ", "two", " ", "possible", " ", "solutions", " ", 
    "for", " ", "p22", " ", "for", " ", "each", " ", "equation"}], " ", 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"firstTwo", "=", 
    RowBox[{"Solve", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"pDot22", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "\[Equal]", "0"}], ",", "p22"}], 
     "]"}]}], " ", "\[IndentingNewLine]", 
   RowBox[{"secondTwo", "=", 
    RowBox[{"Solve", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"pDot22", "[", 
        RowBox[{"[", "2", "]"}], "]"}], "\[Equal]", "0"}], ",", "p22"}], 
     "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.4840641022624426`*^9, 3.484064143075988*^9}, {
  3.4840642761887703`*^9, 3.484064310408227*^9}, {3.4840643718312817`*^9, 
  3.4840643781751533`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"p22", "\[Rule]", 
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"-", "a"}], " ", "b", " ", "r"}], "-", 
       SqrtBox[
        RowBox[{
         RowBox[{
          RowBox[{"-", "2"}], " ", 
          SuperscriptBox["a", "3"], " ", 
          SqrtBox["q"], " ", 
          SuperscriptBox["r", 
           RowBox[{"3", "/", "2"}]]}], "+", 
         RowBox[{
          SuperscriptBox["a", "2"], " ", 
          SuperscriptBox["b", "2"], " ", 
          SuperscriptBox["r", "2"]}]}]]}], 
      SuperscriptBox["a", "3"]]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"p22", "\[Rule]", 
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"-", "a"}], " ", "b", " ", "r"}], "+", 
       SqrtBox[
        RowBox[{
         RowBox[{
          RowBox[{"-", "2"}], " ", 
          SuperscriptBox["a", "3"], " ", 
          SqrtBox["q"], " ", 
          SuperscriptBox["r", 
           RowBox[{"3", "/", "2"}]]}], "+", 
         RowBox[{
          SuperscriptBox["a", "2"], " ", 
          SuperscriptBox["b", "2"], " ", 
          SuperscriptBox["r", "2"]}]}]]}], 
      SuperscriptBox["a", "3"]]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.484064104809383*^9, 3.484064143997886*^9}, {
   3.484064291110982*^9, 3.484064311017614*^9}, 3.48406437851891*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"p22", "\[Rule]", 
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"-", "a"}], " ", "b", " ", "r"}], "-", 
       SqrtBox[
        RowBox[{
         RowBox[{"2", " ", 
          SuperscriptBox["a", "3"], " ", 
          SqrtBox["q"], " ", 
          SuperscriptBox["r", 
           RowBox[{"3", "/", "2"}]]}], "+", 
         RowBox[{
          SuperscriptBox["a", "2"], " ", 
          SuperscriptBox["b", "2"], " ", 
          SuperscriptBox["r", "2"]}]}]]}], 
      SuperscriptBox["a", "3"]]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"p22", "\[Rule]", 
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"-", "a"}], " ", "b", " ", "r"}], "+", 
       SqrtBox[
        RowBox[{
         RowBox[{"2", " ", 
          SuperscriptBox["a", "3"], " ", 
          SqrtBox["q"], " ", 
          SuperscriptBox["r", 
           RowBox[{"3", "/", "2"}]]}], "+", 
         RowBox[{
          SuperscriptBox["a", "2"], " ", 
          SuperscriptBox["b", "2"], " ", 
          SuperscriptBox["r", "2"]}]}]]}], 
      SuperscriptBox["a", "3"]]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.484064104809383*^9, 3.484064143997886*^9}, {
   3.484064291110982*^9, 3.484064311017614*^9}, 3.48406437851891*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"allSol", "=", 
  RowBox[{"Join", "[", 
   RowBox[{"firstTwo", ",", "secondTwo"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4840643806752014`*^9, 3.4840644360981407`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"p22", "\[Rule]", 
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"-", "a"}], " ", "b", " ", "r"}], "-", 
       SqrtBox[
        RowBox[{
         RowBox[{
          RowBox[{"-", "2"}], " ", 
          SuperscriptBox["a", "3"], " ", 
          SqrtBox["q"], " ", 
          SuperscriptBox["r", 
           RowBox[{"3", "/", "2"}]]}], "+", 
         RowBox[{
          SuperscriptBox["a", "2"], " ", 
          SuperscriptBox["b", "2"], " ", 
          SuperscriptBox["r", "2"]}]}]]}], 
      SuperscriptBox["a", "3"]]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"p22", "\[Rule]", 
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"-", "a"}], " ", "b", " ", "r"}], "+", 
       SqrtBox[
        RowBox[{
         RowBox[{
          RowBox[{"-", "2"}], " ", 
          SuperscriptBox["a", "3"], " ", 
          SqrtBox["q"], " ", 
          SuperscriptBox["r", 
           RowBox[{"3", "/", "2"}]]}], "+", 
         RowBox[{
          SuperscriptBox["a", "2"], " ", 
          SuperscriptBox["b", "2"], " ", 
          SuperscriptBox["r", "2"]}]}]]}], 
      SuperscriptBox["a", "3"]]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"p22", "\[Rule]", 
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"-", "a"}], " ", "b", " ", "r"}], "-", 
       SqrtBox[
        RowBox[{
         RowBox[{"2", " ", 
          SuperscriptBox["a", "3"], " ", 
          SqrtBox["q"], " ", 
          SuperscriptBox["r", 
           RowBox[{"3", "/", "2"}]]}], "+", 
         RowBox[{
          SuperscriptBox["a", "2"], " ", 
          SuperscriptBox["b", "2"], " ", 
          SuperscriptBox["r", "2"]}]}]]}], 
      SuperscriptBox["a", "3"]]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"p22", "\[Rule]", 
     FractionBox[
      RowBox[{
       RowBox[{
        RowBox[{"-", "a"}], " ", "b", " ", "r"}], "+", 
       SqrtBox[
        RowBox[{
         RowBox[{"2", " ", 
          SuperscriptBox["a", "3"], " ", 
          SqrtBox["q"], " ", 
          SuperscriptBox["r", 
           RowBox[{"3", "/", "2"}]]}], "+", 
         RowBox[{
          SuperscriptBox["a", "2"], " ", 
          SuperscriptBox["b", "2"], " ", 
          SuperscriptBox["r", "2"]}]}]]}], 
      SuperscriptBox["a", "3"]]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.4840643951598544`*^9, 3.484064436707527*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
  "Now", " ", "p22", " ", "must", " ", "be", " ", "postive", " ", "and", " ", 
   "real", " ", "so", " ", "only", " ", "a", " ", "few", " ", "of", " ", 
   "these", " ", "are", " ", "valid", " ", "solutions"}], " ", 
  "*)"}]], "Input",
 CellChangeTimes->{{3.4840505067449684`*^9, 3.484050538245372*^9}}]
},
WindowSize->{809, 750},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
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, 964, 33, 92, "Input"],
Cell[CellGroupData[{
Cell[1534, 57, 582, 15, 68, "Input"],
Cell[2119, 74, 1489, 42, 84, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3645, 121, 358, 10, 31, "Input"],
Cell[4006, 133, 497, 17, 50, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4540, 155, 324, 9, 31, "Input"],
Cell[4867, 166, 771, 26, 50, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5675, 197, 875, 24, 72, "Input"],
Cell[6553, 223, 1354, 42, 58, "Output"],
Cell[7910, 267, 1302, 40, 58, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9249, 312, 194, 4, 31, "Input"],
Cell[9446, 318, 2397, 77, 108, "Output"]
}, Open  ]],
Cell[11858, 398, 350, 7, 31, "Input"]
}
]
*)

(* End of internal cache information *)