dsolve(diff(y(x),x) = (1+y(x)^4-8*a*x*y(x)^2+16*a^2*x^2+y(x)^6-12*y(x)^4*a*x+48*y(x)^2*a^2*x^2-64*a^3*x^3)/y(x),y(x))