(* 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[     24923,        843]
NotebookOptionsPosition[     22101,        746]
NotebookOutlinePosition[     22541,        763]
CellTagsIndexPosition[     22498,        760]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{"SigmaI", " ", "=", " ", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"5", ",", 
      RowBox[{"-", "1"}], ",", "2"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "1"}], ",", "8", ",", "1"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"2", ",", "1", ",", "4"}], "}"}]}], "}"}], " "}]], "Input",
 CellChangeTimes->{{3.5707317855288568`*^9, 3.570731932293915*^9}, {
   3.57073215737211*^9, 3.5707321609289165`*^9}, {3.570732708068677*^9, 
   3.570732717584694*^9}, 3.5707328146480646`*^9, 3.5707332371592064`*^9, 
   3.5707340016863494`*^9, 3.5707366662038264`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"5", ",", 
     RowBox[{"-", "1"}], ",", "2"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1"}], ",", "8", ",", "1"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"2", ",", "1", ",", "4"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.57073184995697*^9, 3.5707319326371155`*^9}, 
   3.570732161443717*^9, {3.570732708567878*^9, 3.5707327185830956`*^9}, 
   3.570732815022465*^9, 3.570733237393207*^9, 3.5707340020919504`*^9, 
   3.570736402266963*^9, 3.570736707606299*^9}]
}, Open  ]],

Cell["Check that SigmaI is pos def:", "Text",
 CellChangeTimes->{{3.570734571742551*^9, 3.5707345801353655`*^9}, 
   3.570736666219426*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Eigenvalues", "[", 
  RowBox[{"N", "[", "SigmaI", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.5707318897214403`*^9, 3.5707318947134485`*^9}, {
   3.5707327223115025`*^9, 3.570732723668705*^9}, {3.5707328248504825`*^9, 
   3.5707328278768873`*^9}, 3.570736666219426*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "8.349966717646584`", ",", "6.546110732079136`", ",", 
   "2.103922550274274`"}], "}"}]], "Output",
 CellChangeTimes->{{3.5707318954310503`*^9, 3.5707319213738956`*^9}, 
   3.5707327241211057`*^9, {3.570732818860072*^9, 3.5707328281888885`*^9}, 
   3.5707364022825627`*^9, 3.5707367076218987`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"xvec", " ", "=", " ", 
  RowBox[{"{", 
   RowBox[{"x1", ",", "x2", ",", "x3"}], "}"}]}]], "Input",
 CellChangeTimes->{{3.570732729066314*^9, 3.5707327358679256`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"x1", ",", "x2", ",", "x3"}], "}"}]], "Output",
 CellChangeTimes->{3.570732736975528*^9, 3.5707364022825627`*^9, 
  3.570736707637499*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"mu", " ", "=", " ", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"-", "1"}], ",", 
    RowBox[{"-", "1"}], ",", 
    RowBox[{"-", "1"}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.5707327868176155`*^9, 3.570732810311257*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "1"}], ",", 
   RowBox[{"-", "1"}], ",", 
   RowBox[{"-", "1"}]}], "}"}]], "Output",
 CellChangeTimes->{3.570732810716858*^9, 3.5707364022981625`*^9, 
  3.570736707637499*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"expon", " ", "=", " ", 
  RowBox[{"FullSimplify", "[", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"xvec", " ", "-", " ", "mu"}], ")"}], ".", "SigmaI", ".", 
    RowBox[{"(", 
     RowBox[{"xvec", "-", "mu"}], ")"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5707327717323895`*^9, 3.570732782465208*^9}, {
   3.570732833196497*^9, 3.570732843492515*^9}, {3.570733763801532*^9, 
   3.570733793113983*^9}, 3.570736666219426*^9}],

Cell[BoxData[
 RowBox[{"21", "+", 
  RowBox[{"5", " ", 
   SuperscriptBox["x1", "2"]}], "+", 
  RowBox[{"8", " ", 
   SuperscriptBox["x2", "2"]}], "+", 
  RowBox[{"2", " ", "x2", " ", 
   RowBox[{"(", 
    RowBox[{"8", "+", "x3"}], ")"}]}], "+", 
  RowBox[{"2", " ", "x3", " ", 
   RowBox[{"(", 
    RowBox[{"7", "+", 
     RowBox[{"2", " ", "x3"}]}], ")"}]}], "-", 
  RowBox[{"2", " ", "x1", " ", 
   RowBox[{"(", 
    RowBox[{"x2", "-", 
     RowBox[{"2", " ", 
      RowBox[{"(", 
       RowBox[{"3", "+", "x3"}], ")"}]}]}], ")"}]}]}]], "Output",
 CellChangeTimes->{
  3.5707328448809175`*^9, 3.570733242244816*^9, {3.5707337692615414`*^9, 
   3.570733794049985*^9}, 3.570734008363161*^9, 3.5707364023761625`*^9, 
   3.570736707653099*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"exponA", " ", "=", " ", 
  RowBox[{"expon", " ", "/.", " ", 
   RowBox[{"x3", "\[Rule]", "0"}]}]}]], "Input",
 CellChangeTimes->{{3.5707329291210656`*^9, 3.5707329423342886`*^9}}],

