\[ a \left (x y'(x)-y(x)\right )^2-b x^2+x^2 y''(x)=0 \] ✓ Mathematica : cpu = 0.363928 (sec), leaf count = 134
\[\left \{\left \{y(x)\to x \left (\int _1^x\frac {i \sqrt {a} \sqrt {b} Y_1\left (-i \sqrt {a} \sqrt {b} K[2]\right )-i \sqrt {a} \sqrt {b} J_1\left (i \sqrt {a} \sqrt {b} K[2]\right ) c_1}{a \left (Y_0\left (-i \sqrt {a} \sqrt {b} K[2]\right )+J_0\left (i \sqrt {a} \sqrt {b} K[2]\right ) c_1\right ) K[2]}dK[2]+c_2\right )\right \}\right \}\] ✓ Maple : cpu = 1.323 (sec), leaf count = 72
\[\left \{y \left (x \right ) = \left (c_{2}+\int -\frac {\sqrt {-a b}\, \left (c_{1} \BesselY \left (1, \sqrt {-a b}\, x \right )+\BesselJ \left (1, \sqrt {-a b}\, x \right )\right )}{\left (c_{1} \BesselY \left (0, \sqrt {-a b}\, x \right )+\BesselJ \left (0, \sqrt {-a b}\, x \right )\right ) a x}d x \right ) x\right \}\]