y'(x)=\frac{x^3 \sqrt{4 a x-y(x)^2}+x^2 \sqrt{4 a x-y(x)^2}+\sqrt{4 a x-y(x)^2}+2 a}{y(x)}