(* 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[      5920,        202]
NotebookOptionsPosition[      5313,        176]
NotebookOutlinePosition[      5677,        192]
CellTagsIndexPosition[      5634,        189]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{"Integrate", "[", 
  RowBox[{
   FractionBox[
    SuperscriptBox["x", 
     FractionBox["a", "2"]], 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"b", "+", 
       RowBox[{"a", " ", "x"}]}], ")"}], 
     FractionBox[
      RowBox[{"(", 
       RowBox[{"a", "+", "b"}], ")"}], "2"]]], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "0", ",", "\[Infinity]"}], "}"}], ",", 
   RowBox[{"Assumptions", "\[Rule]", 
    RowBox[{
     RowBox[{"Element", "[", 
      RowBox[{"a", ",", "Reals"}], "]"}], "&&", 
     RowBox[{"a", ">", "0"}], "&&", 
     RowBox[{"Element", "[", 
      RowBox[{"b", ",", "Reals"}], "]"}], "&&", 
     RowBox[{"b", ">", "2"}]}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.4803226871198344`*^9, 3.48032270641708*^9}, {
  3.4803227771215625`*^9, 3.4803227986376004`*^9}, {3.480322893467546*^9, 
  3.4803229081865788`*^9}, {3.4803230498455486`*^9, 3.480323053814375*^9}, {
  3.480323161394565*^9, 3.480323168097819*^9}, {3.480323287428235*^9, 
  3.480323287537612*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox["a", 
    RowBox[{
     RowBox[{"-", "1"}], "-", 
     FractionBox["a", "2"]}]], " ", 
   SuperscriptBox["b", 
    RowBox[{"1", "-", 
     FractionBox["b", "2"]}]], " ", 
   RowBox[{"Gamma", "[", 
    RowBox[{"1", "+", 
     FractionBox["a", "2"]}], "]"}], " ", 
   RowBox[{"Gamma", "[", 
    RowBox[{
     RowBox[{"-", "1"}], "+", 
     FractionBox["b", "2"]}], "]"}]}], 
  RowBox[{"Gamma", "[", 
   FractionBox[
    RowBox[{"a", "+", "b"}], "2"], "]"}]]], "Output",
 CellChangeTimes->{3.480322754527378*^9, 3.4803228098253155`*^9, 
  3.480322918671155*^9, 3.480323059095726*^9, 3.4803231713478813`*^9, 
  3.4803232893970227`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Integrate", "[", 
  RowBox[{
   FractionBox[
    SuperscriptBox["x", 
     RowBox[{
      FractionBox["a", "2"], "+", "1"}]], 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{"b", "+", 
       RowBox[{"a", " ", "x"}]}], ")"}], 
     FractionBox[
      RowBox[{"(", 
       RowBox[{"a", "+", "b"}], ")"}], "2"]]], ",", 
   RowBox[{"{", 
    RowBox[{"x", ",", "0", ",", "\[Infinity]"}], "}"}], ",", 
   RowBox[{"Assumptions", "\[Rule]", 
    RowBox[{
     RowBox[{"Element", "[", 
      RowBox[{"a", ",", "Reals"}], "]"}], "&&", 
     RowBox[{"a", ">", "0"}], "&&", 
     RowBox[{"Element", "[", 
      RowBox[{"b", ",", "Reals"}], "]"}], "&&", 
     RowBox[{"b", ">", "4"}]}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.4803226871198344`*^9, 3.48032270641708*^9}, {
   3.4803227771215625`*^9, 3.4803227986376004`*^9}, {3.480322893467546*^9, 
   3.4803229081865788`*^9}, {3.4803230498455486`*^9, 3.480323053814375*^9}, {
   3.480323161394565*^9, 3.480323168097819*^9}, {3.480323287428235*^9, 
   3.480323287537612*^9}, {3.4803246744079895`*^9, 3.480324682345642*^9}, 
   3.4803248321141424`*^9}],

Cell[BoxData[
 FractionBox[
  RowBox[{
   SuperscriptBox["a", 
    RowBox[{
     RowBox[{"-", "2"}], "-", 
     FractionBox["a", "2"]}]], " ", 
   SuperscriptBox["b", 
    RowBox[{"2", "-", 
     FractionBox["b", "2"]}]], " ", 
   RowBox[{"Gamma", "[", 
    RowBox[{"2", "+", 
     FractionBox["a", "2"]}], "]"}], " ", 
   RowBox[{"Gamma", "[", 
    RowBox[{
     RowBox[{"-", "2"}], "+", 
     FractionBox["b", "2"]}], "]"}]}], 
  RowBox[{"Gamma", "[", 
   FractionBox[
    RowBox[{"a", "+", "b"}], "2"], "]"}]]], "Output",
 CellChangeTimes->{3.480322754527378*^9, 3.4803228098253155`*^9, 
  3.480322918671155*^9, 3.480323059095726*^9, 3.4803231713478813`*^9, 
  3.4803232893970227`*^9, 3.4803246864082203`*^9, 3.4803248331141615`*^9}]
}, Open  ]],

Cell["Compute the variance using these expressions", "Text",
 CellChangeTimes->{{3.4803252918729696`*^9, 3.48032530177941*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Simplify", "[", " ", 
  RowBox[{
   RowBox[{
    SuperscriptBox["b", "2"], 
    FractionBox[
     RowBox[{"(", 
      RowBox[{"a", "+", "2"}], ")"}], 
     RowBox[{"a", " ", 
      RowBox[{"(", 
       RowBox[{"b", "-", "2"}], ")"}], " ", 
      RowBox[{"(", 
       RowBox[{"b", "-", "4"}], ")"}]}]]}], "-", 
   SuperscriptBox[
    RowBox[{"(", 
     FractionBox["b", 
      RowBox[{"b", "-", "2"}]], ")"}], "2"]}], "]"}]], "Input",
 CellChangeTimes->{{3.4803252005587163`*^9, 3.4803252372469206`*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{"2", " ", 
   SuperscriptBox["b", "2"], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "2"}], "+", "a", "+", "b"}], ")"}]}], 
  RowBox[{"a", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", "4"}], "+", "b"}], ")"}], " ", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "2"}], "+", "b"}], ")"}], "2"]}]]], "Output",
 CellChangeTimes->{3.480325237684429*^9}]
}, Open  ]]
},
WindowSize->{817, 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, 1015, 27, 66, "Input"],
Cell[1585, 51, 688, 22, 72, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2310, 78, 1122, 29, 66, "Input"],
Cell[3435, 109, 736, 22, 72, "Output"]
}, Open  ]],
Cell[4186, 134, 127, 1, 29, "Text"],
Cell[CellGroupData[{
Cell[4338, 139, 527, 17, 49, "Input"],
Cell[4868, 158, 429, 15, 52, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)