\[ \left (x^2-1\right ) y'(x)+x^2 y''(x)-y(x)=0 \] ✓ Mathematica : cpu = 0.0781562 (sec), leaf count = 40
\[\left \{\left \{y(x)\to c_2 e^{-x} \int _1^xe^{K[1]-\frac {1}{K[1]}}dK[1]+c_1 e^{-x}\right \}\right \}\] ✓ Maple : cpu = 0.9 (sec), leaf count = 51
\[ \left \{ y \left ( x \right ) =\sqrt {x} \left ( {{\rm e}^{-{x}^{-1}}}{\it HeunD} \left ( -4,3,-8,5,{\frac {x-1}{1+x}} \right ) {\it \_C2}+{{\rm e}^{-x}}{\it HeunD} \left ( 4,3,-8,5,{\frac {x-1}{1+x}} \right ) {\it \_C1} \right ) \right \} \]