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