dsolve(diff(y(x),x) = -1/128*I*(32*I*x+64+64*y(x)^4+32*x^2*y(x)^2+4*x^4+64*y(x)^6+48*x^2*y(x)^4+12*x^4*y(x)^2+x^6)/y(x),y(x))