Write the equation as
\[\frac{dy}{dx}=\frac{x^{3}-3xy^{2}}{y^{3}-3x^{2}y}.\]
Numerator and denominator are both homogeneous of degree \(\displaystyle 3\), so the equation is homogeneous of degree zero and the substitution \(\displaystyle y=vx\), \(\displaystyle \frac{dy}{dx}=v+x\frac{dv}{dx}\) applies.
Substituting,
\[v+x\frac{dv}{dx}=\frac{x^{3}-3x\cdot v^{2}x^{2}}{v^{3}x^{3}-3x^{2}\cdot vx}=\frac{x^{3}(1-3v^{2})}{x^{3}(v^{3}-3v)}=\frac{1-3v^{2}}{v^{3}-3v}\qquad(x\neq 0).\]
Hence
\[x\frac{dv}{dx}=\frac{1-3v^{2}-v(v^{3}-3v)}{v^{3}-3v}=\frac{1-3v^{2}-v^{4}+3v^{2}}{v^{3}-3v}=\frac{1-v^{4}}{v^{3}-3v},\]
and the variables separate:
\[\frac{v^{3}-3v}{1-v^{4}}\,dv=\frac{dx}{x}.\]
The left integral is the step to be careful with. The numerator is odd in \(\displaystyle v\), so put \(\displaystyle t=v^{2}\), \(\displaystyle dt=2v\,dv\), and use \(\displaystyle 1-v^{4}=(1-t)(1+t)\):
\[\int\frac{v(v^{2}-3)}{(1-v^{2})(1+v^{2})}dv=\frac{1}{2}\int\frac{t-3}{(1-t)(1+t)}\,dt.\]
Partial fractions: \(\displaystyle \dfrac{t-3}{(1-t)(1+t)}=\dfrac{A}{1-t}+\dfrac{B}{1+t}\) gives \(\displaystyle t-3=A(1+t)+B(1-t)\); putting \(\displaystyle t=1\) gives \(\displaystyle A=-1\) and \(\displaystyle t=-1\) gives \(\displaystyle B=-2\). So
\[\frac{1}{2}\int\left(\frac{-1}{1-t}-\frac{2}{1+t}\right)dt=\frac{1}{2}\log|1-t|-\log|1+t|=\frac{1}{2}\log|1-v^{2}|-\log\left(1+v^{2}\right).\]
Therefore
\[\frac{1}{2}\log|1-v^{2}|-\log\left(1+v^{2}\right)=\log|x|+\log k.\]
Multiply by \(\displaystyle 2\) and combine:
\[\log\frac{|1-v^{2}|}{(1+v^{2})^{2}}=\log\left(Kx^{2}\right)\quad\Longrightarrow\quad \frac{1-v^{2}}{(1+v^{2})^{2}}=Kx^{2},\]
the sign being absorbed into the arbitrary constant \(\displaystyle K\).
Put back \(\displaystyle v=\dfrac{y}{x}\): \(\displaystyle 1-v^{2}=\dfrac{x^{2}-y^{2}}{x^{2}}\) and \(\displaystyle (1+v^{2})^{2}=\dfrac{(x^{2}+y^{2})^{2}}{x^{4}}\), so
\[\frac{x^{2}(x^{2}-y^{2})}{(x^{2}+y^{2})^{2}}=Kx^{2}\quad\Longrightarrow\quad x^{2}-y^{2}=K\left(x^{2}+y^{2}\right)^{2}.\]
Writing \(\displaystyle K=c\), the relation \(\displaystyle x^{2}-y^{2}=c\left(x^{2}+y^{2}\right)^{2}\) satisfies the equation and contains exactly one arbitrary constant, which is the number required by a first-order equation. Hence it is the general solution.