Cell[BoxData[
 RowBox[{"21", "+", 
  RowBox[{"5", " ", 
   SuperscriptBox["x1", "2"]}], "-", 
  RowBox[{"2", " ", "x1", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "6"}], "+", "x2"}], ")"}]}], "+", 
  RowBox[{"16", " ", "x2"}], "+", 
  RowBox[{"8", " ", 
   SuperscriptBox["x2", "2"]}]}]], "Output",
 CellChangeTimes->{3.5707329434106903`*^9, 3.5707332455832214`*^9, 
  3.5707338085112104`*^9, 3.570734012746769*^9, 3.5707364023761625`*^9, 
  3.570736707653099*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"muAsub", " ", "=", " ", 
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"D", "[", 
        RowBox[{"exponA", ",", "x1"}], "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"D", "[", 
        RowBox[{"exponA", ",", "x2"}], "]"}], "\[Equal]", "0"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"x1", ",", "x2"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5707330258568354`*^9, 3.5707330906593494`*^9}, {
  3.570733896854165*^9, 3.57073389933457*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"x1", "\[Rule]", 
     RowBox[{"-", 
      FractionBox["56", "39"]}]}], ",", 
    RowBox[{"x2", "\[Rule]", 
     RowBox[{"-", 
      FractionBox["46", "39"]}]}]}], "}"}], "}"}]], "Output",
 CellChangeTimes->{3.570733043313266*^9, 3.5707330919697514`*^9, 
  3.570733249654828*^9, 3.5707338150008216`*^9, 3.5707338999429708`*^9, 
  3.5707340209679832`*^9, 3.5707364023917627`*^9, 3.5707367076686993`*^9}]
}, Open  ]],

Cell["So this is the mean on the slice.", "Text",
 CellChangeTimes->{{3.5707346570903006`*^9, 3.570734677619937*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"subsA", " ", "=", " ", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"x1", "\[Rule]", 
     RowBox[{"y1", " ", "+", " ", 
      RowBox[{"(", 
       RowBox[{"x1", " ", "/.", " ", "muAsub"}], ")"}]}]}], ",", " ", 
    RowBox[{"x2", "\[Rule]", 
     RowBox[{"y2", "+", 
      RowBox[{"(", 
       RowBox[{"x2", " ", "/.", " ", "muAsub"}], ")"}]}]}]}], 
   "}"}]}]], "Input",
 CellChangeTimes->{{3.570733861379703*^9, 3.5707339474918547`*^9}, 
   3.570734641537073*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"x1", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", 
       FractionBox["56", "39"]}], "+", "y1"}], "}"}]}], ",", 
   RowBox[{"x2", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", 
       FractionBox["46", "39"]}], "+", "y2"}], "}"}]}]}], "}"}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{3.5707339492390575`*^9, 3.570733983902318*^9, 
  3.5707340296259985`*^9, 3.570734641537073*^9, 3.5707364023917627`*^9, 
  3.5707367076686993`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Simplify", "[", 
  RowBox[{"exponA", "  ", "/.", " ", "subsA"}], "]"}]], "Input",
 CellChangeTimes->{{3.5707336045252523`*^9, 3.570733627051691*^9}, {
   3.5707339782239084`*^9, 3.5707339794563103`*^9}, 3.570734641537073*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   FractionBox["115", "39"], "+", 
   RowBox[{"5", " ", 
    SuperscriptBox["y1", "2"]}], "-", 
   RowBox[{"2", " ", "y1", " ", "y2"}], "+", 
   RowBox[{"8", " ", 
    SuperscriptBox["y2", "2"]}]}], "}"}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{{3.5707336179724755`*^9, 3.570733627379292*^9}, {
   3.570733973044699*^9, 3.5707339867415233`*^9}, 3.5707340341968064`*^9, 
   3.5707346415526733`*^9, 3.570736402407363*^9, 3.5707367076842995`*^9}]
}, Open  ]],

