Dado que alguien ya ha demostrado que todas las derivadas direccionales existen, solo argumentaré por qué $f$ no es diferenciable en $0$.
La matriz de Jacobi $A:=Df(0,0)=(1,0)$. Por lo tanto, si $f$ es diferenciable $$\lim_{|\epsilon| \to 0}\frac{f(0+\epsilon)-f(0)-A\epsilon}{|\epsilon|}=0 .$$ Dado que $f(0)=0$ y $A=(1,0)$ esto es equivalente a, $$\lim_{|\epsilon| \to 0}\frac{f(\epsilon)-(1,0)\epsilon}{|\epsilon|}=0 $$ Sea $(x_k)_{k \in \mathbb{N}} \subset \mathbb{R}^2$ una serie con $x_k=(a_k,b_k)$ y $|x_k| \to 0$ para $ (k \to \infty)$. Por lo tanto, si $f$ es diferenciable.
$$0=\lim_{k \to \infty}\frac{f(x_k)-(1,0)x_k}{|x_k|}=\lim_{k \to \infty}\frac{\frac{a_k^3}{a_k^2+b_k^2}-a_k}{\sqrt{a_k^2+b_k^2}}= \lim_{k \to \infty}\frac{-a_k b_k^2}{\sqrt{a_k^2+b_k^2}^3}$$
Si definimos $x_k=(a_k,b_k)= (\frac {1}{k\sqrt{3}},\frac{\sqrt{2}}{k\sqrt{3}})$ entonces $|x_k|=\sqrt{a_k^2+b_k^2}=1/k$ y
$$\lim_{k \to \infty}\frac{-a_k b_k^2}{\sqrt{a_k^2+b_k^2}^3}=\lim_{k \to \infty}\frac{-\frac {1}{k\sqrt{3}} \frac{2}{3k^2}}{(1/k)^3}=\lim_{k \to \infty}-k^3\frac{2}{k^3\sqrt{3}^3}=\frac{-2}{\sqrt{3}^3}\neq 0$$ Obtenemos una contradicción, por lo tanto $f$ no es diferenciable en 0.