As usual in such examples, there is no need to integrate against a test function. It is a simple consequence of One can simply use the fact that if a sequence (or net) of distributions converges in the distributional sense, then so does the sequenceone obtained by differentation also convergesdifferentiating term by term. In particular, this applies ifwhen the original sequence consists of functions which converge in pretty well any sensible classical sense, e.g., locally $L^1$ as isin the case here.
In this example we considerin 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 in the above sense and so and we can differentiate to obtain the required result. Most
Most of the representationsexamples of $\delta$-sequences in the literature can be obtainedverified in this way: consider the terms’ primitives and show that they form a sequence which convergesconverge to the Heaviside function in a suitable sense. The result then follows as above.
The first example (SokhotzkySokhotsky) in the question can be proved in the same way in one line, if you can integrateafter integrating $\dfrac{\epsilon}{x^2+\epsilon^2}.$