Cell["\<\
So SigmaI on the slice is just the restriction of the 3-dim SigmaI: \
{{5,-1},{-1,8}}\
\>", "Text",
 CellChangeTimes->{{3.570734687369954*^9, 3.570734734263636*^9}, {
  3.570736666219426*^9, 3.5707366662350264`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"integral", " ", "=", " ", 
  RowBox[{"Integrate", "[", 
   RowBox[{
    RowBox[{"Exp", "[", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "1"}], "/", "2"}], ")"}], 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"xvec", " ", "-", " ", "mu"}], ")"}], ".", "SigmaI", ".", 
        RowBox[{"(", 
         RowBox[{"xvec", "-", "mu"}], ")"}]}], ")"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"x3", ",", 
      RowBox[{"-", "Infinity"}], ",", "Infinity"}], "}"}], ",", 
    RowBox[{"GenerateConditions", "\[Rule]", "False"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.570734123241763*^9, 3.570734190883482*^9}, {
   3.5707342669180155`*^9, 3.570734268961619*^9}, 3.5707346415526733`*^9, 
   3.5707366662350264`*^9}],

Cell[BoxData[
 RowBox[{
  SuperscriptBox["\[ExponentialE]", 
   RowBox[{
    FractionBox["1", "8"], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "35"}], "-", 
      RowBox[{"16", " ", 
       SuperscriptBox["x1", "2"]}], "-", 
      RowBox[{"50", " ", "x2"}], "-", 
      RowBox[{"31", " ", 
       SuperscriptBox["x2", "2"]}], "+", 
      RowBox[{"4", " ", "x1", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "5"}], "+", 
         RowBox[{"3", " ", "x2"}]}], ")"}]}]}], ")"}]}]], " ", 
  SqrtBox[
   FractionBox["\[Pi]", "2"]]}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{3.57073419568829*^9, 3.5707342733452263`*^9, 
  3.5707346415526733`*^9, 3.5707364066037703`*^9, 3.5707367117715063`*^9}],

