dsolve(diff(y(x),x) = 2*a*x/(-x^3*y(x)+2*a*x^3+2*a*y(x)^4*x^3-16*y(x)^2*a^2*x^2+32*a^3*x+2*a*y(x)^6*x^3-24*y(x)^4*a^2*x^2+96*y(x)^2*x*a^3-128*a^4),y(x))