(* 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[     10462,        340]
NotebookOptionsPosition[      9374,        297]
NotebookOutlinePosition[      9738,        313]
CellTagsIndexPosition[      9695,        310]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"EX1", "=", 
   FractionBox[
    RowBox[{"(", " ", 
     RowBox[{
      RowBox[{"\[Alpha]", " ", "r1"}], " ", "-", " ", "r2"}], " ", ")"}], 
    RowBox[{"\[Alpha]", "-", "1"}]]}], ";"}]], "Input",
 CellChangeTimes->{{3.4789419859254665`*^9, 3.4789419954412127`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"expectation", " ", "of", " ", "X_", 
    RowBox[{"1", "^", "2"}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"EX12", " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"\[Alpha]", " ", 
      FractionBox[
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"r2", "-", "r1"}], ")"}], "2"], 
       RowBox[{"\[Alpha]", "-", "2"}]]}], "-", 
     SuperscriptBox["r2", "2"], "+", 
     RowBox[{"2", " ", "r2", " ", "EX1"}]}], "//", "Simplify"}]}]}]], "Input",
 CellChangeTimes->{{3.478858845531761*^9, 3.478858894766136*^9}, {
  3.478858977875511*^9, 3.478858979500511*^9}, {3.478859073266136*^9, 
  3.478859077766136*^9}, {3.478859270265511*^9, 3.478859282843636*^9}, {
  3.478942036394862*^9, 3.478942037254248*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"2", " ", 
    SuperscriptBox["r2", "2"]}], "-", 
   RowBox[{"2", " ", "r1", " ", "r2", " ", "\[Alpha]"}], "+", 
   RowBox[{
    SuperscriptBox["r1", "2"], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "\[Alpha]"}], ")"}], " ", "\[Alpha]"}]}], 
  RowBox[{"2", "-", 
   RowBox[{"3", " ", "\[Alpha]"}], "+", 
   SuperscriptBox["\[Alpha]", "2"]}]]], "Output",
 CellChangeTimes->{3.478858911203636*^9, 3.478858979906761*^9, 
  3.478859078594261*^9, 3.478859283812386*^9, 3.4789406805227575`*^9, 
  3.4789408671970215`*^9, 3.47894203776988*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"the", " ", "variance", " ", "of", " ", "X_", "1"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"VarX1", "=", 
   RowBox[{"Simplify", "[", " ", 
    RowBox[{"EX12", " ", "-", " ", 
     SuperscriptBox[
      RowBox[{"(", "EX1", ")"}], "2"]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.478859534312386*^9, 3.478859565218636*^9}, {
  3.478859608171761*^9, 3.478859614937386*^9}, {3.4789420050350857`*^9, 
  3.478942006394478*^9}, {3.4789420440512104`*^9, 3.478942044723094*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"r1", "-", "r2"}], ")"}], "2"], " ", "\[Alpha]"}], 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "\[Alpha]"}], ")"}], "2"]}]]], "Output",
 CellChangeTimes->{3.478859565562386*^9, 3.478940680569633*^9, 
  3.4789420070507364`*^9, 3.478942046019985*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"EX2", "=", 
   FractionBox[
    RowBox[{"(", " ", 
     RowBox[{
      RowBox[{"\[Alpha]", " ", "r2"}], " ", "-", " ", "r1"}], " ", ")"}], 
    RowBox[{"\[Alpha]", "-", "1"}]]}], ";"}]], "Input",
 CellChangeTimes->{{3.4789419859254665`*^9, 3.478942020816538*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"the", " ", "expectation", " ", "of", " ", "X_", 
    RowBox[{"2", "^", "2"}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"EX22", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"\[Alpha]", " ", 
      FractionBox[
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"r2", "-", "r1"}], ")"}], "2"], 
       RowBox[{"\[Alpha]", "-", "2"}]]}], "-", 
     SuperscriptBox["r1", "2"], "+", 
     RowBox[{"2", " ", "r1", " ", "EX2"}]}], "//", "Simplify"}]}]}]], "Input",
 CellChangeTimes->{{3.478859880124886*^9, 3.4788599298283863`*^9}, {
  3.4788600851252613`*^9, 3.4788601080158863`*^9}, {3.4789406697257442`*^9, 
  3.4789406738820477`*^9}, {3.4789420518169346`*^9, 3.4789420535044565`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"2", " ", 
    SuperscriptBox["r1", "2"]}], "-", 
   RowBox[{"2", " ", "r1", " ", "r2", " ", "\[Alpha]"}], "+", 
   RowBox[{
    SuperscriptBox["r2", "2"], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "\[Alpha]"}], ")"}], " ", "\[Alpha]"}]}], 
  RowBox[{"2", "-", 
   RowBox[{"3", " ", "\[Alpha]"}], "+", 
   SuperscriptBox["\[Alpha]", "2"]}]]], "Output",
 CellChangeTimes->{{3.4788599225940113`*^9, 3.4788599343127613`*^9}, {
   3.4788600856877613`*^9, 3.4788601083596363`*^9}, 3.4789406806008835`*^9, 
   3.4789420539575872`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"the", " ", "variance", " ", "of", " ", "X_", "2"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"VarX2", "=", 
   RowBox[{"Simplify", "[", " ", 
    RowBox[{"EX22", " ", "-", " ", 
     SuperscriptBox[
      RowBox[{"(", "EX2", ")"}], "2"]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.478859534312386*^9, 3.478859565218636*^9}, {
  3.478859608171761*^9, 3.478859614937386*^9}, {3.4788599497658863`*^9, 
  3.4788599534846363`*^9}, {3.478940961448228*^9, 3.478940967667058*^9}, {
  3.478942025675975*^9, 3.478942026785364*^9}, {3.478942057551383*^9, 
  3.4789420582076416`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"r1", "-", "r2"}], ")"}], "2"], " ", "\[Alpha]"}], 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "\[Alpha]"}], ")"}], "2"]}]]], "Output",
 CellChangeTimes->{3.478859565562386*^9, 3.4788599571408863`*^9, 
  3.478940968385817*^9, 3.478942027363497*^9, 3.4789420587388983`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"compute", " ", "E", 
    RowBox[{"(", 
     RowBox[{"X1", " ", "X2"}], ")"}]}], "  ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"EX1X2", " ", "=", " ", 
   RowBox[{"Simplify", "[", 
    RowBox[{
     FractionBox[
      RowBox[{"-", "1"}], "2"], 
     RowBox[{"(", " ", 
      RowBox[{
       RowBox[{"\[Alpha]", " ", 
        RowBox[{"(", 
         RowBox[{"\[Alpha]", "+", "1"}], ")"}], " ", 
        FractionBox[
         SuperscriptBox[
          RowBox[{"(", 
           RowBox[{"r2", "-", "r1"}], ")"}], "2"], 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"\[Alpha]", "-", "1"}], ")"}], " ", 
          RowBox[{"(", 
           RowBox[{"\[Alpha]", "-", "2"}], ")"}]}]]}], "-", "EX12", "-", 
       "EX22"}], ")"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.478861011231011*^9, 3.478861083090386*^9}, {
  3.478861310746636*^9, 3.478861316356011*^9}, {3.478941262306829*^9, 
  3.4789412716038227`*^9}, {3.478941968831497*^9, 3.478941975753461*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox["r1", "2"], "+", 
   SuperscriptBox["r2", "2"], "-", 
   RowBox[{"r1", " ", "r2", " ", "\[Alpha]"}]}], 
  RowBox[{"2", "-", "\[Alpha]"}]]], "Output",
 CellChangeTimes->{{3.478861065824761*^9, 3.478861083340386*^9}, 
   3.478941060652623*^9, 3.478941272463209*^9, {3.478941971690909*^9, 
   3.4789419765034704`*^9}, 3.478942062348319*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"EX1X2", " ", "-", " ", 
   RowBox[{"EX1", " ", "EX2"}]}], "//", "Simplify"}]], "Input",
 CellChangeTimes->{{3.4789420693796597`*^9, 3.4789420828329563`*^9}}],

