I changed to polar coordinates $x=r \cos \varphi$ and $y =r \sin \varphi$ for
$$f(r,\varphi) = {\rm e}^{(r^2-2 r^2 \cos^2 \varphi)} \left( \sin \varphi \cos \varphi \sin(2015 r \cos \varphi + 2016 r \sin \varphi) + r^2 \sin^3 \varphi \cos \varphi \right) $$
This clearly shows that $f(0,\varphi)=0$ with any divisions by zero.
Now for the continuity question. We have to look at $g(r,\varphi) = \frac{\partial}{\partial r} f(r,\varphi)$.
You can do the math (or use Wolfram) and show that
$$ g(0,\varphi) = 2015 \sin\varphi \cos^2 \varphi + 2016 \sin^2 \varphi \cos\varphi $$
The above is continuous for all $\varphi$ and since $f(0,\varphi)=0$ and $\frac{\partial}{\partial r}f(r,\varphi) |_{r=0} = 0$ it is continuous near $r=0$ also.
The 2nd order taylor expansion of $f(r,\varphi)$ is
$$f(r,\varphi) \approx r \left( (2015 \sin\varphi \cos^2 \varphi + 2016 \sin^2 \varphi \cos \varphi) + r (\sin^3 \varphi \cos \varphi) \right) $$