Indefinite sum

From HandWiki

In discrete calculus the indefinite sum operator (also known as the antidifference operator), denoted by x or Δ1,[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

Δxf(x)=f(x).

More explicitly, if xf(x)=F(x), then

F(x+1)F(x)=f(x).

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: Δ1=1eD1.

Fundamental theorem of discrete calculus

Indefinite sums can be used to calculate definite sums with the formula:[3]

k=abf(k)=Δ1f(b+1)Δ1f(a)

Definitions

Laplace summation formula

xf(x)=0xf(t)dtk=1ckΔk1f(x)k!+C
where ck=01Γ(x+1)Γ(xk+1)dx are the Cauchy numbers of the first kind, also known as the Bernoulli Numbers of the Second Kind.[4]

Newton's formula

xf(x)=k=1(xk)Δk1[f](0)+C=k=1Δk1[f](0)k!(x)k+C
where (x)k=Γ(x+1)Γ(xk+1) is the falling factorial.

Faulhaber's formula

xf(x)=n=1f(n1)(0)n!Bn(x)+C,

provided that the right-hand side of the equation converges.

Mueller's formula

If limx+f(x)=0, then[5]

xf(x)=n=0(f(n)f(n+x))+C.

Euler–Maclaurin formula

xf(x)=0xf(t)dt12f(x)+k=1B2k(2k)!f(2k1)(x)+C

Choice of the constant term

Often the constant C in indefinite sum is fixed from the following condition.

Let

F(x)=xf(x)+C

Then the constant C is fixed from the condition

01F(x)dx=0

or

12F(x)dx=0

Alternatively, Ramanujan's sum can be used:

x1f(x)=f(0)F(0)

or at 1

x1f(x)=F(1)

respectively[6][7]

Summation by parts

Main page: Summation by parts

Indefinite summation by parts:

xf(x)Δg(x)=f(x)g(x)x(g(x)+Δg(x))Δf(x)
xf(x)Δg(x)+xg(x)Δf(x)=f(x)g(x)xΔf(x)Δg(x)

Definite summation by parts:

i=abf(i)Δg(i)=f(b+1)g(b+1)f(a)g(a)i=abg(i+1)Δf(i)

Period rules

If T is a period of function f(x) then

xf(Tx)=xf(Tx)+C

If T is an antiperiod of function f(x), that is f(x+T)=f(x) then

xf(Tx)=12f(Tx)+C

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:

k=1nf(k).

In this case a closed form expression F(k) for the sum is a solution of

F(x+1)F(x)=f(x+1)

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

xa=ax+C
xx=x22x2+C
xxa=Ba+1(x)a+1+C,a
where Ba(x)=aζ(a+1,x), the generalized to real order Bernoulli polynomials.
xxa=(1)a1ψ(a1)(x)Γ(a)+C,a
where ψ(n)(x) is the polygamma function.
x1x=ψ(x)+C
where ψ(x) is the digamma function.
xBa(x)=(x1)Ba(x)aa+1Ba+1(x)+C

Antidifferences of exponential functions

xax=axa1+C

Particularly,

x2x=2x+C

Antidifferences of logarithmic functions

xlogbx=logbΓ(x)+C
xlogbax=logb(ax1Γ(x))+C

Antidifferences of hyperbolic functions

xsinhax=12csch(a2)cosh(a2ax)+C
xcoshax=12csch(a2)sinh(axa2)+C
xtanhax=1aψea(xiπ2a)+1aψea(x+iπ2a)x+C
where ψq(x) is the q-digamma function.

Antidifferences of trigonometric functions

xsinax=12csc(a2)cos(a2ax)+C,a2nπ
xcosax=12csc(a2)sin(axa2)+C,a2nπ
xsin2ax=x2+14csc(a)sin(a2ax)+C,anπ
xcos2ax=x214csc(a)sin(a2ax)+C,anπ
xtanax=ix1aψe2ia(xπ2a)+C,anπ2
where ψq(x) is the q-digamma function.
xtanx=ixψe2i(x+π2)+C=k=1(ψ(kππ2+1x)+ψ(kππ2+x)ψ(kππ2+1)ψ(kππ2))+C
xcotax=ixiψe2ia(x)a+C,anπ2
xsincx=sinc(x1)(12+(x1)(ln(2)+ψ(x12)+ψ(1x2)2ψ(x1)+ψ(1x)2))+C
where sinc(x) is the normalized sinc function.

Antidifferences of inverse hyperbolic functions

xartanhax=12ln(Γ(x+1a)Γ(x1a))+C

Antidifferences of inverse trigonometric functions

xarctanax=i2ln(Γ(x+ia)Γ(xia))+C

Antidifferences of special functions

xψ(x)=(x1)ψ(x)x+C
xΓ(x)=(1)x+1Γ(x)Γ(1x,1)e+C
where Γ(s,x) is the incomplete gamma function.
x(x)a=(x)a+1a+1+C
where (x)a is the falling factorial.
xsexpa(x)=lna(sexpa(x))(lna)x+C
(see super-exponential function)

See also

  • Indefinite product
  • Time scale calculus
  • List of derivatives and integrals in alternative calculi

References

  1. On Computing Closed Forms for Indefinite Summations. Yiu-Kwong Man. J. Symbolic Computation (1993), 16, 355-376[yes|permanent dead link|dead link}}]
  2. "If Y is a function whose first difference is the function y, then Y is called an indefinite sum of y and denoted Δ−1y" Introduction to Difference Equations, Samuel Goldberg
  3. "Handbook of discrete and combinatorial mathematics", Kenneth H. Rosen, John G. Michaels, CRC Press, 1999, ISBN:0-8493-0149-1
  4. Bernoulli numbers of the second kind on Mathworld
  5. Markus Müller. How to Add a Non-Integer Number of Terms, and How to Produce Unusual Infinite Summations (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)
  6. Bruce C. Berndt, Ramanujan's Notebooks , Ramanujan's Theory of Divergent Series, Chapter 6, Springer-Verlag (ed.), (1939), pp. 133–149.
  7. Éric Delabaere, Ramanujan's Summation, Algorithms Seminar 2001–2002, F. Chyzak (ed.), INRIA, (2003), pp. 83–88.
  8. Algorithms for Nonlinear Higher Order Difference Equations, Manuel Kauers

Further reading