Cell[BoxData[
 RowBox[{"-", 
  FractionBox[
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"r1", "-", "r2"}], ")"}], "2"], 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}], " ", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", "\[Alpha]"}], ")"}], "2"]}]]}]], "Output",
 CellChangeTimes->{{3.4789420756141143`*^9, 3.478942083129835*^9}, 
   3.4789422151945477`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Sqrt", "[", 
   RowBox[{"VarX1", " ", "VarX2"}], "]"}], "//", 
  RowBox[{
   RowBox[{"Simplify", "[", 
    RowBox[{"#", ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"\[Alpha]", ">", "2"}], ",", 
       RowBox[{"r2", ">", "r1"}]}], "}"}]}], "]"}], "&"}]}]], "Input",
 CellChangeTimes->{{3.4789420891924133`*^9, 3.4789421362399063`*^9}, {
  3.478942194569152*^9, 3.4789422215384192`*^9}, {3.4789423579160376`*^9, 
  3.478942361619234*^9}, {3.4789423954480085`*^9, 3.4789424759808044`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{"r1", "-", "r2"}], ")"}], "2"], " ", "\[Alpha]"}], 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "2"}], "+", "\[Alpha]"}], ")"}], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "1"}], "+", "\[Alpha]"}], ")"}], "2"]}]]], "Output",
 CellChangeTimes->{{3.478942467996276*^9, 3.478942476308936*^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[545, 20, 300, 8, 47, "Input"],
Cell[CellGroupData[{
Cell[870, 32, 790, 19, 70, "Input"],
Cell[1663, 53, 612, 16, 51, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2312, 74, 541, 12, 54, "Input"],
Cell[2856, 88, 474, 15, 52, "Output"]
}, Open  ]],
Cell[3345, 106, 298, 8, 47, "Input"],
Cell[CellGroupData[{
Cell[3668, 118, 757, 18, 70, "Input"],
Cell[4428, 138, 605, 16, 51, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5070, 159, 641, 14, 54, "Input"],
Cell[5714, 175, 498, 15, 52, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6249, 195, 1028, 28, 75, "Input"],
Cell[7280, 225, 394, 9, 50, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7711, 239, 192, 4, 31, "Input"],
Cell[7906, 245, 453, 15, 52, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8396, 265, 532, 13, 31, "Input"],
Cell[8931, 280, 427, 14, 52, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)