In discrete calculus the indefinite sum operator (also known as the antidifference operator), denoted by
or
,[1][2] is the linear operator, inverse of the forward difference operator
. It relates to the forward difference operator as the indefinite integral relates to the derivative. Thus
![{\displaystyle \Delta \sum _{x}f(x)=f(x)\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/af9bd2f07aabb6d4b5367e283bbac352e74607a7)
More explicitly, if
, then
![{\displaystyle F(x+1)-F(x)=f(x)\,.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0cb5ded3c0403fe535a4dfeeba5dc284efc487d3)
If F(x) is a solution of this functional equation for a given f(x), then so is F(x)+C(x) for any periodic function C(x) with period 1. Therefore, each indefinite sum actually represents a family of functions. However, due to the Carlson's theorem, the solution equal to its Newton series expansion is unique up to an additive constant C. This unique solution can be represented by formal power series form of the antidifference operator:
.
Fundamental theorem of discrete calculus
Indefinite sums can be used to calculate definite sums with the formula:[3]
![{\displaystyle \sum _{k=a}^{b}f(k)=\Delta ^{-1}f(b+1)-\Delta ^{-1}f(a)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2cf851a6acd5515931c322c275b27dff83578c17)
Definitions
The Laplace summation formula allows the indefinite sum to be written as the indefinite integral plus correction terms obtained from iterating the difference operator, although it was originally developed for the reverse process of writing an integral as an indefinite sum plus correction terms. As usual with indefinite sums and indefinite integrals, it is valid up to an arbitrary choice of the constant of integration. Using operator algebra avoids cluttering the formula with repeated copies of the function to be operated on:[4]
In this formula, for instance, the term
represents an operator that divides the given function by two. The coefficients
,
, etc., appearing in this formula are the Gregory coefficients, also called Laplace numbers. The coefficient in the term
is[4]
where the numerator
of the left hand side is called a Cauchy number of the first kind, although this name sometimes applies to the Gregory coefficients themselves.[4]
![{\displaystyle \sum _{x}f(x)=\sum _{k=1}^{\infty }{\binom {x}{k}}\Delta ^{k-1}[f]\left(0\right)+C=\sum _{k=1}^{\infty }{\frac {\Delta ^{k-1}[f](0)}{k!}}(x)_{k}+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6758a97991a0410ddfb47366eb53c4d7597ba5b5)
- where
is the falling factorial.
![{\displaystyle \sum _{x}f(x)=\sum _{n=1}^{\infty }{\frac {f^{(n-1)}(0)}{n!}}B_{n}(x)+C\,,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/03e6a5617ea1d050b57bce325707c3cdb619d3ee)
Faulhaber's formula provides that the right-hand side of the equation converges.
If
then[5]
![{\displaystyle \sum _{x}f(x)=\sum _{n=0}^{\infty }\left(f(n)-f(n+x)\right)+C.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a7585ca679958651dc1ccbaca2f8db35514ff7c8)
![{\displaystyle \sum _{x}f(x)=\int _{0}^{x}f(t)dt-{\frac {1}{2}}f(x)+\sum _{k=1}^{\infty }{\frac {B_{2k}}{(2k)!}}f^{(2k-1)}(x)+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/02e7e47071c024fa40b57b2d120122d362b5bb99)
Choice of the constant term
Often the constant C in indefinite sum is fixed from the following condition.
Let
![{\displaystyle F(x)=\sum _{x}f(x)+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0deda6a6fa4324fbc682482e62510245caa8c47e)
Then the constant C is fixed from the condition
![{\displaystyle \int _{0}^{1}F(x)\,dx=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f66616e00de148b547c3933ea6bb5c59bfea5ff4)
or
![{\displaystyle \int _{1}^{2}F(x)\,dx=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dd34418b71544d24ae8f51d500a7d16959bf99d7)
Alternatively, Ramanujan's sum can be used:
![{\displaystyle \sum _{x\geq 1}^{\Re }f(x)=-f(0)-F(0)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a9a1dcb27f63b109dae35714dc64add00b360705)
or at 1
![{\displaystyle \sum _{x\geq 1}^{\Re }f(x)=-F(1)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/48c1c7c5898a76d3304fe46b819dca27524a436f)
respectively[6][7]
Summation by parts
Indefinite summation by parts:
![{\displaystyle \sum _{x}f(x)\Delta g(x)=f(x)g(x)-\sum _{x}(g(x)+\Delta g(x))\Delta f(x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0c77cafdc86be7d2d92749b5d92df5107dd41b27)
![{\displaystyle \sum _{x}f(x)\Delta g(x)+\sum _{x}g(x)\Delta f(x)=f(x)g(x)-\sum _{x}\Delta f(x)\Delta g(x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8a49348a55afcd42b23bf963e609c4c4123b2d4a)
Definite summation by parts:
![{\displaystyle \sum _{i=a}^{b}f(i)\Delta g(i)=f(b+1)g(b+1)-f(a)g(a)-\sum _{i=a}^{b}g(i+1)\Delta f(i)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c099adf7a192250baab400c763c8d11617be797b)
Period rules
If
is a period of function
then
![{\displaystyle \sum _{x}f(Tx)=xf(Tx)+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/af9d9ed55b94ee534e369385cd77981d7facb070)
If
is an antiperiod of function
, that is
then
![{\displaystyle \sum _{x}f(Tx)=-{\frac {1}{2}}f(Tx)+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/484b172257d17f59e3b6785c0dc335b15a2715aa)
Alternative usage
Some authors use the phrase "indefinite sum" to describe a sum in which the numerical value of the upper limit is not given:
![{\displaystyle \sum _{k=1}^{n}f(k).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/11d24165e6016e88647cb18e8a108c545efaefba)
In this case a closed form expression F(k) for the sum is a solution of
![{\displaystyle F(x+1)-F(x)=f(x+1)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c0b38126b64cb3c5bad9bb3d7116cb68ed2ba095)
which is called the telescoping equation.[8] It is the inverse of the backward difference
operator. It is related to the forward antidifference operator using the fundamental theorem of discrete calculus described earlier.
List of indefinite sums
This is a list of indefinite sums of various functions. Not every function has an indefinite sum that can be expressed in terms of elementary functions.
Antidifferences of rational functions
![{\displaystyle \sum _{x}a=ax+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fa80e06a325da9b183621959fbfd5498ed2fc111)
- From which
can be factored out, leaving 1, with the alternative form
. From that, we have: ![{\displaystyle \sum _{x}x^{0}=\ x}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d5a88c65aa257db38af9798582d5ee0d2aa662d2)
- For the sum below, remember
![{\displaystyle x=x^{1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5435a16e253d6087aaea2e5db93af527003d717e)
![{\displaystyle \sum _{x}x={\frac {x(x+1)}{2}}+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cb2d71f024d53360215282ffc9dd153292e02ff2)
- For positive integer exponents Faulhaber's formula can be used. For negative integer exponents,
![{\displaystyle \sum _{x}{\frac {1}{x^{a}}}={\frac {(-1)^{a+1}\psi ^{(a+1)}(x)}{a!}}+C,\,a\in \mathbb {Z} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/95b59e687b741cf4fda7492e0f83d31e64c59587)
- where
is the polygamma function can be used. - More generally,
![{\displaystyle \sum _{x}x^{a}={\begin{cases}-\zeta (-a,x+1)+C_{1},&{\text{if }}a\neq -1\\\psi (x+1)+C_{2},&{\text{if }}a=-1\end{cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/166abfe57a27ed7ee21519d7b033b70ede5d9bf3)
- where
is the Hurwitz zeta function and
is the Digamma function.
and
are constants which would normally be set to
(where
is the Riemann zeta function) and the Euler–Mascheroni constant respectively. By replacing the variable
with
, this becomes the Generalized harmonic number. For the relation between the Hurwitz zeta and Polygamma functions, refer to Balanced polygamma function and Hurwitz zeta function#Special cases and generalizations. - From this, using
, another form can be obtained: ![{\displaystyle \sum _{x}x^{a}=\int _{0}^{x}-a\zeta (1-a,u+1)du+C,{\text{ if }}a\neq -1}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d205b75acccdb9f387d11c151d7a2a45cc99e4fd)
![{\displaystyle \sum _{x}B_{a}(x)=(x-1)B_{a}(x)-{\frac {a}{a+1}}B_{a+1}(x)+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6d5eb05379f4cf70586cd6cd6e070ad51404959c)
Antidifferences of exponential functions
![{\displaystyle \sum _{x}a^{x}={\frac {a^{x}}{a-1}}+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/846bceb5700f4874f23e3670f011acc1eda2c070)
Particularly,
![{\displaystyle \sum _{x}2^{x}=2^{x}+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/de0b5701a9e59b9c5d83184dd4904cbccc9b0b78)
Antidifferences of logarithmic functions
![{\displaystyle \sum _{x}\log _{b}x=\log _{b}(x!)+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/876b33794a692b307249c47b0329d148902c8dc5)
![{\displaystyle \sum _{x}\log _{b}ax=\log _{b}(x!a^{x})+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7e369f3ce3db5ee0bae0abf7472b72cba47ced64)
Antidifferences of hyperbolic functions
![{\displaystyle \sum _{x}\sinh ax={\frac {1}{2}}\operatorname {csch} \left({\frac {a}{2}}\right)\cosh \left({\frac {a}{2}}-ax\right)+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bb6e83f367e3db30902620e13e706f8683cb9014)
![{\displaystyle \sum _{x}\cosh ax={\frac {1}{2}}\operatorname {csch} \left({\frac {a}{2}}\right)\sinh \left(ax-{\frac {a}{2}}\right)+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fd754181303508dd8646c10f67d92d7a51b7be2c)
![{\displaystyle \sum _{x}\tanh ax={\frac {1}{a}}\psi _{e^{a}}\left(x-{\frac {i\pi }{2a}}\right)+{\frac {1}{a}}\psi _{e^{a}}\left(x+{\frac {i\pi }{2a}}\right)-x+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3e33f983ff967154ef0b6e89b1c9a1da6942c2b0)
- where
is the q-digamma function.
Antidifferences of trigonometric functions
![{\displaystyle \sum _{x}\sin ax=-{\frac {1}{2}}\csc \left({\frac {a}{2}}\right)\cos \left({\frac {a}{2}}-ax\right)+C\,,\,\,a\neq 2n\pi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/b61930bedbb0463086a8f755bc8e240a7d0f19b1)
![{\displaystyle \sum _{x}\cos ax={\frac {1}{2}}\csc \left({\frac {a}{2}}\right)\sin \left(ax-{\frac {a}{2}}\right)+C\,,\,\,a\neq 2n\pi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/6a7d462fcfe49562ab28529104e044c405283f7b)
![{\displaystyle \sum _{x}\sin ^{2}ax={\frac {x}{2}}+{\frac {1}{4}}\csc(a)\sin(a-2ax)+C\,\,,\,\,a\neq n\pi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/17bfdb93a8001d5401f0030b9864b9c669574174)
![{\displaystyle \sum _{x}\cos ^{2}ax={\frac {x}{2}}-{\frac {1}{4}}\csc(a)\sin(a-2ax)+C\,\,,\,\,a\neq n\pi }](https://wikimedia.org/api/rest_v1/media/math/render/svg/852c94a4e59fddf7620abad7dc58f7d83e67dc5a)
![{\displaystyle \sum _{x}\tan ax=ix-{\frac {1}{a}}\psi _{e^{2ia}}\left(x-{\frac {\pi }{2a}}\right)+C\,,\,\,a\neq {\frac {n\pi }{2}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3ef1e8e4a4bef38eae48cb062470371ac007dab2)
- where
is the q-digamma function.
![{\displaystyle \sum _{x}\tan x=ix-\psi _{e^{2i}}\left(x+{\frac {\pi }{2}}\right)+C=-\sum _{k=1}^{\infty }\left(\psi \left(k\pi -{\frac {\pi }{2}}+1-x\right)+\psi \left(k\pi -{\frac {\pi }{2}}+x\right)-\psi \left(k\pi -{\frac {\pi }{2}}+1\right)-\psi \left(k\pi -{\frac {\pi }{2}}\right)\right)+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/db12d63cba3fd1ac9b5b15db58c910ec3a7fb354)
![{\displaystyle \sum _{x}\cot ax=-ix-{\frac {i\psi _{e^{2ia}}(x)}{a}}+C\,,\,\,a\neq {\frac {n\pi }{2}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/80f280ceb381a8018b8229d877df17d61de5e562)
![{\displaystyle \sum _{x}\operatorname {sinc} x=\operatorname {sinc} (x-1)\left({\frac {1}{2}}+(x-1)\left(\ln(2)+{\frac {\psi ({\frac {x-1}{2}})+\psi ({\frac {1-x}{2}})}{2}}-{\frac {\psi (x-1)+\psi (1-x)}{2}}\right)\right)+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/031c299fa06afbfa1d7b1ad75e82a08e91838fa0)
- where
is the normalized sinc function.
Antidifferences of inverse hyperbolic functions
![{\displaystyle \sum _{x}\operatorname {artanh} \,ax={\frac {1}{2}}\ln \left({\frac {\Gamma \left(x+{\frac {1}{a}}\right)}{\Gamma \left(x-{\frac {1}{a}}\right)}}\right)+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4c137f6acedb563932deaad2eb983be6e4cf5edb)
Antidifferences of inverse trigonometric functions
![{\displaystyle \sum _{x}\arctan ax={\frac {i}{2}}\ln \left({\frac {\Gamma (x+{\frac {i}{a}})}{\Gamma (x-{\frac {i}{a}})}}\right)+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d5c52cd5e49b7b03e5d84d2c73e6a23db1ef3fae)
Antidifferences of special functions
![{\displaystyle \sum _{x}\psi (x)=(x-1)\psi (x)-x+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2aad9bf36205bcdbb20b22c05542b46b53f06cc7)
![{\displaystyle \sum _{x}\Gamma (x)=(-1)^{x+1}\Gamma (x){\frac {\Gamma (1-x,-1)}{e}}+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b93931ccb98976fe9028cba00af66b9a16fb00fb)
- where
is the incomplete gamma function.
![{\displaystyle \sum _{x}(x)_{a}={\frac {(x)_{a+1}}{a+1}}+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/30270ff9c716c090bc462e87645d080f2fc32b1a)
- where
is the falling factorial.
![{\displaystyle \sum _{x}\operatorname {sexp} _{a}(x)=\ln _{a}{\frac {(\operatorname {sexp} _{a}(x))'}{(\ln a)^{x}}}+C}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6687c7eb4143eb23825e3346586c45141c274eee)
- (see super-exponential function)
See also
References
- ^ Man, Yiu-Kwong (1993), "On computing closed forms for indefinite summations", Journal of Symbolic Computation, 16 (4): 355–376, doi:10.1006/jsco.1993.1053, MR 1263873
- ^ Goldberg, Samuel (1958), Introduction to difference equations, with illustrative examples from economics, psychology, and sociology, Wiley, New York, and Chapman & Hall, London, p. 41, MR 0094249,
If
is a function whose first difference is the function
, then
is called an indefinite sum of
and denoted by ![{\displaystyle \Delta ^{-1}y}](https://wikimedia.org/api/rest_v1/media/math/render/svg/24a14432fb74f8c1908e2354aa4899770e1828ce)
; reprinted by Dover Books, 1986 - ^ "Handbook of discrete and combinatorial mathematics", Kenneth H. Rosen, John G. Michaels, CRC Press, 1999, ISBN 0-8493-0149-1
- ^ a b c Merlini, Donatella; Sprugnoli, Renzo; Verri, M. Cecilia (2006), "The Cauchy numbers", Discrete Mathematics, 306 (16): 1906–1920, doi:10.1016/j.disc.2006.03.065, MR 2251571
- ^ Markus Müller. How to Add a Non-Integer Number of Terms, and How to Produce Unusual Infinite Summations Archived 2011-06-17 at the Wayback Machine (note that he uses a slightly alternative definition of fractional sum in his work, i.e. inverse to backwards difference, hence 1 as the lower limit in his formula)
- ^ Bruce C. Berndt, Ramanujan's Notebooks Archived 2006-10-12 at the Wayback Machine, Ramanujan's Theory of Divergent Series, Chapter 6, Springer-Verlag (ed.), (1939), pp. 133–149.
- ^ Éric Delabaere, Ramanujan's Summation, Algorithms Seminar 2001–2002, F. Chyzak (ed.), INRIA, (2003), pp. 83–88.
- ^ Algorithms for Nonlinear Higher Order Difference Equations, Manuel Kauers
Further reading
- "Difference Equations: An Introduction with Applications", Walter G. Kelley, Allan C. Peterson, Academic Press, 2001, ISBN 0-12-403330-X
- Markus Müller. How to Add a Non-Integer Number of Terms, and How to Produce Unusual Infinite Summations
- Markus Mueller, Dierk Schleicher. Fractional Sums and Euler-like Identities
- S. P. Polyakov. Indefinite summation of rational functions with additional minimization of the summable part. Programmirovanie, 2008, Vol. 34, No. 2.
- "Finite-Difference Equations And Simulations", Francis B. Hildebrand, Prenctice-Hall, 1968