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