dsolve(diff(y(x),x) = 2*a/x/(-x*y(x)+2*a*x*y(x)^2-8*a^2),y(x))