(* 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[      4900,        169]
NotebookOptionsPosition[      4425,        148]
NotebookOutlinePosition[      4768,        163]
CellTagsIndexPosition[      4725,        160]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"Pminus", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"v", "-", 
        FractionBox[
         SuperscriptBox["v", "2"], 
         RowBox[{"v", "+", "v1"}]]}], ",", "0", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", 
       RowBox[{"v", "-", 
        FractionBox[
         SuperscriptBox["v", "2"], 
         RowBox[{"v", "+", "v2"}]]}], ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "v"}], "}"}]}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"H", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"R", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"v1", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"0", ",", "v3"}], "}"}]}], "}"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.4846432339114623`*^9, 3.4846432704898214`*^9}, 
   3.4846433240526648`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "simplify", " ", "some", " ", "needed", " ", "matrix", " ", 
    "multiplications"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"pt1", "=", 
   RowBox[{"Simplify", "[", 
    RowBox[{"Inverse", "[", 
     RowBox[{
      RowBox[{"H", " ", ".", " ", "Pminus", " ", ".", " ", 
       RowBox[{"Transpose", "[", "H", "]"}]}], "+", "R"}], "]"}], 
    "]"}]}]}]], "Input",
 CellChangeTimes->{{3.484643966208399*^9, 3.484644042802639*^9}, {
  3.4846442059542923`*^9, 3.4846442366888638`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     FractionBox[
      RowBox[{"v", "+", "v1"}], 
      RowBox[{
       RowBox[{"2", " ", "v", " ", "v1"}], "+", 
       SuperscriptBox["v1", "2"]}]], ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0", ",", 
     FractionBox["1", 
      RowBox[{"v", "+", "v3"}]]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.4846439709271793`*^9, 3.484644043193267*^9}, 
   3.4846442375794945`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "compute", " ", "the", " ", "full", " ", "a", " ", "posteriori", " ", 
    "covariance"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"Pplus", " ", "=", 
     RowBox[{"Simplify", "[", 
      RowBox[{"Pminus", "-", 
       RowBox[{"Pminus", " ", ".", " ", 
        RowBox[{"Transpose", "[", "H", "]"}], ".", "pt1", " ", ".", " ", "H", 
        " ", ".", " ", "Pminus"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"MatrixForm", "[", "Pplus", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.48464404920893*^9, 3.4846441181881056`*^9}, {
  3.4846442392513804`*^9, 3.4846442513608327`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {
      FractionBox[
       RowBox[{"v", " ", "v1"}], 
       RowBox[{
        RowBox[{"2", " ", "v"}], "+", "v1"}]], "0", "0"},
     {"0", 
      FractionBox[
       RowBox[{"v", " ", "v2"}], 
       RowBox[{"v", "+", "v2"}]], "0"},
     {"0", "0", 
      FractionBox[
       RowBox[{"v", " ", "v3"}], 
       RowBox[{"v", "+", "v3"}]]}
    },
    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.4846440582089877`*^9, 3.4846441185787334`*^9}, 
   3.484644252032712*^9}]
}, Open  ]]
},
WindowSize->{616, 750},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
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, 1108, 37, 91, "Input"],
Cell[CellGroupData[{
Cell[1678, 61, 549, 14, 52, "Input"],
Cell[2230, 77, 467, 15, 48, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2734, 97, 681, 16, 72, "Input"],
Cell[3418, 115, 991, 30, 104, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)