Some while ago, I showed a slightly dicey proof of the Basel Problem identity, n=11n2=π26, and invited readers to share other proofs with me.

My old friend Jean Reinaud stepped up to the mark with an exercise from his undergraduate textbook:

Question 1.7

The French isn’t that difficult, but just in case:

Determine real numbers α and β such that 0π(αt+βt2)cos(nt)dt=1n2, for any integer n>0.

Hence deduce n=11n2=π26

As usual, I suggest having a go at this yourself before reading on. Below here be spoilers.


Setting up the integral

The first part isn’t too tricky: integrating by parts and noting that sin(nπ)=sin(0)=0 for integer n means the integral evaluates to (a+2πb)cos(πn)an2.

We need the top to be 1, and we split it into two cases: n even and n odd (since cos(πn) takes on different values in those cases.)

In the odd case, cos(πn)=1, so a+2πba=1, giving b=12π.

In the even case, cos(πn)=1, so a2πba=1. However, 2πb=1, so 2a1=1, so a=1.

Our integral is 0π(t22πt)cos(nt)dt.

The tricky bit: the infinite sum

So, to get our result, presumably we just need to sum all of those integrals from 1 to , and boom 1! Out should pop our result.

The brackety bit doesn’t have anything to do with n, so we can write our summed integral as:

0π(t22πt)n=1cos(nt)dt.

Evaluating that sum, we hit a problem: it doesn’t converge for any value of t. For example, for t=π, cos(nt) oscillates between 1 and 1 - and that sum is undefined. Even worse, for t=0, cos(nt) is always 1, so we end up with a sum that goes to infinity. It’s a disaster!

But.

We can engage in a bit of jiggery-pokery with the partial sums, and show that they converge to the right thing.

Let’s start with n=0Ncos(nt). Depending on how interested you are, you can:

  • a) Express each cosine as eint+eint2, add up the resulting geometric series and play about until you get the result;
  • b) Look it up on Wolfram Alpha; or
  • c) Look it up on Wikipedia.

All three, thankfully, give the same result: n=0Ncos(nt)=12+sin((N+12)t)2sin(t2).

Well, that’s a mess. However, it does reduce a pretty substantial sum to something that only depends on N and t.

Evaluating that mess

But how to evaluate that, combined with the quadratic t factor? We’re after:

0π(tt22π)[12+sin((N+12)t)2sin(t2)]dt

Let’s split that hairy square bracket into two parts: the 12 bit will be easy to integrate, but 0π(tt22π)[sin((N+12)t)2sin(t2)]dt… less so. So much less so, I suggest calling it IN.

Jean’s textbook gives us a hint: cot(t)<1t<cosec(t) over the domain we care about. In particular, [sin((N+12)t)2sin(t2)]>[sin((N+12)t)t]

That gives us a lower bound for IN:

IN=0π(tt22π)[sin((N+12)t)2sin(t2)]dt>0π(tt22π)[sin((N+12)t)t]dt

Pushing the t in the denominator back into the quadratic factor, we get:

IN>0π(1t2π)sin((N+12)t)dt.

Stick it in the parts formula, turn the crank, note that sin(Nπ)=0 for integer N and out pops IN>2cos(Nπ)+2(1+2N)π(1+2N)2π.

This is good: the lower bound goes to zero as N gets large.

Great. However, that’s just a lower bound. What about an upper bound?

Can we express it another way?

Of course we can.

Let’s come back to here:

sin((N+12)t)2sin(t2)

… and expand the sine in the numerator. We get sin((N+12)t)=sin(Nt)cos(12t)+cos(Nt)sin(12t).

Of course, there’s a sin(12t) in the denominator as well, so we can simplify to get sin(Nt)cot(12t)+cos(Nt). Suddenly the future is brighter.

… at least momentarily. IN=0π(tt22π)(sin(Nt)cot(12t)+cos(Nt))dt… is not Wolfram-friendly. However, Jean’s textbook comes to the rescue again: we know that cot(12t)<2t.

So, IN<0π(tt22π)(2sin(Nt)t+cos(Nt))dt. If I’ve done my sums right, that’s 12N+Ncos(Nπ)N2 - which gives us an upper bound that goes to 0 as N gets large.

So: IN is no bigger than an upper bound that goes to 0 in the limit, and no smaller than a lower bound that also goes to zero in the limit. By the squeeze theorem, as N gets large, IN goes to 0.

Phew.

Back to the integral

0π(tt22π)[12+sin((N+12)t)2sin(t2)]dt

Everything after the 12+ can be ignored - we just proved that bit of the integral goes to 0 in the limit.

Meanwhile, 120π(tt22π)dt=π26, as required .

* If you’ve enjoyed this post, and want to do something kind: please sign up to be an organ donor. Without someone like you saying “I want to help someone live after I die,” Jean would not be with us today. I’m on the register - I’d be delighted if you joined me.

Footnotes:

1. Avoiding surprise factorial