Cell[BoxData[
 RowBox[{
  SuperscriptBox["\[ExponentialE]", 
   RowBox[{
    FractionBox["1", "8"], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "35"}], "-", 
      RowBox[{"16", " ", 
       SuperscriptBox["x1", "2"]}], "-", 
      RowBox[{"50", " ", "x2"}], "-", 
      RowBox[{"31", " ", 
       SuperscriptBox["x2", "2"]}], "+", 
      RowBox[{"4", " ", "x1", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "5"}], "+", 
         RowBox[{"3", " ", "x2"}]}], ")"}]}]}], ")"}]}]], " ", 
  SqrtBox[
   FractionBox["\[Pi]", "2"]]}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{3.57073419568829*^9, 3.5707342733452263`*^9, 
  3.5707346415526733`*^9, 3.5707364066037703`*^9}],

Cell[BoxData[
 RowBox[{
  SuperscriptBox["\[ExponentialE]", 
   RowBox[{
    FractionBox["1", "8"], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "35"}], "-", 
      RowBox[{"16", " ", 
       SuperscriptBox["x1", "2"]}], "-", 
      RowBox[{"50", " ", "x2"}], "-", 
      RowBox[{"31", " ", 
       SuperscriptBox["x2", "2"]}], "+", 
      RowBox[{"4", " ", "x1", " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "5"}], "+", 
         RowBox[{"3", " ", "x2"}]}], ")"}]}]}], ")"}]}]], " ", 
  SqrtBox[
   FractionBox["\[Pi]", "2"]]}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{3.57073419568829*^9, 3.5707342733452263`*^9, 
  3.5707346415526733`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"exponB", " ", "=", " ", 
  RowBox[{
   RowBox[{"-", "2"}], " ", 
   RowBox[{
    RowBox[{"integral", "[", 
     RowBox[{"[", "1", "]"}], "]"}], "[", 
    RowBox[{"[", "2", "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.570734276028431*^9, 3.5707343262761197`*^9}, 
   3.570734518530857*^9, 3.5707346415526733`*^9}],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "4"], " ", 
  RowBox[{"(", 
   RowBox[{"35", "+", 
    RowBox[{"16", " ", 
     SuperscriptBox["x1", "2"]}], "+", 
    RowBox[{"50", " ", "x2"}], "+", 
    RowBox[{"31", " ", 
     SuperscriptBox["x2", "2"]}], "-", 
    RowBox[{"4", " ", "x1", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "5"}], "+", 
       RowBox[{"3", " ", "x2"}]}], ")"}]}]}], ")"}]}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{{3.5707342784152355`*^9, 3.5707342925644608`*^9}, 
   3.570734328678524*^9, 3.5707345198568597`*^9, 3.5707346415526733`*^9, 
   3.5707364066193705`*^9, 3.5707367117871065`*^9}],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "4"], " ", 
  RowBox[{"(", 
   RowBox[{"35", "+", 
    RowBox[{"16", " ", 
     SuperscriptBox["x1", "2"]}], "+", 
    RowBox[{"50", " ", "x2"}], "+", 
    RowBox[{"31", " ", 
     SuperscriptBox["x2", "2"]}], "-", 
    RowBox[{"4", " ", "x1", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "5"}], "+", 
       RowBox[{"3", " ", "x2"}]}], ")"}]}]}], ")"}]}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{{3.5707342784152355`*^9, 3.5707342925644608`*^9}, 
   3.570734328678524*^9, 3.5707345198568597`*^9, 3.5707346415526733`*^9, 
   3.5707364066193705`*^9}],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "4"], " ", 
  RowBox[{"(", 
   RowBox[{"35", "+", 
    RowBox[{"16", " ", 
     SuperscriptBox["x1", "2"]}], "+", 
    RowBox[{"50", " ", "x2"}], "+", 
    RowBox[{"31", " ", 
     SuperscriptBox["x2", "2"]}], "-", 
    RowBox[{"4", " ", "x1", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "5"}], "+", 
       RowBox[{"3", " ", "x2"}]}], ")"}]}]}], ")"}]}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{{3.5707342784152355`*^9, 3.5707342925644608`*^9}, 
   3.570734328678524*^9, 3.5707345198568597`*^9, 3.5707346415526733`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"muBsub", " ", "=", " ", 
  RowBox[{"Solve", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"D", "[", 
        RowBox[{"exponB", ",", "x1"}], "]"}], "\[Equal]", "0"}], ",", 
      RowBox[{
       RowBox[{"D", "[", 
        RowBox[{"exponB", ",", "x2"}], "]"}], "\[Equal]", "0"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"x1", ",", "x2"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.570734368536594*^9, 3.570734375010605*^9}, 
   3.5707346415526733`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"x1", "\[Rule]", 
     RowBox[{"-", "1"}]}], ",", 
    RowBox[{"x2", "\[Rule]", 
     RowBox[{"-", "1"}]}]}], "}"}], "}"}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{3.5707343780526104`*^9, 3.5707345311980796`*^9, 
  3.5707346415526733`*^9, 3.5707364066349707`*^9, 3.5707367118027067`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"x1", "\[Rule]", 
     RowBox[{"-", "1"}]}], ",", 
    RowBox[{"x2", "\[Rule]", 
     RowBox[{"-", "1"}]}]}], "}"}], "}"}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{3.5707343780526104`*^9, 3.5707345311980796`*^9, 
  3.5707346415526733`*^9, 3.5707364066349707`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"x1", "\[Rule]", 
     RowBox[{"-", "1"}]}], ",", 
    RowBox[{"x2", "\[Rule]", 
     RowBox[{"-", "1"}]}]}], "}"}], "}"}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{3.5707343780526104`*^9, 3.5707345311980796`*^9, 
  3.5707346415526733`*^9}]
}, Open  ]],

Cell["\<\
So the projected mean is of course the restriction of the original mean.\
\>", "Text",
 CellChangeTimes->{{3.570734767243094*^9, 3.57073479348234*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"subsB", " ", "=", " ", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"x1", "\[Rule]", 
     RowBox[{"y1", " ", "+", " ", 
      RowBox[{"(", 
       RowBox[{"x1", " ", "/.", " ", "muBsub"}], ")"}]}]}], ",", " ", 
    RowBox[{"x2", "\[Rule]", 
     RowBox[{"y2", "+", 
      RowBox[{"(", 
       RowBox[{"x2", " ", "/.", " ", "muBsub"}], ")"}]}]}]}], 
   "}"}]}]], "Input",
 CellChangeTimes->{{3.570734419470683*^9, 3.5707344293923006`*^9}, 
   3.5707346415526733`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"x1", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "y1"}], "}"}]}], ",", 
   RowBox[{"x2", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "y2"}], "}"}]}]}], "}"}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{3.5707344298915014`*^9, 3.570734535519287*^9, 
  3.5707346415526733`*^9, 3.5707364066505704`*^9, 3.5707367118183064`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"x1", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "y1"}], "}"}]}], ",", 
   RowBox[{"x2", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "y2"}], "}"}]}]}], "}"}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{3.5707344298915014`*^9, 3.570734535519287*^9, 
  3.5707346415526733`*^9, 3.5707364066505704`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"x1", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "y1"}], "}"}]}], ",", 
   RowBox[{"x2", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "y2"}], "}"}]}]}], "}"}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{3.5707344298915014`*^9, 3.570734535519287*^9, 
  3.5707346415526733`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Simplify", "[", 
  RowBox[{"exponB", "  ", "/.", " ", "subsB"}], "]"}]], "Input",
 CellChangeTimes->{{3.570734450155937*^9, 3.570734456863949*^9}, 
   3.5707346415526733`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"4", " ", 
    SuperscriptBox["y1", "2"]}], "-", 
   RowBox[{"3", " ", "y1", " ", "y2"}], "+", 
   FractionBox[
    RowBox[{"31", " ", 
     SuperscriptBox["y2", "2"]}], "4"]}], "}"}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{3.5707344572539496`*^9, 3.5707345396532946`*^9, 
  3.5707346415526733`*^9, 3.5707364066505704`*^9, 3.570736711833906*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"4", " ", 
    SuperscriptBox["y1", "2"]}], "-", 
   RowBox[{"3", " ", "y1", " ", "y2"}], "+", 
   FractionBox[
    RowBox[{"31", " ", 
     SuperscriptBox["y2", "2"]}], "4"]}], "}"}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{3.5707344572539496`*^9, 3.5707345396532946`*^9, 
  3.5707346415526733`*^9, 3.5707364066505704`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"4", " ", 
    SuperscriptBox["y1", "2"]}], "-", 
   RowBox[{"3", " ", "y1", " ", "y2"}], "+", 
   FractionBox[
    RowBox[{"31", " ", 
     SuperscriptBox["y2", "2"]}], "4"]}], "}"}]], "Output",
 GeneratedCell->False,
 CellAutoOverwrite->False,
 CellChangeTimes->{3.5707344572539496`*^9, 3.5707345396532946`*^9, 
  3.5707346415526733`*^9}]
}, Open  ]],

