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