With the change of variables, we have $\displaystyle\lim_{(x,y) \to (0, a)} \frac{\sin(xy)}x = \displaystyle\lim_{(x,y) \to (0, 0)} \frac{\sin(x(y+a))}x$
Now we perform the exchange of variables $\begin{cases} x = r \cos \theta \\ y = r \sin \theta\end{cases}$ and get
$\displaystyle\lim_{(x,y) \to (0, a)} \frac{\sin(xy)}x = \displaystyle\lim_{r \to 0} \frac{\sin(r^2 \cos \theta \sin \theta + ar \cos\theta)}{r \cos \theta}$.
The Taylor expansion for $\sin(r^2 \cos \theta \sin \theta + ar \cos\theta)$ in $r$ at $r = 0$ is $\sin(r^2 \cos \theta \sin \theta + ar \cos\theta) = ar \cos \theta + \mathcal{O}(r^2).$
If one now tries to take the limit, one will have
$\displaystyle\lim_{r \to 0} \frac{\sin(r^2 \cos \theta \sin \theta + ar \cos\theta)}{r \cos \theta} = \displaystyle\lim_{r \to 0} \frac{ar \cos \theta + \mathcal{O}(r^2)}{r \cos \theta} = \displaystyle\lim_{r \to 0} a + \mathcal{O}(r) = a.$
According to this reasoning, the limit should exist for all $a$, and be equal to $a.$ Is this, however, a valid proof for the limit? I tried an $(\varepsilon, \delta_\varepsilon)$-approach as well, but I didn't get anything "good" to work with.