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