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