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