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