Bourbaki

§ 4. Application to series with positive terms

1. Convergence criteria for series with positive terms

In this section by a series with positive terms we shall understand (by an abuse of language) a series (un)(u_n) of real terms such that un≥0u_n \geq 0 starting from some particular value of n. Everything we shall say about such series will extend immediately by a change of sign to series all of whose terms are ≤0\leq 0 starting from some particular value of n. We have seen (II, p. 64, Example 3) that to every sequence (un)n≥1(\mathbf{u}_n)_{n \geq 1} of points in a normed space E one can associate a step function u\mathbf{u} defined on [1,+∞[[1, +\infty[ by the conditions u(x)=un\mathbf{u}(x) = \mathbf{u}_n for n≤x<n+1n \leq x < n+1: then the series (un)(\mathbf{u}_n) converges if and only if the integral ∫1+∞u(t) dt\int_1^{+\infty} \mathbf{u}(t)\,dt converges.

Let (un)(u_n) and (vn)(v_n) be two series with positive terms, and uu and vv the associated step functions: the relation un≤vnu_n \leq v_n for n≥n0n \geq n_0 is equivalent to u(x)≤v(x)u(x) \leq v(x) for x≥n0x \geq n_0. Thus each of the relations un≼vn,un≪vn,un∼vnu_n \preccurlyeq v_n, u_n \ll v_n, u_n \sim v_n is equivalent to u(x)≼v(x),u(x)≪v(x)u(x) \preccurlyeq v(x), u(x) \ll v(x) and u(x)∼v(x)u(x) \sim v(x) respectively; this remark allows us to translate propositions 1 (V, p. 228) and 6 (V, p. 230) as follows:

Proposition 1

Let* (un)(u_n) and (vn)(v_n) be two series with positive terms. If un≼vnu_n \preccurlyeq v_n and if the series (vn)(v_n) converges, then (un)(u_n) converges; if un≻vnu_n \succ v_n and if ∑n=1∞vn=+∞\sum_{n=1}^{\infty} v_n = +\infty, then ∑n=1∞un=+∞\sum_{n=1}^{\infty} u_n = +\infty.

Proposition 2

Let* (un)(u_n) and (vn)(v_n) be two series with positive terms: 1. If the series (vn)(v_n) converges then the relation un≪vnu_n \ll v_n (resp. un∼vnu_n \sim v_n) implies ∑p=n∞up≪∑p=n∞vp\sum_{p=n}^{\infty} u_p \ll \sum_{p=n}^{\infty} v_p (resp. ∑p=n∞up∼∑p=n∞vp\sum_{p=n}^{\infty} u_p \sim \sum_{p=n}^{\infty} v_p). 2. If ∑n=1∞vn=+∞\sum_{n=1}^{\infty} v_n = +\infty then the relation un≪vnu_n \ll v_n (resp. un∼vnu_n \sim v_n) implies ∑p=1nup≪∑p=1nvp\sum_{p=1}^{n} u_p \ll \sum_{p=1}^{n} v_p (resp. ∑p=1nup∼∑p=1nvp\sum_{p=1}^{n} u_p \sim \sum_{p=1}^{n} v_p).

One obtains convenient convergence criteria by taking for the comparison series (vn)(v_n) in prop. 1 a series whose terms are of the form vn=f(n)v_n = f(n), where ff is a function ≥0\geq 0 defined for every real number x>x0x > x_0 and decreasing on the interval [x0,+∞[[x_0, +\infty[; indeed:

Proposition 3 (Cauchy-Maclaurin criterion)

If* ff is a function ≥0\geq 0 and decreasing on [x0,+∞[[x_0, +\infty[, then the series with general term vn=f(n)v_n = f(n) converges if and only if the integral ∫x0+∞f(t) dt\int_{x_0}^{+\infty} f(t)\,dt converges.

To prove this it is sufficient to note that if vv is the step function associated with the series (vn)(v_n) then v(x+1)≤f(x)≤v(x)v(x+1) \leq f(x) \leq v(x) for every x≥x0x \geq x_0 since ff is decreasing; the proposition is thus a consequence of the comparison principle (II, p. 66, prop. 3).

Since the functions which feature in the logarithmic convergence criteria for integrals (V, p. 229, prop. 2, 3 and 4) are decreasing on an interval [x0,+∞[[x_0, +\infty[, applying prop. 1 and 3 of V, p. 237, gives the following criteria:

Proposition 4 ("logarithmic criterion of order 0")

Let (un)(u_n) be a series with positive terms; if un≼nμu_n \preccurlyeq n^\mu for some μ<−1\mu < -1 then the series (un)(u_n) converges; if un≽nμu_n \succcurlyeq n^\mu for some μ≥−1\mu \geq -1 then the series (un)(u_n) has an infinite sum.

Proposition 5 ("logarithmic criterion of order pp")

Let (un)(u_n) be a series with positive terms. If un≼1nl1(n)l2(n)…lp−1(n)(lp(n))μu_n \preccurlyeq \frac{1}{nl_1(n)l_2(n)\ldots l_{p-1}(n)(l_p(n))^\mu} for some μ>1\mu > 1 then the series (un)(u_n) converges; if un≽1nl1(n)l2(n)…lp−1(n)(lp(n))μu_n \succcurlyeq \frac{1}{nl_1(n)l_2(n)\ldots l_{p-1}(n)(l_p(n))^\mu} for some μ≤1\mu \leq 1 then the series (un)(u_n) has an infinite sum.

If 0≤q<10 \leq q < 1 one has qn≼n−μq^n \preccurlyeq n^{-\mu} for any μ>0\mu > 0; applying the logarithmic criterion of order 0 again proves the convergence of the geometric series ∑n=0∞qn\sum_{n=0}^\infty q^n for ∣q∣<1|q| < 1 (Gen. Top., IV, p. 364). If one applies prop. 1 with vn=qnv_n = q^n one obtains a criterion which may be put in the following form ("Cauchy's criterion"): Let (un)(u_n) be a series with positive terms; if lim sup⁡n→∞(un)1/n<1\limsup_{n \to \infty}(u_n)^{1/n} < 1 then the series (un)(u_n) converges; if lim sup⁡n→∞(un)1/n>1\limsup_{n \to \infty}(u_n)^{1/n} > 1 then the series (un)(u_n) has infinite sum. Indeed, if lim sup⁡n→∞(un)1/n=a<1\limsup_{n \to \infty}(u_n)^{1/n} = a < 1 then un≼qnu_n \preccurlyeq q^n for every qq such that a<q<1a < q < 1. If, on the other hand, lim sup⁡n→∞(un)1/n=a>1\limsup_{n \to \infty}(u_n)^{1/n} = a > 1 then un≽qn>1u_n \succcurlyeq q^n > 1 for infinitely many values of nn, for any qq such that 1<q<a1 < q < a; since unu_n does not tend to 0 one has ∑n=1∞un=+∞\sum_{n=1}^\infty u_n = +\infty.

This criterion is very useful in the theory of entire series, which we shall study later; but it even so it does not permit one to decide the convergence of the series (1/nα)(1/n^\alpha), in other words, its field of application is much more restricted than that of the logarithmic criteria.

2. Asymptotic expansion of the partial sums of a series

For xx real tending to +∞+\infty let E\mathscr{E} be a comparison scale formed by functions each of which is defined on a whole interval [x0,+∞[[x_0, +\infty[ (depending on the function) and is ≥0\geq 0 on this interval. Let (un)(\mathbf{u}_n) be a series whose terms belong to a complete normed space E\mathbf{E}, such that un\mathbf{u}_n admits an asymptotic expansion to precision gαg_\alpha with respect to the scale E′\mathscr{E}' of restrictions to N\mathbf{N} of the functions in E\mathscr{E}:

un=∑λ≤αaλgλ(n)+rα(n). \mathbf{u}_n = \sum_{\lambda \leq \alpha} \mathbf{a}_\lambda g_\lambda(n) + \mathbf{r}_\alpha(n).

Suppose that every partial sum ∑m=1ng(m)\sum_{m=1}^n g(m), where g∈Eg \in \mathscr{E}, admits an asymptotic expansion relative to E′\mathscr{E}'. One can then obtain an asymptotic expansion of the sn=∑m=1nums_n = \sum_{m=1}^n \mathbf{u}_m with respect to E′\mathscr{E}'; again we distinguish two cases:

1° ∑n=1∞gα(n)=+∞\sum_{n=1}^{\infty} g_{\alpha}(n) = +\infty. Then (V, p. 237, prop. 2) one has ∑m=1nrα(m)≪∑m=1ngα(m)\sum_{m=1}^{n} r_{\alpha}(m) \ll \sum_{m=1}^{n} g_{\alpha}(m): by hypothesis one can obtain an asymptotic expansion of

∑λ≤αaλ(∑m=1ng(m)) \sum_{\lambda \leq \alpha} a_{\lambda} \left( \sum_{m=1}^{n} g(m) \right)
(V, p. 222) to a certain precision gσg_{\sigma}; if cgσ(n)c g_{\sigma}(n) is the principal part of ∑m=1ngα(m)\sum_{m=1}^{n} g_{\alpha}(m), one will have an asymptotic expansion of sns_{n} to precision gmin⁡(ρ,σ)g_{\min(\rho,\sigma)}.

2° ∑n=1∞gα(n)\sum_{n=1}^{\infty} g_{\alpha}(n) converges; then let β\beta be the smallest of the indices λ≤α\lambda \leq \alpha such that aλ≠0a_{\lambda} \neq 0 and such that ∑n=1∞gλ(n)\sum_{n=1}^{\infty} g_{\lambda}(n) converges; the series

C=∑n=1∞(un−∑λ<βaλgλ(n)) C = \sum_{n=1}^{\infty} \left( u_{n} - \sum_{\lambda < \beta} a_{\lambda} g_{\lambda}(n) \right)
then converges absolutely, and one can write
sn=∑λ<βaλ(∑m=1ngλ(m))+C−∑β≤λ≤αaλ(∑m=n+1∞gλ(m))−∑m=n+1∞rα(m). s_{n} = \sum_{\lambda < \beta} a_{\lambda} \left( \sum_{m=1}^{n} g_{\lambda}(m) \right) + C - \sum_{\beta \leq \lambda \leq \alpha} a_{\lambda} \left( \sum_{m=n+1}^{\infty} g_{\lambda}(m) \right) - \sum_{m=n+1}^{\infty} r_{\alpha}(m).
Further, ∑m=n+1∞rα(m)≪∑m=n+1∞gα(m)\sum_{m=n+1}^{\infty} r_{\alpha}(m) \ll \sum_{m=n+1}^{\infty} g_{\alpha}(m); if cgσ(n)c g_{\sigma}(n) is the principal part of ∑m=n+1∞gα(m)\sum_{m=n+1}^{\infty} g_{\alpha}(m), and if one has an asymptotic expansion of
∑λ<βaλ(∑m=1ngλ(m))+C−∑β≤λ≤αaλ(∑m=n+1∞gλ(m)) \sum_{\lambda < \beta} a_{\lambda} \left( \sum_{m=1}^{n} g_{\lambda}(m) \right) + C - \sum_{\beta \leq \lambda \leq \alpha} a_{\lambda} \left( \sum_{m=n+1}^{\infty} g_{\lambda}(m) \right)
to precision gρg_{\rho} one thus obtains an asymptotic expansion of sns_{n} to precision gmin⁡(ρ,σ)g_{\min(\rho,\sigma)}.

One is thus led to the particular case of series (g(n))(g(n)) where g∈Eg \in \mathscr{E}. We shall see how, subject to certain conditions, one can straight away obtain a principal part of sn=∑m=1ng(m)s_{n} = \sum_{m=1}^{n} g(m) (when ∑n=1∞g(n)=+∞\sum_{n=1}^{\infty} g(n) = +\infty) or of rn=∑m=n+1∞g(m)r_{n} = \sum_{m=n+1}^{\infty} g(m) (when ∑n=1ng(n)<+∞\sum_{n=1}^{n} g(n) < +\infty).

Proposition 6

Let g be a real function, > 0 and monotone, defined on an interval [x_0, +\infty[ (where x_0 \leq 1 ), and such that log g and x are comparable to order 1.

1° If g is of infinite order relative to e, one has

sn=∑m=1ng(m)∼g(n)if∑n=1∞g(n)=+∞;(1) s_n = \sum_{m=1}^n g(m) \sim g(n) \qquad \text{if} \quad \sum_{n=1}^\infty g(n) = +\infty; \tag{1}
rn=∑m=n+1∞g(m)∼g(n+1)if∑n=1∞g(n)<+∞.(2) r_n = \sum_{m=n+1}^\infty g(m) \sim g(n+1) \qquad \text{if} \quad \sum_{n=1}^\infty g(n) < +\infty. \tag{2}

2° If g is of finite order μ\mu with respect to exe^x one has

sn=∑m=1ng(m)∼μ1−e−μ∫λ0ng(t) dtif∑n=1∞g(n)=+∞;(3) s_n = \sum_{m=1}^n g(m) \sim \frac{\mu}{1 - e^{-\mu}} \int_{\lambda_0}^n g(t) \, dt \qquad \text{if} \quad \sum_{n=1}^\infty g(n) = +\infty; \tag{3}
rn=∑m=n+1∞g(m)∼μ1−e−μ∫n∞g(t) dtif∑n=1∞g(n)<+∞(4) r_n = \sum_{m=n+1}^\infty g(m) \sim \frac{\mu}{1 - e^{-\mu}} \int_n^\infty g(t) \, dt \qquad \text{if} \quad \sum_{n=1}^\infty g(n) < +\infty \tag{4}

(the number μ1−e−μ\frac{\mu}{1 - e^{-\mu}} is to be replaced by 1 in (3) and (4) when μ=0\mu = 0).

1° If g is of order +∞+\infty relative to exe^x one has log⁡g≫x\log g \gg x whence g′/g≫1g'/g \gg 1, or g′≫gg' \gg g, by the hypothesis; g is therefore increasing and tends to +∞+\infty with xx, whence ∑n=1∞g(n)=+∞\sum_{n=1}^\infty g(n) = +\infty. If uu is the step function associated with the series (g(n))(g(n)) (V, p. 237), one has u(x)≤g(x)u(x) \leq g(x) starting from a certain value of xx, so u≪gu \ll g and consequently

sn−1=∫1nu(t) dt≪∫1ng(t) dt≪∫1ng′(t) dt∼g(n); s_{n-1} = \int_1^n u(t) \, dt \ll \int_1^n g(t) \, dt \ll \int_1^n g'(t) \, dt \sim g(n);

since sn=sn−1+g(n)s_n = s_{n-1} + g(n) one has sn∼g(n)s_n \sim g(n). The proof is similar when g is of order −∞-\infty relative to exe^x; we thus obtain formula (2).

2° If g is of finite order μ\mu relative to exe^x one can write g(x)=eμxh(x)g(x) = e^{\mu x} h(x) where h is of order 0 relative to exe^x; further, by hypothesis, log⁡g∼μx\log g \sim \mu x for μ≠0\mu \neq 0 (log⁡g≪x\log g \ll x for μ=0\mu = 0) implies h′≪hh' \ll h. Suppose first that ∑n=1∞g(n)=+∞\sum_{n=1}^\infty g(n) = +\infty (which implies that μ≥0\mu \geq 0; the converse always holds if μ>0\mu > 0, since then g(x)g(x) tends to +∞+\infty with xx); let us evaluate the principal part of ∫n−1ng(t) dt\int_{n-1}^n g(t) \, dt. One can write

∫n−1ng(t) dt=∫n−1neμth(t) dt \int_{n-1}^n g(t) \, dt = \int_{n-1}^n e^{\mu t} h(t) \, dt
=h(n)∫n−1neμt dt+∫n−1neμt(h(t)−h(n)) dt = h(n) \int_{n-1}^n e^{\mu t} \, dt + \int_{n-1}^n e^{\mu t} (h(t) - h(n)) \, dt
=1−e−μμg(n)+∫n−1neμt(h(t)−h(n)) dt. = \frac{1 - e^{-\mu}}{\mu} g(n) + \int_{n-1}^n e^{\mu t} (h(t) - h(n)) \, dt.

Now, the relation h′≪hh' \ll h implies that for every ε>0\varepsilon > 0 there is an n0n_0 such that the relation x≥n0x \geq n_0 implies that ∣h′(x)/h(x)∣≤ε|h'(x)/h(x)| \leq \varepsilon; one deduces from the mean value theorem that −ε≤log⁡∣h(t)/h(n)∣≤ε-\varepsilon \leq \log |h(t)/h(n)| \leq \varepsilon, for n−1≤t≤nn-1 \leq t \leq n, if n≥n0n \geq n_0, whence

∣h(t)−h(n)∣≤(eε−1)h(n) |h(t) - h(n)| \leq (e^\varepsilon - 1)h(n)

and consequently

∣∫n−1neμt(h(t)−h(n)) dt∣≤(eε−1)eμnh(n)=(eε−1)g(n) \left| \int_{n-1}^n e^{\mu t} (h(t) - h(n)) \, dt \right| \leq (e^\varepsilon - 1)e^{\mu n}h(n) = (e^\varepsilon - 1)g(n)

since eμte^{\mu t} is increasing. Since eε−1e^\varepsilon - 1 becomes arbitrarily small with ε\varepsilon one sees that one can write

∫n−1ng(t) dt=1−e−μμg(n)+o(g(n)) \int_{n-1}^n g(t) \, dt = \frac{1 - e^{-\mu}}{\mu} g(n) + o(g(n))

(1−e−μμ being replaced by 1 when μ=0)\left( \frac{1 - e^{-\mu}}{\mu} \text{ being replaced by } 1 \text{ when } \mu = 0 \right). The proposition is then a consequence of prop. 2 of V, p. 237. One argues similarly when ∑n=1∞g(n)\sum_{n=1}^\infty g(n) is finite.

By applying prop. 6 of V, p. 239, repeatedly one can then sometimes obtain an asymptotic expansion for sn=∑m=1ng(m)s_n = \sum_{m=1}^n g(m). First suppose that gg is of order +∞+\infty relative to eλe^\lambda; for every fixed value of pp one can write, by prop. 6,

sn=g(n)+g(n−1)+⋯+g(n−p)+o(g(n−p)) s_n = g(n) + g(n-1) + \cdots + g(n-p) + o(g(n-p))

and it suffices to expand (relative to E′\mathscr{E}') each of the functions g(n−k)g(n-k) (0≤k≤p0 \leq k \leq p), limiting the precision of the expansions to the principal part of g(n−p)g(n-p), to obtain an expansion for the sns_n.

Example

Let g(x)=xλ=exp⁡(λlog⁡x)g(x) = x^\lambda = \exp(\lambda \log x), of order +∞+\infty relative to eλe^\lambda. Taking p=2p = 2 one has

(n−1)log⁡(n−1)=(n−1)log⁡n−1+12n+o(1n) (n-1)\log(n-1) = (n-1)\log n - 1 + \frac{1}{2n} + o\left( \frac{1}{n} \right)

whence (V, p. 225)

(n−1)n−1=1enn−1+12enn−2+o1(nn−2) (n-1)^{n-1} = \frac{1}{e} n^{n-1} + \frac{1}{2e} n^{n-2} + o_1\left( n^{n-2} \right)

and similarly

(n−2)n−2=1e2nn−2+o2(nn−2); (n-2)^{n-2} = \frac{1}{e^2} n^{n-2} + o_2\left( n^{n-2} \right);

consequently

sn=nn+1enn−1+(12e+1e2)nn−2+o3(nn−2). s_n = n^n + \frac{1}{e} n^{n-1} + \left( \frac{1}{2e} + \frac{1}{e^2} \right) n^{n-2} + o_3\left( n^{n-2} \right).

One proceeds similarly (for rnr_n) when gg is of order −∞-\infty relative to eλe^\lambda.

Now if gg is of finite order μ\mu relative to exe^x, and if, for example, ∑n=1∞g(n)=+∞\sum_{n=1}^\infty g(n) = +\infty, one can write

sn=μ1−e−μ∫1ng(t) dt+∑m=1nf1(m) s_n = \frac{\mu}{1 - e^{-\mu}} \int_1^n g(t)\, dt + \sum_{m=1}^n f_1(m)
where f1(n)=g(n)−μ1−e−μ∫1ng(t) dt≪g(n)f_1(n) = g(n) - \frac{\mu}{1 - e^{-\mu}} \int_1^n g(t)\, dt \ll g(n) by prop. 6 of V, p. 239. If one has a principal part cg1(n)c g_1(n) of f1(n)f_1(n) relative to E′\mathscr{E}', and if one can again apply prop. 6 to the function g1g_1 one will obtain a primitive equivalent to ∑m=1nf1(m)\sum_{m=1}^n f_1(m) if ∑n=1∞g1(n)=+∞\sum_{n=1}^\infty g_1(n) = +\infty, and equivalent to ∑m=n+1∞f1(m)\sum_{m=n+1}^\infty f_1(m) in the opposite case (in the latter case one writes ∑m=1nf1(m)=C−∑m=n+1∞f1(m)\sum_{m=1}^n f_1(m) = C - \sum_{m=n+1}^\infty f_1(m), with C=∑n=1∞f1(n)C = \sum_{n=1}^\infty f_1(n)).

Step by step one can thus eventually obtain an expression for sns_n as the sum of a certain number of primitives each of which is negligible with respect to the previous, of a term remaining negligible relative to the last primitive written, and finally a constant (the case where the remainder term tends to 0). It then remains to expand each of the primitives obtained with respect to E′\mathscr{E}' (cf. V, p. 235).

Example

Let g(n)=1ng(n) = \frac{1}{n}; then

sn=∑m=1n1m∼∫1ndtt=log⁡n s_n = \sum_{m=1}^n \frac{1}{m} \sim \int_1^n \frac{dt}{t} = \log n
then
1n−(log⁡n−log⁡(n−1))∼−12n2 \frac{1}{n} - (\log n - \log(n-1)) \sim -\frac{1}{2n^2}
whence
sn=log⁡n+γ+12n+o(1n). s_n = \log n + \gamma + \frac{1}{2n} + o\left(\frac{1}{n}\right).
The constant γ\gamma which appears in this formula plays an important rôle in Analysis (cf. chap. VI and VII); it is known as Euler's constant; one has
γ=0.577 215 664… \gamma = 0.577\,215\,664\ldots
to within 1/1091/10^9.

We shall see in VI, p. 288, how the Euler-Maclaurin summation formula gives an asymptotic expansion of arbitrary order for sns_n (or for rnr_n) in the most important cases.

3. Asymptotic expansion of the partial products of an infinite product

One knows (Gen. Top., V. p. 22 and 23) that for the infinite product with general factor 1+un1 + u_n (un>−1u_n > -1) to be convergent (resp. commutatively convergent) it is necessary and sufficient that the series with general term log⁡(1+un)\log(1 + u_n) should be convergent (resp. commutatively convergent), and that one then has the relation

log⁡∏n=1∞(1+un)=∑n=1∞log⁡(1+un). \log \prod_{n=1}^{\infty} (1 + u_n) = \sum_{n=1}^{\infty} \log(1 + u_n).

When the infinite product converges one knows that unu_n tends to 0; thus log⁡(1+un)∼un\log(1 + u_n) \sim u_n; now one knows that for a series of real numbers to be commutatively convergent it is necessary and sufficient that it should be absolutely convergent (Gen. Top., IV, p. 372, prop. 5); by prop. 1 one thus recovers the fact that the product with general factor 1+un1 + u_n is commutatively convergent if and only if the series with general term unu_n is absolutely convergent (Gen. Top., IV, p. 368, th. 4).

A similar argument applies to an infinite product whose general factor is a complex number 1+un1 + u_n (un≠−1u_n \neq -1). Indeed, for such a product to be commutatively convergent it is necessary and sufficient (Gen. Top., VIII, p. 115, prop. 2) that the infinite product with general factor ∣1+un∣|1 + u_n| should be so, and further, if θn\theta_n is the amplitude of 1+un1 + u_n (taken between −π-\pi and +π+\pi), that the series of the θn\theta_n should be commutatively convergent. Since unu_n tends to 0, log⁡(1+un)\log(1 + u_n) is defined starting from some particular value of nn (III, p. 100) and one has

log⁡(1+un)=log⁡∣1+un∣+iθn; \log(1 + u_n) = \log |1 + u_n| + i \theta_n;

thus, for the product with general factor 1+un1 + u_n to be commutatively convergent it is necessary and sufficient that the series with general term ∣log⁡(1+un)∣|\log(1 + u_n)| be absolutely convergent (Gen Top., VII, p. 84, th. 1); now log⁡(1+un)∼un\log(1 + u_n) \sim u_n (I, p. 18, prop. 5), so one again obtains the condition that the series with general term unu_n should be absolutely convergent (Gen. Top., VIII, p. 116, th. 1).

The relation between infinite products and series of real numbers sometimes allows one to obtain an asymptotic expansion for the partial product pn=∏k=1n(1+uk)p_n = \prod_{k=1}^{n} (1 + u_k); it suffices to have an asymptotic expansion for the partial sum sn=∑k=1nlog⁡(1+uk)s_n = \sum_{k=1}^{n} \log(1 + u_k), then to expand pn=exp⁡(sn)p_n = \exp(s_n); one is thus brought back to the two problems examined earlier (V, p. 238 and p. 226).

Example: Stirling’s formula. Let us seek an asymptotic expansion for n!n!; this leads us to expand sn=∑p=1nlog⁡ps_n = \sum_{p=1}^{n} \log p, and then exp⁡(sn)\exp(s_n). The method of n 2 gives successively

sn=∑p=1nlog⁡p∼∫1nlog⁡t dt=nlog⁡n−n+1 s_n = \sum_{p=1}^{n} \log p \sim \int_{1}^{n} \log t \, dt = n \log n - n + 1

then

log⁡n−∫n−1nlog⁡t dt=log⁡n−(nlog⁡n−(n−1)log⁡(n−1)−1)∼12n \log n - \int_{n-1}^{n} \log t \, dt = \log n - (n \log n - (n-1) \log(n-1) - 1) \sim \frac{1}{2n}
whence
sn=nlog⁡n−n+12log⁡n+o(log⁡n). s_n = n \log n - n + \frac{1}{2} \log n + o(\log n).
Then
log⁡n−∫n−1nlog⁡t dt−12(log⁡n−log⁡(n−1))∼−112n2 \log n - \int_{n-1}^{n} \log t \, dt - \frac{1}{2} (\log n - \log(n-1)) \sim -\frac{1}{12n^2}
whence
sn=nlog⁡n−n+12log⁡n+k+112n+o1(1n)(k constant) s_n = n \log n - n + \frac{1}{2} \log n + k + \frac{1}{12n} + o_1 \left( \frac{1}{n} \right) \quad (k \text{ constant})
and one finally deduces (V, p. 226)
n!=ek nn+1/2 e−n(1+112n+o2(1n)).(5) n! = e^k \, n^{n+1/2} \, e^{-n} \left( 1 + \frac{1}{12n} + o_2 \left( \frac{1}{n} \right) \right). \tag{5}
We shall show in VII, p. 322, that ek=2πe^k = \sqrt{2\pi}. The formula (5) (with this value of kk) is called Stirling's formula. In the same way, for every real number aa not an integer >0> 0, one shows that
(a+1)(a+2)…(a+n)∼K(a) nn+a+12 e−n.(6) (a+1)(a+2)\ldots(a+n) \sim K(a) \, n^{n+a+\frac{1}{2}} \, e^{-n}. \tag{6}
We shall determine the function K(a)K(a) too (VII, p. 18). From formulae (5) and (6) one derives in particular that
(an)∼(−1)nφ(a) n−a−1(7) \binom{a}{n} \sim (-1)^n \varphi(a) \, n^{-a-1} \tag{7}
for every real number aa not an integer >0> 0, where φ(a)\varphi(a) is a function of aa which will be specified in VII, p. 322.

4. Application: convergence criteria of the second kind for series with positive terms

Quite often one meets series (un)(u_n) for which un>0u_n > 0 from a certain point on, and un+1/unu_{n+1}/u_n has an asymptotic expansion which is easy to determine. It is convenient, for such series, to have criteria (called criteria of the second kind) allowing one to determine whether the series is convergent from the form of un+1/unu_{n+1}/u_n alone. The following is such a criterion:

Proposition 7 ("Raabe's test")

Let (un)(u_n) be a series with terms >0> 0 from some point on. If, from a certain stage, un+1/un⩽1−αnu_{n+1}/u_n \leqslant 1 - \frac{\alpha}{n} for some α>1\alpha > 1, then the series (un)(u_n) converges; if, from a certain point on, un+1/un⩾1−1nu_{n+1}/u_n \geqslant 1 - \frac{1}{n}, then the series (un)(u_n) has infinite sum.

Indeed, if un+1/un⩽1−αnu_{n+1}/u_n \leqslant 1 - \frac{\alpha}{n} with α>1\alpha > 1, for all n⩾n0n \geqslant n_0, then one has un≼pn=∏k=n0n(1−αk)u_n \preccurlyeq p_n = \prod_{k=n_0}^{n} \left( 1 - \frac{\alpha}{k} \right). Now, log⁡(1−αn)=−αn−α22n2+o(1n2)\log \left( 1 - \frac{\alpha}{n} \right) = -\frac{\alpha}{n} - \frac{\alpha^2}{2n^2} + o \left( \frac{1}{n^2} \right), whence log⁡pn=\log p_n = −αlog⁡n+k+o(1/n)-\alpha \log n + k + o(1/n) (kk constant), and pn∼ek1nαp_n \sim e^k \frac{1}{n^\alpha}; since α>1\alpha > 1 the logarithmic criterion of order 0 allows one to finish.

If, on the other hand, un+1/un≥1−1nu_{n+1}/u_n \geq 1 - \frac{1}{n} from a certain point, then the same calculation proves that un≥1nu_n \geq \frac{1}{n}, whence the proposition.

One can prove in the same way, using the logarithmic criteria of order >0> 0, the following criterion of the second kind:

Proposition 8

Let (un)(u_n) be a series with terms >0> 0 from a certain point on. If, from a certain point, one has

un+1un≤1−1n−1nl1(n)−…−1nl1(n)l2(n)…lp−1(n)−αnl1(n)l2(n)…lp(n) \frac{u_{n+1}}{u_n} \leq 1 - \frac{1}{n} - \frac{1}{nl_1(n)} - \ldots - \frac{1}{nl_1(n)l_2(n)\ldots l_{p-1}(n)} - \frac{\alpha}{nl_1(n)l_2(n)\ldots l_p(n)}

for some α>1\alpha > 1, then the series (un)(u_n) converges; if, from some point on, one has

un+1un≥1−1n−1nl1(n)−…−1nl1(n)l2(n)…lp(n) \frac{u_{n+1}}{u_n} \geq 1 - \frac{1}{n} - \frac{1}{nl_1(n)} - \ldots - \frac{1}{nl_1(n)l_2(n)\ldots l_p(n)}

then the series (un)(u_n) has an infinite sum.

Example

Consider the hypergeometric series, with general term

un=α(α+1)…(α+n−1) β(β+1)…(β+n−1)1.2…n.γ(γ+1)…(γ+n−1) u_n = \frac{\alpha(\alpha+1)\ldots(\alpha+n-1)\,\beta(\beta+1)\ldots(\beta+n-1)}{1.2\ldots n.\gamma(\gamma+1)\ldots(\gamma+n-1)}

where α,β,γ\alpha, \beta, \gamma are arbitrary real numbers, not integers ≤0\leq 0: it is clear that unu_n is >0> 0 from a certain stage, or <0< 0 from a certain stage on. One has

un+1un=(α+n)(β+n)(n+1)(γ+n)=(1+α+βn+αβn2)(1+γ+1n+γn2)−1=1+α+β−γ−1n+αβ−(α+β)(γ+1)+γ2+γ+1n2+o(1n2). \begin{align*} \frac{u_{n+1}}{u_n} &= \frac{(\alpha+n)(\beta+n)}{(n+1)(\gamma+n)} \\ &= \left(1 + \frac{\alpha+\beta}{n} + \frac{\alpha\beta}{n^2}\right) \left(1 + \frac{\gamma+1}{n} + \frac{\gamma}{n^2}\right)^{-1} \\ &= 1 + \frac{\alpha+\beta-\gamma-1}{n} \\ &\quad + \frac{\alpha\beta-(\alpha+\beta)(\gamma+1)+\gamma^2+\gamma+1}{n^2} + o\left(\frac{1}{n^2}\right). \end{align*}

Raabe's test shows that the series converges for α+β<γ\alpha + \beta < \gamma, and has infinite sum for α+β>γ\alpha + \beta > \gamma; when α+β=γ\alpha + \beta = \gamma the series also has an infinite sum, as is shown by prop. 8.

Remark 1

As a particular instance of Raabe's test one sees that if lim sup⁡n→∞un+1/un<1\limsup_{n \to \infty} u_{n+1}/u_n < 1 the series (un)(u_n) converges; if, on the other hand, lim inf⁡n→∞un+1/un>1\liminf_{n \to \infty} u_{n+1}/u_n > 1, the series (un)(u_n) has an infinite sum (d'Alembert's test).

Remark 2

The criteria of the second kind can be applied only to series whose general term behaves in a very regular way as nn tends to +∞+\infty; in other words, their scope is much more restricted than those of the logarithmic criteria, and it would be a blunder to try to use them beyond the special cases to which they are particularly suited. For example, for the series (un)(u_n) defined by u2m=2−mu_{2m} = 2^{-m}, u2m+1=3−mu_{2m+1} = 3^{-m} one has u2m+1/u2m=(23)mu_{2m+1}/u_{2m} = \left(\frac{2}{3}\right)^m, u2m+2/u2m+1=12(32)mu_{2m+2}/u_{2m+1} = \frac{1}{2} \left(\frac{3}{2}\right)^m; the first of these ratios tends to 0 and the second to +∞+\infty as mm increases indefinitely, so no criterion of the second kind is applicable; nevertheless, since un≼2−n/2u_n \preccurlyeq 2^{-n/2}, it is immediate that the series converges.

Even when un+1/unu_{n+1}/u_n has a simple expression, a direct evaluation of a principal part for unu_n often leads to a result as quickly as the criteria of the second kind. For example, for the hypergeometric series, Stirling’s formula shows immediately that un∼anα+β−γ−1u_n \sim a n^{\alpha + \beta - \gamma - 1}, where aa is a constant ≠0\neq 0, and the logarithmic criterion of order 0 is then applicable.

Exercises

See the exercises for § 4.

Transcribed from the printed volume. Not translated.