 Solution

Choose a standard mollifier $\rho \in C_0^{\infty }({\mathbb{R}}^n)$ with $\rho \ge 0$, $\int \rho =1$ and support in the unit ball, and set ${\rho }_{\epsilon }(z)={\epsilon }^{-n}\rho (z/\epsilon )$. Define the convolutionThe smoothing properties of a mollifier give $u_{\epsilon }\in C^{\infty }$ and ${\partial }_j{u}_{\epsilon }={\rho }_{\epsilon }\ast D_{j}u$. Apply the chain rule to $t\mapsto u_{\epsilon }(y+t(x-y))$ and integrate its derivative by the fundamental theorem of calculus: