Note: Exponential Integration

Steve Trettel

|

Here’s a quite long calculation showing that it’s possible to integrate exponential functions directly from first principles. The length of this calculation alone is a good selling point for the fundamental theorem of calculus, which will make this trivial!

(Exponential Functions) An exponential function is a continuous nonconstant function E:RR satisfying the law of exponents E(x+y)=E(x)E(y) for all x,yR.

This is a functional characterization of exponentials; we don’t specify a formula for how to compute the function but instead we specify how the function ought to behave. There are several facts about exponentials we will need, that follow directly from this definition:

  • Exponentials are always nonzero
  • Exponentials are strictly increasing, or strictly decreasing
  • Exponentials are differentiable everywhere

Our main theorem here is that we compute the Darboux integral of an exponential on an interval, directly from this functional definition. To do so, we make use of minimal facts about integration, so that this can be given early on. One thing we do take advantage of is the ability to compute Darboux integrals with a sequence of partitions (instead of directly from the definition, with infima and suprema)

  • If Pn is a sequence of shrinking partitions of J and limU(f,Pn)=limL(f,Pn) then f is integrable on J and Jf is this common value.

Without further ado, here’s the advertised calculation!

(Integrating Exponentials) Let E be an exponential function, and [a,b] an interval. Then E is Darboux integrable on [a,b] and [a,b]E=E(b)E(a)E(0)

We will show the argument for E an increasing exponential (its base E(1)>1): an identical argument applies to decreasing exponentials (only switching U and L in the computations below).

To show E(x) is integrable, we use @thm-compute-integral-sequence, which assures us it is enough to find a sequence Pn of shrinking partitions where limL(f,Pn)=limU(f,Pn). Indeed - for each n, let Pn denote the evenly spaced partition of [a,b] with widths Δn=(ba)/n Pn=a,a+Δn,a+2Δn,,a+nΔn=b

We will begin by computing the lower sum. Because E is continuous, it achieves a maximum and minimum value on each interval Pi=[ti,ti+1]. And, since E is monotone increasing, this value occurs at the leftmost endpoint. Thus,

L(E,Pn)=0i<ninfPiE(x)|Pi| =0i<nE(ti)Δn =0i<nE(a+iΔn)Δn

Using the law of exponents for E we can simplify this expression somewhat:

E(a+iΔn)=E(a)E(iΔn) =E(a)E(Δn+Δn++Δn) =E(a)E(Δn)E(Δn)E(Δn) =E(a)E(Δn)i

Plugging this back in and factoring out the constants, we see that the summation is actually a partial sum of a geometric series:

0i<nE(a+iΔn)Δn=0i<nE(a)E(Δn)iΔn =E(a)Δn0i<nE(Δn)i

Having previously derived the formula for the partial sums of a geometric series, we can write this in closed form:

0i<nE(Δn)i=1E(Δn)n1E(Δn)

But, we can simplify even further! Using again the laws of exponents we see that E(Δn)n is the same as E(nΔn), and nΔn is nothing other than the width of our entire interval, so ba. Thus the numerator becomes 1E(ba), and putting it all together yields a simple expression for L(E,Pn):

L(E,Pn)=E(a)Δn1E(ba)1E(Δn)

Some algebraic re-arrangement is beneficial: first, note that by the laws of exponents we have

E(a)(1E(ba))=E(a)E(ba)E(a) =E(a)E(b)

Thus for every n we have

L(E,Pn)=(E(a)E(b))Δn1E(Δn)

We are interested in the limit as n: by the limit laws we can pull the constant E(a)E(b) out front, and only concern ourselves with the fraction involving Δn. There’s one final trick: look at the negative reciprocal of this fraction:

1Δn1E(Δn)=E(Δn)1Δn

Because we know E(0)=1 for all exponentials, this latter term is none other than the difference quotient defining the derivative for E! Since we have proven E to be differentiable, we know that evaluating this along any sequence converging to zero yields the derivative at zero. And as Δn0 this implies

limE(Δn)E(0)Δn=E(0)

Thus, our original limit Δn/(1E(Δn)) is the negative reciprocal of this, and

limL(E,Pn)=lim(E(a)E(b))Δn1E(Δn) =(E(a)E(b))limΔn1E(Δn) =(E(a)E(b))1E(0) =E(b)E(a)E(0)

Phew! That was a lot of work! Now we have to tackle the upper sum. But luckily this will not be nearly as bad: we can reuse most of what we’ve done! Since E is monotone increasing, we know that the maximum on any interval occurs at the rightmost endpoint, so

U(E,Pn)=0i<nsupPiE(x)|Pi| =0i<nE(ti+1)Δn =0i<nE(a+(i+1)Δn)Δn

Comparing this with our previous expression for L(E,Pn), we see (unsurprisingly) its identical except for a shift of ii+1. The law of exponents turns this additive shift into a multiplicative one:

U(E,Pn)=0i<nE(a+(i+1)Δn)Δn =0i<nE(Δn)E(a+iΔn)Δn =E(Δn)0i<nE(a+iΔn)Δn =E(Δn)L(E,Pn)

Thus, U(E,Pn)=E(Δn)L(E,Pn) for every n. Since E is continuous, limE(Δn)=E(limΔn)=E(0)=1

And, as L(E,Pn) converges (as we proved above) we can apply the limit theorem for products to get

limU(E,Pn)=lim(E(Δn)L(E,Pn)) =(limE(Δn))(limL(E,Pn)) =limL(E,Pn) =E(b)E(a)E(0)

Thus, the limits of our sequence of upper and lower bounds are equal! And, by the argument at the beginning of this proof, that squeezes L(E) and U(E) to be equal as well. Thus, E is integrable on [a,b] and its value is what we have squeezed:

[a,b]E=E(b)E(a)E(0)

Recalling our definition of exp as the exponential function with exp(0)=1 (in Part II on differentiation), we have

On any interval [a,b] the natural exponential is integrable, and [a,b]exp=exp(b)exp(a)