Solution (source code)

= Solution

Put $Q=1+x^2+y^2>0$. Differentiation of $Q^{-1}$ gives
$$
\boxed{\nabla f(x,y)=\frac{-2}{(1+x^2+y^2)^2}\begin{pmatrix}x\\y\end{pmatrix}.}
$$
The positive denominator also shows that the function and these derivatives are smooth throughout $\mathbb R^2$.