Evaluation of convergent series by using finite parts

T he aim of this article is to present a method for the summation of convergent series by using the finite part of divergent series in the intermediate steps. Using divergent series to find the result of convergent processes is a very old practice as can be seen very clearly by reading chapter 1 of Hardy’s book [1] or by consulting classic texts on series such as Bromwich’s [2]; a great variety of summability methods have been employed throughout the years [1]. This article differs from older methods, however, in that we employ Hadamard finite part limits in our process. Finite part limits are explained in Sections 3 and 4, and the many examples given should convince the reader that the extraction of the finite part of a series that diverges to infinity is a regularization procedure, not a summability method. Hadamard finite parts of limits and integrals have been employed mainly in the areas of partial differential equations, since they were first defined by Hadamard [3] to find fundamental solutions of hyperbolic partial differential equations, and in the theory of distributions and generalized functions [4,5], where they play a rather important part, mostly because of their relationship with pseudofunctions. In Section 5, we find the sum of the convergent series


Introduction
T he aim of this article is to present a method for the summation of convergent series by using the finite part of divergent series in the intermediate steps. Using divergent series to find the result of convergent processes is a very old practice as can be seen very clearly by reading chapter 1 of Hardy's book [1] or by consulting classic texts on series such as Bromwich's [2]; a great variety of summability methods have been employed throughout the years [1]. This article differs from older methods, however, in that we employ Hadamard finite part limits in our process. Finite part limits are explained in Sections 3 and 4, and the many examples given should convince the reader that the extraction of the finite part of a series that diverges to infinity is a regularization procedure, not a summability method. Hadamard finite parts of limits and integrals have been employed mainly in the areas of partial differential equations, since they were first defined by Hadamard [3] to find fundamental solutions of hyperbolic partial differential equations, and in the theory of distributions and generalized functions [4,5], where they play a rather important part, mostly because of their relationship with pseudofunctions.
In Section 5, we find the sum of the convergent series by finding the finite part of the sums of the two divergent series ∞ ∑ n=0 ln n + a + 1 n + a and and substracting those results. Next, in Section 6, we consider a more complicated example, namely the convergent series ∞ ∑ n=0 (n + a) ln n + a + 1 n + a − 1 + 1 2 (n + a) , where we employ standard ideas such as Abel's summation formula and Stirling's asymptotic approximation to find the finite part of the divergent series encountered. Sections 7 and 8 contain a more systematic treatment of a class of convergent series that generalize (1) and (3). The main tool of our analysis is the formula for the finite part of the divergent series ∑ ∞ n=0 (n + a) k ln (n + a) , where H k = ∑ k j=1 (1/j) are the harmonic numbers.
We are then able to find the sum of certain convergent series involving values of the Hurwitz zeta function, namely, those of the form for |b| ≤ |a| , b = −a, k ∈ N. Series of this family have received a lot of attention since the times of Euler, who gave the interesting sums and where γ is Euler's constant. Srivastava [6] considers many of these series in an article that contains a very detailed historical account of the contributions of several authors. Recently, Blagouchine [7] and Coppo [8] have given general formulas for the series for all k, that generalize not only (6) but also the nice Surnayanarayana's formula [9] Interestingly, Surnayanarayana also considers divergent but Cesàro summable sums in [9], as the formula that can also be found by employing exterior Euler summability [10] and is also ammenable of our analysis. Notice that the series (5) diverges to infinity if b = −a. Nevertheless, after reading this article the reader will be able to see that

Preliminaries
In this section, we fix the notation employed and recall several useful well known facts; see [11] and [12] for details.
We shall always consider the principal branches of powers and logarithms, z α and ln z, defined for z ∈ C \ (−∞, 0]. If a ∈ C \ (−∞, 0], the Hurwicz zeta function ζ (s, a) is the analytic continuation to s ∈ C \ {1} of the function The zeta function ζ (s, a) has a simple pole at s = 1, with residue 1. If a = 1 it reduces to the Riemman zeta function, ζ (s, 1) = ζ (s) . We shall also employ the digamma function ψ (a) = Γ (a) /Γ (a) , the logarithmic derivative of the gamma function, whose Mittag-Leffer expansion is the following where γ is Euler's constant. The digamma function is actually the finite part -explained in the next sectionof the negative of the Hurwitz zeta function, −ζ (s, a) , at s = 1, that is for a ∈ C \ (−∞, 0]. Notice also that It is easy to see that for n ≥ 1, and this yields the Taylor expansion which is an absolutely convergent series for |ω| < |a| . The series is divergent for |ω| = |a| but it is Cesàro summable 1 in the this circle except for ω = −a. Sometimes, we shall need to use Abel's summation formula [13,Chp. 1], where A n = ∑ n j=1 a j . We shall also employ Stirling's asymptotic formula [12], as z → ∞ in the region −π < arg z < π.

