\[ y'(x)=\frac {-64 a^3 x^3+48 a^2 x^2 y(x)^2+16 a^2 x^2-12 a x y(x)^4-8 a x y(x)^2+y(x)^6+y(x)^4+1}{y(x)} \] ✓ Mathematica : cpu = 0.315198 (sec), leaf count = 130
\[\text {Solve}\left [2 a \left (x-\frac {1}{2} \text {RootSum}\left [64 \text {$\#$1}^3 a^3-48 \text {$\#$1}^2 a^2 y(x)^2-16 \text {$\#$1}^2 a^2+12 \text {$\#$1} a y(x)^4+8 \text {$\#$1} a y(x)^2+2 a-y(x)^6-y(x)^4-1\& ,\frac {\log (x-\text {$\#$1})}{48 \text {$\#$1}^2 a^2-24 \text {$\#$1} a y(x)^2-8 \text {$\#$1} a+3 y(x)^4+2 y(x)^2}\& \right ]\right )=c_1,y(x)\right ]\] ✗ Maple : cpu = 0. (sec), leaf count = 0 , could not solve
dsolve(diff(y(x),x) = (1+y(x)^4-8*a*x*y(x)^2+16*a^2*x^2+y(x)^6-12*y(x)^4*a*x+48*y(x)^2*a^2*x^2-64*a^3*x^3)/y(x),y(x))