dsolve(diff(y(x),x) = -I*(I*x+x^4+2*x^2*y(x)^2+y(x)^4)/y(x),y(x))