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