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