\[ y(x) (-a-x \tan (x))+x^2 y''(x)-\left (2 x^2 \tan (x)-x\right ) y'(x)=0 \] ✓ Mathematica : cpu = 0.106272 (sec), leaf count = 30
DSolve[(-a - x*Tan[x])*y[x] - (-x + 2*x^2*Tan[x])*Derivative[1][y][x] + x^2*Derivative[2][y][x] == 0,y[x],x]
\[\left \{\left \{y(x)\to c_1 \sec (x) J_{\sqrt {a}}(x)+c_2 \sec (x) Y_{\sqrt {a}}(x)\right \}\right \}\] ✓ Maple : cpu = 0.045 (sec), leaf count = 24
dsolve(x^2*diff(diff(y(x),x),x)-(2*x^2*tan(x)-x)*diff(y(x),x)-(x*tan(x)+a)*y(x)=0,y(x))
\[y \left (x \right ) = \frac {c_{2} \BesselY \left (\sqrt {a}, x\right )+c_{1} \BesselJ \left (\sqrt {a}, x\right )}{\cos \left (x \right )}\]