dsolve(diff(y(x),x) = (2*x^3*y(x)+x^6+x^2*y(x)^2+y(x)^3)/x^4,y(x))