(* 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[      5124,        158]
NotebookOptionsPosition[      4776,        141]
NotebookOutlinePosition[      5140,        157]
CellTagsIndexPosition[      5097,        154]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{"thePSolution", "=", 
  RowBox[{"DSolve", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"p11", "'"}], "[", "t", "]"}], "\[Equal]", 
       RowBox[{"2", " ", 
        RowBox[{"p12", "[", "t", "]"}]}]}], ",", 
      RowBox[{
       RowBox[{
        RowBox[{"p12", "'"}], "[", "t", "]"}], "\[Equal]", 
       RowBox[{
        RowBox[{"p22", "[", "t", "]"}], "-", 
        RowBox[{"b", " ", 
         RowBox[{"p12", "[", "t", "]"}]}]}]}], ",", 
      RowBox[{
       RowBox[{
        RowBox[{"p22", "'"}], "[", "t", "]"}], "\[Equal]", 
       RowBox[{
        RowBox[{
         RowBox[{"-", "2"}], " ", "b", " ", 
         RowBox[{"p22", "[", "t", "]"}]}], "+", 
        SuperscriptBox["s", "2"]}]}], ",", 
      RowBox[{
       RowBox[{"p11", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"p12", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"p22", "[", "0", "]"}], "\[Equal]", 
       SuperscriptBox["s", "2"]}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"p11", ",", "p12", ",", "p22"}], "}"}], ",", 
    RowBox[{"{", "t", "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.4807839068209066`*^9, 3.4807839830411205`*^9}, {
  3.480784042823518*^9, 3.480784058526944*^9}, {3.480784159403881*^9, 
  3.480784160325774*^9}, {3.480784193342033*^9, 3.4807841973264847`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"p11", "\[Rule]", 
     RowBox[{"Function", "[", 
      RowBox[{
       RowBox[{"{", "t", "}"}], ",", 
       FractionBox[
        RowBox[{
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "2"}], " ", "b", " ", "t"}]], " ", 
         SuperscriptBox["s", "2"], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           RowBox[{"2", " ", "b"}], "+", 
           RowBox[{"4", " ", 
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{"b", " ", "t"}]]}], "-", 
           RowBox[{"4", " ", "b", " ", 
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{"b", " ", "t"}]]}], "-", 
           RowBox[{"3", " ", 
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{"2", " ", "b", " ", "t"}]]}], "+", 
           RowBox[{"2", " ", "b", " ", 
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{"2", " ", "b", " ", "t"}]]}], "+", 
           RowBox[{"2", " ", "b", " ", 
            SuperscriptBox["\[ExponentialE]", 
             RowBox[{"2", " ", "b", " ", "t"}]], " ", "t"}]}], ")"}]}], 
        RowBox[{"2", " ", 
         SuperscriptBox["b", "3"]}]]}], "]"}]}], ",", 
    RowBox[{"p12", "\[Rule]", 
     RowBox[{"Function", "[", 
      RowBox[{
       RowBox[{"{", "t", "}"}], ",", 
       FractionBox[
        RowBox[{
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "2"}], " ", "b", " ", "t"}]], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           SuperscriptBox["\[ExponentialE]", 
            RowBox[{"b", " ", "t"}]]}], ")"}], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           RowBox[{"2", " ", "b"}], "+", 
           SuperscriptBox["\[ExponentialE]", 
            RowBox[{"b", " ", "t"}]]}], ")"}], " ", 
         SuperscriptBox["s", "2"]}], 
        RowBox[{"2", " ", 
         SuperscriptBox["b", "2"]}]]}], "]"}]}], ",", 
    RowBox[{"p22", "\[Rule]", 
     RowBox[{"Function", "[", 
      RowBox[{
       RowBox[{"{", "t", "}"}], ",", 
       FractionBox[
        RowBox[{
         SuperscriptBox["\[ExponentialE]", 
          RowBox[{
           RowBox[{"-", "2"}], " ", "b", " ", "t"}]], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1"}], "+", 
           RowBox[{"2", " ", "b"}], "+", 
           SuperscriptBox["\[ExponentialE]", 
            RowBox[{"2", " ", "b", " ", "t"}]]}], ")"}], " ", 
         SuperscriptBox["s", "2"]}], 
        RowBox[{"2", " ", "b"}]]}], "]"}]}]}], "}"}], "}"}]], "Output",
 CellChangeTimes->{3.4807840591832075`*^9, 3.480784161200791*^9, 
  3.4807841978108687`*^9, 3.4808427032146235`*^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, 1394, 38, 98, "Input"],
Cell[1964, 62, 2796, 76, 167, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)