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))