Why does the square integral of a Dirac-Delta-function blow up to infinity? In Griffiths Introduction to Quantum Mechanics  on page 102 it is shown, that the eigenfunctions of the position operator $\hat{x}=x$ are not normalizable.
$$\int_{-\infty}^{\infty}g_{\lambda}\left(x\right)^{*}g_{\lambda}\left(x\right)dx
~=~\left|B_{\lambda}\right|^{2}\int_{-\infty}^{\infty}\delta\left(x-\lambda\right)\delta\left(x-\lambda\right)dx
~=~\left|B_{\lambda}\right|^{2}\delta\left(\lambda-\lambda\right).$$
I understand everything up to this point. But now he concludes that 
$$\left|B_{\lambda}\right|^{2}\delta\left(\lambda-\lambda\right)
~\rightarrow~\infty.$$
How comes?
References:
D.J. Griffiths, Introduction to Quantum Mechanics, (1995) p. 102.
(Link: http://www.fisica.net/quantica/Griffiths%20-%20Introduction%20to%20quantum%20mechanics.pdf)
 A: Let's accept the last step in your "derivation". A squared delta already is an hint of serious problems.
Simply put, the expression $\delta(0)$ makes little sense. Consider the $\delta(x)$ "function" (it is a distribution, but I'll go with the simplest answer possible, just for intuition) and a succession approximating it:
$$
\delta(x-\lambda)=\lim_{\epsilon\to0}\frac{1}{\sqrt{2\pi}\epsilon}e^{\frac{(x-\lambda)^2}{2\epsilon^2}}
$$
Now, you can see that, if you "evaluate" this limit for any $x$ that is non $\lambda$, you get 0: the exponential is always stronger than the power in the denominator, and dictates the limit.
But, if you take $x=\lambda$, as in your case, the exponential disappears: nothing stops the denominator from bringing the whole limit to infinity. So, you can say that $\delta(\lambda-\lambda)=\delta(0)=\infty$. So, the position eigenstate does not have a norm.
A: The way Griffiths demonstrates the problem is purely formalistic; he uses the usual rules defining the behaviour of $\delta$ distribution and applies them to an unusual case: an integral of product of two same delta distributions
$$
\int \delta(x-\lambda)\delta(x-\lambda)\,dx.
$$
The usual rules lead to result
$$
\delta(0)
$$
which is $the~problem$ - the delta distribution is defined through its action on test functions, but $\delta(0)$ is not such delta distribution, because it has no variable argument left. It looks like value of a function $\delta$, only $\delta$ has no values, as it is not a function. So the expression $\delta(0)$ is undefined. It is not true that $\delta(0)$ is infinite by definition, at least not by the usual mathematical definition. In mathematics, $\delta(0)$ is just undefined. This is not a problem, because $\delta$ distribution is meant to be used as linear functional, not as a function to be evaluated.
A: In the traditional Physics literature, we find $\delta$ be defined as
$$\delta(x)=\begin{cases}0, & x\neq 0, \\ \infty, & x = 0\end{cases}$$
and then one "derives" that if $f$ is a function we have
$$\int_{-\infty}^\infty f(x)\delta(x-a)dx=f(a).$$
The problem is that $\delta$ doesn't make sense as a function. Indeed one can show that if one tries to define $\delta$ by this integral, there is no function that will make this equation true.
It turns out that this can be handled by the introduction of the theory of distributions. The idea is: define $\delta$ as something that acts on a function, giving the result of the above integral. In other words: make $\delta$ a linear functional over some functions.
The appropriate setting, at first, to do this, is to consider $\mathcal{D}(\mathbb{R})$ the so-called space of test functions, being those that are $C^\infty$ and vanish outside a closed and bounded interval. Consider now $\mathcal{D}'(\mathbb{R})$ the space of all continuous linear functionals over $\mathcal{D}(\mathbb{R})$. We denote the action of $\xi\in \mathcal{D}'(\mathbb{R})$ on the test function $f$ by $\xi[f]=\langle \xi,f\rangle$.
This means this is the space of all $\xi : \mathcal{D}(\mathbb{R})\to \mathbb{R}$ such that given a sequence of test functions $f_n\in \mathcal{D}(\mathbb{R})$ convergent to $f$, the sequence of numbers $\langle \xi,f_n\rangle$ converges to $\langle \xi,f\rangle$.
Now the Dirac Delta is defined as the distribution $\langle \delta, f\rangle = f(0)$ and the Dirac Delta centered at $a\in \mathbb{R}$ is $\langle \delta_a, f\rangle = f(a)$. It turns out that $\delta$ is what traditional literature calls $\delta(x)$ and $\delta_a$ is what is called $\delta(x-a)$.
The fact is that functions define distributions by means of taking the integral against the function. Also most operations you do on functions like differentiation and taking the Fourier transform, can be lifted to distributions, such that they coincide with the usual notions when restricted to functions. Thus distributions are sometimes called "generalized functions", even though they are not functions.
As you see $\delta(0)$ is not defined. Actually $\delta(x)$ is not defined, because $\delta$ isn't a function in the traditional sense, it doesn't act on numbers. It acts on functions. The point is that $\delta(\lambda - \lambda) = \delta(0)$ is so undefined as $0/0$. It is not infinite, it is just something that doesn't make sense. When you see $\delta(x)$ or $\delta(x-a)$ it is just a notation, to an object that just makes sense when applied to a function.
