(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     52679,       1425]
NotebookOptionsPosition[     50208,       1341]
NotebookOutlinePosition[     50737,       1362]
CellTagsIndexPosition[     50694,       1359]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.6766953708609076`*^9, 3.676695388512907*^9}, 
   3.6766985071433496`*^9}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{
   "Edgeworth", " ", "expansion", " ", "for", " ", "Stirling", " ", "numbers",
     " ", "of", " ", "the", " ", "first", " ", "kind", "\[IndentingNewLine]", 
    "We", " ", "use", " ", "the", " ", "expansion", " ", "and", " ", "the", 
    " ", "notation", " ", "from", " ", "the", " ", "paper", " ", 
    "\[IndentingNewLine]", 
    "\"\<General Edgeworth expansions with applications to profiles of random \
trees\>\"", "\n", " ", "by", " ", "Zakhar", " ", "Kabluchko"}], ",", " ", 
   RowBox[{"Alexander", " ", "Marynych"}], ",", " ", 
   RowBox[{
    RowBox[{
    "Henning", " ", "Sulzbach", "\[IndentingNewLine]", "available", " ", "at",
      " ", 
     RowBox[{"http", ":"}]}], "//", 
    RowBox[{
     RowBox[{
      RowBox[{"arxiv", ".", "org"}], "/", "abs"}], "/", "1606.03920"}]}]}], 
  "\[IndentingNewLine]", "*)"}]], "Input",
 CellChangeTimes->{{3.6766984342063494`*^9, 3.6766984342531495`*^9}, {
   3.676700013434848*^9, 3.676700013918442*^9}, 3.676705964320282*^9, {
   3.6767754215438957`*^9, 3.6767754466576138`*^9}, {3.6767777443720407`*^9, 
   3.676777746056841*^9}, {3.682065494451191*^9, 3.6820655059656954`*^9}, {
   3.682065548297021*^9, 3.6820656098614664`*^9}, {3.682066386102324*^9, 
   3.6820663963514585`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"Parameters", ":", " ", "\[IndentingNewLine]", 
      RowBox[{"phi", 
       RowBox[{"(", "beta", ")"}]}]}], " ", "=", " ", 
     RowBox[{
      RowBox[{"Exp", "[", "beta", "]"}], "-", "1"}]}], ",", 
    "\[IndentingNewLine]", "  ", 
    RowBox[{
     RowBox[{"LogW", "[", "beta", "]"}], " ", "=", " ", 
     RowBox[{
      RowBox[{"LogGamma", "[", "theta", "]"}], "-", 
      RowBox[{"LogGamma", "[", 
       RowBox[{"theta", "*", 
        RowBox[{"Exp", "[", "beta", "]"}]}], "]"}]}]}], ",", 
    "\[IndentingNewLine]", "  ", 
    RowBox[{
     RowBox[{"kappa", "[", 
      RowBox[{"[", "j", "]"}], "]"}], " ", "=", " ", 
     RowBox[{
      RowBox[{"j", "-", 
       RowBox[{"th", " ", "derivative", " ", "of", " ", "phi", 
        RowBox[{"(", "beta", ")"}]}]}], "=", 
      RowBox[{
       RowBox[{
        RowBox[{"Exp", "[", "beta", "]"}], "-", 
        RowBox[{"1", "  ", "at", " ", "beta"}]}], "=", "0"}]}]}], ",", 
    "\[IndentingNewLine]", "  ", 
    RowBox[{
     RowBox[{"chi", "[", 
      RowBox[{"[", "j", "]"}], "]"}], " ", "=", " ", 
     RowBox[{
      RowBox[{"j", "-", 
       RowBox[{
       "th", " ", "derivative", " ", "of", " ", "LogW", " ", "at", " ", 
        "beta"}]}], "=", "0"}]}], ",", "\[IndentingNewLine]", "  ", 
    RowBox[{
     RowBox[{"sigma", "[", 
      RowBox[{"[", "j", "]"}], "]"}], " ", "=", " ", 
     RowBox[{"1", " ", "is", " ", "ignored", " ", 
      RowBox[{"everywhere", "."}]}]}]}], "\[IndentingNewLine]", " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"kappa", "=", 
    RowBox[{"Table", "[", 
     RowBox[{"1", ",", 
      RowBox[{"{", 
       RowBox[{"n", ",", "1", ",", "5"}], "}"}]}], "]"}]}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"LogW", "[", "beta_", "]"}], " ", "=", " ", 
    RowBox[{
     RowBox[{"LogGamma", "[", "theta", "]"}], "-", 
     RowBox[{"LogGamma", "[", 
      RowBox[{"theta", "*", 
       RowBox[{"Exp", "[", "beta", "]"}]}], "]"}]}]}], "\[IndentingNewLine]", 
   RowBox[{"chi", "=", 
    RowBox[{"List", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"D", "[", 
        RowBox[{
         RowBox[{"LogW", "[", "x", "]"}], ",", "x"}], "]"}], "/.", 
       RowBox[{"x", "\[Rule]", "0"}]}], ",", 
      RowBox[{
       RowBox[{"D", "[", 
        RowBox[{
         RowBox[{"D", "[", 
          RowBox[{
           RowBox[{"LogW", "[", "x", "]"}], ",", "x"}], "]"}], ",", "x"}], 
        "]"}], "/.", 
       RowBox[{"x", "\[Rule]", "0"}]}], ",", " ", 
      RowBox[{
       RowBox[{"D", "[", 
        RowBox[{
         RowBox[{"D", "[", 
          RowBox[{
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"LogW", "[", "x", "]"}], ",", "x"}], "]"}], ",", "x"}], 
          "]"}], ",", "x"}], "]"}], "/.", 
       RowBox[{"x", "\[Rule]", "0"}]}]}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.6766990894002905`*^9, 3.6766991123512893`*^9}, {
   3.6766992458466654`*^9, 3.6766993068768744`*^9}, {3.6766993474868207`*^9, 
   3.6766993884047585`*^9}, {3.6767000112936754`*^9, 
   3.6767000327398005`*^9}, {3.6767011128983917`*^9, 3.676701114224349*^9}, {
   3.6767021719437885`*^9, 3.676702198801972*^9}, {3.676702433645372*^9, 
   3.6767024788897715`*^9}, {3.6767025377599716`*^9, 3.676702706112172*^9}, 
   3.6767059666758823`*^9, 3.6767754488868265`*^9, {3.6767777569144406`*^9, 
   3.6767778614126406`*^9}, {3.682065627069825*^9, 3.682065628582987*^9}, {
   3.6820656602713575`*^9, 3.682065727059101*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"1", ",", "1", ",", "1", ",", "1", ",", "1"}], "}"}]], "Output",
 CellChangeTimes->{{3.6766992790240526`*^9, 3.67669930754767*^9}, {
   3.6766993533367834`*^9, 3.676699388919555*^9}, 3.6767000346273766`*^9, 
   3.676701115191518*^9, 3.676702102044732*^9, {3.6767021797592382`*^9, 
   3.676702199940772*^9}, {3.676702437843772*^9, 3.676702479357772*^9}, 
   3.676702580213972*^9, {3.676702627590572*^9, 3.676702669060772*^9}, 
   3.676702706767372*^9, 3.6767029057636547`*^9, {3.676705969842682*^9, 
   3.676705990713882*^9}, 3.6767775566573954`*^9, 3.6767777617972403`*^9, {
   3.6767778371350403`*^9, 3.6767778623038406`*^9}, 3.6820664213912983`*^9, {
   3.6822389825666385`*^9, 3.6822389854106383`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"LogGamma", "[", "theta", "]"}], "-", 
  RowBox[{"LogGamma", "[", 
   RowBox[{
    SuperscriptBox["\[ExponentialE]", "beta"], " ", "theta"}], 
   "]"}]}]], "Output",
 CellChangeTimes->{{3.6766992790240526`*^9, 3.67669930754767*^9}, {
   3.6766993533367834`*^9, 3.676699388919555*^9}, 3.6767000346273766`*^9, 
   3.676701115191518*^9, 3.676702102044732*^9, {3.6767021797592382`*^9, 
   3.676702199940772*^9}, {3.676702437843772*^9, 3.676702479357772*^9}, 
   3.676702580213972*^9, {3.676702627590572*^9, 3.676702669060772*^9}, 
   3.676702706767372*^9, 3.6767029057636547`*^9, {3.676705969842682*^9, 
   3.676705990713882*^9}, 3.6767775566573954`*^9, 3.6767777617972403`*^9, {
   3.6767778371350403`*^9, 3.6767778623038406`*^9}, 3.6820664213912983`*^9, {
   3.6822389825666385`*^9, 3.6822389854116383`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "theta"}], " ", 
    RowBox[{"PolyGamma", "[", 
     RowBox[{"0", ",", "theta"}], "]"}]}], ",", 
   RowBox[{
    RowBox[{
     RowBox[{"-", "theta"}], " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"0", ",", "theta"}], "]"}]}], "-", 
    RowBox[{
     SuperscriptBox["theta", "2"], " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"1", ",", "theta"}], "]"}]}]}], ",", 
   RowBox[{
    RowBox[{
     RowBox[{"-", "theta"}], " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"0", ",", "theta"}], "]"}]}], "-", 
    RowBox[{"3", " ", 
     SuperscriptBox["theta", "2"], " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"1", ",", "theta"}], "]"}]}], "-", 
    RowBox[{
     SuperscriptBox["theta", "3"], " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"2", ",", "theta"}], "]"}]}]}]}], "}"}]], "Output",
 CellChangeTimes->{{3.6766992790240526`*^9, 3.67669930754767*^9}, {
   3.6766993533367834`*^9, 3.676699388919555*^9}, 3.6767000346273766`*^9, 
   3.676701115191518*^9, 3.676702102044732*^9, {3.6767021797592382`*^9, 
   3.676702199940772*^9}, {3.676702437843772*^9, 3.676702479357772*^9}, 
   3.676702580213972*^9, {3.676702627590572*^9, 3.676702669060772*^9}, 
   3.676702706767372*^9, 3.6767029057636547`*^9, {3.676705969842682*^9, 
   3.676705990713882*^9}, 3.6767775566573954`*^9, 3.6767777617972403`*^9, {
   3.6767778371350403`*^9, 3.6767778623038406`*^9}, 3.6820664213912983`*^9, {
   3.6822389825666385`*^9, 3.682238985413638*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"Differential", " ", "operators", " ", "D_", "1"}], ",", 
   RowBox[{"D_", "2"}], ",", 
   RowBox[{"D_", "3.", " ", "Here", " ", "D", " ", "stays", " ", "for", " ", 
    RowBox[{"sigma", "^", 
     RowBox[{"(", 
      RowBox[{"-", "1"}], ")"}]}], " ", 
    RowBox[{"d", "/", "dx"}]}]}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.6767755582456837`*^9, 3.676775569547636*^9}, {
  3.6767778686394405`*^9, 3.6767778789998407`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"D1", "=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"kappa", "[", 
      RowBox[{"[", "3", "]"}], "]"}], "/", "6"}], "*", 
    RowBox[{"D", "^", "3"}]}], " ", "+", " ", 
   RowBox[{
    RowBox[{"chi", "[", 
     RowBox[{"[", "1", "]"}], "]"}], "*", 
    RowBox[{"D", "^", "1"}]}]}]}], "\[IndentingNewLine]", 
 RowBox[{"D2", "=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"kappa", "[", 
      RowBox[{"[", "4", "]"}], "]"}], "/", "12"}], "*", 
    RowBox[{"D", "^", "4"}]}], "+", " ", 
   RowBox[{
    RowBox[{"chi", "[", 
     RowBox[{"[", "2", "]"}], "]"}], "*", 
    RowBox[{"D", "^", "2"}]}]}]}], "\[IndentingNewLine]", 
 RowBox[{"D3", "=", " ", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"kappa", "[", 
      RowBox[{"[", "5", "]"}], "]"}], "/", "20"}], "*", 
    RowBox[{"D", "^", "5"}]}], " ", "+", " ", 
   RowBox[{
    RowBox[{"chi", "[", 
     RowBox[{"[", "3", "]"}], "]"}], "*", 
    RowBox[{"D", "^", "3"}]}]}]}]}], "Input",
 CellChangeTimes->{{3.67669866854315*^9, 3.6766987608397493`*^9}, {
  3.6766988004747496`*^9, 3.676698810304648*^9}, {3.676699007852847*^9, 
  3.676699031917276*^9}, {3.6767011199825644`*^9, 3.676701130453829*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox[
   SuperscriptBox["D", "3"], "6"], "-", 
  RowBox[{"D", " ", "theta", " ", 
   RowBox[{"PolyGamma", "[", 
    RowBox[{"0", ",", "theta"}], "]"}]}]}]], "Output",
 CellChangeTimes->{
  3.6766987227135496`*^9, 3.6766988147857046`*^9, 3.6766990334460273`*^9, 
   3.676699117282858*^9, 3.676699400965678*^9, 3.6767000421932793`*^9, 
   3.676701130921814*^9, 3.676702673370372*^9, 3.6767027138673716`*^9, {
   3.676702902610572*^9, 3.676702908290742*^9}, 3.676705992632682*^9, 
   3.6767775602609262`*^9, 3.6767778646146407`*^9, 3.6820664240432806`*^9, 
   3.6822389868016386`*^9}],

Cell[BoxData[
 RowBox[{
  FractionBox[
   SuperscriptBox["D", "4"], "12"], "+", 
  RowBox[{
   SuperscriptBox["D", "2"], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "theta"}], " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"0", ",", "theta"}], "]"}]}], "-", 
     RowBox[{
      SuperscriptBox["theta", "2"], " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"1", ",", "theta"}], "]"}]}]}], ")"}]}]}]], "Output",
 CellChangeTimes->{
  3.6766987227135496`*^9, 3.6766988147857046`*^9, 3.6766990334460273`*^9, 
   3.676699117282858*^9, 3.676699400965678*^9, 3.6767000421932793`*^9, 
   3.676701130921814*^9, 3.676702673370372*^9, 3.6767027138673716`*^9, {
   3.676702902610572*^9, 3.676702908290742*^9}, 3.676705992632682*^9, 
   3.6767775602609262`*^9, 3.6767778646146407`*^9, 3.6820664240432806`*^9, 
   3.6822389868046384`*^9}],

