4.29.19 \((1-x \cot (x)) y''(x)-x y'(x)+y(x)=0\)

ODE
\[ (1-x \cot (x)) y''(x)-x y'(x)+y(x)=0 \] ODE Classification

[[_2nd_order, _with_linear_symmetries]]

Book solution method
TO DO

Mathematica
cpu = 0.27977 (sec), leaf count = 15

\[\{\{y(x)\to c_1 x+c_2 \sin (x)\}\}\]

Maple
cpu = 3.037 (sec), leaf count = 81

\[\left [y \left (x \right ) = x \textit {\_C2} +\textit {\_C1} x \left (-\frac {i {\mathrm e}^{2 i x}}{2 \sqrt {\frac {{\mathrm e}^{2 i x}}{\left ({\mathrm e}^{2 i x}+1\right )^{2}}}\, \left ({\mathrm e}^{2 i x}+1\right ) x}+\frac {i}{2 \sqrt {\frac {{\mathrm e}^{2 i x}}{\left ({\mathrm e}^{2 i x}+1\right )^{2}}}\, \left ({\mathrm e}^{2 i x}+1\right ) x}\right )\right ]\] Mathematica raw input

DSolve[y[x] - x*y'[x] + (1 - x*Cot[x])*y''[x] == 0,y[x],x]

Mathematica raw output

{{y[x] -> x*C[1] + C[2]*Sin[x]}}

Maple raw input

dsolve((1-x*cot(x))*diff(diff(y(x),x),x)-x*diff(y(x),x)+y(x) = 0, y(x))

Maple raw output

[y(x) = x*_C2+_C1*x*(-1/2*I/(exp(2*I*x)/(exp(2*I*x)+1)^2)^(1/2)/(exp(2*I*x)+1)*e
xp(2*I*x)/x+1/2*I/(exp(2*I*x)/(exp(2*I*x)+1)^2)^(1/2)/(exp(2*I*x)+1)/x)]