(* 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[      5073,        173]
NotebookOptionsPosition[      4614,        152]
NotebookOutlinePosition[      4978,        168]
CellTagsIndexPosition[      4935,        165]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{"Solve", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"0", "\[Equal]", 
      RowBox[{
       RowBox[{"2", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "bf"}], "-", "k"}], ")"}], " ", "pI"}], "+", 
       RowBox[{"2", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "bf"}], "+", "b"}], ")"}], " ", "vI"}], "+", "q", "+", 
       RowBox[{
        SuperscriptBox["k", "2"], "r"}]}]}], ",", 
     RowBox[{"0", "\[Equal]", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "b"}], " ", "vI"}], " ", "+", " ", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "bf"}], "-", "k"}], ")"}], " ", "vI"}], " ", "+", " ", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "bf"}], "+", "b"}], ")"}], "uI"}], "-", "q"}]}], ",", 
     RowBox[{"0", "\[Equal]", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "2"}], " ", "b", " ", "uI"}], "+", "q"}]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"uI", ",", "vI", ",", "pI"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.491133665056736*^9, 3.49113373063612*^9}, 
   3.49113584106889*^9, 3.491135977836525*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"pI", "\[Rule]", 
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        RowBox[{
         RowBox[{"-", "b"}], " ", "bf", " ", "q"}], "-", 
        RowBox[{
         SuperscriptBox["bf", "2"], " ", "q"}], "-", 
        RowBox[{"b", " ", "k", " ", "q"}], "-", 
        RowBox[{
         SuperscriptBox["b", "2"], " ", 
         SuperscriptBox["k", "2"], " ", "r"}], "-", 
        RowBox[{"b", " ", "bf", " ", 
         SuperscriptBox["k", "2"], " ", "r"}], "-", 
        RowBox[{"b", " ", 
         SuperscriptBox["k", "3"], " ", "r"}]}], 
       RowBox[{"2", " ", "b", " ", 
        RowBox[{"(", 
         RowBox[{"bf", "+", "k"}], ")"}], " ", 
        RowBox[{"(", 
         RowBox[{"b", "+", "bf", "+", "k"}], ")"}]}]]}]}], ",", 
    RowBox[{"vI", "\[Rule]", 
     RowBox[{"-", 
      FractionBox[
       RowBox[{
        RowBox[{"b", " ", "q"}], "+", 
        RowBox[{"bf", " ", "q"}]}], 
       RowBox[{"2", " ", "b", " ", 
        RowBox[{"(", 
         RowBox[{"b", "+", "bf", "+", "k"}], ")"}]}]]}]}], ",", 
    RowBox[{"uI", "\[Rule]", 
     FractionBox["q", 
      RowBox[{"2", " ", "b"}]]}]}], "}"}], "}"}]], "Output",
 CellChangeTimes->{3.491133731354884*^9, 3.4911358480846496`*^9, 
  3.4911359788052874`*^9, 3.4911361043537693`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Simplify", "[", 
  RowBox[{"-", 
   FractionBox[
    RowBox[{
     RowBox[{
      RowBox[{"-", "b"}], " ", "bf", " ", "q"}], "-", 
     RowBox[{
      SuperscriptBox["bf", "2"], " ", "q"}], "-", 
     RowBox[{"b", " ", "k", " ", "q"}], "-", 
     RowBox[{
      SuperscriptBox["b", "2"], " ", 
      SuperscriptBox["k", "2"], " ", "r"}], "-", 
     RowBox[{"b", " ", "bf", " ", 
      SuperscriptBox["k", "2"], " ", "r"}], "-", 
     RowBox[{"b", " ", 
      SuperscriptBox["k", "3"], " ", "r"}]}], 
    RowBox[{"2", " ", "b", " ", 
     RowBox[{"(", 
      RowBox[{"bf", "+", "k"}], ")"}], " ", 
     RowBox[{"(", 
      RowBox[{"b", "+", "bf", "+", "k"}], ")"}]}]]}], "]"}]], "Input",
 CellChangeTimes->{{3.491133778590166*^9, 3.491133791684167*^9}, 
   3.491135857709834*^9, 3.491136111197607*^9}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{
    SuperscriptBox["bf", "2"], " ", "q"}], "+", 
   RowBox[{
    SuperscriptBox["b", "2"], " ", 
    SuperscriptBox["k", "2"], " ", "r"}], "+", 
   RowBox[{"b", " ", 
    RowBox[{"(", 
     RowBox[{"bf", "+", "k"}], ")"}], " ", 
    RowBox[{"(", 
     RowBox[{"q", "+", 
      RowBox[{
       SuperscriptBox["k", "2"], " ", "r"}]}], ")"}]}]}], 
  RowBox[{"2", " ", "b", " ", 
   RowBox[{"(", 
    RowBox[{"bf", "+", "k"}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{"b", "+", "bf", "+", "k"}], ")"}]}]]], "Output",
 CellChangeTimes->{3.491133792762313*^9, 3.4911358591942377`*^9, 
  3.491136111775739*^9}]
}, Open  ]]
},
WindowSize->{616, 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[CellGroupData[{
Cell[567, 22, 1188, 36, 77, "Input"],
Cell[1758, 60, 1321, 38, 91, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3116, 103, 824, 23, 52, "Input"],
Cell[3943, 128, 655, 21, 53, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)