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