The Dirac delta function appears in the Sokhotsky formula,
$$\text{Im}\lim_{\epsilon\to 0^+} \frac{1}{x-i\epsilon} = \pi\delta(x),$$
to be understood in the integral sense
$$\text{Im}\lim_{\epsilon\to 0^+} \int \frac{f(y)}{y-x-i\epsilon}dy=\pi f(x),$$
for a real valued function $f(x)$.
I stumbled on an identity that has a similar flavour,
$$\lim_{\epsilon\to 0^+}\int_x^b \frac{\epsilon f(y)}{(y-x)^{1-\epsilon}} dy=f(x).\label{1}\tag{$\ast$}$$
The upper limit $b>x$ of the integral is arbitrary, one may send it to infinity if $f(x)$ has compact support.
A corollary is
$$
\lim_{\epsilon\to 0^+}\int_a^b \frac{\epsilon f(x)}{[(b-x)(x-a)]^{1-\epsilon}}\,dx=\frac{f(a)+f(b)}{b-a}.$$
All of this can be interpreted as a delta function representation in terms of the unit step function $\theta(x)$,
$$\lim_{\epsilon\to 0^+} \frac{\epsilon\theta(x)}{x^{1-\epsilon}}=\delta(x),\tag{$\ast\ast$}$$
acting on compactly supported functions.
Q: One can readily check the formula \eqref{1} for polynomial functions $f(x)$. Is there a more comprehensive derivation? Is this representation of the delta function known?
Best Answer
As usual in such examples, there is no need to integrate against a test function. One can simply use the fact that if a sequence (or net) of distributions converges in the distributional sense, then so does the one obtained by differentiating term by term. In particular, this applies when the sequence consists of functions which converge in pretty well any sensible classical sense, e.g., locally $L^1$ as in the case in point, that of the functions which are defined as $x^{\epsilon}$ on the positive real axis, and $0$ elsewhere. They converge to the Heaviside function and we can differentiate to obtain the required result.
Most of the examples of $\delta$-sequences in the literature can be verified in this way: consider the terms’ primitives and show that they converge to the Heaviside function. The result then follows as above.
The first example (Sokhotsky) in the question can be proved in one line, after integrating $\dfrac{\epsilon}{x^2+\epsilon^2}.$