Modelling Dirac Delta Function and Proving its Area Property I have modeled the Dirac delta function using limiting value of the rectangular pulse and now proving its area property. In the proof, I need to interchange the order of integration and limit, which does not look possible. I looked out at the similar problems on math stackexchange but could not find it anywhere. The details about the problem are as follows:

To model the delta function, first I took the rectangular pulse function as:

\begin{equation}
p_n(t) = \frac{1}{n} \ \textrm{for}\ |t| \le \frac{n}{2}\ \textrm{otherwise}\ 0
\end{equation}

Based on this rectangular pulse, we can define the delta function as follows:

\begin{equation}
\delta(t) = \lim_{n \rightarrow 0}{p_n(t)}
\end{equation}

So the above equation precisely captures the behavior of the delta function. 

Now, I want to prove the area property of the delta function, which is:

\begin{equation}
\int_{-\infty}^{+\infty}{\delta(t) dt} = 1
\end{equation}

By writing the above equation using the definition of the delta function, it results into:

\begin{equation}
\int_{-\infty}^{+\infty}{\lim_{n \rightarrow 0}{p_n(t)}}dt = 1
\end{equation}

Next, if we look at the integrand of the above integral, it is divergent as the function $\lim_{n \rightarrow 0}{p_n(t)}$ has infinite value, so overall its not integrable. Now, I want to swap the limit and integral operator, which is only possible if the limit is convergent (where it is divergent in this case). 

My question is that is there really a way I can swap the order of this limit and integral? If yes then how would I do it?

Note: I also checked out many questions/answers regarding this on math stackexchage, but in all cases they directly take limit outside the integral, i.e., $\lim_{n \rightarrow 0}{\int}_{-\infty}^{+\infty}{p_n(t)} dt$, and it would only be equal to $\int_{-\infty}^{+\infty}{\delta(t) dt}$, if the limits and integral operators can be swapped. Moreover, the integral that I am using is the gauge integral. 

I would highly appreciate your earliest response on it.

Thanks in advance.
Since the support of $\delta $ is a point, the natural interpretation of the gauge integral of $\delta $ is that it should be zero. To get a gauge integral you need to calculate Riemann sums; for an interval $(u_i,u_{i+1}) $, you need to calculate $\delta (t_i)\, (u_{i+1}-u_i) $. You now want to get $\infty\times (u_{i+1}-u_i)=1$, and there's no arithmetic that will give you that.
So, the question about exchanging limit and integral does not really apply, since $\delta $ is not actually a function.
The usual approach to avoid this problem is the following. The $\delta$ is never used on its own. It used in the sense that $$\tag1\int_{\mathbb R} f(t)\,\delta(t)\,dt=f(0)$$ for any (good enough, say) function $f$.
One way to achieve $(1)$ is to consider the left-hand-side as the limit
$$
\lim_n \int_{\mathbb R} f(t) p_n(t)\,dt.
$$
This is still not ideal because you are subjecting yourself to the choice of the $p_n$, and in the end what one cares about is the net result, which is the equality $(1)$. The solution to this is to be more abstract, and consider $\delta$ as a linear functional, that to every $f$ assigns the number $f(0)$.
A second approach, that agrees with the previous one, is to consider $\delta$ as a measure, namely $$\delta(E)=\begin{cases} 1,&amp;\ 0\in E\\ \ \\ 0,&amp;\ 0\not\in E\end{cases}$$ Then one can write
$$\tag2
\int_{\mathbb R} f(t)\,d\delta(t)=f(0)
$$
in the sense of Lebesgue integration.

Edit: why $\delta$ &quot;function&quot;?
On the set $C_0(\mathbb R)$, of continuous functions vanishing at infinity, one considers (bounded) linear functionals, which are linear maps $\varphi:C_0(\mathbb R)\to\mathbb R$. By the Riesz-Markov theorem, for any such $\varphi$ there exists a measure $\mu_\varphi$ such that
$$\tag3
\varphi(f)=\int_{\mathbb R} f(t)\,d\mu_\varphi(t),\ \ \ \ f\in C_0(\mathbb R).
$$
If $\mu_\varphi$ is absolutely continuous with respect to Lebesgue measure, by Radon-Nikodym there exists a function $g_\varphi$ such that
$$\tag4
\varphi(f)=\int_{\mathbb R} f(t)\,g(t)\, dt, \ \ \ \ f\in C_0(\mathbb R). 
$$
Dirac needed to use the functional $\psi(f)=f(0)$. Doing things &quot;the physicist way&quot;, he assumed a $g$ like in $(4)$ would somehow exist, even if not a function, he called it $\delta$ and wrote $(1)$.
Mathematically, what's happening is $(2)$. There is no &quot;$\delta$-function&quot;.
Just because the function under the integral contains an infinite value, it does not mean that the function is not integrable. A simple example would be $f(x)=1/\sqrt{|x|}$. When $x\rightarrow 0$ $f$ goes to infinity. You can however integrate in an interval around zero. You can integrate between $-1$ and $1$, and you get $4$.

Since the integration limits do not change with $n$, you can safely move the limit outside of the integration.