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