\[ \boxed { {\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) =-{\frac { \left ( {x}^{2} \left ( \left ( {x}^{2}-{\it a1} \right ) \left ( {x}^{2}-{\it a2} \right ) + \left ( {x}^{2}-{\it a2} \right ) \left ( {x}^{2}-{\it a3} \right ) + \left ( {x}^{2}-{\it a3} \right ) \left ( {x}^{2}-{\it a1} \right ) \right ) - \left ( {x}^{2}-{\it a1} \right ) \left ( {x}^{2}-{\it a2} \right ) \left ( {x}^{2}-{\it a3} \right ) \right ) {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) }{x \left ( {x}^{2}-{\it a1} \right ) \left ( {x}^{2}-{\it a2} \right ) \left ( {x}^{2}-{\it a3} \right ) }}-{\frac { \left ( A{x}^{2}+B \right ) y \left ( x \right ) }{x \left ( {x}^{2}-{\it a1} \right ) \left ( {x}^{2}-{\it a2} \right ) \left ( {x}^{2}-{\it a3} \right ) }}=0} \]
Mathematica: cpu = 59.615570 (sec), leaf count = 98 \[ \left \{\left \{y(x)\to \text {DifferentialRoot}\left (\{\unicode {f818},\unicode {f817}\}\unicode {f4a1}\left \{\left (A \unicode {f817}^2+B\right ) \unicode {f818}(\unicode {f817})+\left (2 \unicode {f817}^6-\text {a1} \unicode {f817}^4-\text {a2} \unicode {f817}^4-\text {a3} \unicode {f817}^4+\text {a1} \text {a2} \text {a3}\right ) \unicode {f818}'(\unicode {f817})-\unicode {f817} \left (\text {a1}-\unicode {f817}^2\right ) \left (\text {a2}-\unicode {f817}^2\right ) \left (\text {a3}-\unicode {f817}^2\right ) \unicode {f818}''(\unicode {f817})=0,\unicode {f818}(1)=c_1,\unicode {f818}'(1)=c_2\right \}\right )(x)\right \}\right \} \]
Maple: cpu = 7.301 (sec), leaf count = 165 \[ \left \{ y \left ( x \right ) ={\it DESol} \left ( \left \{ {\frac { {\rm d}^{2}}{{\rm d}{x}^{2}}}{\it \_Y} \left ( x \right ) +{\frac { \left ( {x}^{2} \left ( \left ( {x}^{2}-{\it a1} \right ) \left ( {x}^{2 }-{\it a2} \right ) + \left ( {x}^{2}-{\it a2} \right ) \left ( {x}^{2}-{ \it a3} \right ) + \left ( {x}^{2}-{\it a3} \right ) \left ( {x}^{2}-{ \it a1} \right ) \right ) - \left ( {x}^{2}-{\it a1} \right ) \left ( {x} ^{2}-{\it a2} \right ) \left ( {x}^{2}-{\it a3} \right ) \right ) { \frac {\rm d}{{\rm d}x}}{\it \_Y} \left ( x \right ) }{ \left ( {x}^{2}-{ \it a1} \right ) \left ( {x}^{2}-{\it a2} \right ) x \left ( {x}^{2}-{ \it a3} \right ) }}+{\frac { \left ( A{x}^{2}+B \right ) {\it \_Y} \left ( x \right ) }{ \left ( {x}^{2}-{\it a1} \right ) \left ( {x}^{2}-{ \it a2} \right ) x \left ( {x}^{2}-{\it a3} \right ) }} \right \} , \left \{ {\it \_Y} \left ( x \right ) \right \} \right ) \right \} \]