\[ y'(x)=\frac {a^3 x^6}{64}+\frac {3}{32} a^2 b x^5+\frac {3}{16} a^2 x^4 y(x)+\frac {a^2 x^4}{16}+\frac {3}{16} a b^2 x^4+\frac {3}{4} a b x^3 y(x)+\frac {1}{4} a b x^3+\frac {3}{4} a x^2 y(x)^2+\frac {1}{2} a x^2 y(x)-\frac {a x}{2}+\frac {b^3 x^3}{8}+\frac {3}{4} b^2 x^2 y(x)+\frac {b^2 x^2}{4}+\frac {3}{2} b x y(x)^2+b x y(x)+y(x)^3+y(x)^2+1 \] ✓ Mathematica : cpu = 0.282194 (sec), leaf count = 1
\[\{\}\]
✓ Maple : cpu = 0.104 (sec), leaf count = 42
\[ \left \{ y \left ( x \right ) =-{\frac {a{x}^{2}}{4}}-{\frac {bx}{2}}+{\it RootOf} \left ( -x+2\,\int ^{{\it \_Z}}\! \left ( 2\,{{\it \_a}}^{3}+2\,{{\it \_a}}^{2}+b+2 \right ) ^{-1}{d{\it \_a}}+{\it \_C1} \right ) \right \} \]