This article examines the famous Basel problem.


The Basel problem was the problem to evaluate the sum


in other words, \zeta(2) . First, we will see why it must converge here, using the integral test: If the series converge, the integral of the general term must converge. Let's test:

\lim_{x\to\infty}\int x^{-2}dx=\lim_{x\to\infty}\dfrac{x^{-1}}{-1}+C=\lim_{x\to\infty}-\frac{1}{x}+C=C

Hence the series converge. The real problem is, to what number?

Euler's approach

Leonhard Euler was the first one to develop a clever approach to "solve" the problem, though not as rigorously as required. We will first demonstrate this approach and then give a rigorous proof.

The power series for the sine function is


so this holds:


which we can express more clearly as


It is obvious that the coefficient of x^2 is 1/6. We will now obtain this coefficient through a different method. Obviously, the LHS function has zeros at every integer that is a multiple of pi, except 0. Through polynomial factoring, we obtain that we can show this as


Using the difference of two squares formula, we simplify this to


Rearranging factors, we get


Now, with some simple algebra we obtain that the coefficient of x^2 in this expression is


Since if two polynomials are equal, all coefficients must be equal; and this one is equal to the power series expansion, we get


Multiplying thoroughly by \pi^2 , we obtain the famous solution to the Basel problem:


However, as we stated, this proof is not rigorous enough, whilst it is still true: The result is correct.

A rigorous proof

Now, since we know the answer, we are ready to give a rigorous proof. Using Fourier analysis, we know that for a function that is periodic over the interval (-\pi,\pi) with period 2\pi , the Fourier series are (please note that the function must be continous on the open interval and the left hand side aswell as the right hand side derivates must exist and be finite at every point of the open interval if we want equality between the function and its fourier series. If we, in additional, want equality to hold between the function and the fourier series on the whole interval [-\pi,\pi] then it must also hold that f(\pi)=f(-\pi)

f(x)=\dfrac{a_0}{2}+\sum_{r=1}^\infty a_r\cos(rx)+\sum_{r=1}^\infty b_r\sin(rx)


a_0&=\dfrac{1}{\pi}\int\limits_{-\pi}^\pi f(x)dx\\
a_r&=\dfrac{1}{\pi}\int\limits_{-\pi}^\pi f(x)\cos(rx)dx\\
b_r&=\dfrac{1}{\pi}\int\limits_{-\pi}^\pi f(x)\sin(rx)dx

We start out by finding the Fourier series for x^{2n} :

a_0=\dfrac{1}{\pi}\int\limits_{-\pi}^\pi x^{2n}dx=\frac{2\pi^{2n}}{2n+1}

Since the function is even, the sine coefficient will vanish. Hence we just need to evaluate

a_r=\dfrac{1}{\pi}\int\limits_{-\pi}^\pi x^{2n}\cos(rx)dx

First, we will evaluate the indefinite integral. If

F(n)=\int x^{2n}\cos(rx)dx

then, we obtain the functional equation


via integration by parts. By induction or simple observation, this functional equation yields the following formula:


The sine term vanishes in the Fourier series because we will use integer multiples of pi in the sine function, and that yields zero, and the cosine function yields (-1)^r . Hence, we are left with

G_n(r)=\int\limits_{-\pi}^\pi x^{2n}\cos(rx)dx=2(-1)^r\sum_{k=0}^{n-1}\dfrac{(-1)^k(2n)!\pi^{2n-2k-1}}{(2n-2k-1)!r^{2k+2}}

Substituting these into our original Fourier series, we get

x^{2n}=\dfrac{\pi^{2n}}{2n+1}+\frac{1}{\pi}\sum_{r=1}^\infty G_{n}(r)\cos(rx)

We now put n=1 and x=\pi to obtain

\pi^2=\dfrac{\pi^2}{3}+\frac{1}{\pi}\sum_{r=1}^\infty G_1(r)(-1)^r

The general term of the sum simplifies to


Substituting back:


The value we want to evaluate just appeared. We set it to C:


Solving this equation:


and this completes the rigorous proof.