Cell[BoxData[
 RowBox[{
  FractionBox[
   SuperscriptBox["D", "5"], "20"], "+", 
  RowBox[{
   SuperscriptBox["D", "3"], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "theta"}], " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"0", ",", "theta"}], "]"}]}], "-", 
     RowBox[{"3", " ", 
      SuperscriptBox["theta", "2"], " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"1", ",", "theta"}], "]"}]}], "-", 
     RowBox[{
      SuperscriptBox["theta", "3"], " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"2", ",", "theta"}], "]"}]}]}], ")"}]}]}]], "Output",
 CellChangeTimes->{
  3.6766987227135496`*^9, 3.6766988147857046`*^9, 3.6766990334460273`*^9, 
   3.676699117282858*^9, 3.676699400965678*^9, 3.6767000421932793`*^9, 
   3.676701130921814*^9, 3.676702673370372*^9, 3.6767027138673716`*^9, {
   3.676702902610572*^9, 3.676702908290742*^9}, 3.676705992632682*^9, 
   3.6767775602609262`*^9, 3.6767778646146407`*^9, 3.6820664240432806`*^9, 
   3.6822389868086386`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
    "Bell", " ", "polynomials", " ", "are", " ", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"B_", "0"}], "=", "1"}]}], ",", "\[IndentingNewLine]", " ", 
    RowBox[{
     RowBox[{"B_", "1", 
      RowBox[{"(", 
       RowBox[{"z_", "1"}], ")"}]}], " ", "=", " ", 
     RowBox[{"z_", "1"}]}], ",", "\[IndentingNewLine]", " ", 
    RowBox[{
     RowBox[{"B_", "2", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"z_", "1"}], ",", 
        RowBox[{"z_", "2"}]}], ")"}]}], " ", "=", " ", 
     RowBox[{
      RowBox[{"z_", 
       RowBox[{"1", "^", "2"}]}], " ", "+", " ", 
      RowBox[{"z_", "2"}]}]}], ",", "\[IndentingNewLine]", " ", 
    RowBox[{
     RowBox[{"B_", "3", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"z_", "1"}], ",", 
        RowBox[{"z_", "2"}], ",", 
        RowBox[{"z_", "3"}]}], ")"}]}], " ", "=", " ", 
     RowBox[{
      RowBox[{"z_", 
       RowBox[{"1", "^", "3"}]}], " ", "+", " ", 
      RowBox[{"3", "z_", "1", "z_", "2"}], " ", "+", " ", 
      RowBox[{"z_", "3"}]}]}]}], "\[IndentingNewLine]", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{
     "Probabilist", " ", "Hermite", " ", "polynomials", " ", "are", " ", 
      "given", " ", "by", "\[IndentingNewLine]", 
      RowBox[{"He_p", 
       RowBox[{"(", "x", ")"}]}]}], " ", "=", " ", 
     RowBox[{
      RowBox[{"2", "^", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "p"}], "/", "2"}], ")"}]}], "*", 
      RowBox[{"HermiteH", "[", 
       RowBox[{"p", ",", " ", 
        RowBox[{"x", "/", 
         RowBox[{"Sqrt", "[", "2", "]"}]}]}], "]"}], "\[IndentingNewLine]", 
      "To", " ", "check", " ", "this"}]}], ",", " ", 
    RowBox[{"use", ":", " ", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Simplify", "[", 
        RowBox[{
         RowBox[{"2", "^", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"-", "p"}], "/", "2"}], ")"}]}], "*", 
         RowBox[{"HermiteH", "[", 
          RowBox[{"p", ",", " ", 
           RowBox[{"x", "/", 
            RowBox[{"Sqrt", "[", "2", "]"}]}]}], "]"}]}], "]"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"p", ",", "0", ",", "6"}], "}"}]}], "]"}]}]}], 
   "\[IndentingNewLine]", "*)"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "Terms", " ", "in", " ", "the", " ", "Edgeworth", " ", "expansion", " ", 
     "of", " ", "the", " ", "Stirling", " ", "numbers", " ", "are", " ", 
     "denoted", " ", "by", " ", "G_", "0", 
     RowBox[{"(", "x", ")"}]}], ",", " ", 
    RowBox[{"G_", "1", 
     RowBox[{"(", "x", ")"}]}], ",", "..."}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"H0", "=", "1"}], ";"}]}]], "Input",
 CellChangeTimes->{{3.676700088559405*^9, 3.6767001082853546`*^9}, {
   3.676700155608756*^9, 3.6767001576211567`*^9}, {3.676775484212273*^9, 
   3.676775502339847*^9}, {3.676775613182391*^9, 3.6767757137104945`*^9}, {
   3.676775825685906*^9, 3.6767758284152236`*^9}, {3.676775859738592*^9, 
   3.676775951844993*^9}, {3.6767760032739787`*^9, 3.67677603357446*^9}, {
   3.6767775647848387`*^9, 3.6767775718827024`*^9}, 3.6767778907798405`*^9, {
   3.6820657627222443`*^9, 3.6820657910693493`*^9}, {3.6822389722076387`*^9, 
   3.6822389725936384`*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.6767760513432655`*^9, 3.676776070566741*^9}, 
   3.6767775782337804`*^9, 3.6820657805081496`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"H1", "[", "x_", "]"}], "=", 
  RowBox[{"Simplify", "[", " ", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"Expand", "[", "D1", "]"}], "/.", " ", 
      RowBox[{
       RowBox[{"D", "^", "p_"}], " ", "->", " ", 
       RowBox[{
        RowBox[{"2", "^", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "p"}], "/", "2"}], ")"}]}], "*", 
        RowBox[{"HermiteH", "[", 
         RowBox[{"p", ",", 
          RowBox[{"x", "/", 
           RowBox[{"Sqrt", "[", "2", "]"}]}]}], "]"}]}]}]}], ")"}], " ", "/.",
     " ", 
    RowBox[{"D", " ", "->", " ", 
     RowBox[{
      RowBox[{"2", "^", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "1"}], "/", "2"}], ")"}]}], "*", 
      RowBox[{"HermiteH", "[", 
       RowBox[{"1", ",", 
        RowBox[{"x", "/", 
         RowBox[{"Sqrt", "[", "2", "]"}]}]}], "]"}]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{
  3.6766999299007597`*^9, {3.6767000644213943`*^9, 3.676700078681731*^9}, {
   3.6767001540155563`*^9, 3.6767001541735563`*^9}, {3.6767002940303564`*^9, 
   3.6767002987435565`*^9}, {3.676700385926156*^9, 3.6767003940899563`*^9}, {
   3.6767015336427383`*^9, 3.6767015338143377`*^9}, {3.676701651501738*^9, 
   3.676701693275938*^9}, 3.6767755062774453`*^9, {3.676775833156438*^9, 
   3.6767758418162727`*^9}, {3.6767776898830338`*^9, 
   3.6767776917393985`*^9}, {3.682238974708638*^9, 3.6822389750286384`*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "6"], " ", "x", " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"-", "3"}], "+", 
    SuperscriptBox["x", "2"], "-", 
    RowBox[{"6", " ", "theta", " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"0", ",", "theta"}], "]"}]}]}], ")"}]}]], "Output",
 CellChangeTimes->{{3.67670007957092*^9, 3.6767000868122272`*^9}, 
   3.6767001601483564`*^9, 3.6767003099181566`*^9, 3.6767004032999563`*^9, 
   3.6767011708181424`*^9, 3.6767016944615383`*^9, 3.6767026760379715`*^9, 
   3.6767027195737715`*^9, 3.676702911753786*^9, 3.676705996019882*^9, 
   3.676775770288158*^9, 3.6767758430795565`*^9, 3.676777579871749*^9, {
   3.676777680554414*^9, 3.6767776953117294`*^9}, 3.682066429253648*^9, 
   3.682238991826638*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"A11", "=", 
  RowBox[{"Simplify", "[", 
   RowBox[{"Coefficient", "[", 
    RowBox[{
     RowBox[{"H1", "[", "x", "]"}], ",", " ", "x"}], "]"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"A12", "=", 
  RowBox[{"Simplify", "[", 
   RowBox[{"Coefficient", "[", 
    RowBox[{
     RowBox[{"H1", "[", "x", "]"}], ",", " ", 
     RowBox[{"x", "^", "3"}]}], "]"}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.6822391074540024`*^9, 3.6822391627225285`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", 
   FractionBox["1", "2"]}], "-", 
  RowBox[{"theta", " ", 
   RowBox[{"PolyGamma", "[", 
    RowBox[{"0", ",", "theta"}], "]"}]}]}]], "Output",
 CellChangeTimes->{{3.682239135391796*^9, 3.6822391632975864`*^9}}],

Cell[BoxData[
 FractionBox["1", "6"]], "Output",
 CellChangeTimes->{{3.682239135391796*^9, 3.6822391632995863`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"H2", "[", "x_", "]"}], "=", 
  RowBox[{"Simplify", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"1", "/", 
      RowBox[{"(", 
       RowBox[{"2", "!"}], ")"}]}], "*", 
     RowBox[{"Expand", "[", 
      RowBox[{
       RowBox[{"D1", "^", "2"}], " ", "+", " ", "D2"}], "]"}]}], "/.", " ", 
    RowBox[{
     RowBox[{"D", "^", "p_"}], " ", "->", " ", 
     RowBox[{
      RowBox[{"2", "^", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "p"}], "/", "2"}], ")"}]}], "*", 
      RowBox[{"HermiteH", "[", 
       RowBox[{"p", ",", 
        RowBox[{"x", "/", 
         RowBox[{"Sqrt", "[", "2", "]"}]}]}], "]"}]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.6766999313983407`*^9, 3.676699972751169*^9}, {
  3.6767001513927565`*^9, 3.6767001515643563`*^9}, {3.6767003012571564`*^9, 
  3.6767003013351564`*^9}, {3.6767003808073564`*^9, 3.6767003969779563`*^9}, {
  3.676775749611328*^9, 3.6767757634570665`*^9}, {3.676775849020072*^9, 
  3.6767758515974274`*^9}, {3.6767760783660407`*^9, 3.676776079520337*^9}, {
  3.6822389769076385`*^9, 3.6822389772596383`*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "72"], " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"-", "6"}], "+", 
    RowBox[{"27", " ", 
     SuperscriptBox["x", "2"]}], "-", 
    RowBox[{"12", " ", 
     SuperscriptBox["x", "4"]}], "+", 
    SuperscriptBox["x", "6"], "-", 
    RowBox[{"12", " ", "theta", " ", 
     SuperscriptBox["x", "2"], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "3"}], "+", 
       SuperscriptBox["x", "2"]}], ")"}], " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"0", ",", "theta"}], "]"}]}], "+", 
    RowBox[{"36", " ", 
     SuperscriptBox["theta", "2"], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", 
       SuperscriptBox["x", "2"]}], ")"}], " ", 
     SuperscriptBox[
      RowBox[{"PolyGamma", "[", 
       RowBox[{"0", ",", "theta"}], "]"}], "2"]}], "-", 
    RowBox[{"36", " ", 
     SuperscriptBox["theta", "2"], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1"}], "+", 
       SuperscriptBox["x", "2"]}], ")"}], " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"1", ",", "theta"}], "]"}]}]}], ")"}]}]], "Output",
 CellChangeTimes->{3.676699973936754*^9, 3.6767000462708273`*^9, 
  3.6767001621919565`*^9, 3.6767003080149565`*^9, 3.6767004051251564`*^9, 
  3.67670117270973*^9, 3.676701699034338*^9, 3.6767026785183716`*^9, 
  3.6767027215569715`*^9, 3.676702915178033*^9, 3.6767059972210817`*^9, 
  3.676775764564389*^9, 3.6767775821825047`*^9, 3.6767776970900955`*^9, 
  3.6767779014698405`*^9, 3.682066432326828*^9, 3.6822389933576384`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"A21", "=", 
  RowBox[{"Simplify", "[", 
   RowBox[{"H2", "[", "0", "]"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"A22", "=", " ", 
  RowBox[{"Simplify", "[", 
   RowBox[{"Coefficient", "[", 
    RowBox[{
     RowBox[{"H2", "[", "x", "]"}], ",", " ", 
     RowBox[{"x", "^", "2"}]}], "]"}], 
   "]"}]}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.6822391974279985`*^9, 3.682239211122368*^9}, {
  3.6822392745817137`*^9, 3.682239302384494*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "12"], " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"-", "1"}], "-", 
    RowBox[{"6", " ", 
     SuperscriptBox["theta", "2"], " ", 
     SuperscriptBox[
      RowBox[{"PolyGamma", "[", 
       RowBox[{"0", ",", "theta"}], "]"}], "2"]}], "+", 
    RowBox[{"6", " ", 
     SuperscriptBox["theta", "2"], " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"1", ",", "theta"}], "]"}]}]}], ")"}]}]], "Output",
 CellChangeTimes->{
  3.6822392195032063`*^9, {3.682239299298185*^9, 3.68223930274953*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "8"], " ", 
  RowBox[{"(", 
   RowBox[{"3", "+", 
    RowBox[{"4", " ", "theta", " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"0", ",", "theta"}], "]"}]}], "+", 
    RowBox[{"4", " ", 
     SuperscriptBox["theta", "2"], " ", 
     SuperscriptBox[
      RowBox[{"PolyGamma", "[", 
       RowBox[{"0", ",", "theta"}], "]"}], "2"]}], "-", 
    RowBox[{"4", " ", 
     SuperscriptBox["theta", "2"], " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"1", ",", "theta"}], "]"}]}]}], ")"}]}]], "Output",
 CellChangeTimes->{
  3.6822392195032063`*^9, {3.682239299298185*^9, 3.6822393027545304`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"H3", "[", "x_", "]"}], "=", 
  RowBox[{"Simplify", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"1", "/", 
      RowBox[{"(", 
       RowBox[{"3", "!"}], ")"}]}], "*", 
     RowBox[{"Expand", "[", 
      RowBox[{
       RowBox[{"D1", "^", "3"}], " ", "+", " ", 
       RowBox[{"3", " ", "D1", "*", "D2"}], "+", "D3"}], "]"}]}], "/.", " ", 
    RowBox[{
     RowBox[{"D", "^", "p_"}], " ", "->", " ", 
     RowBox[{
      RowBox[{"2", "^", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"-", "p"}], "/", "2"}], ")"}]}], "*", 
      RowBox[{"HermiteH", "[", 
       RowBox[{"p", ",", 
        RowBox[{"x", "/", 
         RowBox[{"Sqrt", "[", "2", "]"}]}]}], "]"}]}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.676699475479638*^9, 3.6766995382428384`*^9}, {
  3.676699596202238*^9, 3.6766996372294383`*^9}, {3.6766997116288385`*^9, 
  3.6766997184168386`*^9}, {3.6766999243900304`*^9, 3.6766999267300005`*^9}, {
  3.6767001482551565`*^9, 3.6767001484891567`*^9}, {3.676700304519556*^9, 
  3.6767003052059565`*^9}, {3.6767003895473566`*^9, 3.6767003998015566`*^9}, {
  3.6767757532920074`*^9, 3.676775755069963*^9}, {3.676776081565754*^9, 
  3.676776093065922*^9}, {3.6822390704233*^9, 3.6822390707903366`*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "6480"], 
  RowBox[{"x", " ", 
   RowBox[{"(", 
    RowBox[{"810", "-", 
     RowBox[{"2115", " ", 
      SuperscriptBox["x", "2"]}], "+", 
     RowBox[{"999", " ", 
      SuperscriptBox["x", "4"]}], "-", 
     RowBox[{"135", " ", 
      SuperscriptBox["x", "6"]}], "+", 
     RowBox[{"5", " ", 
      SuperscriptBox["x", "8"]}], "+", 
     RowBox[{"540", " ", 
      SuperscriptBox["theta", "2"], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "3"}], "-", 
        RowBox[{"4", " ", 
         SuperscriptBox["x", "2"]}], "+", 
        SuperscriptBox["x", "4"]}], ")"}], " ", 
      SuperscriptBox[
       RowBox[{"PolyGamma", "[", 
        RowBox[{"0", ",", "theta"}], "]"}], "2"]}], "-", 
     RowBox[{"1080", " ", 
      SuperscriptBox["theta", "3"], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "3"}], "+", 
        SuperscriptBox["x", "2"]}], ")"}], " ", 
      SuperscriptBox[
       RowBox[{"PolyGamma", "[", 
        RowBox[{"0", ",", "theta"}], "]"}], "3"]}], "-", 
     RowBox[{"540", " ", 
      SuperscriptBox["theta", "2"], " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "3"}], "-", 
        RowBox[{"4", " ", 
         SuperscriptBox["x", "2"]}], "+", 
        SuperscriptBox["x", "4"]}], ")"}], " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"1", ",", "theta"}], "]"}]}], "+", 
     RowBox[{"90", " ", "theta", " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"0", ",", "theta"}], "]"}], " ", 
      RowBox[{"(", 
       RowBox[{"6", "-", 
        RowBox[{"27", " ", 
         SuperscriptBox["x", "2"]}], "+", 
        RowBox[{"12", " ", 
         SuperscriptBox["x", "4"]}], "-", 
        SuperscriptBox["x", "6"], "+", 
        RowBox[{"36", " ", 
         SuperscriptBox["theta", "2"], " ", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "3"}], "+", 
           SuperscriptBox["x", "2"]}], ")"}], " ", 
         RowBox[{"PolyGamma", "[", 
          RowBox[{"1", ",", "theta"}], "]"}]}]}], ")"}]}], "+", 
     RowBox[{"3240", " ", 
      SuperscriptBox["theta", "3"], " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"2", ",", "theta"}], "]"}]}], "-", 
     RowBox[{"1080", " ", 
      SuperscriptBox["theta", "3"], " ", 
      SuperscriptBox["x", "2"], " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"2", ",", "theta"}], "]"}]}]}], ")"}]}]}]], "Output",
 CellChangeTimes->{{3.6766996150030384`*^9, 3.6766996381030383`*^9}, 
   3.676699720710038*^9, 3.6767001648147564`*^9, 3.6767003058923564`*^9, 
   3.6767004119911566`*^9, 3.676701174113721*^9, 3.676701700721138*^9, 
   3.6767026809227715`*^9, 3.6767027262057714`*^9, 3.6767029165663705`*^9, 
   3.676705998422282*^9, 3.6767757731422443`*^9, 3.6767775840076694`*^9, 
   3.676777733352441*^9, 3.676777904012641*^9, 3.682066434604413*^9, 
   3.6822389951207705`*^9, 3.682239071469404*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"A31", "=", " ", 
   RowBox[{"Simplify", "[", 
    RowBox[{"Coefficient", "[", 
     RowBox[{
      RowBox[{"H3", "[", "x", "]"}], ",", "x"}], "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.682239684415635*^9, 3.682239706981635*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "24"], " ", 
  RowBox[{"(", 
   RowBox[{"3", "-", 
    RowBox[{"6", " ", 
     SuperscriptBox["theta", "2"], " ", 
     SuperscriptBox[
      RowBox[{"PolyGamma", "[", 
       RowBox[{"0", ",", "theta"}], "]"}], "2"]}], "+", 
    RowBox[{"12", " ", 
     SuperscriptBox["theta", "3"], " ", 
     SuperscriptBox[
      RowBox[{"PolyGamma", "[", 
       RowBox[{"0", ",", "theta"}], "]"}], "3"]}], "+", 
    RowBox[{"6", " ", 
     SuperscriptBox["theta", "2"], " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"1", ",", "theta"}], "]"}]}], "+", 
    RowBox[{
     RowBox[{"PolyGamma", "[", 
      RowBox[{"0", ",", "theta"}], "]"}], " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"2", " ", "theta"}], "-", 
       RowBox[{"36", " ", 
        SuperscriptBox["theta", "3"], " ", 
        RowBox[{"PolyGamma", "[", 
         RowBox[{"1", ",", "theta"}], "]"}]}]}], ")"}]}], "+", 
    RowBox[{"12", " ", 
     SuperscriptBox["theta", "3"], " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"2", ",", "theta"}], "]"}]}]}], ")"}]}]], "Output",
 CellChangeTimes->{3.682239708402635*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
    "The", " ", "first", " ", "three", " ", "terms", " ", "in", " ", "the", 
     " ", "expansion", " ", "for", " ", "x"}], " ", "=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"a", "/", 
       RowBox[{
        RowBox[{"Sqrt", "[", "w", "]"}], ".", " ", "\[IndentingNewLine]", 
        "The"}]}], " ", "term", " ", 
      RowBox[{"H4", "/", 
       RowBox[{"w", "^", "2"}]}]}], " ", "=", " ", 
     RowBox[{
      RowBox[{"constant", "/", 
       RowBox[{"w", "^", "2"}]}], " ", "+", 
      RowBox[{"o", 
       RowBox[{"(", 
        RowBox[{"1", "/", 
         RowBox[{"w", "^", "2"}]}], ")"}], " ", "is", " ", "ignored", " ", 
       "because", " ", "the", " ", "constant", " ", "does", " ", "not", " ", 
       "depend", " ", "on", " ", 
       RowBox[{"a", ".", "  ", "\[IndentingNewLine]", 
        RowBox[{"Error", ":", " ", 
         RowBox[{"o", 
          RowBox[{"(", 
           RowBox[{"1", "/", 
            RowBox[{"w", "^", "2"}]}], ")"}]}]}]}]}]}]}]}], 
   "\[IndentingNewLine]", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Series", "[", 
   RowBox[{
    RowBox[{"Exp", "[", 
     RowBox[{
      RowBox[{"-", 
       RowBox[{"t", "^", "2"}]}], "/", "2"}], "]"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"t", ",", "0", ",", "4"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.6820658448933496`*^9, 3.6820658458449497`*^9}, {
  3.682239035412799*^9, 3.6822390541566734`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"1", "-", 
   FractionBox[
    SuperscriptBox["t", "2"], "2"], "+", 
   FractionBox[
    SuperscriptBox["t", "4"], "8"], "+", 
   InterpretationBox[
    SuperscriptBox[
     RowBox[{"O", "[", "t", "]"}], "5"],
    SeriesData[$CellContext`t, 0, {}, 0, 5, 1],
    Editable->False]}],
  SeriesData[$CellContext`t, 0, {1, 0, 
    Rational[-1, 2], 0, 
    Rational[1, 8]}, 0, 5, 1],
  Editable->False]], "Output",
 CellChangeTimes->{{3.6767763098622017`*^9, 3.676776312935323*^9}, 
   3.6767775928289003`*^9, 3.6767779058564405`*^9, 3.682066439658781*^9, 
   3.682238997449003*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"A", "=", " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1", "+", 
        RowBox[{
         RowBox[{"H1", "[", 
          RowBox[{"a", "/", 
           RowBox[{"Sqrt", "[", "w", "]"}]}], "]"}], "/", 
         RowBox[{"Sqrt", "[", "w", "]"}]}], " ", "+", " ", 
        RowBox[{
         RowBox[{"H2", "[", 
          RowBox[{"a", "/", 
           RowBox[{"Sqrt", "[", "w", "]"}]}], "]"}], "/", "w"}], " ", "+", 
        RowBox[{
         RowBox[{"H3", "[", 
          RowBox[{"a", "/", 
           RowBox[{"Sqrt", "[", "w", "]"}]}], "]"}], "/", 
         RowBox[{"w", "^", 
          RowBox[{"(", 
           RowBox[{"3", "/", "2"}], ")"}]}]}]}], ")"}], " ", "*", " ", 
      RowBox[{"(", 
       RowBox[{"1", "-", 
        RowBox[{
         RowBox[{"a", "^", "2"}], "/", 
         RowBox[{"(", 
          RowBox[{"2", "w"}], ")"}]}], " ", "+", " ", 
        RowBox[{
         RowBox[{"1", "/", "8"}], " ", "*", 
         RowBox[{
          RowBox[{"a", "^", "4"}], "/", 
          RowBox[{"w", "^", "2"}]}]}]}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"Collect", "[", 
    RowBox[{
     RowBox[{"Coefficient", "[", 
      RowBox[{"A", ",", 
       RowBox[{"1", "/", "w"}]}], "]"}], ",", "a"}], "]"}], 
   "\[IndentingNewLine]", 
   RowBox[{"Collect", "[", 
    RowBox[{
     RowBox[{"Coefficient", "[", 
      RowBox[{"A", ",", 
       RowBox[{"1", "/", 
        RowBox[{"w", "^", 
         RowBox[{"(", 
          RowBox[{"1", "/", "2"}], ")"}]}]}]}], "]"}], ",", "a"}], "]"}], 
   "\[IndentingNewLine]", 
   RowBox[{"Collect", "[", 
    RowBox[{
     RowBox[{"Coefficient", "[", 
      RowBox[{"A", ",", 
       RowBox[{"1", "/", 
        RowBox[{"w", "^", 
         RowBox[{"(", 
          RowBox[{"3", "/", "2"}], ")"}]}]}]}], "]"}], ",", "a"}], "]"}], 
   "\[IndentingNewLine]", 
   RowBox[{"Collect", "[", 
    RowBox[{
     RowBox[{"Coefficient", "[", 
      RowBox[{"A", ",", 
       RowBox[{"1", "/", 
        RowBox[{"w", "^", "2"}]}]}], "]"}], ",", "a"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.6767004199355564`*^9, 3.6767004778757563`*^9}, {
   3.676700512837756*^9, 3.676700579632756*^9}, {3.6767006699363565`*^9, 
   3.676700682797756*^9}, {3.6767008573288865`*^9, 3.676700868047743*^9}, {
   3.676701794723138*^9, 3.676701801965538*^9}, {3.6767018756205378`*^9, 
   3.676701880021738*^9}, {3.676702046207006*^9, 3.676702055054036*^9}, {
   3.676703070769101*^9, 3.676703086526394*^9}, 3.676704129970925*^9, {
   3.676704176268525*^9, 3.676704176642925*^9}, {3.6767761094152546`*^9, 
   3.676776278701801*^9}, {3.6767763160864425`*^9, 3.676776319019167*^9}, 
   3.6820658428341494`*^9, {3.682239058283086*^9, 3.6822390638276405`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"-", 
   FractionBox["1", "12"]}], "-", 
  FractionBox[
   SuperscriptBox["a", "2"], "2"], "-", 
  RowBox[{
   FractionBox["1", "2"], " ", 
   SuperscriptBox["theta", "2"], " ", 
   SuperscriptBox[
    RowBox[{"PolyGamma", "[", 
     RowBox[{"0", ",", "theta"}], "]"}], "2"]}], "+", 
  RowBox[{"a", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", 
      FractionBox["1", "2"]}], "-", 
     RowBox[{"theta", " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"0", ",", "theta"}], "]"}]}]}], ")"}]}], "+", 
  RowBox[{
   FractionBox["1", "2"], " ", 
   SuperscriptBox["theta", "2"], " ", 
   RowBox[{"PolyGamma", "[", 
    RowBox[{"1", ",", "theta"}], "]"}]}]}]], "Output",
 CellChangeTimes->{{3.6767005074361563`*^9, 3.6767005218273563`*^9}, {
   3.6767005555987563`*^9, 3.6767005800851564`*^9}, 3.6767006835777564`*^9, 
   3.6767008689681134`*^9, 3.6767011759389095`*^9, 3.6767017025795383`*^9, 
   3.676701802558338*^9, {3.676701880396138*^9, 3.676701909702629*^9}, 
   3.6767027308751717`*^9, 3.6767029204661956`*^9, {3.6767030776327934`*^9, 
   3.676703086996373*^9}, {3.676704177502925*^9, 3.676704182682125*^9}, 
   3.676706001011882*^9, 3.6767761515590734`*^9, 3.6767762704466124`*^9, 
   3.6767776022979183`*^9, 3.6767779075256405`*^9, 3.6820664427631607`*^9, 
   3.68223907482874*^9}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{{3.6767005074361563`*^9, 3.6767005218273563`*^9}, {
   3.6767005555987563`*^9, 3.6767005800851564`*^9}, 3.6767006835777564`*^9, 
   3.6767008689681134`*^9, 3.6767011759389095`*^9, 3.6767017025795383`*^9, 
   3.676701802558338*^9, {3.676701880396138*^9, 3.676701909702629*^9}, 
   3.6767027308751717`*^9, 3.6767029204661956`*^9, {3.6767030776327934`*^9, 
   3.676703086996373*^9}, {3.676704177502925*^9, 3.676704182682125*^9}, 
   3.676706001011882*^9, 3.6767761515590734`*^9, 3.6767762704466124`*^9, 
   3.6767776022979183`*^9, 3.6767779075256405`*^9, 3.6820664427631607`*^9, 
   3.6822390748317404`*^9}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{{3.6767005074361563`*^9, 3.6767005218273563`*^9}, {
   3.6767005555987563`*^9, 3.6767005800851564`*^9}, 3.6767006835777564`*^9, 
   3.6767008689681134`*^9, 3.6767011759389095`*^9, 3.6767017025795383`*^9, 
   3.676701802558338*^9, {3.676701880396138*^9, 3.676701909702629*^9}, 
   3.6767027308751717`*^9, 3.6767029204661956`*^9, {3.6767030776327934`*^9, 
   3.676703086996373*^9}, {3.676704177502925*^9, 3.676704182682125*^9}, 
   3.676706001011882*^9, 3.6767761515590734`*^9, 3.6767762704466124`*^9, 
   3.6767776022979183`*^9, 3.6767779075256405`*^9, 3.6820664427631607`*^9, 
   3.6822390748337407`*^9}],

Cell[BoxData[
 RowBox[{
  FractionBox[
   SuperscriptBox["a", "4"], "8"], "+", 
  RowBox[{
   SuperscriptBox["a", "3"], " ", 
   RowBox[{"(", 
    RowBox[{
     FractionBox["5", "12"], "+", 
     RowBox[{
      FractionBox["1", "2"], " ", "theta", " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"0", ",", "theta"}], "]"}]}]}], ")"}]}], "+", 
  RowBox[{
   SuperscriptBox["a", "2"], " ", 
   RowBox[{"(", 
    RowBox[{
     FractionBox["5", "12"], "+", 
     RowBox[{
      FractionBox["1", "2"], " ", "theta", " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"0", ",", "theta"}], "]"}]}], "+", 
     RowBox[{
      FractionBox["3", "4"], " ", 
      SuperscriptBox["theta", "2"], " ", 
      SuperscriptBox[
       RowBox[{"PolyGamma", "[", 
        RowBox[{"0", ",", "theta"}], "]"}], "2"]}], "-", 
     RowBox[{
      FractionBox["3", "4"], " ", 
      SuperscriptBox["theta", "2"], " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"1", ",", "theta"}], "]"}]}]}], ")"}]}], "+", 
  RowBox[{"a", " ", 
   RowBox[{"(", 
    RowBox[{
     FractionBox["1", "8"], "+", 
     RowBox[{
      FractionBox["1", "12"], " ", "theta", " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"0", ",", "theta"}], "]"}]}], "-", 
     RowBox[{
      FractionBox["1", "4"], " ", 
      SuperscriptBox["theta", "2"], " ", 
      SuperscriptBox[
       RowBox[{"PolyGamma", "[", 
        RowBox[{"0", ",", "theta"}], "]"}], "2"]}], "+", 
     RowBox[{
      FractionBox["1", "2"], " ", 
      SuperscriptBox["theta", "3"], " ", 
      SuperscriptBox[
       RowBox[{"PolyGamma", "[", 
        RowBox[{"0", ",", "theta"}], "]"}], "3"]}], "+", 
     RowBox[{
      FractionBox["1", "4"], " ", 
      SuperscriptBox["theta", "2"], " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"1", ",", "theta"}], "]"}]}], "-", 
     RowBox[{
      FractionBox["3", "2"], " ", 
      SuperscriptBox["theta", "3"], " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"0", ",", "theta"}], "]"}], " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"1", ",", "theta"}], "]"}]}], "+", 
     RowBox[{
      FractionBox["1", "2"], " ", 
      SuperscriptBox["theta", "3"], " ", 
      RowBox[{"PolyGamma", "[", 
       RowBox[{"2", ",", "theta"}], "]"}]}]}], ")"}]}]}]], "Output",
 CellChangeTimes->{{3.6767005074361563`*^9, 3.6767005218273563`*^9}, {
   3.6767005555987563`*^9, 3.6767005800851564`*^9}, 3.6767006835777564`*^9, 
   3.6767008689681134`*^9, 3.6767011759389095`*^9, 3.6767017025795383`*^9, 
   3.676701802558338*^9, {3.676701880396138*^9, 3.676701909702629*^9}, 
   3.6767027308751717`*^9, 3.6767029204661956`*^9, {3.6767030776327934`*^9, 
   3.676703086996373*^9}, {3.676704177502925*^9, 3.676704182682125*^9}, 
   3.676706001011882*^9, 3.6767761515590734`*^9, 3.6767762704466124`*^9, 
   3.6767776022979183`*^9, 3.6767779075256405`*^9, 3.6820664427631607`*^9, 
   3.682239074835741*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"Here", " ", "we", " ", "insert", " ", "a"}], " ", "=", " ", 
     RowBox[{"gam", " ", "+", " ", 
      RowBox[{"chi", "[", 
       RowBox[{"[", "1", "]"}], "]"}], "-", 
      RowBox[{"1", "/", "2"}]}]}], ",", " ", 
    RowBox[{
    "where", " ", "gam", " ", "is", " ", "the", " ", "new", " ", 
     "variable"}]}], "*)"}], 
  RowBox[{
   RowBox[{"Poly", "[", "gam_", "]"}], "=", " ", 
   RowBox[{"Simplify", "[", 
    RowBox[{
     RowBox[{"Coefficient", "[", 
      RowBox[{"A", ",", 
       RowBox[{"1", "/", 
        RowBox[{"w", "^", "2"}]}]}], "]"}], " ", "/.", " ", 
     RowBox[{"a", "\[Rule]", 
      RowBox[{"gam", "+", 
       RowBox[{"chi", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "-", 
       RowBox[{"1", "/", "2"}]}]}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.67670206004994*^9, 3.6767020710829277`*^9}, {
   3.676703163839727*^9, 3.6767032422780623`*^9}, 3.6767032979887486`*^9, {
   3.6767779473126407`*^9, 3.676777969643241*^9}, {3.6820658700737495`*^9, 
   3.68206587090055*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "384"], " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"-", "1"}], "+", 
    RowBox[{"2", " ", "gam"}], "-", 
    RowBox[{"2", " ", "theta", " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"0", ",", "theta"}], "]"}]}]}], ")"}], " ", 
  RowBox[{"(", 
   RowBox[{"1", "+", 
    RowBox[{"18", " ", "gam"}], "+", 
    RowBox[{"44", " ", 
     SuperscriptBox["gam", "2"]}], "+", 
    RowBox[{"24", " ", 
     SuperscriptBox["gam", "3"]}], "+", 
    RowBox[{"4", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "19"}], "+", 
       RowBox[{"6", " ", "gam"}]}], ")"}], " ", 
     SuperscriptBox["theta", "2"], " ", 
     SuperscriptBox[
      RowBox[{"PolyGamma", "[", 
       RowBox[{"0", ",", "theta"}], "]"}], "2"]}], "+", 
    RowBox[{"24", " ", 
     SuperscriptBox["theta", "3"], " ", 
     SuperscriptBox[
      RowBox[{"PolyGamma", "[", 
       RowBox[{"0", ",", "theta"}], "]"}], "3"]}], "-", 
    RowBox[{"24", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "5"}], "+", 
       RowBox[{"6", " ", "gam"}]}], ")"}], " ", 
     SuperscriptBox["theta", "2"], " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"1", ",", "theta"}], "]"}]}], "-", 
    RowBox[{"2", " ", "theta", " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"0", ",", "theta"}], "]"}], " ", 
     RowBox[{"(", 
      RowBox[{"13", "+", 
       RowBox[{"44", " ", "gam"}], "-", 
       RowBox[{"12", " ", 
        SuperscriptBox["gam", "2"]}], "+", 
       RowBox[{"72", " ", 
        SuperscriptBox["theta", "2"], " ", 
        RowBox[{"PolyGamma", "[", 
         RowBox[{"1", ",", "theta"}], "]"}]}]}], ")"}]}], "+", 
    RowBox[{"96", " ", 
     SuperscriptBox["theta", "3"], " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"2", ",", "theta"}], "]"}]}]}], ")"}]}]], "Output",
 CellChangeTimes->{{3.6767032428864546`*^9, 3.676703298799938*^9}, 
   3.676704186192125*^9, 3.676705883709882*^9, 3.6767060040538816`*^9, 
   3.6767776084441996`*^9, 3.676777970516841*^9, 3.682066445651142*^9, 
   3.6822390808613434`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", "\[IndentingNewLine]", 
   RowBox[{
   "Now", " ", "we", " ", "can", " ", "compute", " ", "the", " ", "quantity", 
    " ", "denoted", " ", "by", " ", 
    RowBox[{"s", "^"}], "*", 
    RowBox[{"(", "theta", ")"}], " ", "in", " ", "the", " ", "paper", " ", 
    "\"\<Asymptotic expansions for the Ewens distribution and the Stirling \
numbers of the first kind\>\""}], "\[IndentingNewLine]", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"Simplify", "[", 
   RowBox[{
    RowBox[{"Poly", "[", 
     RowBox[{"1", "/", "2"}], "]"}], " ", "-", 
    RowBox[{"Poly", "[", 
     RowBox[{
      RowBox[{"-", "1"}], "/", "2"}], "]"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.6767021062410517`*^9, 3.6767021168236485`*^9}, 
   3.676702743548372*^9, {3.676705724406682*^9, 3.6767057268422823`*^9}, {
   3.6820659162049494`*^9, 3.6820659954277496`*^9}, {3.6820660896888237`*^9, 
   3.682066130578162*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox["1", "2"], " ", 
  SuperscriptBox["theta", "2"], " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"2", " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"1", ",", "theta"}], "]"}]}], "+", 
    RowBox[{"theta", " ", 
     RowBox[{"PolyGamma", "[", 
      RowBox[{"2", ",", "theta"}], "]"}]}]}], ")"}]}]], "Output",
 CellChangeTimes->{{3.676702108583007*^9, 3.6767021173852377`*^9}, 
   3.676702451435372*^9, 3.676702744905572*^9, 3.6767029386785784`*^9, 
   3.6767030960041695`*^9, 3.676703245495621*^9, {3.676703276177828*^9, 
   3.6767033015923023`*^9}, 3.6767033875248013`*^9, 3.676704188910525*^9, 
   3.6767057275286818`*^9, 3.676705893401482*^9, 3.676706005910282*^9, 
   3.6767776097565746`*^9, 3.6767779723732405`*^9, 3.682066448412325*^9, 
   3.6822390846007175`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "The", " ", "resulting", " ", "expression", " ", "is", " ", "positive"}], 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"Poly", "[", 
      RowBox[{"1", "/", "2"}], "]"}], " ", "-", 
     RowBox[{"Poly", "[", 
      RowBox[{
       RowBox[{"-", "1"}], "/", "2"}], "]"}]}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"theta", ",", "0", ",", "1"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.676702750677572*^9, 3.6767027707949715`*^9}, {
  3.676704207515725*^9, 3.676704219623325*^9}, {3.6767043307533245`*^9, 
  3.6767043318453245`*^9}, {3.676705760086282*^9, 3.676705770169882*^9}, {
  3.676705829746282*^9, 3.676705832429482*^9}, {3.676706021911282*^9, 
  3.676706025877882*^9}, {3.676706663199882*^9, 3.676706668619082*^9}, {
  3.6767779805788403`*^9, 3.676777983028041*^9}, {3.6820658974537497`*^9, 
  3.6820659138961496`*^9}, {3.682239089407198*^9, 3.682239090166274*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], LineBox[CompressedData["
1:eJwVxXc81QsDB2CzjLJSsrLdBie5Qsr5cpGi7DIyi0RuSIko2YVE3EP2zogj
DhKVOOk1s4o6J+TK5lfXrCve9/3j+Twy570t3NhYWFge/s//N3Kb6m2aPkNW
Z7XYs5ER/pJDT8xnRFQN0lmmYcdeR5B/z07NqBbVRw71N/6UFwXkcQn50BuK
JnA+2E0omdeRjbsGqzJFrWAW1bxkXksny8dJFm5K2eLuwtp67oEe8rrRhRQX
RUdo2nzx2mYwSJ6JGKC7UFwwamfU+bJohGztuvrioagrik96/NpaPE6m64nV
tWRcxG5W+fQD0jNkFTntyiUpDwTPd25pkPxGzmBzLlHIvwwzu30lbcGLZK6x
sLyzilfA5S7YIaK8Sn4+2tfVYuUDCVJYWhDpX7KPmauRPcUXa1827DvubJI/
H4zWuy96DbsrbOTosRxIzt7dpHDuOjbGRDR6lbfCiL/02MsMf6QGZ+uzrnGj
dqHzMCEViIQrjRNH2QTg5ehYFe1yE/sd2ZOM8gUh202QpPODIPTJ6ZGuzQ7E
lwvtNVe8jQC37b2vfu2CgWRBwbR7CMQid2SrLuzG7dY0gddWd6CZn8Di9o8Y
1NSV7ttSQlE8oLEyclgKOTtNo2JEI0AL5RHjsZJHUFcK/ezzCFjZfrhNS1PA
2chRNrlzkbhHgveuWUVsW/YNaciIQu+pk9zzxfsQMJAcMCd1D5WrQ3Z2qQdh
Gfe59lnTPShwNkuLU1RA0ldcjnCJQZC2s3Cv6yGMV9f5SubHgrOa9/dYgd9h
nvTJ00QxHim31+WWctShZCxXIvY2HtrsyaOfEzTAxe41OeH+ADSepoPBkZpw
qt94+soqAXE6lVOdwVrYuS+BZENJBGjTLUH5ZIRyVf92TzQZwfrf629N6mG/
YeLdKp9kbHN3oQwl6qMv0nua8TYZMQG8y6NkA8izK5WR/P/CxbIWs6Wc42j9
VUAa6KOAQzr61HrESWxbpByWinuExELRUqMOU9Qcuk458fcjmClf4p80NYOj
j+Wqr1YaVj40zvg1mIE6z19Pn0qDxMel/V4Uc1hM3T3maZABFpM3/zafs0Tq
55t6NRtZUM22j/eSs4auhG3B8JlsXC4ndfLmWWPGToOTqzwbhg46azMyNtgz
F3ihXDkHi+Fr9hcVbOHtoa7AwsxBIne2NvnoOQi5UosLNfPA8+6v/k+5TrCx
zaV++14Ai7L+LQJlriAYVjx5YoXg7ZL888iaK6Icudws9QrhLEoS/mrghuoL
3mK1yYXgG7QTbx5zA98V7cggjSI0hiUNxMq6403YR7sttx5DXOdnjvxzD6g+
EdwiwVWK9NEq72JDb7SR3jh3qZQiZG1yV0aBN5yfBjTcti3Fg0DvdQ5WH9yv
HfH5UlqK6celz/c5+WCqqZxRdLoMEQeunTTe44us9yefHkp6goi5J9/6H18F
72aoveEeKvLafHr3Tl5Hqoe1o9JxKrJkPbzXxf0hP6DkLHiFijiZHapbzf1B
Lhm8wHhBxdzFxGevG/xRX0KpqSNVgvTnYN/bpBuYt5xrNO6uhGboracsJoE4
U5zS4be9Cv7/OJdengmGgvnCVEssDXRjaSO6YxjEQgy13HNp0NbQXfUPC4NA
eU4sTx0NrqynLY8XheEnlyXJfIwGIUFThtZCGHqa6vyGNWtgcoQ7PDAkHDdV
QjfWvtYg134Ln87jCHTzCwuTdOqgwzdXE8MbjRvdR5C6VI9d4vcDV+TiIEIT
HGZlfYWO4dnICusk1I7pdb471ozAgIJ89xeP0OyZKW47R0eCQztLrlkOlhW1
+cv/fIt0u2C+Ono+2H6UCYgrtEOj8l1R/MvHUN90rrWs6MRSAAs2C8vQGaHW
3br/Ha6XJyxW5lEhlTD7o2ShBzERXZke/FWYXR0Zfq/Sh5TDtBNbNGkomGlo
54vvB73MbZtHdi34JB8T4+MDuOrfOXyEpx5xHjUXc+Q/wDQ9vH35bANWdMrG
L4UM4qfnrblXoy9wTZVO29UyBD9tZUlz+SaMVNBje2U+QXYz4niU/2v86zmQ
vu8mA1XchYZ5rc04oLt142IDEzoPM9nc2ekwzJNxfyoyjEzt77PfHd7ghqST
7Bf7ERhFTyQIF7RCae4M7+DpUUw9FHFY+voWO/i/cdF/jqJOWqpPz6ANHK/S
qXIZX8Cd7Ukdjm/HZbWGj9zGY+gLNb/RvtyBspPnc16ujOFmRpxKiXUXlmfF
BCYof6M4NqHsa143aiwN2ZP1x7EuXDiRwdmDkBOR2bSJcUiIUKJ0Knog89lo
KD/yK2TtGsK/W/ZCSyMiWVR1Aieynn/ZzdUH+z3vfrD2T2CKq+hDYGUfltvv
j/mETqLYY/mOoXs/Ujkcr1nIT8HRypL/kNAAdgzX5LHRp1A6PrB3rn4AnLdC
rY74TkO05rUoh897sJQ4UseEZrBFUzBop9gHXE2j3hdqnkGrxYJrUN8HvDhw
l0vwyiwouiytqdGDEDS6qhTGN4eMxZiEVJUhMJV7+gMb5/AjV8NwdHIIlVmQ
H3ecxxVtzqssSR/hQFUXlNyYx+42tXctWp9wJzI3lFm6gDmHmiD6t08on5JM
5ztFgDVY4j9JaQzczYx0HjYhUCm5uXgvkwFXi3mFCnMC/YN6UndyGJBobKw8
bU1A1f+bn1cRAzEPzrXGnScwnNy006CKgUvqj77zBBLwjXygtdLGgGyE8Imt
RQR+ulLYbX4wsKEZvH2wmMCQt46YyToDH+f/7isqI5Btp6aiv8lAgnW1w/Gn
BDTHVmxVOJnY3G/hF9lI4A+7vgIuQSaYvQ+y2PsJiGzVFanfy8SzqNUL/e8J
MMVUZakHmEg+6rQvf4hASrSLUiGJCeNCEu2PYQLPKAHkRDUmngd0tYVOE3B7
72frocMERfnwA/M5AtRFPUcnPSZ8xzKsZAgCFR2nzp85zsTeU14jTUsEREc4
L+meYoKDdaAgYZXA9mc0Dw1TJkZrjno6/yTQE/7IU9mCiUbP/IMqvwh0JlZ7
yp1hIkWKd3lzk0DAITZPURsm/gsI911a
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  ImageSize->{808., Automatic},
  Method->{},
  PlotRange->{{0, 1}, {-1.461301188145163*^-7, 0.4428771529426494}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], 
    Scaled[0.02]}]], "Output",
 CellChangeTimes->{{3.676702759181572*^9, 3.6767027712161717`*^9}, 
   3.676702943487163*^9, {3.676703256859276*^9, 3.676703303589077*^9}, 
   3.67670339155355*^9, {3.6767041915041246`*^9, 3.676704219966525*^9}, 
   3.676704332750125*^9, 3.676705562113282*^9, 3.676705736130282*^9, 
   3.676705770778282*^9, 3.676705832944282*^9, 3.676705912720682*^9, {
   3.676706008830482*^9, 3.676706052469282*^9}, {3.6767066646838818`*^9, 
   3.676706669087082*^9}, 3.6767776125177217`*^9, {3.6767779747912407`*^9, 
   3.6767779852140408`*^9}, {3.6822390873919964`*^9, 3.682239090675325*^9}}]
}, Open  ]]
},
WindowSize->{1898, 855},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
PrintingCopies->1,
PrintingPageRange->{32000, 32000},
PrintingOptions->{"Magnification"->1.,
"PaperOrientation"->"Portrait",
"PaperSize"->{595.1999999999999, 841.8}},
Magnification->2.,
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 122, 2, 60, "Input"],
Cell[682, 24, 1290, 26, 251, "Input"],
Cell[CellGroupData[{
Cell[1997, 54, 3514, 93, 403, "Input"],
Cell[5514, 149, 744, 11, 59, "Output"],
Cell[6261, 162, 840, 15, 65, "Output"],
Cell[7104, 179, 1524, 37, 106, "Output"]
}, Open  ]],
Cell[8643, 219, 482, 11, 60, "Input"],
Cell[CellGroupData[{
Cell[9150, 234, 1209, 36, 137, "Input"],
Cell[10362, 272, 614, 13, 85, "Output"],
Cell[10979, 287, 863, 22, 85, "Output"],
Cell[11845, 311, 1013, 26, 85, "Output"]
}, Open  ]],
Cell[12873, 340, 3398, 89, 593, "Input"],
Cell[16274, 431, 146, 2, 60, "Input"],
Cell[CellGroupData[{
Cell[16445, 437, 1446, 37, 99, "Input"],
Cell[17894, 476, 756, 16, 106, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[18687, 497, 477, 13, 99, "Input"],
Cell[19167, 512, 251, 7, 106, "Output"],
Cell[19421, 521, 116, 2, 106, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19574, 528, 1107, 28, 60, "Input"],
Cell[20684, 558, 1546, 41, 149, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22267, 604, 483, 12, 137, "Input"],
Cell[22753, 618, 542, 16, 106, "Output"],
Cell[23298, 636, 643, 18, 106, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23978, 659, 1252, 30, 60, "Input"],
Cell[25233, 691, 2894, 77, 233, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[28164, 773, 293, 7, 99, "Input"],
Cell[28460, 782, 1134, 33, 149, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29631, 820, 1489, 40, 251, "Input"],
Cell[31123, 862, 620, 18, 110, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[31780, 885, 2762, 74, 289, "Input"],
Cell[34545, 961, 1335, 33, 110, "Output"],
Cell[35883, 996, 652, 9, 84, "Output"],
Cell[36538, 1007, 652, 9, 84, "Output"],
Cell[37193, 1018, 2887, 78, 288, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[40117, 1101, 1088, 29, 99, "Input"],
Cell[41208, 1132, 2063, 58, 227, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[43308, 1195, 936, 21, 213, "Input"],
Cell[44247, 1218, 812, 18, 106, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[45096, 1241, 991, 22, 99, "Input"],
Cell[46090, 1265, 4102, 73, 1062, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)