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