\[ y'(x) \left (\sin (\alpha ) \left (y(x)^2-x^2\right )-2 x \cos (\alpha ) y(x)+\sqrt {x^2+y(x)^2} y(x)\right )+\cos (\alpha ) \left (y(x)^2-x^2\right )+2 x \sin (\alpha ) y(x)+x \sqrt {x^2+y(x)^2}=0 \] ✓ Mathematica : cpu = 61.416 (sec), leaf count = 17681
DSolve[2*x*Sin[alpha]*y[x] + Cos[alpha]*(-x^2 + y[x]^2) + x*Sqrt[x^2 + y[x]^2] + (-2*x*Cos[alpha]*y[x] + Sin[alpha]*(-x^2 + y[x]^2) + y[x]*Sqrt[x^2 + y[x]^2])*Derivative[1][y][x] == 0,y[x],x]
\[ \text {Too large to display} \] ✓ Maple : cpu = 0.857 (sec), leaf count = 127
dsolve((y(x)*(y(x)^2+x^2)^(1/2)+(y(x)^2-x^2)*sin(alpha)-2*x*y(x)*cos(alpha))*diff(y(x),x)+x*(y(x)^2+x^2)^(1/2)+2*x*y(x)*sin(alpha)+(y(x)^2-x^2)*cos(alpha) = 0,y(x))
\[y \left (x \right ) = \RootOf \left (-\ln \left (x \right )+\int _{}^{\textit {\_Z}}-\frac {\textit {\_a}^{3} \cos \left (2 \alpha \right )+3 \textit {\_a}^{2} \sin \left (2 \alpha \right )+\textit {\_a}^{3}-3 \textit {\_a} \cos \left (2 \alpha \right )-\sin \left (2 \alpha \right )+\sqrt {2}\, \sqrt {\left (\textit {\_a}^{2}+1\right ) \left (\textit {\_a}^{2} \cos \left (2 \alpha \right )+2 \textit {\_a} \sin \left (2 \alpha \right )+\textit {\_a}^{2}-\cos \left (2 \alpha \right )+1\right )}+\textit {\_a}}{\left (\textit {\_a}^{2}+1\right ) \left (\textit {\_a}^{2} \cos \left (2 \alpha \right )+2 \textit {\_a} \sin \left (2 \alpha \right )+\textit {\_a}^{2}-\cos \left (2 \alpha \right )+1\right )}d \textit {\_a} +c_{1}\right ) x\]