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