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