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