ODE
\[ y(x) \left (\text {a0}+4 \text {a1} \sinh ^2(x)-\text {a2} \text {csch}^2(x)\right )+y''(x)+\coth (x) y'(x)=0 \] ODE Classification
[[_2nd_order, _with_linear_symmetries]]
Book solution method
TO DO
Mathematica ✗
cpu = 9.17248 (sec), leaf count = 0 , could not solve
DSolve[(a0 - a2*Csch[x]^2 + 4*a1*Sinh[x]^2)*y[x] + Coth[x]*Derivative[1][y][x] + Derivative[2][y][x] == 0, y[x], x]
Maple ✓
cpu = 2.958 (sec), leaf count = 71
\[\left [y \left (x \right ) = \textit {\_C1} \left (\sinh ^{\sqrt {\mathit {a2}}}\left (x \right )\right ) \HeunC \left (0, \sqrt {\mathit {a2}}, -\frac {1}{2}, -\mathit {a1} , \frac {1}{4}+\frac {\mathit {a0}}{4}+\frac {\mathit {a2}}{4}, -\left (\sinh ^{2}\left (x \right )\right )\right )+\textit {\_C2} \left (\sinh ^{-\sqrt {\mathit {a2}}}\left (x \right )\right ) \HeunC \left (0, -\sqrt {\mathit {a2}}, -\frac {1}{2}, -\mathit {a1} , \frac {1}{4}+\frac {\mathit {a0}}{4}+\frac {\mathit {a2}}{4}, -\left (\sinh ^{2}\left (x \right )\right )\right )\right ]\] Mathematica raw input
DSolve[(a0 - a2*Csch[x]^2 + 4*a1*Sinh[x]^2)*y[x] + Coth[x]*y'[x] + y''[x] == 0,y[x],x]
Mathematica raw output
DSolve[(a0 - a2*Csch[x]^2 + 4*a1*Sinh[x]^2)*y[x] + Coth[x]*Derivative[1][y][x] + Derivative[2][y][x] == 0, y[x], x]
Maple raw input
dsolve(diff(diff(y(x),x),x)+diff(y(x),x)*coth(x)+(a0+4*a1*sinh(x)^2-a2*csch(x)^2)*y(x) = 0, y(x))
Maple raw output
[y(x) = _C1*sinh(x)^(a2^(1/2))*HeunC(0,a2^(1/2),-1/2,-a1,1/4+1/4*a0+1/4*a2,-sinh(x)^2)+_C2*sinh(x)^(-a2^(1/2))*HeunC(0,-a2^(1/2),-1/2,-a1,1/4+1/4*a0+1/4*a2,-sinh(x)^2)]