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