\[ y'(x)=\frac {a^3 x^3 y(x)^3+3 a^2 x^2 y(x)^2+a^2 x+3 a x y(x)+1}{a^3 x^3} \] ✓ Mathematica : cpu = 0.146897 (sec), leaf count = 49
\[\left \{\left \{y(x)\to -\frac {1}{a x}-\frac {1}{\sqrt {-2 x+c_1}}\right \},\left \{y(x)\to -\frac {1}{a x}+\frac {1}{\sqrt {-2 x+c_1}}\right \}\right \}\] ✓ Maple : cpu = 0.046 (sec), leaf count = 41
\[ \left \{ y \left ( x \right ) =-{\frac {1}{\sqrt {{\it \_C1}-2\,x}}}-{\frac {1}{ax}},y \left ( x \right ) ={\frac {1}{\sqrt {{\it \_C1}-2\,x}}}-{\frac {1}{ax}} \right \} \]