Lambdia

Both Term Tests Return One, and the Sum Still Stops Below Two

The ratio and root tests both return 1 on the sum of e to the minus root n, which settles nothing, and the usual write-up of the problem then quotes 4 over e as the answer. That number is the floor rather than the cap: the sum is 1.6704068, which is 13.52 percent above it, and the usable bound comes from integrating from 0 instead of from 1, giving exactly 2. Where each bar of width one sits relative to its index is the single step that decides which way the inequality points.

Add up ene^{-\sqrt{n}} for every whole number n1n \ge 1. The terms shrink, but slowly: the hundredth is e104.5×105e^{-10} \approx 4.5 \times 10^{-5} and the ten-thousandth is still e100e^{-100}, which is small but arrived at unhurriedly. Does the total stay finite?

It does, and it lands at 1.6704068. The interesting part is that the two tests most people reach for first both refuse to answer, and that the standard write-up of this problem gets the bound backwards.

Two tests, no verdict

The ratio test looks at consecutive terms:

an+1an=enn+1=exp ⁣(1n+1+n)1\frac{a_{n+1}}{a_n} = e^{\sqrt{n} - \sqrt{n+1}} = \exp\!\left(\frac{-1}{\sqrt{n+1} + \sqrt{n}}\right) \longrightarrow 1
(1)

The exponent goes to zero like 1/(2n)1/(2\sqrt{n}), so the ratio approaches 1 from below and never gets there. The root test does no better: the nn-th root of the nn-th term is e1/ne^{-1/\sqrt{n}}, which also tends to 1. A limit of 1 is the one value at which both tests are silent by construction, and silence is not divergence.

Nor does the shrinking of the terms settle anything. The harmonic series has terms going to zero and diverges anyway, so an0a_n \to 0 is a requirement rather than an argument. At this point the instruments have run out, and the move is to stop weighing individual terms and start measuring area instead.

Standing the terms up as bars

Consider the function f(x)=exf(x) = e^{-\sqrt{x}} on the positive axis. It is decreasing everywhere, which the derivative confirms:

f(x)=ex2x<0(x>0)f'(x) = -\frac{e^{-\sqrt{x}}}{2\sqrt{x}} < 0 \qquad (x > 0)
(2)

Now draw a bar for each term, one unit wide, with height f(n)=enf(n) = e^{-\sqrt{n}}. Its area is the term itself, so the sum of the series is the total area of the bars. The question is where to put each bar, and this is the step that decides everything.

Put the bar for term nn on the interval [n1,n][n-1, n]. Its height is the value of ff at the right end of that interval, and since ff is decreasing, the right end is where ff is smallest there. So the bar sits entirely under the curve, every bar does, and the bars do not overlap:

n1en    0exdx=2\sum_{n \ge 1} e^{-\sqrt{n}} \;\le\; \int_0^{\infty} e^{-\sqrt{x}}\,dx = 2
(3)

Put the same bar on [n,n+1][n, n+1] instead and the height becomes the value at the left end, which is the largest value on that interval, so each bar now covers the curve rather than hiding under it. That comparison runs the other way:

n1en    1exdx=4e=1.4715178\sum_{n \ge 1} e^{-\sqrt{n}} \;\ge\; \int_1^{\infty} e^{-\sqrt{x}}\,dx = \frac{4}{e} = 1.4715178
(4)
Fig. 1 — The same bars, moved one unit. On the left they are bounded by the area from zero, on the right they bound the area from one. The two inequalities point in opposite directions.

This is where the usual presentation of the problem goes wrong. It computes the integral from 1, gets 4/e4/e, and reports that the sum is just under that. The sum is not under 4/e4/e. It is 13.52 percent above it, because the integral from 1 is a floor. The cap you want is the integral from 0, and it happens to be a rounder number.

The integral test, with both directions kept

Let ff be positive and decreasing on [0,)[0, \infty). Then the bars on [n1,n][n-1,n] give n1f(n)0f\sum_{n \ge 1} f(n) \le \int_0^{\infty} f, and the bars on [n,n+1][n, n+1] give n1f(n)1f\sum_{n \ge 1} f(n) \ge \int_1^{\infty} f. Both inequalities need monotonicity and neither one is reversible. Convergence of the integral and convergence of the sum stand or fall together, but the numerical bounds are different objects and sit on opposite sides of the answer.

