\[ a x y(x)^2+\left (x^2-1\right ) y'(x)+x y(x)=0 \] ✓ Mathematica : cpu = 0.0366478 (sec), leaf count = 31
\[\left \{\left \{y(x)\to -\frac {e^{c_1}}{a e^{c_1}-\sqrt {x^2-1}}\right \}\right \}\]
✓ Maple : cpu = 0.022 (sec), leaf count = 22
\[ \left \{ y \left ( x \right ) = \left ( \sqrt {x-1}\sqrt {1+x}{\it \_C1}-a \right ) ^{-1} \right \} \]