ODE
\[ -\left (4-x^2\right ) y'(x)+x y''(x)+2 x y(x)=0 \] ODE Classification
[[_2nd_order, _exact, _linear, _homogeneous]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.208334 (sec), leaf count = 68
\[\left \{\left \{y(x)\to \frac {1}{30} \left (\sqrt {2 \pi } c_2 e^{-\frac {x^2}{2}} x^5 \text {erfi}\left (\frac {x}{\sqrt {2}}\right )+30 c_1 e^{-\frac {x^2}{2}} x^5-2 c_2 \left (x^4+x^2+3\right )\right )\right \}\right \}\]
Maple ✓
cpu = 0.22 (sec), leaf count = 55
\[\left [y \left (x \right ) = \textit {\_C1} \,{\mathrm e}^{-\frac {x^{2}}{2}} x^{5}+\textit {\_C2} \left (-i \sqrt {\pi }\, \sqrt {2}\, \left (x^{4}+x^{2}+3\right )+\erf \left (\frac {i \sqrt {2}\, x}{2}\right ) \pi \,{\mathrm e}^{-\frac {x^{2}}{2}} x^{5}\right )\right ]\] Mathematica raw input
DSolve[2*x*y[x] - (4 - x^2)*y'[x] + x*y''[x] == 0,y[x],x]
Mathematica raw output
{{y[x] -> ((30*x^5*C[1])/E^(x^2/2) - 2*(3 + x^2 + x^4)*C[2] + (Sqrt[2*Pi]*x^5*C[2]*Erfi[x/Sqrt[2]])/E^(x^2/2))/30}}
Maple raw input
dsolve(x*diff(diff(y(x),x),x)-(-x^2+4)*diff(y(x),x)+2*x*y(x) = 0, y(x))
Maple raw output
[y(x) = _C1*exp(-1/2*x^2)*x^5+_C2*(-I*Pi^(1/2)*2^(1/2)*(x^4+x^2+3)+erf(1/2*I*2^(1/2)*x)*Pi*exp(-1/2*x^2)*x^5)]