Cell["\<\
But the projected SigmaI is different: {{4,-3/2},{-3/2,31/4}}
\
\>", "Text",
 CellChangeTimes->{{3.5707348030607567`*^9, 3.570734841951625*^9}, 
   3.5707366662350264`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Sigma", " ", "=", " ", 
  RowBox[{"Inverse", "[", "SigmaI", "]"}]}]], "Input",
 CellChangeTimes->{{3.5707367477295694`*^9, 3.5707367790232244`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     FractionBox["31", "115"], ",", 
     FractionBox["6", "115"], ",", 
     RowBox[{"-", 
      FractionBox["17", "115"]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["6", "115"], ",", 
     FractionBox["16", "115"], ",", 
     RowBox[{"-", 
      FractionBox["7", "115"]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", 
      FractionBox["17", "115"]}], ",", 
     RowBox[{"-", 
      FractionBox["7", "115"]}], ",", 
     FractionBox["39", "115"]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5707367794912252`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"redSigma", " ", "=", " ", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"Sigma", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}], ",", 
      RowBox[{
       RowBox[{"Sigma", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "[", 
       RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{
       RowBox[{"Sigma", "[", 
        RowBox[{"[", "2", "]"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}], ",", 
      RowBox[{
       RowBox[{"Sigma", "[", 
        RowBox[{"[", "2", "]"}], "]"}], "[", 
       RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.570736801580864*^9, 3.570736848739747*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     FractionBox["31", "115"], ",", 
     FractionBox["6", "115"]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["6", "115"], ",", 
     FractionBox["16", "115"]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.570736849160948*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"redSigmaI", " ", "=", " ", 
  RowBox[{"Inverse", "[", "redSigma", "]"}]}]], "Input",
 CellChangeTimes->{{3.570736854464957*^9, 3.5707368627641716`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"4", ",", 
     RowBox[{"-", 
      FractionBox["3", "2"]}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", 
      FractionBox["3", "2"]}], ",", 
     FractionBox["31", "4"]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.570736863138572*^9}]
}, Open  ]]
},
WindowSize->{807, 848},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
Magnification:>FEPrivate`If[
  FEPrivate`Equal[FEPrivate`$VersionNumber, 6.], 2., 2. Inherited],
FrontEndVersion->"7.0 for Microsoft Windows (64-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, 622, 15, 57, "Input"],
Cell[1192, 39, 563, 14, 82, "Output"]
}, Open  ]],
Cell[1770, 56, 139, 2, 55, "Text"],
Cell[CellGroupData[{
Cell[1934, 62, 295, 5, 57, "Input"],
Cell[2232, 69, 339, 7, 82, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2608, 81, 190, 4, 57, "Input"],
Cell[2801, 87, 177, 4, 82, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3015, 96, 245, 7, 57, "Input"],
Cell[3263, 105, 231, 7, 82, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3531, 117, 451, 10, 92, "Input"],
Cell[3985, 129, 742, 22, 119, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4764, 156, 203, 4, 57, "Input"],
Cell[4970, 162, 475, 13, 84, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5482, 180, 546, 15, 126, "Input"],
Cell[6031, 197, 471, 12, 104, "Output"]
}, Open  ]],
Cell[6517, 212, 117, 1, 55, "Text"],
Cell[CellGroupData[{
Cell[6659, 217, 485, 14, 92, "Input"],
Cell[7147, 233, 555, 17, 104, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7739, 255, 249, 4, 57, "Input"],
Cell[7991, 261, 524, 13, 104, "Output"]
}, Open  ]],
Cell[8530, 277, 226, 5, 88, "Text"],
Cell[CellGroupData[{
Cell[8781, 286, 794, 21, 228, "Input"],
Cell[9578, 309, 763, 23, 127, "Output"],
Cell[10344, 334, 739, 23, 102, "Output"],
Cell[11086, 359, 715, 23, 102, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11838, 387, 341, 9, 57, "Input"],
Cell[12182, 398, 671, 19, 104, "Output"],
Cell[12856, 419, 647, 19, 79, "Output"],
Cell[13506, 440, 619, 18, 79, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14162, 463, 522, 15, 126, "Input"],
Cell[14687, 480, 396, 11, 82, "Output"],
Cell[15086, 493, 372, 11, 57, "Output"],
Cell[15461, 506, 348, 11, 57, "Output"]
}, Open  ]],
Cell[15824, 520, 161, 3, 88, "Text"],
Cell[CellGroupData[{
Cell[16010, 527, 487, 14, 92, "Input"],
Cell[16500, 543, 474, 14, 82, "Output"],
Cell[16977, 559, 450, 14, 57, "Output"],
Cell[17430, 575, 426, 14, 57, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17893, 594, 199, 4, 57, "Input"],
Cell[18095, 600, 437, 12, 109, "Output"],
Cell[18535, 614, 415, 12, 84, "Output"],
Cell[18953, 628, 391, 12, 84, "Output"]
}, Open  ]],
Cell[19359, 643, 182, 5, 88, "Text"],
Cell[CellGroupData[{
Cell[19566, 652, 172, 3, 57, "Input"],
Cell[19741, 657, 631, 22, 162, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[20409, 684, 808, 24, 130, "Input"],
Cell[21220, 710, 319, 11, 104, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21576, 726, 176, 3, 57, "Input"],
Cell[21755, 731, 330, 12, 104, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
