(***********************************************************************

                    Mathematica-Compatible Notebook

This notebook can be used on any computer system with Mathematica 3.0,
MathReader 3.0, or any compatible application. The data for the notebook 
starts with the line containing stars above.

To get the notebook into a Mathematica-compatible application, do one of 
the following:

* Save the data starting with the line of stars above into a file
  with a name ending in .nb, then open the file inside the application;

* Copy the data starting with the line of stars above to the
  clipboard, then use the Paste menu command inside the application.

Data for notebooks contains only printable 7-bit ASCII and can be
sent directly in email or through ftp in text mode.  Newlines can be
CR, LF or CRLF (Unix, Macintosh or MS-DOS style).

NOTE: If you modify the data for this notebook not in a Mathematica-
compatible application, you must delete the line below containing the 
word CacheID, otherwise Mathematica-compatible applications may try to 
use invalid cache data.

For more information on notebooks and Mathematica-compatible 
applications, contact Wolfram Research:
  web: http://www.wolfram.com
  email: info@wolfram.com
  phone: +1-217-398-0700 (U.S.)

Notebook reader applications are available free of charge from 
Wolfram Research.
***********************************************************************)

(*CacheID: 232*)


(*NotebookFileLineBreakTest
NotebookFileLineBreakTest*)
(*NotebookOptionsPosition[     21127,        719]*)
(*NotebookOutlinePosition[     21889,        746]*)
(*  CellTagsIndexPosition[     21845,        742]*)
(*WindowFrame->Normal*)



