This is documentation for Mathematica 3, which was
based on an earlier version of the Wolfram Language.
View current documentation (Version 11.2)
 Documentation / Mathematica / The Mathematica Book / Advanced Mathematics / Series, Limits and Residues  /

3.6.7 Summation of Series


Evaluating sums.




  • Mathematica recognizes this as the power series expansion of


    .
  • In[1]:= Sum[x^n/n!, {n, 0, Infinity}]

    Out[1]=

  • This sum comes out in terms of a Bessel function.
  • In[2]:= Sum[x^n/(n!^2), {n, 0, Infinity}]

    Out[2]=

  • Here is another sum that can be done in terms of common special functions.
  • In[3]:= Sum[n! x^n/(2n)!, {n, 1, Infinity}]

    Out[3]=

  • Generalized hypergeometric functions are not uncommon in sums.
  • In[4]:= Sum[x^n/(n!^4), {n, 0, Infinity}]

    Out[4]=

    There are many analogies between sums and integrals. And just as it is possible to have indefinite integrals, so indefinite sums can be set up by using symbolic variables as upper limits.

  • This is effectively an indefinite sum.
  • In[5]:= Sum[k, {k, 0, n}]

    Out[5]=

  • This sum comes out in terms of incomplete gamma functions.
  • In[6]:= Sum[x^k/k!, {k, 0, n}]

    Out[6]=

  • This sum involves polygamma functions.
  • In[7]:= Sum[1/k^4, {k, 1, n}]

    Out[7]=




  • Taking the difference between results for successive values of


    gives back the original summand.
  • In[8]:= FullSimplify[ % - (% /. n->n-1) ]

    Out[8]=

    Mathematica can do essentially all sums that are found in books of tables. Just as with indefinite integrals, indefinite sums of expressions involving simple functions tend to give answers that involve more complicated functions. Definite sums, like definite integrals, often, however, come out in terms of simpler functions.

  • This indefinite sum gives a quite complicated result.
  • In[9]:= Sum[Binomial[2k, k]/3^(2k), {k, 0, n}]

    Out[9]=

  • The definite form is much simpler.
  • In[10]:= Sum[Binomial[2k, k]/3^(2k), {k, 0, Infinity}]

    Out[10]=

  • Here is a slightly more complicated definite sum.
  • In[11]:= Sum[PolyGamma[k]/k^2, {k, 1, Infinity}]

    Out[11]=