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