\[ y'(x)=x^2 \sqrt {a^2 x^2+2 a b x+4 a y(x)+b^2-4 c}-\frac {a x}{2}-\frac {b}{2} \] ✓ Mathematica : cpu = 0.632958 (sec), leaf count = 164
\[\text {Solve}\left [-\frac {\sqrt {a^2 x^2+2 a b x+4 a y(x)+b^2-4 c}}{2 a}-\frac {b \log \left (\sqrt {a^2 x^2+2 a b x+4 a y(x)+b^2-4 c}+a x+b\right )}{2 a}+\frac {b \tanh ^{-1}\left (\frac {2 a^2 x+2 a b}{2 a \sqrt {a^2 x^2+2 a b x+4 a y(x)+b^2-4 c}}\right )}{2 a}+\frac {b \log (4 a y(x)-4 c)}{4 a}+\frac {x^3}{3}=c_1,y(x)\right ]\] ✓ Maple : cpu = 0.492 (sec), leaf count = 39
\[ \left \{ {\it \_C1}+{\frac {2\,a{x}^{3}}{3}}-\sqrt {{a}^{2}{x}^{2}+2\,abx+{b}^{2}+4\,ay \left ( x \right ) -4\,c}=0 \right \} \]