A possibility is to represent the point $(x,y)$ by a complex number:
$$z = x + j y = \rho e^{j \theta}$$
then
$$ xy = \rho^2 \frac{\sin 2\theta}{2} $$
$$ x^2 - y^2 = \rho^2 (2\cos^2 \theta - 1) $$
$$ x^2 + y^2 = \rho^2$$
And finally:
$$f(x, y) = \frac{xy\,(x^2-y^2)}{x^2+y^2} = f(z) = \rho^2 \frac{\sin2\theta \, (2\cos^2 \theta - 1)}{2} $$
When the point $(x, y)$ tends to $0$, $\rho$ necessarily tends to $0$ and therefore, whatever the value of $\theta$, $f(z) = f(x, y)$ tends to $0$
This method may look overskill for this simple function. The sandwich method proposed by OP works. However, OP was asking for a possible general method. This one often gives good results even in more difficult use cases.