Notebook[{

Cell[CellGroupData[{
Cell["The solution of the problem 96-1 from the SIAM Review", "Title",
  TextAlignment->Center,
  FontSize->24],

Cell["by Victor Adamchik", "Subsubtitle",
  TextAlignment->Center],

Cell[CellGroupData[{

Cell["The problem 96-1.", "Section"],

Cell[TextData[{
  "Define, for integer ",
  Cell[BoxData[
      \(TraditionalForm\`n \[GreaterEqual] 2\)]],
  ", the integral"
}], "Text"],

Cell[BoxData[
    \(TraditionalForm
    \`I\_n = \[Integral]\_0\%1
          \( x\^\(n - 1\)\/\((\((1 - x)\)\^n + x\^n)\)\^2\) 
          \[DifferentialD]x\)], "Text",
  TextAlignment->Center],

Cell[TextData[{
  "(1) Evaluate ",
  Cell[BoxData[
      \(TraditionalForm\`I\_n\)]],
  "as a finite sum of elementary functions."
}], "Text"],

Cell[TextData[{
  "(2) How does ",
  Cell[BoxData[
      \(TraditionalForm\`\(I\_n\ \)\)]],
  "behave for ",
  StyleBox["n",
    FontSlant->"Italic"],
  " large?"
}], "Text"],

Cell[TextData[{
  "\tThe problem arose in the study of the expected number of edges in a \
random sphere-of-influence graph. A random sphere-of-influence graph is \
created as follows. Begin by throwing ",
  StyleBox["m",
    FontSlant->"Italic"],
  " points randomly onto the unit cube in ",
  StyleBox["n",
    FontSlant->"Italic"],
  " dimentions. For each point ",
  Cell[BoxData[
      \(TraditionalForm\`x\_i\)]],
  " create a sphere of influence by drawing an open ball about ",
  Cell[BoxData[
      \(TraditionalForm\`x\_i\)]],
  " with radius equal to the distance to ",
  Cell[BoxData[
      \(TraditionalForm\`x\_i\)]],
  "'s nearest neighbor. Then draw an edge between two points if their spheres \
intersect. The integral ",
  Cell[BoxData[
      \(TraditionalForm\`I\_n\)]],
  " arose in calculating the probability that there exists an edge between \
two points."
}], "Text"]
}, Closed]],

Cell[CellGroupData[{

Cell["The solution of 1).", "Section"],

Cell[TextData[{
  "Setting ",
  Cell[BoxData[
      \(TraditionalForm\`x = z\/\(z + 1\)\)]],
  ", we find"
}], "Text"],

Cell[BoxData[
    \(TraditionalForm
    \`I\_n = \(\[Integral]\_0\%1
            \( x\^\(n - 1\)\/\((\((1 - x)\)\^n + x\^n)\)\^2\) 
            \[DifferentialD]x\  = \ 
        \[Integral]\_0
            \%\[Infinity]\(\(\( z\^\(n - 1\)\)(z + 1)\)\^\(n - 1\)\/\((
                    z\^n\  + \ 1)\)\^2\) \[DifferentialD]x\)\)], "Text",
  TextAlignment->Center],

Cell[TextData[{
  "Expanding ",
  Cell[BoxData[
      \(TraditionalForm\`\((z + 1)\)\^\(n - 1\)\)]],
  " into the series"
}], "Text"],

Cell[BoxData[
    FormBox[
      RowBox[{\(\((z + 1)\)\^\(n - 1\)\), "=", 
        RowBox[{\(\[Sum]\+\(k = 0\)\%\(n - 1\)\), 
          RowBox[{
            TagBox[
              RowBox[{"(", GridBox[{
                    {\(n - 1\)},
                    {"k"}
                    }], ")"}],
              Binomial,
              Editable->False], " ", \(z\^k\)}]}]}], TraditionalForm]], "Text",\

  TextAlignment->Center],

Cell["and changing the order of summation and integration, we get", "Text"],

Cell[BoxData[
    FormBox[
      RowBox[{\(I\_n\), " ", "=", " ", 
        RowBox[{\(\[Sum]\+\(k = 0\)\%\(n - 1\)\), 
          RowBox[{
            TagBox[
              RowBox[{"(", GridBox[{
                    {\(n - 1\)},
                    {"k"}
                    }], ")"}],
              Binomial,
              Editable->False], " ", 
            \(\[Integral]\_0
                \%\[Infinity]\( z\^\(\ k\  + \ n - 1\)\/\((z\^n\  + \ 1)\)\^2
                      \) \[DifferentialD]z\)}]}]}], TraditionalForm]], "Text",\

  TextAlignment->Center],

Cell["Evaluating the inner integral ", "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(\[Integral]\_0\%\[Infinity]\( z\^\(k + n - 1\)\/\((z\^n + 1)\)\^2\) 
        \[DifferentialD]z\)], "Input"],

Cell[BoxData[
    FormBox[
      RowBox[{"If", "[", 
        RowBox[{
        \(Re(n) > 0 \[And] Re(n - k) > 0 \[And] Re(k + n) > 0\), ",", 
          \(\(k\ \[Pi]\ \(csc(\(k\ \[Pi]\)\/n)\)\)\/n\^2\), ",", 
          RowBox[{
            SubsuperscriptBox["\[Integral]", "0", 
              InterpretationBox["\[Infinity]",
                DirectedInfinity[ 1]]], 
            \(\(z\^\(k + n - 1\)\/\((z\^n + 1)\)\^2\) 
              \[DifferentialD]z\)}]}], "]"}], TraditionalForm]], "Output"]
}, Open  ]],

Cell["we finally obtain ", "Text"],

Cell[BoxData[
    FormBox[
      RowBox[{\(I\_n\), "=", 
        RowBox[{\(1\/n\), "+", 
          RowBox[{\(\(\[Pi]\ \)\/n\^2\), 
            RowBox[{\(\[Sum]\+\(k = 1\)\%\(n - 1\)\), 
              RowBox[{
                TagBox[
                  RowBox[{"(", GridBox[{
                        {\(n - 1\)},
                        {"k"}
                        }], ")"}],
                  Binomial,
                  Editable->False], " ", "k", " ", 
                \(csc(\(\[Pi]\ k\)\/n)\)}]}]}]}]}], TraditionalForm]], "Text",\

  TextAlignment->Center]
}, Closed]],

Cell[CellGroupData[{

Cell["The solution of 2).", "Section"],

Cell[TextData[{
  "\tThe asymptotic expansion of the given integral for ",
  StyleBox["n",
    FontSlant->"Italic"],
  " large can be derived by using the above finite sum and the power series \
of ",
  Cell[BoxData[
      \(TraditionalForm\`csc(z)\)]],
  ":"
}], "Text"],

Cell[BoxData[
    FormBox[
      RowBox[{\(\(\[Pi]\ k\)\/n\ \(csc(\(\[Pi]\ k\)\/n)\)\), " ", "=", " ", 
        RowBox[{"1", " ", "+", 
          RowBox[{\(\[Sum]\+\(j = 1\)\%\[Infinity]\), 
            RowBox[{\((2 - 2\^\(2 - 2\ j\))\), \(k\^\(2\ j\)\/n\^\(2\ j\)\), 
              TagBox[
                RowBox[{"\[Zeta]", "(", 
                  TagBox[\(2\ j\),
                    (Editable -> True)], ")"}],
                InterpretTemplate[ Zeta[ #]&]]}]}]}]}], TraditionalForm]], 
  "Text",
  TextAlignment->Center],

Cell[TextData[{
  "under the assumption ",
  Cell[BoxData[
      \(TraditionalForm\`\[LeftBracketingBar]k\[RightBracketingBar] < n\)]],
  "."
}], "Text"],

Cell[TextData[{
  "Substituting this series into the finite sum for ",
  Cell[BoxData[
      FormBox[
        FormBox[\(I\_n\),
          "TraditionalForm"], TraditionalForm]]],
  " "
}], "Text"],

Cell[BoxData[
    FormBox[
      RowBox[{\(I\_n\), "=", 
        RowBox[{\(1\/n\), "+", 
          RowBox[{\(1\/n\), 
            RowBox[{\(\[Sum]\+\(k = 1\)\%\(n - 1\)\), 
              RowBox[{
                TagBox[
                  RowBox[{"(", GridBox[{
                        {\(n - 1\)},
                        {"k"}
                        }], ")"}],
                  Binomial,
                  Editable->False], " ", 
                RowBox[{"(", 
                  RowBox[{"1", " ", "+", 
                    RowBox[{\(\[Sum]\+\(j = 1\)\%\[Infinity]\), 
                      RowBox[{
                      \((2 - 2\^\(2 - 2\ j\))\), 
                        \(k\^\(2\ j\)\/n\^\(2\ j\)\), 
                        TagBox[
                          RowBox[{"\[Zeta]", "(", 
                            TagBox[\(2\ j\),
                              (Editable -> True)], ")"}],
                          InterpretTemplate[ Zeta[ #]&]]}]}]}], 
                  ")"}]}]}]}]}]}], TraditionalForm]], "Text",
  TextAlignment->Center],

Cell["and changing the order of summation, we find that", "Text"],

Cell[BoxData[
    FormBox[
      RowBox[{\(I\_n\), "=", 
        RowBox[{\(1\/n\), "+", 
          RowBox[{\(1\/n\), 
            RowBox[{\(\[Sum]\+\(k = 1\)\%\(n - 1\)\), 
              TagBox[
                RowBox[{"(", GridBox[{
                      {\(n - 1\)},
                      {"k"}
                      }], ")"}],
                Binomial,
                Editable->False]}]}], "  ", "+", " ", 
          RowBox[{\(1\/n\), 
            RowBox[{\(\[Sum]\+\(j = 1\)\%\[Infinity]\), 
              RowBox[{\(\((2 - 2\^\(2 - 2\ j\))\)\/n\^\(2\ j\)\), 
                TagBox[
                  RowBox[{"\[Zeta]", "(", 
                    TagBox[\(2\ j\),
                      (Editable -> True)], ")"}],
                  InterpretTemplate[ Zeta[ #]&]], 
                RowBox[{"(", 
                  RowBox[{\(\[Sum]\+\(k = 1\)\%\(n - 1\)\), 
                    RowBox[{
                      TagBox[
                        RowBox[{"(", GridBox[{
                              {\(n - 1\)},
                              {"k"}
                              }], ")"}],
                        Binomial,
                        Editable->False], \(k\^\(2  j\)\)}]}], 
                  ")"}]}]}]}]}]}], TraditionalForm]], "Text",
  TextAlignment->Center],

Cell["Taking into account that", "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(1\/n + \[Sum]\+\(k = 1\)\%\(n - 1\)Binomial[n - 1, k]\/n\)], "Input"],

Cell[BoxData[
    \(TraditionalForm\`\(\(-2\) + 2\^n\)\/\(2\ n\) + 1\/n\)], "Output"]
}, Open  ]],

Cell["we arrive at", "Text"],

Cell[BoxData[
    FormBox[
      RowBox[{\(I\_n\), "=", 
        RowBox[{\(2\^\(n - 1\)\/n\), " ", "+", " ", 
          RowBox[{\(1\/n\), 
            RowBox[{\(\[Sum]\+\(j = 1\)\%\[Infinity]\), 
              RowBox[{\(\((2 - 2\^\(2 - 2\ j\))\)\/n\^\(2\ j\)\), 
                TagBox[
                  RowBox[{"\[Zeta]", "(", 
                    TagBox[\(2\ j\),
                      (Editable -> True)], ")"}],
                  InterpretTemplate[ Zeta[ #]&]], 
                RowBox[{"(", 
                  RowBox[{\(\[Sum]\+\(k = 1\)\%\(n - 1\)\), 
                    RowBox[{
                      TagBox[
                        RowBox[{"(", GridBox[{
                              {\(n - 1\)},
                              {"k"}
                              }], ")"}],
                        Binomial,
                        Editable->False], \(k\^\(2  j\)\)}]}], 
                  ")"}]}]}]}]}]}], TraditionalForm]], "Text",
  TextAlignment->Center],

Cell[TextData[{
  "To find the asymptotic expansion of  ",
  Cell[BoxData[
      \(TraditionalForm\`\(I\_n, \)\)]],
  " when ",
  Cell[BoxData[
      \(TraditionalForm\`n \[Rule] \[Infinity]\)]],
  " we need to find an asymptotic expantion of the inner sum. First of all we \
observe that"
}], "Text"],

Cell[BoxData[
    FormBox[
      RowBox[{
        RowBox[{\(\[Sum]\+\(k = 1\)\%\(n - 1\)\), 
          RowBox[{
            TagBox[
              RowBox[{"(", GridBox[{
                    {\(n - 1\)},
                    {"k"}
                    }], ")"}],
              Binomial,
              Editable->False], \(k\^\(2  j\)\)}]}], " ", "=", 
        \(lim\+\(z \[Rule] 1\)\[ThinSpace]\((
            \(\((z  d\/dz)\)\^\(2  j\)\) \((z + 1)\)\^\(n - 1\))\)\)}], 
      TraditionalForm]], "Text",
  TextAlignment->Center],

Cell[TextData[{
  "Now let us find several values of this sum when ",
  StyleBox["j",
    FontSlant->"Italic"],
  " is changing form 1 to 4"
}], "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(Table[
        Factor[Nest[z*D[#1, \ z]\  &\ , \ \((z\  + \ 1)\)^\((n\  - \ 1)\), \ 
              \n\ \ \ \ 2*j] /. \ z\  -> \ 1], \ {j, \ 1, \ 4}] // TableForm
      \)], "Input"],

Cell[BoxData[
    FormBox[
      InterpretationBox[GridBox[{
            {\(2\^\(n - 3\)\ \((n - 1)\)\ n\)},
            {\(2\^\(n - 5\)\ \((n - 1)\)\ n\ \((n\^2 + 3\ n - 6)\)\)},
            {
              \(2\^\(n - 7\)\ \((n - 1)\)\ n\ 
                \((n\^4 + 10\ n\^3 - 5\ n\^2 - 70\ n + 80)\)\)},
            {
              \(2\^\(n - 9\)\ \((n - 1)\)\ n\ 
                \((n\^6 + 21\ n\^5 + 63\ n\^4 - 385\ n\^3 - 280\ n\^2 + 
                    2212\ n - 1904)\)\)}
            },
          RowSpacings->1,
          ColumnSpacings->3,
          RowAlignments->Baseline,
          ColumnAlignments->{Left}],
        TableForm[ {
          Times[ 
            Power[ 2, 
              Plus[ -3, n]], 
            Plus[ -1, n], n], 
          Times[ 
            Power[ 2, 
              Plus[ -5, n]], 
            Plus[ -1, n], n, 
            Plus[ -6, 
              Times[ 3, n], 
              Power[ n, 2]]], 
          Times[ 
            Power[ 2, 
              Plus[ -7, n]], 
            Plus[ -1, n], n, 
            Plus[ 80, 
              Times[ -70, n], 
              Times[ -5, 
                Power[ n, 2]], 
              Times[ 10, 
                Power[ n, 3]], 
              Power[ n, 4]]], 
          Times[ 
            Power[ 2, 
              Plus[ -9, n]], 
            Plus[ -1, n], n, 
            Plus[ -1904, 
              Times[ 2212, n], 
              Times[ -280, 
                Power[ n, 2]], 
              Times[ -385, 
                Power[ n, 3]], 
              Times[ 63, 
                Power[ n, 4]], 
              Times[ 21, 
                Power[ n, 5]], 
              Power[ n, 6]]]}]], TraditionalForm]], "Output"]
}, Open  ]],

Cell["\<\
It is easy to see that elements of this list can be written as\
\>",
   "Text"],

Cell[BoxData[
    \(TraditionalForm
    \`2\^\(n\  - \ 2  j - 1\)\ 
      \(\(n\^\(2\ j\)\)(1\  + e\_1\/n + e\_2\/n\^2 + \  ... )\)\)], "Text",
  TextAlignment->Center],

Cell[TextData[{
  "where ",
  Cell[BoxData[
      \(TraditionalForm\`e\_p\)]],
  " are coefficients depend on ",
  StyleBox["j ",
    FontSlant->"Italic"],
  "only. If we divide the above list by the common multiplicator we obtain"
}], "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(TableForm[
      Table[Collect[
            \((Nest[z*D[#1, \ z]\  &\ , \ \((z\  + \ 1)\)^\((n\  - \ 1)\), \ 
                    \n\ \ \ \ 2*j] /. \ z\  -> \ 1)\)\ \ 
              1/\((2^\((n - 2  j - 1)\)\ n^\((2  j)\))\), \ n], 
          \ {j, \ 1, \ 5}] /. \ n^s_\  :> \ 0\  /; \ s\  < \(-2\)]\)], "Input"],

Cell[BoxData[
    FormBox[
      InterpretationBox[GridBox[{
            {\(1 - 1\/n\)},
            {\(1 + 2\/n - 9\/n\^2\)},
            {\(1 + 9\/n - 15\/n\^2\)},
            {\(1 + 20\/n + 42\/n\^2\)},
            {\(1 + 35\/n + 270\/n\^2\)}
            },
          RowSpacings->1,
          ColumnSpacings->3,
          RowAlignments->Baseline,
          ColumnAlignments->{Left}],
        TableForm[ {
          Plus[ 1, 
            Times[ -1, 
              Power[ n, -1]]], 
          Plus[ 1, 
            Times[ -9, 
              Power[ n, -2]], 
            Times[ 2, 
              Power[ n, -1]]], 
          Plus[ 1, 
            Times[ -15, 
              Power[ n, -2]], 
            Times[ 9, 
              Power[ n, -1]]], 
          Plus[ 1, 
            Times[ 42, 
              Power[ n, -2]], 
            Times[ 20, 
              Power[ n, -1]]], 
          Plus[ 1, 
            Times[ 270, 
              Power[ n, -2]], 
            Times[ 35, 
              Power[ n, -1]]]}]], TraditionalForm]], "Output"]
}, Open  ]],

Cell["We can guess that", "Text"],

Cell[BoxData[
    \(e\_1[j_] := a\_1\ j\^2 + a\_2\ j + a\_3\)], "Input"],

Cell[BoxData[
    \(e\_2[j_] := b\_1\ j\^4 + b\_2\ j\^3 + b\_3\ j\^2 + b\_4\ j + b\_5\)], 
  "Input"],

Cell[TextData[{
  "when ",
  Cell[BoxData[
      \(TraditionalForm\`j\  = 1, \ 2, \ \( ...  . \)\)]],
  "  Then from the second raw it follows"
}], "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(Solve[{e\_1[1] == \(-1\), e\_1[2] == 2, e\_1[3] == 9}, {a\_1, a\_2, 
        a\_3}]\)], "Input"],

Cell[BoxData[
    \(TraditionalForm
    \`{{a\_1 \[Rule] 2, a\_2 \[Rule] \(-3\), a\_3 \[Rule] 0}}\)], "Output"]
}, Open  ]],

Cell["Thus", "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(Factor[e\_1[j] /. %\[LeftDoubleBracket]1\[RightDoubleBracket]]\)], 
  "Input"],

Cell[BoxData[
    \(TraditionalForm\`j\ \((2\ j - 3)\)\)], "Output"]
}, Open  ]],

Cell[TextData[{
  "Similar for ",
  Cell[BoxData[
      \(TraditionalForm\`e\_\(2 : \)\)]]
}], "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(Solve[{e\_2[1] == 0, e\_2[2] == \(-9\), e\_2[3] == \(-15\), 
        e\_2[4] == 42, e\_2[5] == 270}, {b\_1, b\_2, b\_3, b\_4, b\_5}]\)], 
  "Input"],

Cell[BoxData[
    \(TraditionalForm
    \`{{b\_1 \[Rule] 2, b\_2 \[Rule] \(-10\), b\_3 \[Rule] 23\/2, 
        b\_4 \[Rule] \(-\(7\/2\)\), b\_5 \[Rule] 0}}\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(Factor[e\_2[j] /. %\[LeftDoubleBracket]1\[RightDoubleBracket]]\)], 
  "Input"],

Cell[BoxData[
    \(TraditionalForm
    \`1\/2\ \((j - 1)\)\ j\ \((2\ j - 7)\)\ \((2\ j - 1)\)\)], "Output"]
}, Open  ]],

Cell["\<\
Finally, we derive the following asymptoric expansion for the inner \
sum\
\>", "Text"],

Cell[BoxData[
    FormBox[
      RowBox[{
        RowBox[{\(\[Sum]\+\(k = 1\)\%\(n - 1\)\), 
          RowBox[{
            TagBox[
              RowBox[{"(", GridBox[{
                    {\(n - 1\)},
                    {"k"}
                    }], ")"}],
              Binomial,
              Editable->False], \(k\^\(2  j\)\)}]}], " ", "\[TildeEqual]", 
        \(2\^\(n\  - \ 2  j - 1\)\ 
          \(\(n\^\(2\ j\)\)(
            1\  + \(j\ \((2\ j - 3)\)\)\/n + 
                \(\((j - 1)\)\ j\ \((2\ j - 7)\)\ 
                    \((2\ j - 1)\)\)\/\(2  n\^2\) + \  ... )\)\)}], 
      TraditionalForm]], "Text",
  TextAlignment->Center],

Cell[TextData[{
  "Thus, the integral ",
  Cell[BoxData[
      \(TraditionalForm\`I\_n\)]],
  " can be rewritten as"
}], "Text"],

Cell[BoxData[
    FormBox[
      RowBox[{\(I\_n\), "\[TildeEqual]", 
        RowBox[{\(2\^\(n - 1\)\/n\), " ", "+", " ", 
          RowBox[{\(1\/n\), 
            RowBox[{\(\[Sum]\+\(j = 1\)\%\[Infinity]\), 
              RowBox[{\(\((2 - 2\^\(2 - 2\ j\))\)\/n\^\(2\ j\)\), 
                TagBox[
                  RowBox[{
                    RowBox[{"\[Zeta]", "(", 
                      TagBox[\(2\ j\),
                        (Editable -> True)], ")"}], 
                    \(2\^\(n\  - \ 2  j - 1\)\), " ", \(n\^\(2\ j\)\)}],
                  InterpretTemplate[ Zeta[ #]&]], 
                \((1\  + \(j\ \((2\ j - 3)\)\)\/n + 
                    \(\((j - 1)\)\ j\ \((2\ j - 7)\)\ 
                        \((2\ j - 1)\)\)\/\(2  n\^2\))\)}]}]}]}]}], 
      TraditionalForm]], "Text",
  TextAlignment->Center],

Cell["The infinite sum is doable:", "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(2\^\(n - 1\)\/n + 
      \(\(2\^\(n - 1\)\ \)\/n\) 
        \(\[Sum]\+\(j = 1\)\%\[Infinity]\((
            \((2 - 2\^\(2 - 2\ j\))\)\/2\^\(2\ j\)\ Zeta[2\ j]\ 
              \((1 + \(j\ \((2\ j - 3)\)\)\/n + 
                  \(\((j - 1)\)\ j\ \((2\ j - 7)\)\ \((2\ j - 1)\)\)\/\(2\ 
                      n\^2\))\))\)\)\)], "Input"],

Cell[BoxData[
    \(TraditionalForm
    \`\(2\^\(n - 9\)\ 
          \((128\ \[Pi]\ n\^2 - 256\ n\^2 + 16\ \[Pi]\^3\ n - 128\ \[Pi]\ n + 
              5\ \[Pi]\^5 - 48\ \[Pi]\^3)\)\)\/n\^3 + 2\^\(n - 1\)\/n\)], 
  "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(2^n\ Collect[%/2^n, \ n, \ Factor]\)], "Input"],

Cell[BoxData[
    \(TraditionalForm
    \`2\^n\ \((
        \[Pi]\/\(4\ n\) + \(\[Pi]\ \((\(-8\) + \[Pi]\^2)\)\)\/\(32\ n\^2\) + 
          \(\[Pi]\^3\ \((\(-48\) + 5\ \[Pi]\^2)\)\)\/\(512\ n\^3\))\)\)], 
  "Output"]
}, Open  ]],

Cell[TextData[{
  "This provides first three asymptotic terms of ",
  Cell[BoxData[
      \(TraditionalForm\`I\_n\)]],
  ", when ",
  Cell[BoxData[
      \(TraditionalForm\`n \[Rule] \[Infinity]\)]],
  ":"
}], "Text"],

Cell[BoxData[
    \(TraditionalForm
    \`I\_n \[TildeEqual] 
      2\^n\ \((\[Pi]\/\(4\ n\) + 
            \(\[Pi]\ \((\(-8\) + \[Pi]\^2)\)\)\/\(32\ n\^2\) + 
            \(\[Pi]\^3\ \((\(-48\) + 5\ \[Pi]\^2)\)\)\/\(512\ n\^3\))\)\)], 
  "Text",
  TextAlignment->Center]
}, Closed]]
}, Open  ]]
},
FrontEndVersion->"NeXT 3.1",
ScreenRectangle->{{0, 1053}, {0, 832}},
WindowToolbars->{"RulerBar", "EditBar"},
Evaluator->"aconitum",
WindowSize->{564, 600},
WindowMargins->{{Automatic, 138}, {Automatic, 24}},
Magnification->1,
StyleDefinitions -> "ArticleClassic.nb"
]


(***********************************************************************
Cached data follows.  If you edit this Notebook file directly, not using
Mathematica, you must remove the line containing CacheID at the top of 
the file.  The cache data will then be recreated when you save this file 
from within Mathematica.
***********************************************************************)

(*CellTagsOutline
CellTagsIndex->{}
*)

(*CellTagsIndex
CellTagsIndex->{}
*)

(*NotebookFileOutline
Notebook[{

Cell[CellGroupData[{
Cell[1739, 51, 111, 2, 83, "Title"],
Cell[1853, 55, 66, 1, 51, "Subsubtitle"],

Cell[CellGroupData[{
Cell[1944, 60, 36, 0, 53, "Section"],
Cell[1983, 62, 138, 5, 26, "Text"],
Cell[2124, 69, 192, 5, 42, "Text"],
Cell[2319, 76, 142, 5, 26, "Text"],
Cell[2464, 83, 174, 8, 26, "Text"],
Cell[2641, 93, 891, 24, 127, "Text"]
}, Closed]],

Cell[CellGroupData[{
Cell[3569, 122, 38, 0, 27, "Section"],
Cell[3610, 124, 118, 5, 29, "Text"],
Cell[3731, 131, 361, 8, 43, "Text"],
Cell[4095, 141, 133, 5, 26, "Text"],
Cell[4231, 148, 422, 13, 50, "Text"],
Cell[4656, 163, 75, 0, 26, "Text"],
Cell[4734, 165, 558, 16, 58, "Text"],
Cell[5295, 183, 46, 0, 26, "Text"],

Cell[CellGroupData[{
Cell[5366, 187, 127, 2, 49, "Input"],
Cell[5496, 191, 494, 11, 85, "Output"]
}, Open  ]],
Cell[6005, 205, 34, 0, 26, "Text"],
Cell[6042, 207, 561, 16, 50, "Text"]
}, Closed]],

Cell[CellGroupData[{
Cell[6640, 228, 38, 0, 27, "Section"],
Cell[6681, 230, 271, 9, 43, "Text"],
Cell[6955, 241, 525, 12, 58, "Text"],
Cell[7483, 255, 153, 5, 26, "Text"],
Cell[7639, 262, 195, 7, 26, "Text"],
Cell[7837, 271, 1043, 26, 59, "Text"],
Cell[8883, 299, 65, 0, 26, "Text"],
Cell[8951, 301, 1270, 32, 58, "Text"],
Cell[10224, 335, 40, 0, 26, "Text"],

Cell[CellGroupData[{
Cell[10289, 339, 89, 1, 58, "Input"],
Cell[10381, 342, 85, 1, 41, "Output"]
}, Open  ]],
Cell[10481, 346, 28, 0, 26, "Text"],
Cell[10512, 348, 969, 23, 58, "Text"],
Cell[11484, 373, 301, 9, 43, "Text"],
Cell[11788, 384, 523, 15, 50, "Text"],
Cell[12314, 401, 152, 5, 26, "Text"],

Cell[CellGroupData[{
Cell[12491, 410, 201, 4, 80, "Input"],
Cell[12695, 416, 1688, 54, 86, "Output"]
}, Open  ]],
Cell[14398, 473, 89, 3, 26, "Text"],
Cell[14490, 478, 168, 4, 35, "Text"],
Cell[14661, 484, 243, 8, 43, "Text"],

Cell[CellGroupData[{
Cell[14929, 496, 332, 6, 107, "Input"],
Cell[15264, 504, 1039, 36, 125, "Output"]
}, Open  ]],
Cell[16318, 543, 33, 0, 26, "Text"],
Cell[16354, 545, 72, 1, 28, "Input"],
Cell[16429, 548, 101, 2, 28, "Input"],
Cell[16533, 552, 155, 5, 26, "Text"],

Cell[CellGroupData[{
Cell[16713, 561, 116, 2, 39, "Input"],
Cell[16832, 565, 111, 2, 27, "Output"]
}, Open  ]],
Cell[16958, 570, 20, 0, 26, "Text"],

Cell[CellGroupData[{
Cell[17003, 574, 98, 2, 25, "Input"],
Cell[17104, 578, 68, 1, 27, "Output"]
}, Open  ]],
Cell[17187, 582, 102, 4, 26, "Text"],

Cell[CellGroupData[{
Cell[17314, 590, 168, 3, 53, "Input"],
Cell[17485, 595, 169, 3, 41, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[17691, 603, 98, 2, 25, "Input"],
Cell[17792, 607, 108, 2, 41, "Output"]
}, Open  ]],
Cell[17915, 612, 97, 3, 26, "Text"],
Cell[18015, 617, 647, 18, 83, "Text"],
Cell[18665, 637, 128, 5, 26, "Text"],
Cell[18796, 644, 821, 18, 91, "Text"],
Cell[19620, 664, 43, 0, 26, "Text"],

Cell[CellGroupData[{
Cell[19688, 668, 355, 7, 102, "Input"],
Cell[20046, 677, 224, 5, 43, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[20307, 687, 67, 1, 25, "Input"],
Cell[20377, 690, 216, 5, 44, "Output"]
}, Open  ]],
Cell[20608, 698, 217, 8, 26, "Text"],
Cell[20828, 708, 271, 7, 41, "Text"]
}, Closed]]
}, Open  ]]
}
]
*)




(***********************************************************************
End of Mathematica Notebook file.
***********************************************************************)