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