Evaluating both integrals

Substitute u=xu = \sqrt{x}, so x=u2x = u^2 and dx=2ududx = 2u\,du. The improper integral becomes an ordinary one:

0exdx=02ueudu=2Γ(2)=2\int_0^{\infty} e^{-\sqrt{x}}\,dx = \int_0^{\infty} 2u\,e^{-u}\,du = 2\,\Gamma(2) = 2
(5)

One integration by parts gives ueudu=(u+1)eu\int u e^{-u} du = -(u+1)e^{-u}, and the same antiderivative evaluated from u=1u = 1 supplies the other integral, 22/e=4/e2 \cdot 2/e = 4/e. It also gives the tail in closed form, which is what makes a numerical check honest:

Nexdx=2(1+N)eN\int_N^{\infty} e^{-\sqrt{x}}\,dx = 2\left(1 + \sqrt{N}\right)e^{-\sqrt{N}}
(6)

At N=4,000,000N = 4{,}000{,}000 that tail is below 1030010^{-300}, so the partial sum to four million terms, 1.6704068180, brackets the true value to a width no arithmetic will ever notice. Sharpening (4) by adding the first term to it gives the classical sandwich 4/eSe1+4/e4/e \le S \le e^{-1} + 4/e, that is 1.4715S1.83941.4715 \le S \le 1.8394, which is tighter than the cap of 2 and less quotable.

Fig. 2 — The sum and its bounds. The number the standard write-up quotes as the answer is the floor, and it is 13.52 percent below the truth.

The same argument on a whole family

Nothing in the derivation used the square root beyond its monotonicity, so replace it with any positive power. For α>0\alpha > 0, the substitution u=xαu = x^{\alpha} turns the integral into a gamma function:

0exαdx=1αΓ ⁣(1α)=Γ ⁣(1+1α)\int_0^{\infty} e^{-x^{\alpha}}\,dx = \frac{1}{\alpha}\,\Gamma\!\left(\frac{1}{\alpha}\right) = \Gamma\!\left(1 + \frac{1}{\alpha}\right)
(7)

So enα\sum e^{-n^{\alpha}} converges for every positive α\alpha, with cap Γ(1+1/α)\Gamma(1 + 1/\alpha). At α=1/2\alpha = 1/2 that is Γ(3)=2\Gamma(3) = 2, which is where the 2 in (3) comes from. Both term tests still return 1 for every α<1\alpha < 1, so the whole family is invisible to them and visible to one integral.

Where the method needs care

The monotonicity in (2) is not decoration, and it is easy to see what it is holding up. Take g(x)=excos20(πx)g(x) = e^{-\sqrt{x}}\cos^{20}(\pi x). Since cos(πn)=±1\cos(\pi n) = \pm 1, this agrees with our terms at every single integer, so the sum of g(n)g(n) is still 1.6704068. But gg dips almost to zero between the integers, and its integral from 0 to infinity is about 0.366, well below the sum. A bound read off that curve would be flatly false. The sampled values do not control the area on their own; a decreasing curve is what ties the two together, because it attains its extremes at the ends of each interval.

The cap from 0 also requires the integrand to behave at the origin. Here f(0)=1f(0) = 1 and there is nothing to worry about, but for a function with a singularity at 0, say f(x)=ex/xf(x) = e^{-\sqrt{x}}/x, the integral from 0 diverges and the only usable bound is f(1)+1ff(1) + \int_1^{\infty} f. The version of the integral test that survives that case is the sandwich, not the clean cap.

One alternative is often suggested and worth pricing honestly. Comparing with 1/n21/n^2 does work, because en1/n2e^{-\sqrt{n}} \le 1/n^2 exactly when n2lnn\sqrt{n} \ge 2\ln n, and the difference n2lnn\sqrt{n} - 2\ln n is increasing past n=16n = 16 and first turns positive at n=75n = 75. So the comparison settles convergence, and it says nothing at all about the first seventy-four terms. You have to add those up by hand before the argument produces a number, which is why the integral is the route worth remembering.

Sources and further reading

Both bar placements were checked term by term over the first two hundred indices and against fine Riemann sums of both integrals, the two term-test limits were evaluated out to n=1012n = 10^{12}, and the same machinery was run on the harmonic series as a control, where the integral diverges and the partial sums track lnn\ln n as they should.

Commentaires · 0

Sois le premier à commenter.