The finite part of limits
Let us now recall the notion of the finite part of a limit [13,Section 2.4]. Let X be a topological space, and let x 0 ∈ X. Suppose F, the basic functions, is a family of strictly positive functions defined for x ∈ V \ {x 0 } , where V is a neighborhood of x 0 , such that all of them tend to infinity at x 0 and such that, given two different The finite part of the limit of G (x) as x → x 0 with respect to F exists and equals A if we can write 2 G (x) = G 1 (x) + G 2 (x) , where G 1 , the infinite part, is a linear combination of the basic functions and where G 2 , the finite part, has the property that the limit A = lim x→x 0 G 2 (x) exists. We then employ the notation We shall be interested in the Hadamard finite parts at infinity in this article. It is the case when F is the family of powers and logarithms. Definition 2. The Hadamard finite part limit corresponds to the case when x 0 = 0 and F is the family of functions x −α |ln x| β , where α > 0 and β ≥ 0 or where α = 0 and β > 0 or when x 0 = ∞ and F is the family of functions x α |ln x| β , where α > 0 and β ≥ 0 or where α = 0 and β > 0. We then use the simpler notations Applying the notion of finite part of limits we can define the finite part of integrals at singularities [13, Section 2.4], but in this article, we will be interested in the finite part of series, 1 Details on summability can be found in [1,13] 2 Such a decomposition, if it exists, is unique since any finite number of elements of F has to be linearly independent.
Observe that the finite part of a series gives a finite sum to a series that diverges to infinity; it does not work for oscillatory series. This makes the finite part process very different from summability methods like Abel summability of Cesàro summability. The finite part process is more of a regularization procedure than a summability method [13,Section 2.4].
Needless to say, many of the usual operations that are valid for convergent series may not work for the finite parts. For example, F.p. ∑ ∞ n=n 0 +1 a n−1 does not have to be equal to F.p. ∑ ∞ n=n 0 a n . The notion of finite part integrals and its name were introduced by Hadamard [3], who used them in his study of fundamental solutions of partial differential equations. The finite part method has been a very important part of the theory of distributions since its beginnings [5] since it is closely related to the notion of pseudofunction.

Some finite part limits
In this section, for future reference, we will give the formulas for several useful finite part limits. Our definition says that F.p. lim N→∞ N β = 0 for all β ∈ R. While it is true that for all a ∈ C, we have F.p. lim In a similar fashion, while F.p. lim N→∞ N β ln N = 0 for all β ∈ R, by definition, when we replace N by N + a we have the ensuing formulas.
where H k is the harmonic sum Observe now that which yields (25).
Next, we would like to consider the finite part limit of several finite sums. We have that ∑ ∞ n=1 n −β = ζ (β) , a convergent series if e β > 1. We now consider the finite part limit F.p. lim N→∞ ∑ N n=1 n −β for e β ≤ 1; do we obtain ζ (β)? The answer is yes almost all the time, but not always. Proof. Formula (29) is trivial, while (30) follows from the well known asymptotic relation In the other cases, we just employ the asymptotic expansion that comes from the Euler-Maclaurin formula.
Similar formulas are obtained if we replace n by n + a in the previous lemma.
Proof. The proof is basically the same as that of the previous lemma if we use the asymptotic expansion (56) instead of (32) and the formula [14] that corresponds to (31).
It is important to point out that while ζ (β) is obtained by putting a = 1 in ζ (β, a) , the formulas in the Lemmas 2 and 3, are not exactly the same: the fact that we obtain a k in (33) and not 0 as in (29) is explained because one sum is from n = 0 to N while the other from n = 1 to N.

A simple example
Let us consider the convergent series where a is a constant, a = 0, −1, −2, . . . . We can easily evaluate it by using finite parts, since we have a telescoping sum, which, when we notice that F.p. lim N→∞ ln (N + a + 1) = 0, gives the finite part limit There is also a simple alternative evaluation of σ 0 (a) , without the use of finite parts, since (13) gives We also observe that from (17), If a = 1, we obtain the well known formula that goes back to Euler,

Another example
Our next task is to find the sum of the convergent series First, we approximate the sum n ln n + a + 1 n + a = N ∑ n=0 n (ln (n + a + 1) − ln (n + a)) (46) by employing Abel's summation formula (18), and then Stirling's formula (19), to obtain and consequently If we now recall (39), we obtain F.p. and Formulas (49) so that ν 1 (a) = σ 1 (a) , that is, Notice that if a = 1 we obtain Surnayanarana's formula [9] while a = 2 yields the result [6]

The functions ϕ k (a)
The Euler-Maclaurin formula yields the asymptotic development as N → ∞. Differentiation of (56) gives the expansion Let us define We may evaluate ϕ 0 (a) in two ways, one using Stirlings formula (19), the other using (57) that gives ϕ 0 (a) = −ζ (0, a) , so that Next we consider ϕ 1 (a) . We have from (57) that ϕ 1 (a) is given as the finite part of the limit of as N → ∞. Hence using the Lemma 1 we obtain More generally, we have the following formula.
where H k are the harmonic numbers (26).
We also have because of (25). Consequently, P k+1 is the a polynomial solution of the equation Actually it is the only polynomial solution that satisfies the initial condition since P k+1 (0) is the constant term in the polynomial P k+1 (a) , and (57) tell us that that constant term is precisely [11]. Hence H k B k+1 (a + 1) / (k + 1) is another polynomial solution of the initial value problem for the difference equation. Therefore Formula (62) is obtained since as is well known [11].
When a = 1 it is convenient to consider a variation of the function ϕ k (1) .
Proof. Indeed, using the Lemma 1, we obtain 8. The series ∑ ∞ n=2 (−1) n ζ (n, a) b n+k /(n + k) We now employ finite part limits to find a formula for the sum of some general convergent series.
Proof. Indeed, from (17), we obtain Observe now that If we notice now that and substitute, we obtain (69).