An electrodynamics lecture asks me to prove that $$ \nabla \times \left( \frac{\vec{M} \times \vec{x}}{ |\vec{x}|^3} \right) = \frac{8 \pi}{3} \vec{M} \delta^3(\vec{x})- \frac{\vec{M}}{|\vec{x}|^3}+ \frac{3\vec{M} \cdot \vec{x}}{|\vec{x}|^5}\vec{x} $$ where $\vec{M}$ is some constant vector. I have already shown that for $\vec{x} \neq 0$ $$ \nabla \times \left( \frac{\vec{M} \times \vec{x}}{ |\vec{x}|^3} \right) = - \frac{\vec{M}}{|\vec{x}|^3}+ \frac{3\vec{M} \cdot \vec{x}}{|\vec{x}|^5}\vec{x} $$ To determine the value at the singularity I was thinking to consider the integral over any volume V
$$ \int \int \int_V \nabla \times \left( \frac{\vec{M} \times \vec{x}}{ |\vec{x}|^3} \right) d^3x =- \int \int _{dV} \left( \frac{\vec{M} \times \vec{x}}{ |\vec{x}|^3} \right) \times \vec{da} \ \ \ \ \ \ \ \ (*) $$ where the equality is due to a special case of the Gaussian theorem. I think I have to show $$ \int \int \int_V \nabla \times \left( \frac{\vec{M} \times \vec{x}}{ |\vec{x}|^3} \right) d^3x = \frac{8\pi}{3} \vec{M} + \int \int \int_V \left( {- \frac{\vec{M}}{|\vec{x}|^3}+ \frac{3\vec{M} \cdot \vec{x}}{|\vec{x}|^5}\vec{x} } \right )d^3x $$ for all volumes that include the origin.
Is this approach correct? I am stuck at evaluating the integral (*). Can someone give me some pointers? I did not know what this method of integrating singularities by using the delta function is called and I couldn't find anything helpful in my literature or on the internet. Thanks!