\[ y''(x)=-\frac {b x y'(x)}{a \left (x^2-1\right )}-\frac {y(x) \left (c x^2+d x+e\right )}{a \left (x^2-1\right )^2} \] ✓ Mathematica : cpu = 50.3879 (sec), leaf count = 1763961
\[ \text {Too large to display} \] ✓ Maple : cpu = 0.181 (sec), leaf count = 561
\[y \relax (x ) = \left (-\frac {1}{2}+\frac {x}{2}\right )^{\frac {2 a +\sqrt {4 a^{2}+\left (-4 b -4 c -4 d -4 e \right ) a +b^{2}}}{4 a}} \left (x^{2}-1\right )^{-\frac {b}{4 a}} \left (\hypergeom \left (\left [-\frac {-\sqrt {4 a^{2}+\left (-4 b -4 c -4 d -4 e \right ) a +b^{2}}+2 \sqrt {a^{2}+\left (-2 b -4 c \right ) a +b^{2}}+\sqrt {4 a^{2}+\left (-4 b -4 c +4 d -4 e \right ) a +b^{2}}-2 a}{4 a}, -\frac {-\sqrt {4 a^{2}+\left (-4 b -4 c -4 d -4 e \right ) a +b^{2}}-2 \sqrt {a^{2}+\left (-2 b -4 c \right ) a +b^{2}}+\sqrt {4 a^{2}+\left (-4 b -4 c +4 d -4 e \right ) a +b^{2}}-2 a}{4 a}\right ], \left [-\frac {-2 a +\sqrt {4 a^{2}+\left (-4 b -4 c +4 d -4 e \right ) a +b^{2}}}{2 a}\right ], \frac {1}{2}+\frac {x}{2}\right ) \left (\frac {1}{2}+\frac {x}{2}\right )^{\frac {2 a -\sqrt {4 a^{2}+\left (-4 b -4 c +4 d -4 e \right ) a +b^{2}}}{4 a}} c_{1}+\hypergeom \left (\left [\frac {\sqrt {4 a^{2}+\left (-4 b -4 c -4 d -4 e \right ) a +b^{2}}-2 \sqrt {a^{2}+\left (-2 b -4 c \right ) a +b^{2}}+\sqrt {4 a^{2}+\left (-4 b -4 c +4 d -4 e \right ) a +b^{2}}+2 a}{4 a}, \frac {\sqrt {4 a^{2}+\left (-4 b -4 c -4 d -4 e \right ) a +b^{2}}+2 \sqrt {a^{2}+\left (-2 b -4 c \right ) a +b^{2}}+\sqrt {4 a^{2}+\left (-4 b -4 c +4 d -4 e \right ) a +b^{2}}+2 a}{4 a}\right ], \left [\frac {2 a +\sqrt {4 a^{2}+\left (-4 b -4 c +4 d -4 e \right ) a +b^{2}}}{2 a}\right ], \frac {1}{2}+\frac {x}{2}\right ) \left (\frac {1}{2}+\frac {x}{2}\right )^{\frac {2 a +\sqrt {4 a^{2}+\left (-4 b -4 c +4 d -4 e \right ) a +b^{2}}}{4 a}} c_{2}\right )\]