\[ x^2 \left (-y'(x)\right )+y''(x)+x y(x)=0 \] ✓ Mathematica : cpu = 0.0397665 (sec), leaf count = 41
DSolve[x*y[x] - x^2*Derivative[1][y][x] + Derivative[2][y][x] == 0,y[x],x]
\[\left \{\left \{y(x)\to c_1 x-\frac {c_2 \sqrt [3]{-x^3} \Gamma \left (-\frac {1}{3},-\frac {x^3}{3}\right )}{3 \sqrt [3]{3}}\right \}\right \}\] ✓ Maple : cpu = 0.097 (sec), leaf count = 48
dsolve(diff(diff(y(x),x),x)-x^2*diff(y(x),x)+x*y(x)=0,y(x))
\[y \left (x \right ) = \frac {-\left (-x^{3}\right )^{\frac {2}{3}} 3^{\frac {1}{3}} {\mathrm e}^{\frac {x^{3}}{3}} c_{2}+x^{3} \left (c_{2} \Gamma \left (\frac {2}{3}\right )-c_{2} \Gamma \left (\frac {2}{3}, -\frac {x^{3}}{3}\right )+c_{1}\right )}{x^{2}}\]