The integral is an example of an* improper* integral: The *integrand* is not defined at 0, and the interval of integration is unbounded. The issue at 0 is not really serious, since

Since the interval is unbounded, this integral is actually defined as a limit,

In this post I want to show how to evaluate this limit. The problem is that the **sine integral**

does not admit an *elementary* expression, i.e., it cannot be expressed in terms of by composing exponentials, logarithms, trigonometric functions, their inverses, and polynomials.

(To get a sense of the difficulty in finding an expression for , it may be useful to play a little trying integration by parts or similar tricks.)

This is a classical integral and appears in applications with some frequency, so there are several standard methods of computing its value. Usually, these methods involve techniques from complex analysis (either the Fourier transform or contour integrals), although there are some elementary approaches as well, such as the one I illustrate here. “Elementary” is to be taken with a grain of salt; although the argument below is easy to follow with some patience, it is difficult to come up with. There are good references where the integral is computed, but it usually takes me a while to remember where to look, so I decided to write this post and to provide a reference here as well.

An excellent book where the computation can be found more or less in complete detail (except for what I call Claims 1 and 2 below) is** Introduction to Calculus and Analysis, vol. I**, by Richard Courant and Fritz John, Interscience Publishers (1965), reprinted by Springer (1989). The argument I present here is essentially the same they use, although I have arranged some of the details in a different way.

I will show that

Let’s begin by showing that the integral converges. Although we do not need this, I’ll also show that the convergence is *conditional *but not *absolute*, i.e., the integral diverges.

**Fact 1.** converges.

**Proof.** Let . Notice that .

Notice also that , since (as )

This proves that as any is in for a unique , and

So .

Now observe that for all , since for any , and therefore

so

Also, the alternate in sign (being positive for even values of and negative for odd values). It follows (by Leibniz’s test) that converges, and therefore converges as well.

Continuing with the argument just shown, we also see that, since the terms alternate, then (for any )

showing that the integral converges relatively fast, which helps explain the shape of the graph displayed above.

**Fact 2.** diverges.

**Proof.** It suffices to see that diverges. For this we now look for an estimate bounding from below. For this, notice that for all , so

and , which diverges.

From the argument above, we also have that for any ,

so, although divergent, the function grows as at a similar speed as the harmonic series (which grows very slowly).

Having established both the convergence of the integral and found a bound for how quickly it converges, we now proceed to find its value.

For this, note that for any constant , the change of variables shows that

The key trick that allows us to compute the integral is the observation that

This is not particularly hard, but the algebra is a tad messy, so I will postpone its proof until the end.

Let’s assume for now that this is the case, so

To compute this limit, we make use of a nice trigonometric identity:

**Fact 3.** Let . Then

for , where for the right hand side of the displayed equation is understood as its limit as .

**Proo****f.** Let’s begin by verifying that the result holds for . In this case, the definition of gives . Also, using l’Hôpital’s rule,

as required. Let’s now argue that the result also holds for .

When , the right hand side equals , so the result holds. We argue by induction, showing that if is as stated, then the corresponding formula also holds for .

First, notice that which (upon expanding) equals , so

(1)

Also, or

(2)

Adding (1) and (2) we get

that (dividing both sides by ) gives us

which, by the assumption that the result holds for , gives us that satisfies the desired relation as well.

**Remark.** For there is a much more natural proof of Fact 3 using complex numbers: Simply notice that , where and denotes the real part of the complex number Hence,

since the expression being added is a (finite) geometric series. Now: , so

Also, which we can expand as , so

and the result follows.

It follows from Fact 3 that

since for all we have

With this computation, we can finally conclude that

All that remains is to verify the following:

**Fact 4.** .

**Proof.** Let for . I first argue that can be defined at in a way that makes it continuous at this point (and therefore it is continuous in the whole interval ). I then argue that (extending this way) is also defined and continuous. I use this to find an expression for the integral by means of integration by parts, and the limit of this expression is easy to compute, giving the desired result.

**Claim 1.** With as above,

- , so if we set , the new function is continuous in all of .
- exists and equals .

**Proof.** We use l’Hôpital’s rule (repeatedly), and to keep the argument readable, write to indicate that .

We could prove items 1 and 2 separately, but for brevity’s sake I will prove both of them at the same time. It suffices to show that , because (in particular) this limit can only exist if , which proves 1, but it also proves 2, since then the limit can be rewritten (by writing instead of ) as

which is precisely the definition of .

Now,

This completes the proof of Claim 1.

There are of course other ways of proving Claim 1. For example, one could compare the power series expansions about 0 of and .

**Claim 2.** , so is continuous in all of .

**Proof**. For , Using the notation as in the previous claim, we see that

This completes the proof of Claim 2.

(As with Claim 1, we could have argued here by considering power series expansions.)

The point of Claims 1 and 2 is twofold: They allow us to use integration by parts to find an expression for

and they allow us to show easily that the limit of this expression (as ) is 0 as required.

Since both and are continuous on , we can integrate by parts the expression above, setting and , obtaining because .

Now: Since is continuous on , is bounded by some constant , so and therefore

as , as wanted.

This completes the proof of Fact 4.

And with this, we are done.

Nice explanation! Thanks a lot!