4.33.16 \(y(x) (a+b x)+2 (1-x) x y''(x)+(1-2 x) y'(x)=0\)

ODE
\[ y(x) (a+b x)+2 (1-x) x y''(x)+(1-2 x) y'(x)=0 \] ODE Classification

[_Jacobi]

Book solution method
TO DO

Mathematica
cpu = 2.29214 (sec), leaf count = 0 , DifferentialRoot result

\[\left \{\left \{y(x)\to \text {DifferentialRoot}\left (\{\unicode {f818},\unicode {f817}\}\unicode {f4a1}\left \{(-a-\unicode {f817} b) \unicode {f818}(\unicode {f817})+(2 \unicode {f817}-1) \unicode {f818}'(\unicode {f817})+2 (\unicode {f817}-1) \unicode {f817} \unicode {f818}''(\unicode {f817})=0,\unicode {f818}(2)=c_1,\unicode {f818}'(2)=c_2\right \}\right )(x)\right \}\right \}\]

Maple
cpu = 0.146 (sec), leaf count = 35

\[ \left \{ y \left ( x \right ) ={\it \_C1}\,{\it MathieuC} \left ( b+2\,a,-{\frac {b}{2}},\arccos \left ( \sqrt {x} \right ) \right ) +{\it \_C2}\,{\it MathieuS} \left ( b+2\,a,-{\frac {b}{2}},\arccos \left ( \sqrt {x} \right ) \right ) \right \} \] Mathematica raw input

DSolve[(a + b*x)*y[x] + (1 - 2*x)*y'[x] + 2*(1 - x)*x*y''[x] == 0,y[x],x]

Mathematica raw output

{{y[x] -> DifferentialRoot[Function[{\[FormalY], \[FormalX]}, {(-a - \[FormalX]*
b)*\[FormalY][\[FormalX]] + (-1 + 2*\[FormalX])*Derivative[1][\[FormalY]][\[Form
alX]] + 2*(-1 + \[FormalX])*\[FormalX]*Derivative[2][\[FormalY]][\[FormalX]] == 
0, \[FormalY][2] == C[1], Derivative[1][\[FormalY]][2] == C[2]}]][x]}}

Maple raw input

dsolve(2*x*(1-x)*diff(diff(y(x),x),x)+(1-2*x)*diff(y(x),x)+(b*x+a)*y(x) = 0, y(x),'implicit')

Maple raw output

y(x) = _C1*MathieuC(b+2*a,-1/2*b,arccos(x^(1/2)))+_C2*MathieuS(b+2*a,-1/2*b,arcc
os(x^(1/2)))